按频率抽取基2-快速傅里叶逆变换算法_MATLAB代码

《按频率抽取基2-快速傅里叶逆变换算法_MATLAB代码》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、functionx=MyIFFT_FB(y)%MyIFFT_TB:MyInverseFastFourierTransformTimeBased%按频率抽取基2-傅里叶逆变换算法%input:%y--傅里叶正变换结果,1N的向量%output:%x--逆变换结果，1N的向量%参考文献：%http://wenku.baidu.com/view/fea1e985b9d528ea81c779ee.htmlN=length(y);x=conj(y);%求共轭x=MyFFT_FB(x);%求FFTx=conj(x);%求共轭x=x./N;%除以Nend%%内嵌函数===============

2、=======================================functiony=MyFFT_FB(x,n)%MYFFT_TB:MyFastFourierTransformFrequencyBased%按频率抽取基2-fft算法%input:%x--输入的一维样本%n--变换长度，缺省时n=length(x)当n小于x数据长度时，x数据被截断到第n个数据%当n大于时，x数据在尾部补0直到x含n个数据%output:%y--1n的向量，快速傅里叶变换结果%variabledefine:%N--一维数据x的长度%xtem--临时储存x数据用%m,M--对N进行分解N=2

3、^mM,M为不能被2整除的整数%two_m--2^m%adr--变址，1N的向量%l--当前蝶形运算的级数%W--长为N/2的向量，记录W(0,N),W(1,N),...W(N/2-1,N)%d--蝶形运算两点间距离%t--第l级蝶形运算含有的奇偶数组的个数%mul--标量，乘数%ind1,ind2--标量，下标%tem--标量，用于临时储存%参考文献：%http://wenku.baidu.com/view/fea1e985b9d528ea81c779ee.html%%输入参数个数检查msg=nargchk(1,2,nargin);error(msg);%%输入数据截断或加0N=

4、length(x);ifnargin==2ifN

5、;forr=0:d-1tem=y(ind1)+y(ind2);y(ind2)=y(ind1)-y(ind2);y(ind1)=tem;ind1=ind1+1;ind2=ind2+1;endend%乘forr=0:d-1mul=W(rt+1);forii=0:t-1y(ii2d+d+1+r)=y(ii2d+d+1+r)mul;endendd=d/2;t=t2;end%%直接傅里叶变换ifM~=1%N分解含有非2因数M时，对y中每M个数据做直接傅里叶变换forii=1:two_my((ii-1)M+1:iiM)=DDFT(y((ii-1)M+1:iiM));endend%%变址输出y=

6、y(adr+1);else%%如果N不能被2整除y=DDFT(x);endendfunctiony=DDFT(x)%%直接离散傅里叶变换%input:%x--样本数据，N维向量%output:%y--N维向量%参考文献：%结构动力学，克拉夫，P82%variabledefine%s--sum，用于求和N=length(x);y=zeros(size(x));forn=1:Ns=0;form=1:Ns=s+x(m)exp(-i2pi(m-1)(n-1)/N);endy(n)=s;endendfunction[m,M]=factorize(N)%%对N分解m=0;whiletrueif

7、mod(N,2)==0m=m+1;N=N/2;elseM=N;break;endendendfunctionadr=address(m,M,two_m)%%变址%b--2^mm的矩阵，用来存储二进制数据%ds--数，公差adr=zeros(two_m,M);b=de2bi(0:two_m-1,m);%转换为2进制注：matlab中二进制[011]=6b=b(:,end:-1:1);%逆序adr(:,1)=bi2de(b);%2进制转换为10进制ifM~=1ds=two