a)输出原始数据.docVIP

#include stdio.h #include math.h #define PI 3.1415926 #define N 8 void yuanshi int a[][8] //输出原始数据 int i,j; for i 0;i 8;i++ for j 0;j 8;j++ printf %d ,a[i][j] ; if j 7 printf \n ; void DCT int f[8][8],float m[8][8] //DCT变换子程序 int x,y,i,j,l; float coff[8][8],doff[8][8],D[8][8] 0 ; for y 0;y N;y++ //计算第一行的coff coff[0][y] sqrt 2.0 /sqrt N * sqrt 1.0 /sqrt 2.0 ; for x 1;x N;x++ //计算后面N-1行的coff for y 0;y N;y++ coff[x][y] sqrt 2.0 /sqrt N *cos 2* y+1 -1 *x*PI/ 2*N ; printf COFF: \n ; //输出矩阵coff for i 0;i N;i++ for j 0;j N;j++ printf %f ,coff[i][j] ; if j 7 printf \n ; for x 0;x N;x++ //求coff的转置矩阵doff for y 0;y N;y++ doff[x][y] coff[y][x]; for i 0;i N;i++ //矩阵doff和f矩阵相乘，结果为矩阵D for l 0;l N;l++ D[i][l] 0; for j 0;j N;j++ D[i][l] D[i][l]+doff[i][j]* 1.0*f[j][l] ; printf \n ; for i 0;i N;i++ //矩阵D和矩阵coff相乘,结果为矩阵m，即经过DCT变换//得到的矩阵 for l 0;l N;l++ m[i][l] 0; for j 0;j N;j++ m[i][l] m[i][l]+D[i][j]*coff[j][l]; printf 输出DCT矩阵\n ; for i 0;i N;i++ //输出经过DCT变换的结果矩阵m for l 0;l N;l++ printf %f ,m[i][l] ; printf \n ; void sch int te[8][8],float m[8][8],int n[8][8] //量化子程序 int x,y; int tr[8][8] ; for x 0;x N;x++ for y 0;y N;y++ tr[x][y] int m[x][y]/te[x][y] ; n[x][y] tr[x][y]; printf %d ,tr[x][y] ; if y 7 printf \n ; int bianma int n[8][8],int c[] //一维行程编码 int k[64],a[128],i,j,t 0,e 1,d 0,l; for i 0;i N;i++ //把二维的编码数组转换成一维数组 for j 0;j N;j++ k[t++] n[i][j]; for i 0;i 64;i++ //计算编码后的数据,将其存在数组a中 l i+1; if k[i] k[l] e++; else a[d++] e; a[d++] k[i]; e 1; printf 编码输出:\n ; for i 0;i d;i++ //输出编码后的数据 printf %d ,a[i] ; c[i] a[i];//把得到的编码后的数组a的值复制到数组c中 printf \n ; return d;//返回数组c的长度 void main /*8*8图像块数据如下*/ float m[8][8],b[8][8]; int Y 128,i 0,j 0,s,number,n[8][8],c[52],data; int a[][8] 130,130,130,129,134,133,129,130 , 130,130,130,129,134,133,130,130 , 130,130,130,129,132,132,130,130 , 129,130,130,129,130,130,129,129 , 127,128,127,129,131,129,131,130 , 127,128,127,128,127,128,132,132 , 125,126,129,129,127,129,133,132 , 127,125