精品自拍影视三级_在线天堂网www资源AV_制服丝袜综合另类中文字幕_中国日产乱码无线码_十八禁啪啪污污网站免费下载_2021国产福利三级大片_新婚少妇在线观看一区_欧美人与动人物交欧美精品_长篇色翁荡熄小说全文_久草91视频三级亚洲aⅴ

1 線性插值

線性插值是指插值函數(shù)為一次多項(xiàng)式的插值方式，其在插值節(jié)點(diǎn)上的插值誤差為零。

線性插值相比其他插值方式，如拋物線插值，具有簡(jiǎn)單、方便的特點(diǎn)。

線性插值的幾何意義即為概述圖中利用過(guò)A點(diǎn)和B點(diǎn)的直線來(lái)近似表示原函數(shù)。

線性插值可以用來(lái)近似代替原函數(shù)，也可以用來(lái)計(jì)算得到查表過(guò)程中表中沒有的數(shù)值。

1.2 基礎(chǔ)知識(shí)

已知函數(shù)在區(qū)間上個(gè)互異點(diǎn)上的函數(shù)值，若存在一簡(jiǎn)單函數(shù)，使

并要求誤差

的絕對(duì)值在整個(gè)區(qū)間上比較小。這樣的問(wèn)題稱為插值問(wèn)題。

其中：

：插值節(jié)點(diǎn)

：被插值函數(shù)

：插值函數(shù)

：插值區(qū)間

如果在插值區(qū)間內(nèi)部用代替則稱為內(nèi)插；在插值區(qū)間以外，用代替則稱為外插。

1.3 簡(jiǎn)介

線性插值是一種較為簡(jiǎn)單的插值方法，其插值函數(shù)為一次多項(xiàng)式。線性插值，在各插值節(jié)點(diǎn)上插值的誤差為0。

設(shè)函數(shù)在兩點(diǎn)，上的值分別為，，求多項(xiàng)式

使?jié)M足

由解析幾何可知

稱為在處的一階均差，記以。于是，得

如果按照整理，則

以上插值多項(xiàng)式為一次多項(xiàng)式，這種插值稱為線性插值。

1.4 幾何意義

線性插值的幾何意義如圖1所示，即為利用過(guò)點(diǎn)和的直線來(lái)近似原函數(shù)。

1.5 應(yīng)用

1）線性插值在一定允許誤差下，可以近似代替原來(lái)函數(shù)；

2）在查詢各種數(shù)值表時(shí)，可通過(guò)線性插值來(lái)得到表中沒有的數(shù)值。

2一維線性插值仿真實(shí)例

首先，在Baltamulink中，添加輸出正弦波模塊、一維插值模塊、信號(hào)合并模塊、輸出模塊，建立一維插值仿真模型，每個(gè)模塊參數(shù)都設(shè)置為默認(rèn)值；模型如下圖所示；

設(shè)置仿真參數(shù)：

仿真時(shí)長(zhǎng)：10s；步長(zhǎng)0.01s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：

綠色代表原正弦波數(shù)據(jù)；

橙色表示一維插值后的正弦波數(shù)據(jù)。

      已知質(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)的北太真元模型，并運(yùn)行仿真。彈簧-質(zhì)量-阻尼系統(tǒng)如下圖所示：

建立理論數(shù)學(xué)模型。對(duì)于無(wú)外力的系統(tǒng)，根據(jù)牛頓定理可以得到：

mx’’ + bx’ + kx = 0

代入數(shù)值并整理得：

x’’ = -3x’ - 90x

在北太真元建立模型如下圖所示：

設(shè)置仿真參數(shù)：

仿真時(shí)長(zhǎng)：10s；步長(zhǎng)0.01s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：

      已知二自由度質(zhì)量-彈簧-阻尼串聯(lián)系統(tǒng)模型如下圖所示：

      這個(gè)系統(tǒng)由兩個(gè)質(zhì)量塊（小車）和三組彈簧阻尼器組成，假設(shè)地面是光滑的，這樣系統(tǒng)中沒有摩擦作用。 u1(t) 和 u2(t) 分別是兩個(gè)質(zhì)量塊所受的外力， x1(t) 和 x2(t) 分別是兩個(gè)質(zhì)量塊的位移。 m1,2、k1,2,3 和 b1,2,3 分別對(duì)應(yīng)圖中的質(zhì)量、彈簧剛度和阻尼系數(shù)。是對(duì)于整個(gè)系統(tǒng)而言，輸入兩個(gè)外力，輸出兩個(gè)位移，因此這是一個(gè)多輸入多輸出系統(tǒng)。兩個(gè)小車都只能沿橫向左右運(yùn)動(dòng)，因此為二自由度系統(tǒng)。

