【信號與系統】筆記(2)——連續系統的時域分析

Author:AXYZdong
自動化專業 工科男
有一點思考,有一點想法,有一點理性!

前言

以下內容是關於連續系統的時域分析,重點難點,考試常考,考前複習。

一、系統的微分方程及其響應

1、LTI系統的微分方程

描述線性時不變(LTI)系統的輸入–輸出特性的是常係數線性微分方程

時域分析方法:從系統的模型(微分方程)出發,在時域研究輸入信號通過系統響應後的變化規律,是研究系統時域特性的重要方法。

對於電系統,建立其微分方程的依據:
KCL:i(t)=0KVL:u(t)=0VCR:UR(t)=Ri(t),uL(t)=Ldi(t)dt,iC(t)=Cdu(t)dt KCL:\sum i(t)=0\\ KVL: \sum u(t)=0\\ VCR:U_R(t)=R \cdot i(t),u_L(t)=L\cdot \frac{di(t)}{dt},i_C(t)=C \cdot \frac{du(t)}{dt}

在這裏插入圖片描述
對於圖(a)中RC電路
RCduC(t)dt+uC(t)=uS(t):uC(t)+1RCuC(t)=1RCuS(t) RC \cdot \frac{du_C(t)}{dt} +u_C(t)=u_S(t)\\ 即:u ' _C(t) + \frac{1}{RC} \cdot u_C(t) = \frac{1}{RC} \cdot u_S(t)
對於圖(b)中RL電路
LdiL(t)Rdt+iL(t)=iS(t):iL(t)+RLiL(t)=RLiS(t) \frac{L \cdot di_L(t)}{R \cdot dt} +i_L(t)=i_S(t)\\ 即:i ' _L(t) + \frac{R}{L} \cdot i_L(t) = \frac{R}{L} \cdot i_S(t)
以上可以總結出一般形式:y(t)+ay(t)=χ(t)y '(t)+ay(t)=\chi(t)

2、系統的響應

2.1、零輸入響應(儲能響應)

從觀察的初始時刻起不再施加輸入信號,僅由該時刻系統本身的起始儲能狀態引起的響應稱爲零輸入響應(ZIR)

2.2、零狀態響應(受激響應)

當系統的儲能狀態爲零時,由外加激勵信號(輸入)產生的響應稱爲零狀態響應(ZSR)

對於一階系統方程
y(t)+ay(t)=χ(t) y '(t)+ay(t)=\chi(t)
零狀態響應:
yZS(t)=eat0tχ(t)eaτdτ,(t0) y_{ZS}(t)=\rm e^{-at} \int_{0-}^t \chi(t) e^{a \tau} d \tau,(t \geq0)

2.3、完全響應

y(t)=yZI(t)+yZS(t),[=+] y(t)=y_{ZI}(t)+y_{ZS}(t),[全響應=零輸入+零狀態]

