| 查看: 1111 | 回复: 0 | |||
[资源]
【分享】二维傅里叶变换c程序
|
|
#include #include #include #include"FFT.h" //申请二维动态数组的函数 void MySpace(float ***p, int m, int n) { int i; *p=(float **)calloc(m,sizeof(float *)); for(i=0; i } //释放申请的二维动态数组 void FreeMySpace(float ***p, int m) { int i; for(i=0; i free(*p); } //矩阵转置a[m][n]->b[n][m] void Zhuan(float ***a,float ***b,int m,int n) { int i,j; for(i=0; i } //2D-FFT tag=1,正变换;tag=-1,反变换 void FFT2D(float **xr,float **xi,int m,int n,int tag) { int i,j,k1,k2; float **xtr, **xti; MySpace(&xtr,n,m); MySpace(&xti,n,m); k1=(int)(log(n)/log(2)+0.5); k2=(int)(log(m)/log(2)+0.5); for(i=0; i Zhuan(&xr,&xtr,m,n); Zhuan(&xi,&xti,m,n); for(i=0; i Zhuan(&xtr,&xr,n,m); Zhuan(&xti,&xi,n,m); FreeMySpace(&xtr,n); FreeMySpace(&xti,n); } void main() { int i,j; const int M=64; const int N=256; float **xr, **xi; FILE *fp1, *fp2; MySpace(&xr,M,N); MySpace(&xi,M,N); fp1=fopen("2D-test-64-256.dat","rb" ); fp2=fopen("2D-64-256.dat","wb" ); for(i=0; i FFT2D(xr,xi,M,N,1); for(i=0; i FreeMySpace(&xr,M); FreeMySpace(&xi,M); fclose(fp1); fclose(fp2); } 原始图形：  变换结果：  [ Last edited by baobiao007 on 2011-2-13 at 17:54 ] |
回复此楼» 猜你喜欢
诚邀
已经有0人回复
-
考研求调剂
已经有4人回复
-
地球物理学和空间物理学论文润色/翻译怎么收费?
已经有285人回复
-
考研求调剂
已经有2人回复
-
东北石油大学地球科学学院地质资源与地质工程专业接收调剂
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
【分享】一天征服傅里叶变换【已搜索无重复】
已经有242人回复
- 准晶：开除、胡说与擦肩而过
已经有18人回复
- 【分享】《傅里叶变换红外光谱分析》第二版【已搜索无重复】
已经有60人回复
- 求助：能不能帮忙写一个用matlab求一个函数的傅里叶变换的程序
已经有16人回复
- matlab 傅里叶变换 频谱分析
已经有12人回复
- 【求助】紧急求助：三角函数的傅里叶变换
已经有14人回复
- 【求助】原位傅里叶变换红外光谱仪
已经有5人回复
- 求助C语言有限元一维、二维、矩形网格以及三角形网格程序
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
100