武科大matlab仿真第七章系統(tǒng)時(shí)間響應(yīng)及其仿真(存儲(chǔ)版)
【正文】 12121,22,2,324213121???????????????????????????????????????????????????????????????????????????KKyhthfKKKyhthfKKyhthfKythfKmmmmmmmm)19(61211312113161 43211 KKKKyy mm ????????? ?????????? ?????第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 仿真算法 數(shù)值積分方法的選擇 在選擇積分方法時(shí)應(yīng)考慮以下幾個(gè)問(wèn)題 。 ? 固定步長(zhǎng)與變步長(zhǎng) ? 固定步長(zhǎng) ： 在整個(gè)仿真計(jì)算過(guò)程中 ， 積分步長(zhǎng) h始終不變 。 若系統(tǒng)的數(shù)學(xué)模型為高階微分方程 ， 則應(yīng)將高階微分方程轉(zhuǎn)化成一階微分方程組 。 % 21 yy ?? 12212 )1( yyyy ????第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) （ 3）編寫 MATLAB主程序 [T,Y]=ode45(39。Solution39。y239。這種調(diào)用格式不繪制仿真曲線圖。 其調(diào)用格式與 STEP函數(shù)相同 。 [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)曲線。,5,10)。\leftarrow輸出 39。 。 :06:1808:06:18February 4, 2023 ? 1意志 堅(jiān) 強(qiáng)的人能把世界放在手中像泥 塊 一 樣 任意揉捏。 :06:1808:06Feb234Feb23 ? 1越是無(wú)能的人，越喜 歡 挑剔 別 人的 錯(cuò)兒 。 2023年 2月 4日星期六 8時(shí) 6分 18秒 08:06:184 February 2023 ? 1一 個(gè) 人即使已登上 頂 峰，也仍要自強(qiáng)不息。 , February 4, 2023 ? 閱讀 一切好 書 如同和 過(guò) 去最杰出的人 談話 。 :06:1808:06Feb234Feb23 ? 1世 間 成事，不求其 絕對(duì)圓滿 ，留一 份 不足，可得無(wú)限完美。 :06:1808:06:18February 4, 2023 ? 1他 鄉(xiāng) 生白 發(fā) ， 舊國(guó)見青 山。輸入 \rightarrow39。 [u,t]=gensig(39。 Tf為信號(hào)的時(shí)間區(qū)間 。gx39。 第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 【 例 4】 已知系統(tǒng)模型 ，求其單位階躍響應(yīng)。y139。b:39。 % 生成 2行 1列的零陣 dy(1)= y(2)。 第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 系統(tǒng)仿真的 MATLAB函數(shù) 數(shù)值積分方法的 MATLAB函數(shù) ? 【 說(shuō)明 2】 ? 對(duì)于剛性微分方程 （ 特征值數(shù)值相差較大 ） ， 可用 ode15s， 其調(diào)用格式與 ode45相同 。 截?cái)嗾`差 舍入誤差 總誤差 e h 誤差與積分步長(zhǎng) 顯然選擇一個(gè)合適的積分步長(zhǎng)可使總誤差達(dá)到最小 。 一般與系統(tǒng)中具有最大時(shí)間常數(shù) Tmax的環(huán)節(jié)有關(guān) ， 要求計(jì)算步長(zhǎng) h取得很大 。 Ki的個(gè)數(shù)與 yk+1泰勒展開式所取的項(xiàng)數(shù)有關(guān)(即 RK算法的階數(shù) )， 同時(shí)還與計(jì)算區(qū)間內(nèi)所取導(dǎo)數(shù)值的點(diǎn)數(shù)有關(guān) 。 3)3(!312)2(!211 ??????? kkkkk hyhyhyyy第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 仿真算法 數(shù)值積分的幾個(gè)算法 ? 龍格 庫(kù)塔法 ? 二階龍格 庫(kù)塔 （ RK） 法 取 (3)式的前三項(xiàng) ， 則有 )3(39。 3)3(!312)2(!211 ??????? kkkkk hyhyhyyy第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 仿真算法 數(shù)值積分的幾個(gè)算法 ? 梯度法 梯度法是歐拉法的改進(jìn) 。 ),2,1,0()3(39。 對(duì)連續(xù)時(shí)間系統(tǒng)進(jìn)行數(shù)字動(dòng)態(tài)仿真 ， 主要是兩種方法： ① 基于數(shù)值積分的仿真方法 。 第七章 系 統(tǒng)時(shí)間響應(yīng) 及其仿 真 仿真技術(shù) 第七章 系統(tǒng)時(shí)間響應(yīng)及其仿真 仿真算法 ? 對(duì)系統(tǒng)的時(shí)間響應(yīng)進(jìn)行動(dòng)態(tài)仿真 ， 采用什么樣的仿真算法是一個(gè)至關(guān)重要的問(wèn)題 。 積分值與實(shí)際微分方程解的誤差取決于步長(zhǎng) h和計(jì)算所用的階數(shù) ， 它是數(shù)值積分的基礎(chǔ) 。1 mmm hyyy ???)5(,(),(),(1111200001 00????????????????? hytfyyhytfyyhytfyhyyymmmmty?)3(39。 )3(39。 若設(shè) a1=a2，則 ??? ?? ?? 1 2/12121bbaa即二階 RK法公式為 )12()(2),(),(211121?????
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1