快速傅立叶变换程序(更新版)..doc

1．加图标ID_FFT，类向导在视类添加相应函数 2．BMP文件中增加： BITMAPINFO* lpDIB = NULL; //原时域图像 BITMAPINFO* lpDIB_FFT = NULL; //频谱图像 static DWORD LineBytes = 0; //图像每行像素所占字节数 static LONG size = 0; //图像实际数据所占字节数 注意：在LoadBmpFile函数中 //删除LONG size= sizeof(BITMAPINFOHEADER) + NumColors * sizeof(RGBQUAD) + ImgSize;中的LONG //删除 DWORD LineBytes = (bi.biWidth * bi.biBitCount + 31)/32 * 4;中的DWORD //在“if (NULL == (lpBitsInfo = (LPBITMAPINFO) malloc(size))) return FALSE;”后面给lpDIB分配内存“if (NULL == (lpDIB= (LPBITMAPINFO) malloc(size))) return FALSE; “ //在“fread((char*)lpBitsInfo, size, 1, fp);”后面添加“lpDIB = lpBitsInfo” 在视类OnDraw前中增加： extern BITMAPINFO *lpDIB_FFT; extern BITMAPINFO *lpDIB; OnDraw中增加： if (lpDIB_FFT) { lpBits = (LPVOID)lpDIB_FFT-bmiColors[lpDIB_FFT-bmiHeader.biClrUsed]; StretchDIBits( pDC-GetSafeHdc(), 600,0, lpDIB_FFT-bmiHeader.biWidth, lpDIB_FFT-bmiHeader.biHeight, 0,0, lpDIB_FFT-bmiHeader.biWidth, lpDIB_FFT-bmiHeader.biHeight, lpBits, lpDIB_FFT, // bitmap data DIB_RGB_COLORS, SRCCOPY); } 3．在视类添加 void Fourier(); void CFFTView::OnFft() { if (lpDIB_FFT) free(lpDIB_FFT); Fourier(); Invalidate(); } void CFFTView::OnUpdateFft(CCmdUI* pCmdUI) { pCmdUI-Enable(lpDIB != NULL); } 4．BMP文件前面添加 #include complex using namespace std; （PPT中Copy） #define PI 3.1415926535 4．BMP文件中添加函数void Fourier(){} 5．BMP文件中添加函数一维傅里叶变换函数FFT() (红色加粗是需要自己写的，其他的为PPT中Copy) void FFT(complexdouble * TD, complexdouble * FD, int r) { // 计算付立叶变换点数 LONG count = 1 r; // 计算加权系数 int i; double angle; complexdouble* W = new complexdouble[count / 2]; for(i = 0; i count / 2; i++) { angle = -i * PI * 2 / count; W[i] = complexdouble (cos(angle), sin(angle)); } // 将时域点写入X1 complexdouble* X1 = new complexdouble[count]; memcpy(X1, TD, sizeof(complexdouble) * count); // 采用蝶形算法进行快速付立叶变换，输出为频域值X2 complexdouble* X2 = new complexdouble[count]; int k,j,p,size; complexdouble* temp; for (k = 0; k r; k++) { for (j = 0; j 1 k; j++) { size = 1 (r-