概述
处女作
为什么写
从一个脚本坑坑洼洼地补,最终成功运行,有些经验,有些被忽视的小问题其实很致命,我找不到人说。。。有意向成为一名程序员,目前大三,正在学习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实现为什么写第一份代码分享所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
发表评论 取消回复