exa3_2.m

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

%  exa3-2_casfilter.m, for example 3-2
% to test casfilter.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];
% 
[b0,B,A]=dir2cas(b,a) %IIR系统结构的直接型参数转换为级联型参数
[b,a]=cas2dir(b0,B,A) %系统结构的级联型参数转换为直接型参数
y=casfiltr(b0,B,A,x); %利用级联型系统结构求系统的单位阶跃响应
%y=filter(b,a,x);
% 求所给系统的输出，本例实际上是求所给系统的阶跃响应；
plot(t,x,'r.',t,y,'k-');
grid on;
ylabel('x(n) 与 y(n)')
xlabel('n')