%Samp4_2
b=0.15;a=[1 -0.8];    %濾波器分子和分母多項式系數
N=100;dt=1;n=0:N-1;t=n*dt;%時間序列
x=2*sin(0.05*pi*t)+0.2*randn(1,N);  %濾波器的輸入信號
imp=[1,zeros(1,N-1)];   %脈沖序列
h=filter(b,a,imp);    %輸入脈沖序列得到脈沖響應
yc=conv(h,x);    %脈沖響應與輸入信號卷積
y=yc(1:N);      %由于長度與輸入信號長度不相等，截斷使之相等
y1=filter(b,a,x);   %運用濾波器函數求得輸入x后系統的輸出
subplot(3,1,1),
plot(t,x);ylabel('x');title('輸入信號');  %繪制輸入信號
subplot(3,1,2),plot(t,y); %繪制卷積輸出信號
title('卷積得到的輸出信號')
subplot(3,1,3), plot(t,y1);  %繪制濾波器輸出信號
title('filter函數得到的輸出信號');xlabel('時間/s')