注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

being23

写给未来的自己

 
 
 

日志

 
 
关于我

真正的坚定,就是找到力量去做自己喜欢的事情,并为之努力,这样才会觉得生活是幸福的。

网易考拉推荐

FFT  

2009-04-06 23:33:14|  分类: 准备工作 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

%fft of x
%x is the input vector
%N is the number of the points
%at 024009 2009-4-3
function y = fft_custom(x,N)

 

%逆序操作
temp_reverse=[0:N-1];
LH=N/2;
J=LH;
N1=log2(N);
for I=temp_reverse(2):temp_reverse(N-1)
    if I<J
        T=temp_reverse(I+1);
        temp_reverse(I+1)=temp_reverse(J+1);
        temp_reverse(J+1)=T;
    end
    K=LH;

    while J>=K
        J=J-K;
        K=K/2;
    end
    J=J+K;
end
temp_reverse=temp_reverse+1;
for i=1:N
    y(i)=x(temp_reverse(i));
end

%蝶形计算
M=log2(N);
Wn=exp^(-j*2*pi/N);
for L=1:M                       %表示蝶形运算的级数
    B=2^(L-1);                  %每个蝶形输入数据相距的点数,和不同蝶形数目   
    for J=0:B-1;                %每级参与计算的蝶形数
        P=2^(M-L)*J;           %当前蝶形中旋转因子的指数          
        for k=J:2^L:N-1;        
            T=y(k+1)+y(k+1+B)*Wn^P;            
            y(k+1+B)=y(k+1)-y(k+1+B)*Wn^P;              
            y(k+1)=T;
        end
    end
end

 

   

  评论这张
 
阅读(291)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017