www.gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第5章/li5_6.m

     clear;
n=5; %仿真计算的时间序列点数
f=[1,zeros(1,n-1)];%输入：单位数字冲激信号
x=zeros(2,n+1);%状态变量存储矩阵初始化
x(:,1)=[0;0]; %初始状态赋值
for i=1:n
    x(1,i+1)=-0.5.*(x(1,i)+f(i)); %状态方程1
    x(2,i+1)=f(n); %状态方程2
    y(n)=x(1,i)+x(2,i)+f(i);%输出方程
end
t=0:n-1; %得到序列对应的离散时间点并作出波形
subplot(411);
stem(t,f);%输入信号波形
axis([-1 n 0 1.5]);
subplot(412);
stem(t,x(1,1:n));%状态1的波形
axis([-1 n -0.6 0.6]);
subplot(413);
stem(t,x(2,1:n));%状态2的波形
axis([-1 n 0 1.5]);
subplot(414);
stem(t,y);%输入信号波形
axis([-1 n -0.5 1.2]);%输出信号波形