快速傅里叶变换
分别用直接卷积和快速卷积两种方法求以下两个序列的卷积。
tic、toc,利用matlab计时,tic为起始时间,toc为结束时间,最终显示时间=toc-tic
直接卷积:conv,快速卷积:利用fft,利用时域的卷积等于频域的乘积
x(n)=0.2n*R19(n) , h(n)=(0.8)n*R15(n)。
%直接卷积与快速卷积所花时间的比较
% x(n)=0.2n*R19(n) , h(n)=(0.8)^n*R15(n)。Rx(n)长度为x+1
xn=0.2*(1:20);%生成xn序列,长度为20
hn=0.8.^(1:16);%生成hn序列,长度为16
tic;%利用matlab进行计时,tic表示计时开始,toc表示计时结束,显示时间=toc-tic
yn=conv(xn,hn);%利用直接卷积计算xn与hn的卷积
toc;%计时结束,输出直接卷积所花费的时间
L=pow2(nextpow2(16+20-1));%取快速卷积的长度L≥M+N-1,且保证L为2的整次幂
tic;%同理,获取快速卷积的时间
Xk=fft(xn,L);%用fft计算x(n)的离散傅里叶变换,长度为L,不足补零
Hk=fft(hn,L);%用fft计算h(n)的离散傅里叶变换
Yk=Xk.*Hk;%依次将Xk与Hk中元素相乘,则输出结果为yn的频域形式
yn=ifft(Yk,L);%利用ifft计算yn
toc;%计时结束
subplot(221),stem(xn,'.');ylabel(