%本程序用來求兩個離散信號的卷積y(n)=x(n)*h(n)，并且繪制出x(n)、h(n)和y(n)的波形圖
%
%執行本程序時，程序提示你輸入x(n)、h(n)的表達式，形式如下：
%      x(n)=[1,4,3,5,1,2,3,5]
%      nx=-4:3  nx表示序列  x(n)中各樣值對應的橫坐標n          
%      
%      h(n)=[4,2,4,0,4,2]
%      nh=-4:2  nh表示序列  h(n)中各樣值對應的橫坐標n 
%
%輸入兩個信號
x=input('請輸入序列x(n)=')
nx=input('請給定x(n)各樣值對應的橫坐標范圍nx=')
h=input('請輸入序列h(n)=')
nh=input('請給定h(n)各樣值對應的橫坐標范圍nh=')

%求卷積
[y,ny]=convm(x,nx,h,nh)

a=min(nx);
b=max(nx);
c=min(ny);
d=max(ny);
e=min(nh);
f=max(nh);
x1=zeros(1,a-c);
x2=zeros(1,d-b);
x3=[x1 x x2]
h1=zeros(1,e-c);
h2=zeros(1,d-f);
h3=[h1 h h2]

subplot(3,1,1)
stem(ny,x3)
title('輸入x(n)');
subplot(3,1,2)
stem(ny,h3)
title('系統沖激響應h(n)');
subplot(3,1,3)
stem(ny,y)
title('輸出y(n)');