原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Contents
DCT变换编码C语言
MPEG4中逆DCT变换
DCT变换探究
快速DCT变换
DCT变换编码C语言
#include memory.h #include stdio.h #include math.h #include time.h #define PI 3.1415926 #define CLK_TCK CLOCKS_PER_SEC int N; void DCT(double *f,double *F) { ??int n,m,x; ?double *dTemp = new double[N*N];//中间矩阵 ?double *coff = new double[N*N];//变换系数 ?coff[0] = 1/sqrt(N); ?for( m=1; mN; m++ ) ??coff[m] = sqrt(2)/sqrt(N); ?memset( dTemp, 0, sizeof(double)*N*N ); ?memset( F, 0, sizeof(double)*N*N );
?//一维变换 ?for(n=0;nN;n++) ??for(m=0;mN;m++) ???for(x=0;xN;x++) ????dTemp[m*N+n] += f[x*N+n] * coff[m] * cos( (2*x+1) * PI * m/(2*N) );
?//第二次一维变换 ?for(m=0;mN;m++)c ??for(n=0;nN;n++) ???for(x=0;xN;x++) ????F[m*N+n] += dTemp[m*N+x] * coff[n] * cos( (2*x+1) * PI * n/(2*N) );
?delete []dTemp; ?delete []coff; }
void iDCT(double *f,double *F) { ?int m,y,x; ?double *dTemp=new double[N*N];//中间矩阵 ?double *coff=new double[N*N];//变换系数 ?coff[0]=1/sqrt(N); ?for(m=1;mN;m++) ??coff[m]=sqrt(2)/sqrt(N); ?memset(dTemp,0,sizeof(double)*N*N); ?memset(F,0,sizeof(double)*N*N);
?//一维变换 ?for(x=0;xN;x++) ??for(y=0;yN;y++) ???for(m=0;mN;m++) ????dTemp[x*N+y]+=F[x*N+m]*coff[m]*cos((2*y+1)*PI*m/(2*N));
?//第二次一维变换 ?for(y=0;yN;y++) ??for(x=0;xN;x++) ???for(m=0;mN;m++) ????F[x*N+y]+=dTemp[m*N+y]*coff[m]*cos((2*x+1)*PI*m/(2*N));
?delete []dTemp; ?delete []coff;
}
int main() {? ??clock_t start,end; ?start=clock(); ?int i; ?long L; ? ?printf(变换维数:); ?scanf(%d,N); ? ?double *f=new double[N*N];//初始矩阵 ?double *F=new double[N*N];//变换后输出矩阵?
memset(F,0,sizeof(double)*N*N);//初始化为0
?for(i=0;iN*N;i++) ?{ ??printf(f[%d][%d]:,i/N,i%N); ??scanf(%lf,f[i]); ?}
?printf(循环次数:); ?scanf(%d,L);
//输出初始矩阵 ?printf(变换前:\n); ?for(i=1;i=N*N;i++) ?{ ??printf(%f\t,f[i-1]); ??if(i%N==0) ???printf(\n); ?}
?for(i=0;iL;i++) ??DCT(f,F);//变换 //输出变换后矩阵 ?printf(变换后:\n);
?for(i=1;i=N*N;i++) ?{ ??printf(%f\t,F[i-1]); ??if(i%N==0) ???printf(\n); ?}?
for(i=0;iL;i++) ?iDCT(f,F); ?//输出反变换后矩阵 ?printf(反变换后:\n); ?for(i=1;i=N*N;i++) ?{??printf(%f\t,f[i-1]); ??if(i%N==0) ???printf(\n); ?} ?//pr
文档评论（0）