基于matlab的零輸入和零狀態(tài)響應(yīng)

一、 課程設(shè)計(jì)題目： 信號(hào)系統(tǒng)的時(shí)域分析 二、課程設(shè)計(jì)目的： 1、學(xué)習(xí)MATLAB 軟件的使用。 2、使學(xué)生掌握利用工具軟件來(lái)實(shí)現(xiàn)信號(hào)系統(tǒng)基本概念、基本原 理的方法。 3、通過(guò)編程對(duì) matlab 軟件的具體應(yīng)用有了更好的了解，進(jìn)一步 加強(qiáng)了對(duì)函數(shù)卷積，零輸入，零響應(yīng)狀態(tài)這三種函數(shù)狀態(tài)的理解。 二、 基本原理 1、卷積 信號(hào)的卷積是數(shù)學(xué)上的一種積分運(yùn)算，兩個(gè)信號(hào)的卷積定義為： 信號(hào)的卷積運(yùn)算在系統(tǒng)分析中主要用于求解系統(tǒng)的零狀態(tài)響應(yīng)。一 般情況，卷積積分的運(yùn)算比較困難，但在 MATLAB 中則變得十分 簡(jiǎn)單，MATLAB 中是利用 conv 函數(shù)來(lái)實(shí)現(xiàn)卷積的。 Conv 函數(shù)功能：實(shí)現(xiàn)二個(gè)函數(shù) 和 的卷積。 格式： 說(shuō)明： 表示二個(gè)函數(shù)， 表示兩個(gè)函數(shù)的卷積結(jié)果。 例：已知兩信號(hào)f1(t)= ε(t+1)+ε(t-2) f2(t)=ε(t-3)+ε(t-8) 求卷積MATLAB 程序如下： t1=-1:0.01:2; f1=ones(size(t1)); %高度為一的門函數(shù)，時(shí)間從t=-1到t=2 t2=3:0.01:8; f2=ones(size(t2)); %高度為一的門函數(shù)，時(shí)間從t=3到t=8 g=conv(f1,f2); %對(duì)f1和f2進(jìn)行卷積 t3=2:0.01:10; subplot(3,1,1),plot(t1,f1),title(‘f1’);grid %畫(huà)f1的波 形 subplot(3,1,2),plot(t2,f2) ,title(‘f2’);grid %畫(huà)f2的波形 subplot(3,1,3),plot(t3,g) ,title(‘f1*f2’);grid %畫(huà)g的波形2、零輸入，零狀態(tài)響應(yīng) 一般的連續(xù)時(shí)間系統(tǒng)分析有以下幾個(gè)步驟: ①求解系統(tǒng)的零輸 入響應(yīng); ②求解系統(tǒng)的零狀態(tài)響應(yīng); ③求解系統(tǒng)的全響應(yīng); ④分析 系統(tǒng)的卷積；⑤畫(huà)出它們的圖形. 下面以具體的微分方程為例說(shuō)明 利用MATLAB 軟件分析系統(tǒng)的具體方法. 2.1．連續(xù)時(shí)間系統(tǒng)的零輸入響應(yīng)描述n階線性時(shí)不變(LTI)連續(xù)系統(tǒng)的微分方程為：已知y及各階導(dǎo)數(shù)的初始值為y(0),y (1) (0)，… y (n-1) (0)， 求系 統(tǒng)的零輸入響應(yīng)。 當(dāng)LIT系統(tǒng)的輸入為零時(shí)，其零輸入響應(yīng)為微分方程的其次解 (即令微分方程的等號(hào)右端為零) ，其形式為(設(shè)特征根均為單根) 1 1 2 1 1 1 1 n n m n n m m n n m d y d y dy d u du a a a a y b b b u dt dt dt dt dt ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?其中 p 1 ,p 2 ,…,p n 是特征方程 a 1 λ n +a 2 λ n-1 +…+a n λ+a n =0 的根，它們 可以用root(a)語(yǔ)句求得。各系數(shù) 由y及其各階導(dǎo)數(shù)的初始值來(lái) 確定。對(duì)此有 寫(xiě)成矩陣形式為: P 1 n-1 C 1 + P 2 n-1 C 2 +…+ P n n-1 C n =D n-1 y 0 1 0 1 2 2 0 1 1 1 1 1 2 0 1 1 1 n n n n n n n C y p p p C Dy p p p C D y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? A A A A A A 即 V?C=Y 0 其解為：C=V\Y 0式中 V為范德蒙矩陣，在matlab的特殊矩陣庫(kù)中有vander。 以下面式子為例： y″(t)+3y′(t)+6y(t)=6f′(t)-8f′(t) 初始條件為y(0_)=0,y′(0_)=10; MATLAB程序： a=[1,3,6]; 1 2 1 2 ( ) n p t pt p t n y t Ce C e C e ? ? ? ? ? ? ? ? 1 2 0 n C C C y ? ? ? ? ? ? ? ? 1 1 2 2 0 n n pC p C p C Dy ? ? ? ? ? ? ? ? 1 1 1 1 1 1 2 2 0 n n n n n n p C p C p C D y ? ? ? ? ? ? ? ? ? ? ? ? 1 2 1 1 1 1 2 1 1 1 n n n n n p p p V p p p ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? A A A A ? ? 1 2 n C C C C ? ? ? ? 1 0 0 0 n C y Dy D y ? ? ? ? ? ? ? ? ?n=length(a)-1;Y0=[0,10]; p=roots(a);V=rot90(vander(p));c=V\Y0 ; dt=0.002;te=9; t=0:dt:te;y=zeros(1,length(t)); for k=1:n y=y+c(k)*exp(p(k)*t);end plot(t,y);grid xlabel( t ) ;ylabel( y ); title( 零輸入響應(yīng) ); 2.2連續(xù)時(shí)間系統(tǒng)零狀態(tài)響應(yīng)的數(shù)值計(jì)算我們知道，LTI連續(xù)系統(tǒng)可用如下所示的線性常系數(shù)微分方程來(lái) 描述， ( ) ( ) 0 0 ( ) ( ) N M i j i j i j a y t b f t ? ? ? ? ? 例如，對(duì)于以下方程: 3 2 1 0 3 2 1 0 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) a y t a y t a y t a y t b f t b f t b f t b f t ? ? ? ? ? ? ? 可用 輸入函數(shù) ，得出它的 3 2 1 0 3 2 1 0 [ , , , ], [ , , , ], a a a a a b b b b b ? ? ( ) u f t ? 沖激響應(yīng)h ，再根據(jù)LTI系統(tǒng)的零狀態(tài)響應(yīng)y(t)是激勵(lì)u(t) 與沖激響應(yīng)h(t)的卷積積分。 注意，如果微分方程的左端或右端表達(dá)式中有缺項(xiàng)，則其向量 a 或 b中的對(duì)應(yīng)元素應(yīng)為零，不能省略不寫(xiě)，否則出錯(cuò)。 求函數(shù)的零狀態(tài)響應(yīng)y″(t)+3y′(t)+6y(t)=6f′(t)-8f′(t) 及初始狀態(tài) 。輸入函數(shù) 。 (0 ) (0 ) 0 zs zs y y ? ? ? ? ) * 5 cos( ) * 2 sin( t t y ? ? 先求出系統(tǒng)的沖激響應(yīng)，寫(xiě)出其特征方程0 3 2 2 ? ? ? ? ? 求出其特征根為p和p，及相應(yīng)的留數(shù)r，r；則沖激響應(yīng)為1 2 1 2 ( ) pt p t h t re re ? ? 輸入y(t)可用輸入u(t)與沖激響應(yīng)h(t)的卷積求得。 MATLAB程序： a=[1,3,6]; n=length(a)-1; Y0=[0,10]; b=[1,6,8]; dt=0.001;te=9;

總結(jié)

以上是生活随笔為你收集整理的matlab求解零输入响应,基于matlab的零输入和零状态响应的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。