DCT.rar_C++ dct图像处理_DCT图像处理_dct 图像处理

#include<stdio.h> #include<math.h> #define PI 3.141592653 float Construct(float y[][8]); float *Turnover(float a[][8], int, int ); float Multiplication(float a[][8], float b[][8], int, int, float c[][8] ); void main() { int i,j,m,n; float a[][8] = {{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8}}; float b[][8] = {{1,1,1,1,1,1,1,1},{2,2,2,2,2,2,2,2},{3,3,3,3,3,3,3,3},{4,4,4,4,4,4,4,4},{5,5,5,5,5,5,5,5},{6,6,6,6,6,6,6,6},{7,7,7,7,7,7,7,7},{8,8,8,8,8,8,8,8}}; float c[8][8] = {{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0}}; float d[8][8]; float y[8][8]; float y1[8][8]; float M[8][8] = {{1,-2,3,4,5,6,7,8},{1,2,-3,4,5,6,7,8},{1,2,3,-4,5,6,7,8},{1,2,3,4,-5,6,7,8},{1,2,3,4,5,-6,7,8},{1,2,3,4,5,6,-7,8},{1,2,3,4,5,6,7,-8},{-1,2,3,4,5,6,7,8}}; for( i=0;i<8;i++) { for( j=0;j<8;j++) { printf("%7.4f ",M[i][j]); } printf("\n"); } printf("\n\n"); m = 8; n = 8; Construct(y); //Turnover( a, d, m, n); y1 = *Turnover( y, m, n); Multiplication( y, M, m, n, c); Multiplication( c, y1, m, n, c); //---------------------------打印结果-------------------------- for( i=0;i<8;i++) { for( j=0;j<8;j++) { printf("%7.4f ",y[i][j]); } printf("\n"); } printf("\n\n"); // for(i=0;i<n;i++) { for( j=0;j<m;j++) printf("%7.4f ",y1[i][j]); printf("\n"); } printf("\n\n"); // for(i=0;i<m;i++) { for( j=0;j<m;j++) printf("%7.4f ",c[i][j]); printf("\n"); } printf("\n\n"); } //------------------------------构建数组------------------------------ float Construct(float y[][8]) { int i,j; int x = 8; for( i=0;i<8;i++) { for( j=0;j<8;j++) { y[i][j] = cos(2*PI*i*j/(2*x-1)); y[i][j] = y[i][j]/2; } y[0][i] = y[i][0]/sqrt(2); y[i][0] = y[0][i]; } return y[8][8]; } //-----------------------矩阵的转置---------------------------------- float *Turnover(float y[][8], int m,int n) { int i,j; float y1[8][8]; for(i=0;i<n;i++) { for( j=0;j<m;j++) y1[i][j] = y[j][i]; } return y1[]; } //-----------------------矩阵的乘法---------------------------------- float Multiplication(float a[][8],float b[][8],int m,int n, float c[][8]) { int i,j,k; for( i=0;i<m;i++) for( j=0;j<m;j++) for( k=0;k<n;k++) c[i][j] = c[i][j]+a[i][k]*b[k][i]; return c[8][8]; } /* #include<stdio.h> #include<math.h> #define PI 3.141592653 float Construct(float y[][8]); float Turnover(float a[][8], float b[][8], int, int ); float Multiplication(float a[][8], float b[][8], int, int, float c[][8] ); void main() { int i,j,m,n; float a[][8] = {{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8},{1,2,3,4,5,6,7,8}}; float b[][8] = {{1,1,1,1,1,1,1,1},{2,2,2,2,2,2,2,2},{3,3,3,3,3,3,3,3},{4,4,4,4,4,4,4,4},{5,5,5,5,5,5,5,5},{6,6,6,6,6,6,6,6},{7,7,7,7,7,7,7,7},{8,8,8,8,8,8,8,8}}; float c[8][8] = {{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0},{0,0,0,0,0,0,0,0}}; float d[8][8]; float y[8][8]; float y1[8][8]; float M[8][8] = {{1,-2,3,4,5,6,7,8},{1,2,-3,4,5,6,7,8},{1,2,3,-4,5,6,7,8},{1,2,3,4,-5,6,7,8},{1,2,3,4,5,-6,7,8},{1,2,3,4,5,6,-7,8},{1,2,3,4,5,6,7,-8},{-1,2,3,4,5,6,7,8}}; for( i=0;i<8;i++) { for( j=0;j<8;j++) { printf("%7.4f ",M[i][j]); } printf("\n"); } printf("\n\n"); m = 8; n = 8; Construct(y); //Turnover( a, d, m, n); Turnover( y, y1, m, n); Multiplication( y, M, m, n, c); Multiplication( c, y1, m, n, c); //---------------------------打印结果-------------------------- for( i=0;i<8;i++) { for( j=0;j<8;j++) { printf("%7.4f ",y[i][j]); } printf("\n"); } printf("\n\n"); // for(i=0;i<n;i++) { for( j=0;j<m;j++) printf("%7.4f ",y1[i][j]); printf("\n"); } printf("\n\n"); // for(i=0;i<m;i++) { for( j=0;j<m;j++) printf("%7.4f ",c[i][j]); printf("\n"); } printf("\n\n"); } //------------------------------构建数组------------------------------ float Construct(float y[][8]) { int i,j; int x = 8; for( i=0;i<8;i++) { for( j=0;j<8;j++) { y[i][j] = cos(2*PI*i*j/(2*x-1)); y[i][j] = y[i][j]/2; } y[0][i] = y[i][0]/sqrt(2); y[i][0] = y[0][i]; } return y[8][8]; } //-----------------------矩阵的转置---------------------------------- float Turnover(float y[][8],float y1[][8], int m,int n) { int i,j; for(i=0;i<n;i++) { for( j=0;j<m;j++) y1[i][j] = y[j][i]; } return y1[8][8]; } //-----------------------矩阵的乘法---------------------------------- float Multiplication(float a[][8],float b[][8],int m,int n, float c[][8]) { int i,j,k; for( i=0;i<m;i++) for( j=0;j<m;j++) for( k=0;k<n;k++) c[i][j] = c[i][j]+a[i][k]*b[k][i]; return c[8][8]; } */