根據(jù)牛頓第二定律，物體所受合力等于物體的慣性力，而慣性力是物體質(zhì)量與加速度的乘積?？梢岳斫鉃槲矬w受力后產(chǎn)生加速度，而有加速度存在就會(huì)產(chǎn)生運(yùn)動(dòng)趨勢(shì)，造成物體運(yùn)動(dòng)。

對(duì)兩個(gè)質(zhì)量塊分別進(jìn)行受力分析：

首先需定義系統(tǒng)的狀態(tài)變量，在這種情況下一般用物體的位移和速度作為狀態(tài)變量。設(shè)狀態(tài)向量 Z分別對(duì)應(yīng)質(zhì)量塊1的位移和速度、質(zhì)量塊2的質(zhì)量和速度:

根據(jù)上面對(duì)模型的數(shù)學(xué)推導(dǎo)，當(dāng)將各狀態(tài)向量 z 對(duì)時(shí)間求導(dǎo)（微分），可以將系統(tǒng)整理為各狀態(tài)量的一階微分方程組：

因?yàn)檫@是一個(gè)線性系統(tǒng)，因此系統(tǒng)的狀態(tài)可以表示為矩陣形式：

式中，A 為系統(tǒng)狀態(tài)矩陣；B 為輸入矩陣；u 為輸入（控制）向量：

      有了系統(tǒng)的狀態(tài)空間方程，接下來(lái)考慮系統(tǒng)的輸出，也就是我們希望得到的量。我們假設(shè)對(duì)兩個(gè)質(zhì)量塊的位移感興趣，即 z1 和 z3 ，那么就將它兩個(gè)狀態(tài)作為系統(tǒng)的輸出，則輸出方程的矩陣形式為：y = Cz;

式中，y 為輸出向量；C 為輸出矩陣：

以上，使用狀態(tài)空間模型對(duì)系統(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ù)：

兩個(gè)質(zhì)量塊所受的外力u1是常量模塊 = 2 N； u2是階躍型號(hào)=1N。

仿真時(shí)長(zhǎng)：10s；步長(zhǎng)0.01s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：

求解二階微分方程：x’’(t) + 0.4x’(t) + 0.9x(t) = 0.7u(t) 的解，其中u(t)是脈沖信號(hào)。

在北太真元建立模型如下圖所示：

設(shè)置仿真參數(shù)：

脈沖信號(hào)u(t)為方波信號(hào)模塊，參數(shù)振幅 = 1；周期 = 2.5；脈沖寬度 = 50；相位 = 0；

仿真時(shí)長(zhǎng)：10s；步長(zhǎng)0.01s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：

力——質(zhì)量系統(tǒng)，要拉動(dòng)一個(gè)箱子（拉力f=1N），箱子質(zhì)量為M（1kg），箱子與地面的摩擦力為[(b=0.4N.m/s)]，其大小與車子的速度成正比。如下圖所示：

其運(yùn)動(dòng)方程式為：

F - bx’ = Mx’’

拉力作用時(shí)間為2s。

在北太真元建立模型如下圖所示：

設(shè)置仿真參數(shù)：

設(shè)置stepInputOne模塊的階躍時(shí)間為0，表示摩擦力作用時(shí)間；

設(shè)置stepInputTwo模塊的階躍時(shí)間為2，表示拉力作用時(shí)間；

設(shè)置gainTwo模塊的增益值為0.4表示摩擦力；

仿真時(shí)長(zhǎng)：10s；步長(zhǎng)0.01s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：

因?yàn)橛心Σ亮Υ嬖?，箱子最終會(huì)停止前進(jìn)。

零極點(diǎn)增益模型實(shí)際上是傳遞函數(shù)的另一種表現(xiàn)形式，其原理是分別對(duì)源系統(tǒng)傳遞函數(shù)的分子、分母進(jìn)行因式分解處理，以獲得系統(tǒng)零點(diǎn)和幾點(diǎn)的表示形式。

