[原创]自己编写的FFT程序运行速度问题

limsnudt · 发表于 2006-3-29 12:45

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

这里是小弟自己编写的二维快速Fourier变换的程序（Matlab）,一维速度很快，但二维就很慢，请教各位大侠如何改进！！！（这是我的作业），谢谢！ 可以给我发邮件 limsnudt@sohu.com

[此贴子已经被aspen于2006-3-30 9:32:21编辑过]

limsnudt · 发表于 2006-3-29 17:54

function y=daoxu(a,n) %对输入的2的n次幂个点进行倒序运算 t=[1:2^n]; b=dec2bin(t-1,n); if rem(n,2)==0 m=n*0.5; else m=(n+1)*0.5; end b1=b.';%转置之后来交换行，看看速度会不会快起来 for k=1:m d=b1(k,:); b1(k,:)=b1(n+1-k,:);%换行 b1(n+1-k,:)=d; end b=b1.'; for l=1:2^n a1(l)=a(bin2dec(b(l,:))+1); end y=a1;
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=ft1(a) %已经进行倒序之后的一维快速Fourier变换，b为经过倒序之后的输入序列，N为凑为2的次幂的幂指数 %此处的Fourier正变换前面不带任何系数 M=length(a); N=log2(M); b=daoxu(a,N); for k=1:2:(2^N-1) c(k)=b(k)+b(k+1); c(k+1)=b(k)-b(k+1); end %以上计算的是最底层的单个点的Fourier变换 %具体的计算过程采用从下到上分成各个小部分进行计算的方法 for n=2:N for t=1:2^(N-n) m=2^(n-1); p=t*2^n-(2^n-1); u=0; while (p<=t*2^n-2^(n-1))&(u<=(m-1)) d(p)=c(p)+c(p+2^(n-1))*((exp(-j*2*pi/(2*m)))^u); d(p+m)=c(p)-c(p+2^(n-1))*((exp(-j*2*pi/(2*m)))^u); p=p+1; u=u+1; end end c=d; end d=[]; y=c;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=haofft(a) tic; %二维快速Fourier变换,即对一维的情况进行循环 n=length(a); b=zeros(n); c=zeros(n); for k=1:n b(k,:)=ft1(a(k,:)); end a=[]; for l=1:n c(:,l)=ft1(b(:,l)); end y=c; toc;

happy · 发表于 2006-3-29 19:26

matlab不是自带了吗？

账号		自动登录	找回密码
密码			我要加入

[原创]自己编写的FFT程序运行速度问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

不好意思，程序没有上来

回复：（limsnudt）[原创]请大侠赐教，小弟谢谢你！