%===================================================================
%                   程序名：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');
%==========================================================================