式中，k為系統(tǒng)增益；zi(i = 1,2,3,...,m)為分子多項(xiàng)式的根，稱為系統(tǒng)的零點(diǎn)；pj(j = 1,2,...,n)是分母多項(xiàng)式的根，稱為系統(tǒng)的極點(diǎn)。

傳遞函數(shù)的分母多項(xiàng)式就是它的特征多項(xiàng)式，它等于零的方程就是傳遞函數(shù)的特征方程，特征方程的根也就是傳遞函數(shù)的極點(diǎn)。傳遞函數(shù)的極點(diǎn)決定了所描述系統(tǒng)的自由運(yùn)動(dòng)狀態(tài)；零點(diǎn)影響系統(tǒng)各模態(tài)在系統(tǒng)響應(yīng)中的比重。

零點(diǎn)增益模型的命令格式如下：

ZPG = zpk(z, p, k)

其中ZPG是建立的零極點(diǎn)增益模型；z、p、k分別是系統(tǒng)的零點(diǎn)向量、極點(diǎn)向量和增益。

例：利用Baltamulink建立系統(tǒng)

G(s) = 18(s + 2) / (s + 0.4)(s + 15)(s + 25)

的零點(diǎn)增益模型，進(jìn)行系統(tǒng)仿真。

將上面零點(diǎn)增益模型進(jìn)行轉(zhuǎn)換得到模型的傳遞函數(shù)如下：

G(s) = 18(s + 2)/(s^3 + 40.4s^2 + 391s + 150)

根據(jù)該傳遞函數(shù)模型，在北太真元建立模型如下圖所示：

設(shè)置仿真參數(shù)：

仿真時(shí)長(zhǎng)：20s；步長(zhǎng)0.01s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：

針對(duì)標(biāo)準(zhǔn)二階系統(tǒng)傳遞函數(shù)：

改變?chǔ)疲ㄗ枘岜龋┖挺豱（自由振蕩頻率）的參數(shù)設(shè)置，觀察對(duì)系統(tǒng)輸出的影響。

在二階系統(tǒng)自由振蕩頻率ωn不變的情況下，改變阻尼系數(shù)ζ為無(wú)阻尼（ζ= 0）、欠阻尼（0<ζ< 1）、臨界阻尼（ζ= 1）和過(guò)阻尼（ζ> 1）的4中狀態(tài)，分別取ζ= 0，ζ= 0.5，ζ= 1，ζ= 2帶入二階系統(tǒng)傳遞函數(shù)  中，搭建4個(gè)不同的仿真模型，輸出結(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]；

階躍信號(hào)模塊的階躍時(shí)間為0；

仿真時(shí)長(zhǎng)：6s；步長(zhǎng)0.01s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：

結(jié)論：從仿真結(jié)果可以看出，改變阻尼比，系統(tǒng)的超調(diào)量也在變化，系統(tǒng)達(dá)到穩(wěn)態(tài)的時(shí)間也發(fā)生變化。當(dāng)ζ= 0無(wú)阻尼時(shí)，出現(xiàn)等幅振蕩曲線，超調(diào)量為100%，穩(wěn)態(tài)時(shí)間為無(wú)窮大；當(dāng)ζ< 1時(shí)，信號(hào)曲線衰減振蕩，有超調(diào)量；當(dāng)ζ≥ 1時(shí)，沒有超調(diào)量，隨著ζ增大，達(dá)到穩(wěn)態(tài)的時(shí)間也增大。

已知開環(huán)系統(tǒng)傳遞函數(shù)：

首先將G(s)化為唯一標(biāo)準(zhǔn)形式：

此系統(tǒng)由比例環(huán)節(jié)、積分環(huán)節(jié)、慣性環(huán)節(jié)、一階微分環(huán)節(jié)和振蕩環(huán)節(jié)共5個(gè)環(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]；

仿真時(shí)長(zhǎng)：10s；步長(zhǎng)0.01s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：

已知某系統(tǒng)的開環(huán)對(duì)數(shù)頻率特性如下圖所示，試確定其開環(huán)傳遞函數(shù)。

