1、int i,j;int aaMN;for(i=0;iM;i+)for(j=0;jN;j+)aaij=(int)(aij);/printf(%d,aaij);n);void Fast(float aNM, float bMN)/aji=bij;void Mult(float aMN,float bMN,float cMN)/DCT任int i,j,k;float sum;sum=0;for(k=0;kk+)sum+=aik*bkj;cij=sum;void IMult(float aMN,float bMN,float cMN)/DCT任sum+=akj*bik;void LH(float aM
2、N,float bMN,float cMN)/cij=2*aij/bij;void FLH(float aMN,float bMN,float cMN)/cij=aij*bij/2;void main()/float DCTMN,TDCTNM;float tempMN,FMN;float Ck;float qMN;/float fMN;/DCT任float ffMN;/DCT任float lhMN;/float flhMN;/FILE *fp=fopen(YuanShi.txt,r),*fp2=fopen(Q.txtif(fp=NULL)/error!fscanf(fp,%f,&fij);fc
3、lose(fp);if(fp2=NULL)/fscanf(fp2,qij);fclose(fp2);if(i=0) Ck=1.000000/sqrt(2);else Ck=1;DCTij=sqrt(2.000000/N)*Ck*cos(2*j+1)*i*P/(2*N);/DCTFast(TDCT,DCT);/Mult(DCT,f,temp);/DCT任Mult(temp,TDCT,F);LH(F,q,lh);/FLH(lh,q,flh);/IMult(DCT,flh,temp);/DCT任IMult(temp,TDCT,ff);freopen(DCT.txtw,stdout);/DCT任FDCT.txtPrint(F);fclose(stdout);LH.txt/lhLH.txtPrint(lh);FLH.txt/flhFLH.txtPrint(flh);IDCT.txt/DCT任ffIDCT.txtPrint(ff);