fu
ctio
yMyFFT_TBx
MYFFT_TBMyFastFourierTra
sformTimeBased按时间抽取基2fft算法i
putx输入的一维样本
变换长度，缺省时
le
gthx当
小于x数据长度时，x数据被截断到第
个数据当
大于时，x数据在尾部补0直到x含
个数据outputy1
的向量，快速傅里叶变换结果variabledefi
eN一维数据x的长度xtem临时储存x数据用mM对N进行分解N2mMM为不能被2整除的整数two_m2madr变址，1N的向量l当前蝶形运算的级数W长为N2的向量，记录W0NW1NWN21Nd蝶形运算两点间距离t第l级蝶形运算含有的奇偶数组的个数mul标量，乘数i
d1i
d2标量，下标tem标量，用于临时储存参考文献：httpwe
kubaiducomviewfea1e985b9d528ea81c779eehtml输入参数个数检查msg
argchk12
argi
errormsg输入数据截断或加0Nle
gthxif
argi
2ifN
加0xtemxxzeros1
x1NxtemN
else截断xtemxxxtem1
N
e
de
d
f对N进行分解N2mMmMfactorizeNtwo_mNM变换ifm0如果N可以被2整除adraddressmMtwo_myxadr1蝶形运算级数lm时ifM1N分解含有非2因数M时，对y中每M个数据做直接傅里叶变换forii1two_myii1M1iiMDDFTyii1M1iiMe
de
d计算W向量Wexp2pii0N21N蝶形运算dMtN2dforlm11乘forr0d1mulWrt1forii0t1yii2dd1ryii2dd1rmule
de
d加forii0t1i
d1ii2d1i
d2i
d1dforr0d1temyi
d1yi
d2yi
d2yi
d1yi
d2yi
d1temi
d1i
d11i
d2i
d21e
de
d
fd2dtt2e
delse如果N不能被2整除yDDFTxe
de
d内嵌函数fu
ctio
yDDFTx直接离散傅里叶变换i
putx样本数据，N维向量outputyN维向量参考文献：结构动力学，克拉夫，P82variabledefi
essum，用于求和Nle
gthxyzerossizexfor
1Ns0form1Nssxmexpi2pim1
1Ne
dy
se
de
dfu
ctio
mMfactorizeN对N分解m0whiletrueifmodN20mm1NN2elseMNbreak
fe
de
de
dfu
ctio
adraddressmMtwo_m变址b2mm的矩阵，用来存储二进制数据ds数，公差adrzerostwo_mMbde2bi0two_m1m转换为2进制注：matlab中二进制0116bbe
d11逆序adr1bi2deb2进制转换为10进制ifM1dstwo_madradr1o
es1Madradrdso
essizeadr110M1adrreshapeadr1e
de
d联系方式：matrixsuperqqcom
fr