根據(jù)對(duì)數(shù)幅頻特性曲線，可以寫出開環(huán)傳遞函數(shù)的表達(dá)形式如下：

根據(jù)對(duì)數(shù)頻率特性的坐標(biāo)特點(diǎn)有，可以確定開環(huán)增益。

根據(jù)相頻特性的變化趨勢(shì)（-270°-> -90°），可以判定系統(tǒng)為非最小相角系統(tǒng)。

G(s)中一階復(fù)合微分環(huán)節(jié)和慣性環(huán)節(jié)至少有一個(gè)是“非最小相角”的，將系統(tǒng)可能的開環(huán)零點(diǎn)極點(diǎn)分布畫出來(lái)，如下表所示：

分析相角的變化趨勢(shì)，可見，只有當(dāng)慣性環(huán)節(jié)極點(diǎn)在右半s平面，一階復(fù)合微分環(huán)節(jié)零點(diǎn)在左半s平面是，相角才符合從-270°到-90°的變化規(guī)律。因此可以確定系統(tǒng)的開環(huán)傳遞函數(shù)為：

對(duì)于最小相角系統(tǒng)，對(duì)數(shù)幅頻特性與對(duì)數(shù)相頻特性之間存在唯一確定的對(duì)應(yīng)關(guān)系，根據(jù)對(duì)數(shù)幅頻特性就完全可以確定相應(yīng)的對(duì)數(shù)相頻特性和傳遞函數(shù)，反之亦然。由于對(duì)數(shù)幅頻特性容易繪制，所以在分析最小相角系統(tǒng)時(shí)，通常只畫其對(duì)數(shù)幅頻特性，對(duì)數(shù)相頻特性則只需概略畫出，或者不畫。

根據(jù)該傳遞函數(shù)模型，在北太真元搭建最小相角系統(tǒng)模型如下圖所示：

設(shè)置仿真參數(shù)：

從上到下，傳遞函數(shù)參數(shù)依次為：

num = [1]；den = [1 0]；

num = [1 1]；den = [1 -1]；

仿真時(shí)長(zhǎng)：10s；步長(zhǎng)0.01s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：

控制系統(tǒng)穩(wěn)定與否是絕對(duì)穩(wěn)定性的概念。而對(duì)一個(gè)穩(wěn)定的系統(tǒng)而言，還有一個(gè)穩(wěn)定的程度，即相對(duì)穩(wěn)定性的概念。相對(duì)穩(wěn)定性與系統(tǒng)的動(dòng)態(tài)性能指標(biāo)有著密切的關(guān)系。在設(shè)計(jì)一個(gè)控制系統(tǒng)時(shí)，不僅要求它必須是絕對(duì)穩(wěn)定的，而且還應(yīng)保證系統(tǒng)具有一定的穩(wěn)定程度。只有這樣，才能不致因系統(tǒng)參數(shù)變化而導(dǎo)致系統(tǒng)性能變差甚至不穩(wěn)定。

對(duì)于一個(gè)最小相角系統(tǒng)而言，曲線越靠近點(diǎn)，系統(tǒng)階躍響應(yīng)的振蕩就越強(qiáng)烈，系統(tǒng)的相對(duì)穩(wěn)定性就越差。因此，可用曲線對(duì)點(diǎn)的接近程度來(lái)表示系統(tǒng)的相對(duì)穩(wěn)定性。通常，這種接近程度是以相角裕度和幅值裕度來(lái)表示的。

要計(jì)算相角裕度，首先要知道截止頻率。求較方便的方法是先由繪制曲線，由與線的交點(diǎn)確定。而求幅值裕度首先要知道相角交界頻率，對(duì)于階數(shù)不太高的系統(tǒng)，直接解三角方程是求較方便的方法。通常是將寫成虛部和實(shí)部， 令虛部為零而解得。

【例6-10】 某單位反饋系統(tǒng)的開環(huán)傳遞函數(shù)為：

試求時(shí)系統(tǒng)的相角裕度和幅值裕度。將該開環(huán)傳遞函數(shù)變換為：

在實(shí)際工程設(shè)計(jì)中，只要繪出曲線即可。

根據(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]；

仿真時(shí)長(zhǎng)：10s；步長(zhǎng)0.1s；求解器：ode4

得到的仿真結(jié)果，如下圖所示：