长沙做网站智投未来,浏览器网站免费进入,成都 网站设计公司,深圳网页设计兴田德润i简介function yMyFFT_FB(x,n)%MYFFT_TB:My Fast Fourier Transform Frequency Based%按频率抽取基2-fft算法%input:% x -- 输入的一维样本% n -- 变换长度#xff0c;缺省时 nlength(x) 当n小于x数据长度时#xff0c;x数据被截断到第n个数据 % 当n大于时#xff0c;x数据在尾部…function yMyFFT_FB(x,n)%MYFFT_TB:My Fast Fourier Transform Frequency Based%按频率抽取基2-fft算法%input:% x -- 输入的一维样本% n -- 变换长度缺省时 nlength(x) 当n小于x数据长度时x数据被截断到第n个数据 % 当n大于时x数据在尾部补0直到 x 含n个数据%output:% y -- 1*n的向量快速傅里叶变换结果%variable define:% N -- 一维数据x的长度% xtem -- 临时储存x数据用% m,M -- 对N进行分解 N2^m*M,M为不能被2整除的整数% two_m -- 2^m% adr -- 变址1*N的向量% l -- 当前蝶形运算的级数% W -- 长为 N/2的向量记录 W(0,N),W(1,N),...W(N/2-1,N)% d -- 蝶形运算两点间距离% t -- 第l级蝶形运算含有的奇偶数组的个数% mul -- 标量乘数% ind1,ind2 -- 标量下标% tem -- 标量用于临时储存%参考文献% http://www.wendangwang.com/view/fea1e985b9d528ea81c779ee.html%% 输入参数个数检查msgnargchk(1,2,nargin);error(msg);%% 输入数据截断或加0Nlength(x);if nargin2if Nxtemx;xzeros(1,n);x(1:N)xtem;Nn;else % 截断xtemx;xxtem(1:n);Nn;endend
