samp4_2.m

MATLAB7.x数字信号处理 光盘内容

%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')