respon.m

经典《信号与系统》教程的matlab例程,对深入理解信号与系统相关概念有很大帮助

%===================================================================
%                   程序名：respon.m
% 本程序是用状态变量法求解系统响应，并验证手工计算的结果的通用程序
%===================================================================
clear;
t=0:.001:5;
a=input('请输入矩阵A：');
b=input('请输入矩阵B：');
c=input('请输入矩阵C：');
d=input('请输入矩阵D：');
x0=input('请输入状态变量的初始条件矩阵X（0）：');
et=input('请输入激励信号表达式：');
%=======================================================================
%以下程序段求解系统的零输入响应、零状态响应和完全响应，并绘制出响应曲线
%=======================================================================
clf;
figure(1);
[y1,x]=initial(a,b,c,d,x0,t);    %计算零输入响应
plot(t,y1,'r');hold on;
[y2,x]=lsim(a,b,c,d,et,t);       %计算零状态响应
plot(t,y2,'b');hold on;
[y3,x]=lsim(a,b,c,d,et,t,x0);    %计算完全响应
plot(t,y3,'g');hold on;grid on;
title('红——零输入响应rzi，蓝——零状态响应rzs，绿——完全相应r=rzi+rzs');
disp('按任意键继续');
pause;%暂停，按任意键继续执行程序
%=========================================================================
%以下程序段绘制手工求解的系统的零输入响应、零状态响应和完全响应曲线，供实
%验者验证自己计算的结果是否正确：
%=========================================================================
figure(2);
f1=input('请输入手工计算而得到零输入响应表达式：');
plot(t,f1,'r');hold on;
f2=input('请输入手工计算而得到零状态响应表达式：');
plot(t,f2,'b');hold on;
f3=input('请输入手工计算而得到完全响应表达式：');
plot(t,f3,'g');hold on;grid on;
title('红——零输入响应rzi，绿——零状态响应rzs，蓝——完全相应r=rzi+rzs');
%==========================================================================