DCT算法.docVIP

VOID ? WINAPI ? DCT(double ? *dpf, ? double ? *dpF, ? int ? r) ? ? { ? ? double ? PI ? = ? 3.1415926; ? ? // ? 离散余弦变换点数 ? ? LONG lNum; ?? ? // ? 循环变量 ? int i; ?? ? // ? 中间变量 ?? double dTemp; ? ? complexdouble ? *comX; ?? ?// ? 离散余弦变换点数 ?? lNum ? = ? 1r; ? // ? 分配内存 ?? comX ? = ? new ? complexdouble[lNum*2]; ? // ? 赋初值0 ?? memset(comX, ? 0, ? sizeof(complexdouble) ? * ? lNum ? * ? 2); ?? ? // ? 将时域点转化为复数形式 ?? for(i=0;ilNum;i++) ? ? { ?? comX[i] ? = ? complexdouble ? (dpf[i], ? 0); ? ? } ?? ? // ? 调用快速付立叶变换 ?? FFT_1D(comX,comX,r+1); ? ? ?? // ? 调整系数 ?? dTemp ? = ? 1/sqrt(lNum); ?? ? ? // ? 求dpF[0] ?? dpF[0] ? = ? comX[0].real() ? * ? dTemp; ?? ? ? dTemp ? *= ? sqrt(2); ? ? ?? // ? 求dpF[u] ?? for(i ? = ? 1; ? i ? ? lNum; ? i++) ? ? { ? ? dpF[i]=(comX[i].real() ? * ? cos(i*PI/(lNum*2)) ? ? + ? comX[i].imag() ? * ? sin(i*PI/(lNum*2))) ? * ? dTemp; ?? } ?? ? ? // ? 释放内存 ?? ? ? ? ? delete ? comX; ? ? } ? ? ??? VOID ? WINAPI ? FFT_1D(complexdouble ? * ? pCTData, ? complexdouble ? * ? pCFData, ? int ? nLevel) ? ? { ? ? int i; ? ? int ? j; ? ? int? k; ? ? ?? double ? PI ? = ? 3.1415926; ? ? ?? int nCount ? =0 ? ; ?? ? nCount ? =(int)pow(2,nLevel) ? ; ?? ? int nBtFlyLen; ? ? nBtFlyLen ? = ? 0 ? ; ?? ? double dAngle; ?? ? complexdouble ? *pCW ? ; ?? ? pCW ? ? = ? new ? complexdouble[nCount ? / ? 2];? ? ? for(i ? = ? 0; ? i ? ? nCount ? / ? 2; ? i++) ? ? { ? ? dAngle ? = ? -2 ? * ? PI ? * ? i ? / ? nCount; ? ? pCW[i] ? = ? complexdouble ? ( ? cos(dAngle), ? sin(dAngle) ? ); ? ? } ? ?? complexdouble ? *pCWork1,*pCWork2; ? ? ? ?? pCWork1 ? = ? new ? complexdouble[nCount]; ?? ? pCWork2 ? = ? new ? complexdouble[nCount]; ?? ? complexdouble ? *pCTmp; ?? ? memcpy(pCWork1, ? pCTData, ? sizeof(complexdouble) ? * ? nCount); ? ? ?? int ? nInter; ? ? ? nInter ? = ? 0; ?? ? ? for(k ? = ? 0; ? k ? ? nLevel; ? k++)? { ?? for(j ? = ? 0; ? j ? ? (int)pow(2,k); ? j++) ? ? { ?? ? nBtFlyLen ? = ? (int)pow( ? 2,(nLevel-k) ? );? ? ? for(i ? = ? 0; ? i ? ? nBtFlyLen/2; ? i++) ? ? { ? ?