exa3_1.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 16 行

%  exa3-1_filter.m, for example 3-1
% to test filter.m and to obtain the step response.

clear;
x=ones(1,100);
t=1:100;
b=[0.00183,0.00734,0.01101,0.00737,0.00183];
a=[1,-3.054,3.829,-2.292,.551];
% 
y=filter(b,a,x);
% 求所给系统的输出，本例实际上是求所给系统的阶跃响应；
plot(t,x,'r.',t,y,'k-');
grid on;
ylabel('x(n) 与 y(n)')
xlabel('n')