2.4、例一(y(t)+6y(t)+8y(t)=f(t),t>0y ''(t) + 6y '(t) + 8y(t) =f(t),t>0)

已知某二階限行連續時間系統的動態方程:
y(t)+6y(t)+8y(t)=f(t),t>0y(0)=1,y(0)=2, y''(t) + 6y'(t) + 8y(t) =f(t),t>0\\ 初始條件:y({0_-})=1,y' ({0_-})=2,求系統的零輸入響應
解:

s2+6s+8=0,s1=2s2=4yX(t)=k1e2t+k2e4t,t>0yX(0+)=yX(0)=yX(0)=k1+k2=1y(0)=yX(0)=2k14k2=2yX(0)y(0)k1=3k2=2系統的特徵方程:s^2+6s+8=0,特徵根:s_1=-2,s_2=-4 \\ 故系統的零輸入響應:y_X(t)=k_1e^{-2t} + k_2e^{-4t},t>0 \\ y_X({0_+})=y_X({0_-}) = y_X({0}) = k_1+k_2=1\\ y'(0)=y_X'({0_-}) = -2k_1-4k_2=2\\ y_X({0})和y'(0)代入 可解出:k_1=3,k_2=-2

可得零輸入響應:yX(t)=3e2t2e4t,t0y_X(t)=3e^{-2t}-2e^{-4t},t \ge0

2.5、例二(y(t)+3y(t)+2y(t)=2f(t)+6f(t),t>0y ''(t) + 3y'(t) + 2y(t) =2f'(t)+6f(t),t>0)

y(t)+3y(t)+2y(t)=2f(t)+6f(t),t>0y(0)=2,y(0)=0,f(t)=ϵ(t), y ''(t) + 3y'(t) + 2y(t) =2f'(t)+6f(t),t>0\\ 初始條件:y({0_-})=2,y \prime ({0_-})=0,f(t)= \epsilon(t), \\ 求系統的零輸入響應和零狀態響應\\
解:
(1)yX(t)0yX(t)滿(1)零輸入響應y_X(t)激勵爲0,故y_X(t)滿足

yX(t)+3yX(t)+2yX(t)=0s2+3s+2=0,s1=2s2=1yX(t)=k1e2t+k2et,t>0yX(0+)=yX(0)=yX(0)=2y(0)=yX(0)=0yX(0)y(0)k1=4k2=2y _X''(t) + 3y_X'(t) + 2y_X(t) =0\\ 系統的特徵方程:s^2+3s+2=0,特徵根:s_1=-2,s_2=-1 \\ 故系統的零輸入響應:y_X(t)=k_1e^{-2t} + k_2e^{-t},t>0 \\ y_X({0_+})=y_X({0_-}) = y_X({0}) =2\\ y'(0)=y_X'({0_-}) = 0\\ y_X({0})和y'(0)代入 可解出:k_1=4,k_2=-2

可得零輸入響應:yX(t)=4et2e2t,t>0y_X(t)=4e^{-t}-2e^{-2t},t >0

(2)yf(t)滿(2)零狀態響應y_f(t)滿足
yf(t)+3yf(t)+2yf(t)=2δ(t)+6ϵ(t),yf(0)=yf(0+)=0 y _f''(t) + 3y_f'(t) + 2y_f(t) =2 \delta (t)+6 \epsilon(t),並有y_f({0_-}) = y_f({0_+}) = 0\\
δ(t),yf(t)δ(t)yf(t),yf(0+)yf(0)yf(t)t=0yf(0)=yf(0+)=0[yf(0+)yf(0)]+3[yf(0)yf(0+)]+200+yf(t)dt=2+600+ϵ(t)dt由於上式等號右端含有 \delta(t) , 故y _f''(t)含有 \delta(t),從而y_f'(t)躍變 ,即y'_f({0_+}) \neq y'_f({0_-})\\ 而y_f(t) 在t=0處連續,即y_f({0_-}) = y_f({0_+})=0,上式兩邊積分有:\\ [y'_f({0_+}) - y'_f({0_-})]+3[y_f({0_-}) -y_f({0_+})]+2\int_{0-}^{0+}y_f(t) dt=2+6\int_{0-}^{0+} \epsilon(t)dt

整理得:yf(0+)=2+yf(0)=2y'_f({0_+}) =2+ y'_f({0_-})=2

t>0t>0時,有yf(t)+3yf(t)+2yf(t)=6y _f''(t) + 3y_f'(t) + 2y_f(t) =6

不難求出其齊次解爲:Cf1et+Cf2e2t,C_{f1}e^{-t}+C_{f2}e^{-2t},其特解爲常數33

yf(t)=Cf1et+Cf2e2t+3\therefore y_f(t)=C_{f1}e^{-t}+C_{f2}e^{-2t}+3

代入初值求得:yf(t)=4et+e2t+3,t0y_f(t)=-4e^{-t}+e^{-2t}+3,t \geq0

2.6、響應的分類

分類標準 對應響應
響應的不同起因 儲能響應、受激響應
系統的性質和輸入信號的性質 自由響應(取決於系統性質,即特徵根)、強迫響應(取決於輸入信號的形式)
響應的變化形式 瞬態響應(t無限增大,響應趨於零)、穩態響應(響應恆定或爲某個穩態函數)

二、階躍響應

1、定義

LTI系統在零狀態下,由單位階躍信號引起的響應稱爲單位階躍響應,簡稱階躍響應,記爲s(t)s(t)
在這裏插入圖片描述

2、一階系統方程的階躍響應

對於一階系統方程
y(t)+ay(t)=bϵ(t) y '(t)+ay(t)=b\epsilon(t)
則零狀態響應:
yZS(t)=eat0tχ(t)eaτdτ,(t0) y_{ZS}(t)=\rm e^{-at} \int_{0-}^t \chi(t) e^{a \tau} d \tau,(t \geq0)
則階躍響應:
y(t)=s(t)=eat0tbϵ(t)eaτdτ=ba(1eat),(t0) y (t)=s(t)=\rm e^{-at} \int_{0-}^t b\epsilon(t) e^{a \tau} d \tau = \frac {b}{a} \cdot(1-e^{-at}),(t \geq0)

3、階躍響應的測量

在這裏插入圖片描述

三、衝激響應

1、定義

儲能狀態爲零的系統,在單位衝激信號作用下產生的零狀態響應稱爲衝激響應,記爲h(t)h(t)

2、一階系統的衝激響應

對於一階系統方程
y(t)+ay(t)=bδ(t) y '(t)+ay(t)=b\delta(t)
則衝激響應:
y(t)=h(t)=eat0tbδ(t)eaτdτ=beatϵ(t) y (t)=h(t)=\rm e^{-at} \int_{0-}^t b\delta(t) e^{a \tau} d \tau = b \cdot e^{-at} \cdot \epsilon(t)

3、例一(y(t)+5y(t)+6y(t)=f(t)y''(t)+5y'(t)+6y(t)=f(t)

描述某系統的微分方程爲y(t)+5y(t)+6y(t)=f(t)y''(t)+5y'(t)+6y(t)=f(t),求其衝激響應h(t)h(t)

解:根據h(t)h(t)定義有h(t)+5h(t)+6h(t)=δ(t),h(0)=h(0)=0h(0+)h(0+)h''(t)+5h'(t)+6h(t)= \delta(t),\\ 並且有 h'({0_-}) = h({0_-}) = 0,先求h'({0_+}) 和 h({0_+})

δ(t),h(t)δ(t)h(t),h(0+)h(0)h(t)t=0h(0)=h(0+)=0[h(0+)h(0)]+5[h(0)h(0+)]+600+h(t)dt=1由於上式等號右端含有 \delta(t) , 故h''(t)含有 \delta(t),從而h'(t)躍變 ,即h' ({0_+}) \neq h'({0_-})\\ 而h(t) 在t=0處連續,即h({0_-}) = h({0_+})=0,上式兩邊積分有:\\ [h'({0_+}) - h'({0_-})]+5[h({0_-}) -h({0_+})]+6\int_{0-}^{0+}h(t) dt=1

整理得:h(0+)=1+h(0)=1h'({0_+}) =1+ h'({0_-})=1

t>0t>0時,有h(t)+5h(t)+6h(t)=0h''(t)+5h'(t)+6h(t)=0

不難求出其齊次解爲:C1e2t+C2e3t,C_{1}e^{-2t}+C_{2}e^{-3t},

h(t)=(C1e2t+C2e3t)ϵ(t)\therefore h(t)=(C_{1}e^{-2t}+C_{2}e^{-3t}) \cdot \epsilon(t)

代入初值求得:h(t)=(e2te3t)ϵ(t)h(t)=(e^{-2t}-e^{-3t}) \cdot \epsilon(t)

4、例二(y(t)+5y(t)+6y(t)=f(t)+2f(t)+3f(t)y''(t)+5y'(t)+6y(t)=f''(t)+2f'(t)+3f(t)

描述某系統的微分方程爲y(t)+5y(t)+6y(t)=f(t)+2f(t)+3f(t)y''(t)+5y'(t)+6y(t)=f''(t)+2f'(t)+3f(t),求其衝激響應h(t)h(t)

解:根據h(t)h(t)定義有
h(t)+5h(t)+6h(t)=δ(t)+2δ(t)+3δ(t)(1)h''(t)+5h'(t)+6h(t)=\delta''(t) + 2\delta'(t) +3\delta(t) \\ \tag{1}
h(0)=h(0)=0h(0+)h(0+)並且有 h'({0_-}) = h({0_-}) = 0,先求h'({0_+}) 和 h({0_+})
由方程可知:h(t)h(t)中含有δ(t)\delta(t)
故令:
{h(t)=aδ(t)+P1(t),[Pi(t)δ(t)]h(t)=aδ(t)+bδ(t)+P2(t)h(t)=aδ(t)+bδ(t)+cδ(t)+P3(t)\begin{cases} h(t)=a \delta(t) +P_1(t) , [P_i(t)中爲不含有\delta(t的函數) ]\\ h'(t)=a \delta'(t) + b \delta(t) + P_2(t)\\ h''(t)=a \delta''(t) + b \delta'(t) + c \delta(t) + P_3(t)\\ \end{cases}
代入式(1)整理得:
aδ(t)+(b+5a)δ(t)+(6a+5b+c)δ(t)+P1(t)+P2(t)+P3(t)=δ(t)+2δ(t)+3δ(t) a \delta''(t) + (b+5a) \delta'(t) + (6a+5b+c) \delta(t) + P_1(t)+P_2(t)+ P_3(t)=\delta''(t) + 2\delta'(t) +3\delta(t)
利用δ(t)\delta(t)係數匹配,得a=1,b=3,c=12a=1,b=-3,c=12
h(t)=δ(t)+P1(t)(2)h(t)= \delta(t) +P_1(t) \tag{2}
h(t)=δ(t)3δ(t)+P2(t)(3)h'(t)= \delta'(t) -3 \delta(t) + P_2(t) \tag{3}
h(t)=δ(t)3δ(t)+12δ(t)+P3(t)(4)h''(t)= \delta''(t) -3 \delta'(t) + 12 \delta(t) + P_3(t) \tag{4}
300+h(0+)h(0)=3400+h(0+)h(0)=12對式(3)從0_{-}到0_{+}積分h(0_{+})-h(0_{-})=-3 \\ 對式(4)從0_{-}到0_{+}積分h'(0_{+})-h'(0_{-})=12

t>0t>0時,有h(t)+5h(t)+6h(t)=0,23h''(t)+5h'(t)+6h(t)=0,特徵根-2,-3

不難求出其齊次解爲:C1e2t+C2e3t,C_{1}e^{-2t}+C_{2}e^{-3t},

h(t)=(C1e2t+C2e3t)ϵ(t)\therefore h(t)=(C_{1}e^{-2t}+C_{2}e^{-3t}) \cdot \epsilon(t)

代入初始條件h(0+)=3,h(0+)=12h(0_{+})=-3,h'(0_{+}) =12求得:h(t)=(3e2t6e3t)ϵ(t)h(t)=(3e^{-2t}-6e^{-3t}) \cdot \epsilon(t)

2h(t)=δ(t)+(3e2t6e3t)ϵ(t)結合式(2),h(t)= \delta(t)+(3e^{-2t}-6e^{-3t}) \cdot \epsilon(t)

5、階躍響應與衝激響應的關係

{h(t)=ds(t)dts(t)=th(τ)dτ \begin{cases} h(t)= \frac{ds(t)}{dt}\\ \\ s(t)= \int _{-\infty}^{t}h( \tau ) d \tau \end{cases}

四、卷積及其應用

1、信號的時域分解與卷積積分

1.1任意信號的分解

在這裏插入圖片描述
0"“0"號脈衝高度f(0)f(0),寬度Δ\Delta,用p(t)p(t)表示爲 f(0)Δp(t)f(0) \Delta p(t)

1"“1"號脈衝高度f(Δ)f(\Delta),寬度Δ\Delta,用p(tΔ)p(t-\Delta)表示爲 f(Δ)Δp(tΔ)f(\Delta) \Delta p(t-\Delta)

1"“-1"號脈衝高度f(Δ)f(-\Delta),寬度Δ\Delta,用p(t+Δ)p(t+\Delta)表示爲 f(Δ)Δp(t+Δ)f(-\Delta) \Delta p(t+\Delta)

f(t)^=n=f(nΔ)Δp(tnΔ)\hat {f(t)}=\sum_{n=-\infty}^{\infty}f(n\Delta) \Delta p(t-n\Delta)

limΔ0f(t)^=f(t)=f(τ)δ(tτ)dτ \lim_{\Delta \to 0} \hat{f(t)} =f(t) = \int _{-\infty}^{\infty} f(\tau) \delta(t- \tau) d \tau

1.2任意信號作用下的零輸入響應

在這裏插入圖片描述
根據h(t)h(t)定義:
δ(t)    h(t) \delta(t) \implies h(t)
由時不變性:
δ(tτ)    h(tτ) \delta(t - \tau) \implies h(t - \tau)
由齊次性:
f(τ)δ(tτ)    f(τ)h(tτ) f(\tau) \delta(t - \tau) \implies f(\tau) h (t - \tau)
由疊加性:
+f(τ)δ(tτ)dτ    +f(τ)h(tτ)dτ \int_{- \infty}^{+\infty} f(\tau) \delta(t - \tau) d\tau \implies \int_{- \infty}^{+\infty} f(\tau) h (t - \tau) d\tau
yf(t)=+f(τ)h(tτ)dτ, y_f(t)=\int_{- \infty}^{+\infty} f(\tau) h(t - \tau) d\tau ,卷積積分

1.3卷積積分的定義

已知定義在區間(,+)(-\infty,+\infty)上的兩個函數f1(t)f_1(t)f2(t)f_2(t),則定義積分f(t)=+f1(τ)f2(tτ)dτf(t)=\int_{- \infty}^{+\infty} f_1(\tau) f_2(t - \tau) d\tau
f1(t)f_1(t)f2(t)f_2(t)的卷積積分,簡稱卷積;記爲
f(t)=f1(t)f2(t) f(t)=f_1(t) *f_2(t)
溫馨提醒:τ\tau 爲積分變量,積分後的結果爲關於 tt 的函數
yf(t)=+f(τ)h(tτ)dτ=f1(t)h(t) y_f(t)=\int_{- \infty}^{+\infty} f(\tau) h(t - \tau) d\tau =f_1(t) *h(t)

2、卷積的圖解法

f1(t)f2(t)=+f1(τ)f2(tτ)dτ f_1(t) *f_2(t) = \int_{- \infty}^{+\infty} f_1(\tau) f_2(t - \tau) d\tau
四步圖解法:
(1)換元:tτf1(τ),f2(τ)t 換爲 \tau \to f_1(\tau), f_2(\tau)

(2)反轉平移(摺疊平移):f2(τ)f2(τ),tf2((τt))=f2(tτ)f_2(\tau)反轉\to f_2(-\tau),再右移t\to f_2(-(\tau -t))=f_2(t- \tau)【左加右減在τ\tau的裏面】

(3)乘積:f1(τ)f2(tτ)f_1(\tau) f_2(t - \tau)

(4)積分:τ\tau-\infty++\infty對乘積項積分

在這裏插入圖片描述
總結:圖解法步驟比較繁雜,但是按照四步法“戰略”就可以一步步把題目搞定。不過,圖解法對於求某一時刻的卷積值還是比較方便的,對於簡單的信號,通過畫圖就可以直觀的求出某一時刻的卷積值。

五、卷積積分的性質

1、奇異(衝激)函數的卷積特性

1、f(t)δ(t)=δ(t)f(t)=f(t)f(t)*\delta(t) = \delta(t)* f(t) = f(t)

2、f(t)δ(t)=f(t),f(t)δ(n)(t)=f(n)(t)f(t)* \delta'(t) = f'(t) , f(t)* \delta^{(n)}(t) = f^{(n)}(t)

3、f(t)ϵ(t)=+f(τ)ϵ(tτ)dτ=tf(τ)dτf(t)*\epsilon(t) = \int_{-\infty}^{+\infty} f(\tau) \epsilon(t- \tau)d\tau = \int_{-\infty}^{t} f(\tau) d\tau

ϵ(t)ϵ(t)=tϵ(t)\epsilon(t)*\epsilon(t) = t\epsilon(t)【重點關注】

2、卷積的微積分性質

在這裏插入圖片描述

3、卷積的時移特性

卷積的時移特性說白了就是:一個卷積積分,時間 tt 無論左移還是右移,其積分值等於相應函數左移或右移後的函數值。下面通過一個公式來說明:
若:f(t)=f1(t)f2(t),f(t)=f_1(t)*f_2(t),
則:f1(tt1)f2(tt2)=f1(tt1t2)f2(t)=f1(t)f2(tt1t2)=f(tt1t2)f_1(t-t_1)*f_2(t-t_2) \\ = f_1(t-t_1-t_2)*f_2(t)\\ = f_1(t)*f_2(t-t_1-t_2)\\ =f (t-t_1-t_2)

總結

連續系統的時域分析,重點難點內容,多看、多做、多動手。

零輸入、零狀態、階躍、衝激
卷積積分【重點關注】

上一篇筆記承蒙各位小夥伴的厚愛,訪問量蹭蹭上漲,然而我不食言,繼續更,希望喜歡。
如果文章中有哪些知識點寫錯的地方,歡迎私信我
如果文章中哪些知識點不理解的地方,歡迎私信我
如果文章中還有哪些需要補充的地方,歡迎私信我

碼字不易,大家的支持就是我堅持下去的動力。點贊後不要忘了關注我哦!

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章