武科大matlab仿真第七章系統(tǒng)時(shí)間響應(yīng)及其仿真(文件)
【正文】 總誤差 e h 誤差與積分步長(zhǎng) 顯然選擇一個(gè)合適的積分步長(zhǎng)可使總誤差達(dá)到最小 。 ? 此外 ， h還應(yīng)與模型的信號(hào)響應(yīng)情況有關(guān) ， 例如在穩(wěn)態(tài)時(shí) ， 可取較大的步長(zhǎng) ， 見(jiàn)上圖 。 第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 系統(tǒng)仿真的 MATLAB函數(shù) 數(shù)值積分方法的 MATLAB函數(shù) ? 【 說(shuō)明 2】 ? 對(duì)于剛性微分方程 （ 特征值數(shù)值相差較大 ） ， 可用 ode15s， 其調(diào)用格式與 ode45相同 。 第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 【 例 2】 已知二階 微分方程 求時(shí)間區(qū)間 t=[0， 20]微分方程的解。 % 生成 2行 1列的零陣 dy(1)= y(2)。,[0 20],[0,1])。b:39。time s39。y139。 其中 y1(紅線 )為微分方程的解 。 第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 【 例 4】 已知系統(tǒng)模型 ，求其單位階躍響應(yīng)。r39。gx39。 ? 信號(hào)發(fā)生器函數(shù) GENSIG ? GENSIG可為系統(tǒng)時(shí)間響應(yīng)產(chǎn)生周期輸入信號(hào)，其調(diào)用格式為： 其中 ,Type為信號(hào)類型： ‘sin’— 正弦波 ‘square’— 方波 39。 Tf為信號(hào)的時(shí)間區(qū)間 。Ys為響應(yīng)值向量 。 [u,t]=gensig(39。r39。輸入 \rightarrow39。 , February 4, 2023 ? 雨中 黃 葉 樹 ，燈下白 頭 人。 :06:1808:06:18February 4, 2023 ? 1他 鄉(xiāng) 生白 發(fā) ， 舊國(guó)見(jiàn)青 山。 2023年 2月 4日星期六 8時(shí) 6分 18秒 08:06:184 February 2023 ? 1做前，能 夠環(huán)視 四周；做 時(shí) ， 你 只能或者最好沿著以腳 為 起點(diǎn)的射 線 向前。 :06:1808:06Feb234Feb23 ? 1世 間 成事，不求其 絕對(duì)圓滿 ，留一 份 不足，可得無(wú)限完美。 。 , February 4, 2023 ? 閱讀 一切好 書 如同和 過(guò) 去最杰出的人 談話 。 勝 人者有力，自 勝 者強(qiáng)。 2023年 2月 4日星期六 8時(shí) 6分 18秒 08:06:184 February 2023 ? 1一 個(gè) 人即使已登上 頂 峰，也仍要自強(qiáng)不息。 2023年 2月 4日星期六 上午 8時(shí) 6分 18秒 08:06: ? 1最具挑 戰(zhàn) 性的挑 戰(zhàn) 莫 過(guò) 于提升自我。 :06:1808:06Feb234Feb23 ? 1越是無(wú)能的人，越喜 歡 挑剔 別 人的 錯(cuò)兒 。 2023年 2月 4日星期六 8時(shí) 6分 18秒 08:06:184 February 2023 ? 1空山新雨后，天 氣晚來(lái) 秋。 :06:1808:06:18February 4, 2023 ? 1意志 堅(jiān) 強(qiáng)的人能把世界放在手中像泥 塊 一 樣 任意揉捏。 , February 4, 2023 ? 很多事情努力了未必有 結(jié) 果，但是不努力卻什 么 改 變 也 沒(méi) 有。 。 :06:1808:06Feb234Feb23 ? 1故人江海 別 ，幾度隔山川。\leftarrow輸出 39。.39。,5,10)。 lsim(sys,U,T) %基本調(diào)用格式 lsim(sys1,sys2,…,U,T) %繪制多個(gè)系統(tǒng)對(duì)同一個(gè)任意輸入響應(yīng)曲線 [Ys,Ts]=lsim(sys,U,T) %該格式不繪制響應(yīng)曲線 第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 【 例 7】 已知系統(tǒng)模型 ， 計(jì)算系統(tǒng)在周期為 5s的方波信號(hào)作用下的響應(yīng) 。 [U,T]=gensig(Type,Tau) [U,T]=gensig(Type,Tau,Tf,Ts) 第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 時(shí)間響應(yīng)仿真的 MATLAB函數(shù) ? 信號(hào)發(fā)生器和任意輸入響應(yīng)函數(shù) ? 任意輸入響應(yīng)函數(shù) LSIM ? LSIM用來(lái)仿真系統(tǒng)對(duì)任意輸入的時(shí)間響應(yīng)，并繪制響應(yīng)曲線。— 周期脈沖波 Tau為信號(hào)周期； U為信號(hào)值向量 。 其調(diào)用格式與 STEP函數(shù)相同 。y39。這種調(diào)用格式不繪制仿真曲線圖。 ? 階躍響應(yīng)仿真函數(shù) （ STEP） ? 基本調(diào)用格式 ? 對(duì)于 LTI連續(xù) (或離散 )時(shí)間系統(tǒng)，以下調(diào)用格式可用于繪制系統(tǒng)單位階躍響應(yīng)曲線。y239。Position Y39。Solution39。r39。 % 21 yy ?? 12212 )1( yyyy ????第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) （ 3）編寫 MATLAB主程序 [T,Y]=ode45(39。 dy為 2*1數(shù)組 ， 其維數(shù)等于微分方程的階數(shù) 。 若系統(tǒng)的數(shù)學(xué)模型為高階微分方程 ， 則應(yīng)將高階微分方程轉(zhuǎn)化成一階微分方程組 。 其目的是在保證一定計(jì)算精度的前提下 ， 盡可能選擇較大步長(zhǎng) 。 ? 固定步長(zhǎng)與變步長(zhǎng) ? 固定步長(zhǎng) ： 在整個(gè)仿真計(jì)算過(guò)程中 ， 積分步長(zhǎng) h始終不變 。 ? 局部截?cái)嗾`差 ： 由積分方法和階次的限制而引起的誤差 。 第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 仿真算
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1