MATLAB周期时域曲线,MATLAB系列之连续时间系统的时域分析
輸入-輸出描述法:
使用傳輸函數(shù)tf描述系統(tǒng)。
sys?= tf(b,a)
其中b表示微分方程右側(cè)(即激勵)系數(shù)的行矢量,a表示微分方程左側(cè)(即響應(yīng))系數(shù)的行矢量。
狀態(tài)變量分析法:
使用狀態(tài)空間函數(shù)ss描述系統(tǒng)。
sys = ss(a,b,c,d)
其中矩陣a表示狀態(tài)方程中由狀態(tài)組成的向量的系數(shù)矩陣;矩陣b表示狀態(tài)方程中由激勵組成的向量的系數(shù)矩陣;矩陣c表示輸出方程中由狀態(tài)組成的向量的系數(shù)矩陣;矩陣d表示輸出方程中由激勵組成的向量的系數(shù)矩陣。
微分方程的齊次解:
齊次解即系統(tǒng)特征方程的根,用roots函數(shù)計算。
a = roots(p) ,其中p為特征方程的系數(shù)由高至低排列構(gòu)成的行矢量,返回值a是特征根組成的列向量。
微分方程的特解:
特解即系統(tǒng)函數(shù)(微分方程)在給定激勵信號作用下的輸出,使用lsim函數(shù)對LTI系統(tǒng)進行仿真。
y =
lsim(sys,u,t),其中sys表示LTI系統(tǒng)(即上文求出的系統(tǒng)模型),矢量u和t分別表示激勵信號的抽樣值和抽樣時間,返回值y為對應(yīng)于上述抽樣時間的系統(tǒng)響應(yīng)值。一般情況,t設(shè)定為列向量。
矩陣u的行數(shù)要和向量t的長度(即采樣點數(shù))一樣,u的每一列即為一個輸入信號。
零輸入響應(yīng)和零狀態(tài)響應(yīng):
y = lsim(sys,u,t,x0) ,
其中sys表示LTI系統(tǒng),矢量u和t分別表示激勵信號的抽樣值和抽樣時間,矢量x0表示該系統(tǒng)的初始狀態(tài),返回值y表示系統(tǒng)響應(yīng)值。注意:該用法僅適于用狀態(tài)方程描述的系統(tǒng)。
[y,t,x] = lsim(sys,u,t,x0),y為系統(tǒng)輸出,t為系統(tǒng)抽樣時間,x為系統(tǒng)狀態(tài)。
沖激響應(yīng)和階躍響應(yīng):
(1)分別用沖激信號和階躍信號作為激勵,lsim函數(shù)可仿真出沖激響應(yīng)和階躍響應(yīng)。
(2)impulse(sys)和step(sys)直接產(chǎn)生LTI系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)。
卷積:
(1)離散卷積:w = conv(u,v)其中u和v是兩個要做卷積的序列。
(2)連續(xù)卷積:自己定義函數(shù)。
function [w,tw] = conv1(u,tu,v,tv)
%u,v表示兩個序列;tu,tu表示它們的抽樣時間
%w,wt分別表示卷積結(jié)果及其抽樣時間
T = tu(2) - tu(1);
w = T*conv(u,v);
tw = tu(1) + tv(1) + T*[0:length(u) + length(v) - 2];
return
總結(jié)
以上是生活随笔為你收集整理的MATLAB周期时域曲线,MATLAB系列之连续时间系统的时域分析的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 水滴筹怎么捐款
- 下一篇: 新手炒股需要了解股票基础知识,基础是重中