北太真元多域動態(tài)系統(tǒng)建模仿真平臺基于北太天元數(shù)值計算通用軟件研發(fā)，是建立在北太天元之上的一個重量級工具箱，同時它也可以與北太天元其他工具箱聯(lián)動。北太真元主要是解決復(fù)雜系統(tǒng)的圖形化、模塊化建模、仿真與驗證，可應(yīng)用于汽車、航空航天、船舶、電子電力等行業(yè)。
剛剛開始接觸北太真元的仿真，準備先跑一下example試一試水。結(jié)果跑標準案例的時候出現(xiàn)了報錯：%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%錯誤: 文件Comput_FixedorMaxStep第82行的表達式: 錯誤: 文件Comput_FixedorMaxStep第82行的表達式: max是腳本，調(diào)用時不能帶有返回值或者輸入?yún)?shù)由下面文件調(diào)用:1: Comput_FixedorMaxStep的第1行;2: 道路路面不平順模擬仿真_Baltamulink2023ForAll的第20行;3: 道路路面不平順模擬仿真_Baltamulink2023ForAll的第19行;函數(shù)執(zhí)行中顯示有錯誤信息，請反饋給開發(fā)團隊。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%這是啥情況？我用simulink的時候也妹有出現(xiàn)這樣的情況？有木有大佬幫忙看一看？
嘗試復(fù)現(xiàn)基于十字路口的交通流仿真，顯示drivingScenario沒有被賦值。請問，目前是否不支持該功能？
已知完全可控開環(huán)系統(tǒng)狀態(tài)方程，其中：A = [-6 -5 -2；1 0 0;0 1 0]，B = [1; 0; 0]，C = [0 0 1]，D = 0增益矩陣K = [8 43 78]；增益L = 80。在Baltamulink中構(gòu)建如下圖所示的仿真模型：設(shè)置仿真參數(shù)：階躍信號模塊：階躍時間為0；Sum1模塊：符號為 ++++Sum2模塊：符號為 +-Sum3模塊：符號為 ++++Sum4模塊：符號為 +++Gain1模塊 = 80Gain2模塊 = -6Gain3模塊 = -5Gain4模塊 = -2Gain5模塊 = -6Gain6模塊 = -5Gain7模塊 = -2Gain8模塊 = 8Gain9模塊 = 43Gain10模塊 = 78仿真時長：12；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：藍色為原閉環(huán)系統(tǒng)的階躍響應(yīng)曲線；紅色為狀態(tài)反饋控制后的階躍響應(yīng)曲線。
串聯(lián)超前校正是在系統(tǒng)中串聯(lián)一個校正環(huán)節(jié)形成閉環(huán)系統(tǒng)。當時為串聯(lián)超前校正。T1和T2的值是根據(jù)控制指標及被控系統(tǒng)的相位裕度得出的。已知被控對象開環(huán)傳遞函數(shù)為，使用超前校正環(huán)節(jié)進行仿真，比較校正前后系統(tǒng)的動態(tài)特性參數(shù)。在Baltamulink中構(gòu)建如下圖所示的仿真模型：設(shè)置仿真參數(shù)：階躍信號模塊：階躍時間為1；Sum模塊：符號為 +-仿真時長：7；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：藍色為校正前；紅色為校正后。從仿真結(jié)果中可以看出，校正前系統(tǒng)超調(diào)量達到80%，在穩(wěn)態(tài)誤差為2%的情況下，穩(wěn)態(tài)時間為4.9秒左右；校正后系統(tǒng)超調(diào)量為27%，穩(wěn)態(tài)時間為0.1s。由此可以看出，在系統(tǒng)出現(xiàn)超調(diào)前完成了校正，校正速度較快，穩(wěn)態(tài)時間減少為1/50，超調(diào)量明顯降低。
利用狀態(tài)方程模塊建模：若令：，那么微分方程： 可寫成：寫成狀態(tài)方程為： 式中， 在Baltamulink中構(gòu)建求解微分方程的模型并仿真，根據(jù)系統(tǒng)狀態(tài)方程構(gòu)建如下圖所示的仿真模型：模型中各個模塊說明如下。（1） u(t) 模塊：設(shè)置階躍時間為 0。（2） stateSpace模塊：A、B、C、D 系數(shù)依次為 [0,1;-0.4,-0.2]、[0;0.2]、[1,0] 和 0。仿真時長：20s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
利用傳遞函數(shù)模塊建模：對方程： 兩邊取拉普拉斯變換，得：經(jīng)整理得傳遞函數(shù)：在Baltamulink中構(gòu)建求解微分方程的模型并仿真，根據(jù)系統(tǒng)傳遞函數(shù)構(gòu)建如下圖所示的仿真模型：模型中各個模塊說明如下。（1） u(t) 模塊：設(shè)置階躍時間為 0。（2） transferFunc 模塊：分子多項式系數(shù) [0.2]；分母多項式的系數(shù) [1,0.2,0.4]。仿真時長：20s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
控制系統(tǒng)穩(wěn)定與否是絕對穩(wěn)定性的概念。而對一個穩(wěn)定的系統(tǒng)而言，還有一個穩(wěn)定的程度，即相對穩(wěn)定性的概念。相對穩(wěn)定性與系統(tǒng)的動態(tài)性能指標有著密切的關(guān)系。在設(shè)計一個控制系統(tǒng)時，不僅要求它必須是絕對穩(wěn)定的，而且還應(yīng)保證系統(tǒng)具有一定的穩(wěn)定程度。只有這樣，才能不致因系統(tǒng)參數(shù)變化而導(dǎo)致系統(tǒng)性能變差甚至不穩(wěn)定。對于一個最小相角系統(tǒng)而言，曲線越靠近點，系統(tǒng)階躍響應(yīng)的振蕩就越強烈，系統(tǒng)的相對穩(wěn)定性就越差。因此，可用曲線對點的接近程度來表示系統(tǒng)的相對穩(wěn)定性。通常，這種接近程度是以相角裕度和幅值裕度來表示的。要計算相角裕度，首先要知道截止頻率。求較方便的方法是先由繪制曲線，由與線的交點確定。而求幅值裕度首先要知道相角交界頻率，對于階數(shù)不太高的系統(tǒng)，直接解三角方程是求較方便的方法。通常是將寫成虛部和實部， 令虛部為零而解得?！纠?-10】 某單位反饋系統(tǒng)的開環(huán)傳遞函數(shù)為：試求時系統(tǒng)的相角裕度和幅值裕度。將該開環(huán)傳遞函數(shù)變換為：在實際工程設(shè)計中，只要繪出曲線即可。根據(jù)該傳遞函數(shù)模型，在北太真元搭建穩(wěn)定裕度系統(tǒng)模型如下圖所示： 設(shè)置仿真參數(shù)：從上到下，傳遞函數(shù)參數(shù)依次為：num = [52]；den = [1 0]；num = [1]；den = [1 1]；num = [1]；den = [1 5]；仿真時長：10s；步長0.1s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
已知某系統(tǒng)的開環(huán)對數(shù)頻率特性如下圖所示，試確定其開環(huán)傳遞函數(shù)。 根據(jù)對數(shù)幅頻特性曲線，可以寫出開環(huán)傳遞函數(shù)的表達形式如下： 根據(jù)對數(shù)頻率特性的坐標特點有，可以確定開環(huán)增益。根據(jù)相頻特性的變化趨勢（-270°-> -90°），可以判定系統(tǒng)為非最小相角系統(tǒng)。G(s)中一階復(fù)合微分環(huán)節(jié)和慣性環(huán)節(jié)至少有一個是“非最小相角”的，將系統(tǒng)可能的開環(huán)零點極點分布畫出來，如下表所示： 分析相角的變化趨勢，可見，只有當慣性環(huán)節(jié)極點在右半s平面，一階復(fù)合微分環(huán)節(jié)零點在左半s平面是，相角才符合從-270°到-90°的變化規(guī)律。因此可以確定系統(tǒng)的開環(huán)傳遞函數(shù)為：對于最小相角系統(tǒng)，對數(shù)幅頻特性與對數(shù)相頻特性之間存在唯一確定的對應(yīng)關(guān)系，根據(jù)對數(shù)幅頻特性就完全可以確定相應(yīng)的對數(shù)相頻特性和傳遞函數(shù)，反之亦然。由于對數(shù)幅頻特性容易繪制，所以在分析最小相角系統(tǒng)時，通常只畫其對數(shù)幅頻特性，對數(shù)相頻特性則只需概略畫出，或者不畫。根據(jù)該傳遞函數(shù)模型，在北太真元搭建最小相角系統(tǒng)模型如下圖所示： 設(shè)置仿真參數(shù)：從上到下，傳遞函數(shù)參數(shù)依次為：num = [1]；den = [1 0]；num = [1 1]；den = [1 -1]；仿真時長：10s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
已知開環(huán)系統(tǒng)傳遞函數(shù)：首先將G(s)化為唯一標準形式：此系統(tǒng)由比例環(huán)節(jié)、積分環(huán)節(jié)、慣性環(huán)節(jié)、一階微分環(huán)節(jié)和振蕩環(huán)節(jié)共5個環(huán)節(jié)組成。慣性環(huán)節(jié)轉(zhuǎn)折頻率：一階復(fù)合微分環(huán)節(jié)轉(zhuǎn)折頻率：振蕩環(huán)節(jié)轉(zhuǎn)折頻率：開環(huán)增益：K=4，積分環(huán)節(jié)數(shù)v=1，低頻起始段由K/s=4/s決定。根據(jù)該傳遞函數(shù)模型，在北太真元建立模型如下圖所示： 設(shè)置仿真參數(shù)：從上到下，傳遞函數(shù)參數(shù)依次為：num = [1]；den = [1 0]；num = [1]；den = [1 0.5]；num = [64 128]；den = [1 3.2 64]；仿真時長：10s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
針對標準二階系統(tǒng)傳遞函數(shù)：改變ζ（阻尼比）和ωn（自由振蕩頻率）的參數(shù)設(shè)置，觀察對系統(tǒng)輸出的影響。在二階系統(tǒng)自由振蕩頻率ωn不變的情況下，改變阻尼系數(shù)ζ為無阻尼（ζ= 0）、欠阻尼（0<ζ< 1）、臨界阻尼（ζ= 1）和過阻尼（ζ> 1）的4中狀態(tài)，分別取ζ= 0，ζ= 0.5，ζ= 1，ζ= 2帶入二階系統(tǒng)傳遞函數(shù) 中，搭建4個不同的仿真模型，輸出結(jié)果，觀察仿真結(jié)果得出結(jié)論。根據(jù)該傳遞函數(shù)模型，在北太真元建立模型如下圖所示： 設(shè)置仿真參數(shù)：從上到下，傳遞函數(shù)參數(shù)依次為：num = [4]；den = [1 0 4]；num = [4]；den = [1 2 4]；num = [4]；den = [1 4 4]；num = [4]；den = [1 8 4]；階躍信號模塊的階躍時間為0；仿真時長：6s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示： 結(jié)論：從仿真結(jié)果可以看出，改變阻尼比，系統(tǒng)的超調(diào)量也在變化，系統(tǒng)達到穩(wěn)態(tài)的時間也發(fā)生變化。當ζ= 0無阻尼時，出現(xiàn)等幅振蕩曲線，超調(diào)量為100%，穩(wěn)態(tài)時間為無窮大；當ζ< 1時，信號曲線衰減振蕩，有超調(diào)量；當ζ≥ 1時，沒有超調(diào)量，隨著ζ增大，達到穩(wěn)態(tài)的時間也增大。
零極點增益模型實際上是傳遞函數(shù)的另一種表現(xiàn)形式，其原理是分別對源系統(tǒng)傳遞函數(shù)的分子、分母進行因式分解處理，以獲得系統(tǒng)零點和幾點的表示形式。 式中，k為系統(tǒng)增益；zi(i = 1,2,3,...,m)為分子多項式的根，稱為系統(tǒng)的零點；pj(j = 1,2,...,n)是分母多項式的根，稱為系統(tǒng)的極點。傳遞函數(shù)的分母多項式就是它的特征多項式，它等于零的方程就是傳遞函數(shù)的特征方程，特征方程的根也就是傳遞函數(shù)的極點。傳遞函數(shù)的極點決定了所描述系統(tǒng)的自由運動狀態(tài)；零點影響系統(tǒng)各模態(tài)在系統(tǒng)響應(yīng)中的比重。零點增益模型的命令格式如下：ZPG = zpk(z, p, k)其中ZPG是建立的零極點增益模型；z、p、k分別是系統(tǒng)的零點向量、極點向量和增益。例：利用Baltamulink建立系統(tǒng)G(s) = 18(s + 2) / (s + 0.4)(s + 15)(s + 25)的零點增益模型，進行系統(tǒng)仿真。將上面零點增益模型進行轉(zhuǎn)換得到模型的傳遞函數(shù)如下：G(s) = 18(s + 2)/(s^3 + 40.4s^2 + 391s + 150)根據(jù)該傳遞函數(shù)模型，在北太真元建立模型如下圖所示： 設(shè)置仿真參數(shù)：仿真時長：20s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
力——質(zhì)量系統(tǒng)，要拉動一個箱子（拉力f=1N），箱子質(zhì)量為M（1kg），箱子與地面的摩擦力為[(b=0.4N.m/s)]，其大小與車子的速度成正比。如下圖所示： 其運動方程式為：F - bx’ = Mx’’拉力作用時間為2s。在北太真元建立模型如下圖所示： 設(shè)置仿真參數(shù)：設(shè)置stepInputOne模塊的階躍時間為0，表示摩擦力作用時間；設(shè)置stepInputTwo模塊的階躍時間為2，表示拉力作用時間；設(shè)置gainTwo模塊的增益值為0.4表示摩擦力；仿真時長：10s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示： 因為有摩擦力存在，箱子最終會停止前進。
求解二階微分方程：x’’(t) + 0.4x’(t) + 0.9x(t) = 0.7u(t) 的解，其中u(t)是脈沖信號。 在北太真元建立模型如下圖所示： 設(shè)置仿真參數(shù)：脈沖信號u(t)為方波信號模塊，參數(shù)振幅 = 1；周期 = 2.5；脈沖寬度 = 50；相位 = 0；仿真時長：10s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
已知二自由度質(zhì)量-彈簧-阻尼串聯(lián)系統(tǒng)模型如下圖所示： 這個系統(tǒng)由兩個質(zhì)量塊（小車）和三組彈簧阻尼器組成，假設(shè)地面是光滑的，這樣系統(tǒng)中沒有摩擦作用。 u1(t) 和 u2(t) 分別是兩個質(zhì)量塊所受的外力， x1(t) 和 x2(t) 分別是兩個質(zhì)量塊的位移。 m1,2、k1,2,3 和 b1,2,3 分別對應(yīng)圖中的質(zhì)量、彈簧剛度和阻尼系數(shù)。是對于整個系統(tǒng)而言，輸入兩個外力，輸出兩個位移，因此這是一個多輸入多輸出系統(tǒng)。兩個小車都只能沿橫向左右運動，因此為二自由度系統(tǒng)。根據(jù)牛頓第二定律，物體所受合力等于物體的慣性力，而慣性力是物體質(zhì)量與加速度的乘積?？梢岳斫鉃槲矬w受力后產(chǎn)生加速度，而有加速度存在就會產(chǎn)生運動趨勢，造成物體運動。對兩個質(zhì)量塊分別進行受力分析： 首先需定義系統(tǒng)的狀態(tài)變量，在這種情況下一般用物體的位移和速度作為狀態(tài)變量。設(shè)狀態(tài)向量 Z分別對應(yīng)質(zhì)量塊1的位移和速度、質(zhì)量塊2的質(zhì)量和速度: 根據(jù)上面對模型的數(shù)學(xué)推導(dǎo)，當將各狀態(tài)向量 z 對時間求導(dǎo)（微分），可以將系統(tǒng)整理為各狀態(tài)量的一階微分方程組： 因為這是一個線性系統(tǒng)，因此系統(tǒng)的狀態(tài)可以表示為矩陣形式： 式中，A 為系統(tǒng)狀態(tài)矩陣；B 為輸入矩陣；u 為輸入（控制）向量： 有了系統(tǒng)的狀態(tài)空間方程，接下來考慮系統(tǒng)的輸出，也就是我們希望得到的量。我們假設(shè)對兩個質(zhì)量塊的位移感興趣，即 z1 和 z3 ，那么就將它兩個狀態(tài)作為系統(tǒng)的輸出，則輸出方程的矩陣形式為：y = Cz;式中，y 為輸出向量；C 為輸出矩陣： 以上，使用狀態(tài)空間模型對系統(tǒng)完成描述。令，質(zhì)量：m = 1kg；彈簧剛度： k = 1N/m；阻尼系數(shù)：b = 1N.s/m。則，狀態(tài)空間方程系數(shù)如下：A = [0 1 0 0;-2 -2 1 1;0 0 0 1;1 1 -2 -2];B = [0 0;1 0;0 0;0 1];C = [1 0 0 0;0 0 1 0];D = [0 0;0 0]; 在北太真元建立模型如下圖所示： 設(shè)置仿真參數(shù)：兩個質(zhì)量塊所受的外力u1是常量模塊 = 2 N； u2是階躍型號=1N。仿真時長：10s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
已知質(zhì)量塊質(zhì)量 m = 1kg，阻尼 b = 3 N.s/m，彈簧系數(shù) k = 90 N/m，且物塊的初始位移 x(0) = 0.04m，其初始速度為x’(0) = 0.01 m/s。創(chuàng)建該系統(tǒng)的北太真元模型，并運行仿真。彈簧-質(zhì)量-阻尼系統(tǒng)如下圖所示： 建立理論數(shù)學(xué)模型。對于無外力的系統(tǒng)，根據(jù)牛頓定理可以得到：mx’’ + bx’ + kx = 0代入數(shù)值并整理得：x’’ = -3x’ - 90x在北太真元建立模型如下圖所示： 設(shè)置仿真參數(shù)：仿真時長：10s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
1 線性插值線性插值是指插值函數(shù)為一次多項式的插值方式，其在插值節(jié)點上的插值誤差為零。線性插值相比其他插值方式，如拋物線插值，具有簡單、方便的特點。線性插值的幾何意義即為概述圖中利用過A點和B點的直線來近似表示原函數(shù)。線性插值可以用來近似代替原函數(shù)，也可以用來計算得到查表過程中表中沒有的數(shù)值。1.2 基礎(chǔ)知識已知函數(shù)在區(qū)間上個互異點上的函數(shù)值，若存在一簡單函數(shù)，使 并要求誤差 的絕對值在整個區(qū)間上比較小。這樣的問題稱為插值問題。其中：：插值節(jié)點：被插值函數(shù)：插值函數(shù)：插值區(qū)間如果在插值區(qū)間內(nèi)部用代替則稱為內(nèi)插；在插值區(qū)間以外，用代替則稱為外插。1.3 簡介線性插值是一種較為簡單的插值方法，其插值函數(shù)為一次多項式。線性插值，在各插值節(jié)點上插值的誤差為0。設(shè)函數(shù)在兩點，上的值分別為，，求多項式 使?jié)M足 由解析幾何可知 稱為在處的一階均差，記以。于是，得 如果按照整理，則 以上插值多項式為一次多項式，這種插值稱為線性插值。1.4 幾何意義線性插值的幾何意義如圖1所示，即為利用過點和的直線來近似原函數(shù)。 1.5 應(yīng)用1）線性插值在一定允許誤差下，可以近似代替原來函數(shù)；2）在查詢各種數(shù)值表時，可通過線性插值來得到表中沒有的數(shù)值。2一維線性插值仿真實例首先，在Baltamulink中，添加輸出正弦波模塊、一維插值模塊、信號合并模塊、輸出模塊，建立一維插值仿真模型，每個模塊參數(shù)都設(shè)置為默認值；模型如下圖所示； 設(shè)置仿真參數(shù)：仿真時長：10s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示： 綠色代表原正弦波數(shù)據(jù)；橙色表示一維插值后的正弦波數(shù)據(jù)。
基于Baltamulink的衰減曲線法整定參數(shù)使用4:1衰減曲線法設(shè)計下列被控傳遞函數(shù)的PI控制器，分別計算P控制、PI控制的參數(shù)值，并繪制控制前后系統(tǒng)的單位階躍響應(yīng)曲線。4:1衰減法控制參數(shù)計算公式如下表所示：4:1衰減法控制被控傳遞函數(shù)方程如下：Gp(s) =1 / 100^3 + 80^s + 17s + -1;調(diào)節(jié)參數(shù)時，比例系數(shù)由小變大，并增加擾動觀察響應(yīng)過程，知道響應(yīng)曲線峰值衰減比為4:1，記錄此時的比例系數(shù)Kp為Ks，兩個峰值之間的時間周期為周期Ts。假設(shè)：響應(yīng)曲線峰值衰減比為4:1時的比例系數(shù) Ks= 4.74; 兩個峰值之間的時間周期 Ts = 21.9967；則，按照上面4:1衰減法控制表中的計算公式可得：P控制：比例系數(shù) Kp = Ks = 4.74；PI控制：比例系數(shù) Kp = Ks/1.2 = 3.95； 積分時間常數(shù) Ti = 0.5 * Ts = 10.9984; 積分系數(shù) Ki = Kp / Ti = 3.95 / 10.9984 = 0.3591;P（比例）控制首先，把控制器設(shè)置成純比例控制，即令積分系數(shù)Ki和微分系數(shù)Kd為零，在北太真元建立模型，形成比例控制系統(tǒng)，結(jié)構(gòu)如下圖所示； 設(shè)置參數(shù)：仿真時長：30s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示： PI（比例積分）控制首先，在純比例控制系統(tǒng)的基礎(chǔ)上增加積分系數(shù)Ki，令微分系數(shù)Kd為零，在北太真元建立模型，形成比例控制系統(tǒng)，結(jié)構(gòu)如下圖所示； 設(shè)置參數(shù)：仿真時長：30s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示：
基于Baltamulink狀態(tài)空間模型的汽車時域特性仿真問題：利用汽車橫擺角速度傳遞函數(shù)和質(zhì)心側(cè)偏角傳遞函數(shù)，對汽車時域響應(yīng)進行仿真，繪制汽車橫擺角速度和質(zhì)心偏側(cè)角的時域特性曲線。汽車時域響應(yīng)仿真所需參數(shù)見下表。 取狀態(tài)向量為X = [β ωr]’，輸入向量U = [δ1]，輸出向量為Y = [β ωr]’，狀態(tài)空間方程為： 式中，A = [(K1+K2)/mu, (aK1 - BK2)/mu2-1; (aK1-bK2)/Iz, (a2K1+b2K2)/Iz*u] 稱為系統(tǒng)矩陣；B = [-K1/mu; -aK1/Iz] 稱為控制矩陣；C = [1 0; 0 1] 稱為輸出矩陣；D = [0; 0] 稱為傳遞矩陣。汽車速度分別選取20m/s、30m/s、40m/s；在仿真時間0s時給前輪一個階躍信號，使前輪轉(zhuǎn)角從0°轉(zhuǎn)到10°，并保持不變。根據(jù)汽車狀態(tài)空間模型，建立模型，繪制不同車速下的汽車橫擺角速度和質(zhì)心側(cè)偏角的時域特性曲線。 首先：通過北太天元計算汽車狀態(tài)空間方程的系統(tǒng)矩陣和控制矩陣，在北太天元依次輸入下面語句；>> m=2050;Iz=5600;a=1.5;b=1.8;L=3.3;>> k1=-38900;k2=-39200;>> u= [20 30 40];>> a11 = (k1 + k2)/m./u;>> a12 = (a*k1 - b*k2 -m.*u.^2)/m./u.^2;>> a21 = (a*k1 - b*k2)/Iz;>> a22=(a^2*k1 + b^2*k2)/Iz./u;>> b11 = -k1/m./u;>> b21 = -a*k1/Iz;得到結(jié)果如下圖1所示 ；圖1將命令行窗口，和工作區(qū)窗口放大后如圖2、圖3所示；圖2 圖3 因為，汽車狀態(tài)空間方程的系統(tǒng)矩陣為：A = [a11, a12; a21, a22]；控制矩陣為：B = [ b11; b21]； 所以，從圖3紅色框中可以得到各項系數(shù)如下：當汽車速度 s = 20 m/s 時，系統(tǒng)矩陣：A = [-1.90488,-0.98511;2.18036,-1.91547]; 控制矩陣：B = [0.94878; 10.4196];當汽車速度 s = 30 m/s 時，系統(tǒng)矩陣：A = [-1.26992,-0.993382;2.180436,-1.27698]; 控制矩陣：B = [0.63252; 10.4196];當汽車速度 s = 40 m/s 時，系統(tǒng)矩陣：A = [-0.952439,-0.996277;2.180436,-0.957737]; 控制矩陣：B = [0.47439; 10.4196];狀態(tài)方程輸出矩陣C = [1 0; 0 1]；傳遞矩陣D = [0; 0]。又因為，在仿真時間0s時給前輪一個階躍信號，使前輪轉(zhuǎn)角從0°轉(zhuǎn)到10°；所以模型還需一個階躍信號模塊，階躍時間=0；且，還需一個增益模塊，增益= pi*10/180 = 0.1745。通過北太真元建立汽車狀態(tài)空間模型，如下圖所示： 設(shè)置參數(shù)：仿真時長：10s；步長0.01s；求解器：ode4得到的仿真結(jié)果，如下圖所示： 上半部分代表汽車橫擺角速度時域特性曲線；即： 墨綠色代表速度20m/s時的特性曲線；綠色代表速度30m/s時的特性曲線；紅色代表速度40m/s時的特性曲線。 下半部分代表汽車質(zhì)心側(cè)偏角時域特性曲線；即： 紫色代表速度20m/s時的特性曲線；橙色代表速度30m/s時的特性曲線；藍色代表速度40m/s時的特性曲線。
基于Baltamulink傳遞函數(shù)分析阻尼系數(shù)和固有頻率對性能的影響實例：已知傳遞函數(shù)G(s) = ω2/ s2+ 2ζωs + ω2，分析阻尼系數(shù)和固有頻率對性能的影響。（1）假設(shè) ω = 1, ζ = 0, 0.8, 1.5;（2）假設(shè) ζ = 1, ω = 1, 2, 3; 從（1）可得，阻尼系數(shù)傳遞函數(shù)的系數(shù)可以是：G(s1) = [1; 1 0 1];G(s2) = [1; 1 1.6 1];G(s3) = [1; 1 3 1]; 通過北太真元建立“阻尼系數(shù)對系統(tǒng)性能的影響”模型，如下圖所示： 設(shè)置參數(shù)：仿真時長：20s；步長0.1s；求解器：ode4得到的仿真結(jié)果，如下圖所示： 墨綠色代表阻尼系數(shù)=0時的特性曲線；淡綠色代表阻尼系數(shù)=0.8時的特性曲線；紅色代表阻尼系數(shù)=1.5時的特性曲線。從結(jié)果圖中可以看出，阻尼系數(shù)決定了系統(tǒng)的振蕩幅度，阻尼系數(shù)越小，振蕩幅度越大。從（2）可得，固有頻率傳遞函數(shù)的系數(shù)可以是：G(s1) = [1; 1 0.5 1];G(s2) = [4; 1 1 4];G(s3) = [9; 1 1.5 9];通過北太真元建立“固有頻率對系統(tǒng)性能的影響”模型，同上模型；設(shè)置參數(shù)：仿真時長：20s；步長0.1s；求解器：ode4得到的仿真結(jié)果，如下圖所示： 墨綠色代表固有頻率=1時的特性曲線；淡綠色代表固有頻率=2時的特性曲線；紫色代表固有頻率=3時的特性曲線。從結(jié)果圖中可以看出，固有頻率決定了系統(tǒng)的振蕩頻率，固有頻率越大，系統(tǒng)的振蕩越高，響應(yīng)速度也越快。
公眾號
版權(quán)所有 ©北太振寰（重慶）科技有限公司 | 渝ICP備2022008375號-4 法律聲明 隱私政策
