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