信号处理之FFT与CZT变换

FFT

用FFT算法计算序列x(n)=[2，1，3，2，1，5，1]与h(n)=[1，2，-1，-3]的线性卷积，画出输入、输出序列的波形图。

x=[2 1 3 2 1 5 1];  
h=[1 2 -1 -3];  
N=length(x)+length(h)-1;
n=0:N-1;
x=[x,zeros(1,N-length(x))];
h=[h,zeros(1,N-length(h))];
X=fft(x);
H=fft(h);
Y=X.*H;
y=ifft(Y);
subplot(1,3,1);
stem(n,x,'.');
xlabel('n');
ylabel('x(n)');
title(