卷积的matlab程序.txt

对学好matlab很有用

function juanji(a,b)
%一般情况下输入a,b都是实数序列，这里忽略a,b为复序列的情况
clf;
subplot(2,2,1);
k=0:length(a)-1;
stem(k,a);
xlabel('n');ylabel('a');
grid;
subplot(2,2,2);
k=0:length(b)-1;
stem(k,b);
xlabel('n');ylabel('b');
grid;
N=max(length(a),length(b));
c=fft(a,N).*fft(b,N);
yuanjuan=ifft(c)
L=length(a)+length(b)-1;
c=fft(a,L).*fft(b,L);
xianjuan=ifft(c)
subplot(2,2,3);
k=0:N-1;
stem(k,yuanjuan);
xlabel('n');ylabel('圆周卷积(循环卷积)');
grid;
subplot(2,2,4);
k=0:L-1;
stem(k,xianjuan);
xlabel('n');ylabel('线性卷积');
grid;