基于matlab的信號(hào)仿真分析的設(shè)計(jì)___課設(shè)說(shuō)明書(shū)(編輯修改稿)
2024-11-12 07:12 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 MATLAB 卷積積分法求解系統(tǒng)的零狀態(tài)響應(yīng)。 一般的連續(xù)時(shí)間系統(tǒng)分析有以下幾個(gè)步驟 : 唐山學(xué)院課程實(shí)踐 8 ① 求解系統(tǒng)的零輸入響應(yīng) 。 ② 求解系統(tǒng)的零狀態(tài)響應(yīng) 。 ③ 求解系統(tǒng)的全響應(yīng) 。 ④ 分析系統(tǒng)的卷積； ⑤ 畫(huà)出它們的圖形。 下面以具體的微分方程為例說(shuō)明利 用 MATLAB 軟件分析系統(tǒng)的具體方法 我們知道， LTI 連續(xù)系統(tǒng)可用如下所示的線性常系數(shù)微分方程來(lái)描述， 例如，對(duì)于以下方程 : 39。39。39。 39。39。 39。 39。39。39。 39。39。 39。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? ? ? ? ? ? ?可用 3 2 1 0 3 2 1 0[ , , , ] , [ , , , ] ,a a a a a b b b b b??輸入函數(shù) ()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) 及初始狀態(tài) 39。( 0 ) ( 0 ) 0zs zsyy????。輸入函數(shù) ( ) sin ( 3 * ) c o s( 2 * )f t t t??。 建模 先求出系統(tǒng)的沖擊響應(yīng)，寫(xiě)出其特征方程 2 5 4 0??? ? ? 求出其特征根為 p 和 p，及相應(yīng)的留數(shù) r， r；則沖擊響應(yīng)為 1212() p t p th t r e r e?? 輸入 y（ t）可用輸入 u（ t） 與沖擊響應(yīng) h（ t）的卷積求得。 卷積的計(jì)算 連續(xù)時(shí)間信號(hào) 1()ft和 2()ft的卷積運(yùn)算可用信號(hào)的分段求和來(lái)實(shí)現(xiàn)，即： 如果只求當(dāng) t =n（ n為整數(shù)）時(shí) f (t)的值 f (n ，則上式可得： 39。39。 39。 39。( ) 5 ( ) 4 ( ) 2 ( ) 4 ( )y t y t y t f t f t? ? ? ?1 2 1 2 1 20( ) ( ) * ( ) ( ) ( ) l im ( ) ( )kf t f t f t f t f t d f k f t k?????? ?? ? ? ?? ? ? ? ? ? ? ? ???( ) ( )00( ) ( )NMijijija y t b f t?????唐山學(xué)院課程實(shí)踐 9 式中的 實(shí)際上就是連續(xù)時(shí)間信號(hào) 1()ft和 2()ft經(jīng)等時(shí)間間隔 均勻抽樣的離 散序列 1()fk? 和 2()fk? 的卷積和。當(dāng) 足夠小時(shí)， ()fn?就是卷積積分的結(jié)果 ——連續(xù)時(shí)間信號(hào) f (t)的較好數(shù)值近似。 下面是利用 MATLAB 實(shí)現(xiàn)連續(xù)信號(hào)卷積的通用程序 conv()，該程序在計(jì)算出卷積積分的數(shù)值近似的同時(shí)，還繪制出 f (t)的時(shí)域波形圖。應(yīng)注意，程序中是如何設(shè)定 f (t)的時(shí)間長(zhǎng)度。 MATLAB程序： f1=input(39。輸入函數(shù) f1=39。)。 f2=input(39。輸入函數(shù) f2=39。)。 dt=input(39。dt=39。)。 y=conv(f1,f2)。 plot(dt*([1:length(y)]1),y)。grid on。 title(39。卷積 39。)。 xlabel(39。t39。)。 ylabel(39。f1*f239。) 程序運(yùn)行結(jié)果： 輸入以下數(shù)據(jù)： f1=sin(3*t) f2=cos(3*t+2) dt= 得出圖形如圖 3—5所示： 12( ) [( ) ]k f k f n k???? ? ? ? ??1 2 1 2( ) ( ) ( ) ( ) [ ( ) ]kkf n f k f t k f k f n k??? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ? ???0 2 4 6 8 10 12 3 0 0 2 0 0 1 0 00100200300400卷積tf1*f2圖 3— 5 唐山學(xué)院課程實(shí)踐 10 在 MATLAB 中，對(duì)于連續(xù) LTI 系統(tǒng)的沖激響 應(yīng)和階躍響應(yīng)的數(shù)值解，可分別用控制系統(tǒng)工具箱提供的函數(shù) impulse 和 step 來(lái)求解。其語(yǔ)句格式分別為： y = impulse(sys,t)。 y = step(sys,t)。 其中， t 表示計(jì)算機(jī)系統(tǒng)響應(yīng)的時(shí)間抽樣點(diǎn)向量， sys 表示 LTI 系統(tǒng)模型。 下面我們用幾個(gè)例子來(lái)驗(yàn)證以上知識(shí)點(diǎn) （ 1）求零輸入相應(yīng)：求齊次微分方程在給定初始條件下的零輸入相應(yīng)，y39。39。(t)+4y(t)=0,y(0_)=1,y39。(0_)=1。 解 在 MATLAB 命令窗口鍵入如下語(yǔ)句： eq=39。D2y+4*y=039。 con=39。y(0)=1,Dy(0)=139。 y=dsolve(eq,con)。y=simplify(y) y = 1/2*sin(2*t)+cos(2*t) ezplot(y),grid 零輸入響應(yīng)圖如圖 3—6 所示 圖 3—6 （ 2）求零狀態(tài)響應(yīng)：已知系統(tǒng) y39。39。(t)+4y(t)=x(t),輸入 x(t)為 cos(t)u(t),求系統(tǒng)的零狀態(tài)響應(yīng)。 解 在 MATLAB 命令窗口鍵入如下語(yǔ)句 eq=39。D2y+4*y=cos(t)*heaviside(t)39。 con=39。y()=0,Dy()=039。 y=dsolve(eq,con)。y=simplify(y) y = 唐山學(xué)院課程實(shí)踐 11 1/3*heaviside(t)*(cos(t)cos(2*t)) ezplot(y,[0,10]） 零狀態(tài)響應(yīng)圖如圖 3—7 所示 圖 3—7 （ 3 ） 求沖 擊響 應(yīng) 和階 躍 響應(yīng) ： 已知 某 LTI 系 統(tǒng) 的 微分 方程 為y’’(t)+2y’(t)+32y(t)=f’(t)+16f”(t),試用 MATLAB命令繪出 0≦ t≦ 4范圍內(nèi)系統(tǒng)的沖激響應(yīng) h（ t）和階躍響應(yīng) g（ t）。 解 在 MATLAB 命令窗口 鍵入如下語(yǔ)句 t=0::4。 sys=tf([1,16],[1,2,32])。 h=impulse(sys,t)。 g=step(sys,t)。 subplot(211) plot(t,h),grid on xlabel(39。Time(sec)39。),ylabel(39。h(t)39。) title(39。沖激響應(yīng) 39。) 沖激響應(yīng)圖如圖 3—8（ 1）所示 唐山學(xué)院課程實(shí)踐 12 圖 3—8（ 1） subplot(212) plot(t,g),grid on xlabel(39。Time(sec)39。),ylabel(39。g(t)39。) title(39。階躍響應(yīng) 39。) 階躍響應(yīng)圖如圖 3—8（ 2）所示 圖 3—8（ 2） 唐山學(xué)院課程實(shí)踐 13 （ 4）求卷積：求信號(hào) f1(t)=t 和 f2(t)=sin(t)的卷積的解
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號(hào)-1