信號系統matlab實例2——matlab求解LTI系統響應、matlab求解LTI衝激響應和階躍響應

求解系統響應

求解系統零狀態響應函數lsim

y=lsim（sys，f，t）

t：表示座標範圍及單位長度。
f：表示輸入信號的函數。
sys：表示系統的微分方程。

獲取系統的微分方程tf

sys=tf（b，a）

b：微分方程等號右端的係數向量。
a：表示微分方程等號左邊的係數向量。
例子：

sys=tf([1],[1,2,77]);%tf函數獲得系統模型sys
t=0:0.01:5;%對t進行離散抽樣（0~5單位長度爲0.01）
f=10*sin(2*pi*t);%輸入函數
y=lsim(sys,f,t);%求系統的零狀態響應
plot(t,y);%畫圖
xlabel('Time(sec)');
ylabel('y(t)');

求解衝激響應和階躍響應

衝激響應函數impulse

impulse（b，a）

階躍響應函數step

step（b，a）

a，b表示方程左右兩端的係數向量

a=[7 4 6];%構造係數向量
b=[1 1];
subplot(2,1,1)
impulse(b,a);%求系統的衝激響應並做圖
subplot(2,1,2)
step(b,a);%求系統的階躍響應並做圖