我是靠谱客的博主 会撒娇向日葵,最近开发中收集的这篇文章主要介绍基2频率抽样FF的matlab实现为什么写第一份代码分享,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

处女作

为什么写

从一个脚本坑坑洼洼地补,最终成功运行,有些经验,有些被忽视的小问题其实很致命,我找不到人说。。。有意向成为一名程序员,目前大三,正在学习pythonmatlabc,本科专业为电子信息科学与技术,专业课在搞通信,所以matlab可能用的频次更高。以后也可能分享FPGA开发过程中的经验,拭目以待吧!

第一份代码分享

在学习快速傅里叶变换中,FFT算法听的头昏脑胀,不知所云,一时难以接受。通过研究程序搞明白了,程序总结了其迭代规律,得其精髓;
下面给出频率抽样FFT matlab代码。ps:fft在matlab属于内部函数,help也不告诉你

function y=DIF_FFT(x)
%输入自然序,输出倒位序
m=ceil(log2(length(x)));
N=2^m;
if length(x)<N
    x=[x,zeros(1,N-length(x))];
end
for mm=1:m  %级数控制
    Nmr=2^(m-mm); %与时域抽取不同,指数部分为:总级数-当前级
    u=1;
    N=exp(-i*2*pi/(2*Nmr));
    for j=1:Nmr %控制不同的叠加系数
     for k=j :2*Nmr:N %控制相同的系数有几次,每隔2*Nmr出现一个相同系数
            kp=k+Nmr;
            t=(x(k)-x(kp))*u;%一定要保留一个中间变量,否则修改xk与xkp
            %蝶形的递推由上级的两个值确定
            x(k)=x(k)+x(kp);
            x(kp)=t;
             end;
        u=u*WN;
    end
    end
nxd=bin2dec(fliplr(dec2bin([1:N]-1,m)))+1;  
y=x(nxd);%将输出次序修正为自然序

# 简单地结束语
写文档是门手艺,我在排版上有点差,第一次也就不折腾了,码农嘛,算法和思想才关键。
留个QQ吧,寻找码农共同进修啊
My:2991447145

最后

以上就是会撒娇向日葵为你收集整理的基2频率抽样FF的matlab实现为什么写第一份代码分享的全部内容,希望文章能够帮你解决基2频率抽样FF的matlab实现为什么写第一份代码分享所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(53)

评论列表共有 0 条评论

立即
投稿
返回
顶部