fftjuanji.m

来自「数字信号处理程序——利用FFT实现快速傅立叶卷积的matlab程序」· M 代码 · 共 20 行

%用FFT实现序列 x(n) = 0.8^n*R12(n) 和h(n) = R6(n)的线性卷积；
%实现程序：
n = [0:1:11];
m = [0:1:5];
N1 = length(n);
N2 = length(m);
%生成xn
xn = 0.8.^n;
%生成hn
hn = ones(1,N2);
N = N1 + N2 -1;
XK = fft(hn,N);
HK = fft(hn,N);
YK = XK.*HK;
yn = ifft(YK,N);
if all(imag(xn) ==0)&(all(imag(xn) ==0))
   yn = real(yn);
end
x = 0:N - 1;
stem(x,yn,'.')