ch2example7prg1.m

一些有用的MATLAB通信信号处理的源代码程序

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