基于matlab的FFT算法实现实序列线性卷积方法一-新建 1文本文档.txt

h=input('请输入给定序列h=');
x=input('请输入序列x=');
lenx=length(x);
lenh=length(h);
N=lenx+lenh-1;
x=[x zeros(1,N-lenx)];
subplot(221);
stem(x);
h=[h zeros(1,N-lenh)];
subplot(222);
stem(h);
N=2^(ceil(log2(N)));
H=fft1(h);
X=fft1(x);
Y=H.*X;
y=real(ifft1(Y));
subplot(223);
stem(y);
z=conv(x,h);
subplot(224);
stem(z)
end