查看“︁訊號與系統/由系統微分方程式求系統脈衝響應”︁的源代码

== '''線性非時變系統''' ==

 ●連續時間線性非時變系統一般可以線性常微分方程式表示 : 
  
<math>\sum_{n=0}^N a_n{d^ny(t) \over dt^n} = \sum_{m=0}^M bm{d^mx(t) \over dt^m}</math>

其中<math>\mathbf y(t)</math>為系統輸出，<math>\mathbf x(t)</math>為系統輸入。

 ●此一系統對於任意輸入訊號<math>\mathbf x(t)</math>的零狀態響應可用旋積運算求得 : 

<math>\mathbf y(t)=x(t)*h(t)</math>   
  
其中<math>\mathbf h(t)</math>為系統的單位脈衝響應。

 ●本附錄將介紹如何由描述系統的線性常微分方程式找出系統的單位脈衝響應。



----


== '''LIT系統的表示''' ==






----


== '''尋找系統單位脈衝響應的第一種方法''' ==

 步驟一 : 將系統的微分方程式中，等號右邊直接以<math>\mathbf x(t)</math>代替(即子系統<math>\mathbf \alpha</math>的輸入<math>\mathbf x(t)=\delta(t)</math> ，此時的輸出即為子系統<math>\mathbf \alpha</math>之單位脈衝響應<math>\mathbf y_1(t)=h_a(t)</math>)

 步驟二 : 將形成新方程式取積分，積分範圍由<math>\mathbf t=0^-</math>到<math>t=0^+</math>。觀察此一新方程式等號左、右兩邊的特性，可求得一組在<math>\mathbf t=0^+</math>時之初始條件。

 步驟三 : 利用步驟二的初始條件，可求得子系統<math>\mathbf \alpha</math>之單位脈衝響應<math>\mathbf h_a(t)</math>。

 步驟四 : 系統的單位脈衝響應<math>\mathbf h(t)</math>可由<math>\mathbf h_a(t)</math>及其微分之線性組合獲得。也就是 : 

<math>\mathbf h(t)=\sum_{m=0}^M bm{d^mx(t) \over dt^m}</math>




----


== '''範例3.23''' ==

 試求LTI系統<math>\tau_0{dy(t) \over dt}+y(t)={d \over dt}x(t)+x(t)</math> 之單位脈衝響應<math>h(t)</math>。

【解】
        
 步驟一 :<math>\tau_0{dh_a(t) \over dt}+h_a(t)=\delta(t)</math>

 步驟二 :<math>\int_{0^-}^{0^+} \tau{dh_0(t) \over dt}\, dt+\int_{0^-}^{0^+} h_0(t)\, dt=\int_{0^-}^{0^+} \delta(t)\, dt</math>

      
觀察步驟一的方程式，等號左邊包含<math>h_a(t)</math>及<math>h_a(t)</math>的一切微分，而等式右邊僅含有<math>\delta(t)</math>。故可推論，<math>h_a(t)</math>不會包含<math>\delta(t)</math>(否則等號左邊會有<math>\delta(t)</math>的一次微分)但在<math>t=0</math>處是一不連續點，即<math>h_a(t^+)-h_a(t^-)=k<\infty</math>。

<math>\Rightarrow \tau_0[h_a(0^+)-h_a(0^-)]+0=1</math>

 又由於系統單位脈衝響應是假設無初始條件即<math>h_a(0^-)=0</math>

<math>\Rightarrow h_a(0^+)={1 \over \tau_0}</math>

 步驟三 : 對於<math>t>0</math>，步驟一的方程式可寫成

<math>\tau_0{dh_a(t) \over dt}+h_a(t)=0</math>       <math>t>0</math> 

<math>\Rightarrow h_a(t)=Ae^{-t/ \delta_0}</math>       <math>t>0</math> 

 由步驟二知

<math>h_a(0^+)={1 \over \tau_0}</math>

<math>\Rightarrow h_a(t)={1 \over \tau_0}e^{-t/ \tau_0}</math>

 步驟四 : 系統單位脈衝響應

<math>h(t)={d \over dt}h_a(t)</math>

<math>={d \over dt}[{1 \over \tau_0}e^{-t/ \tau_0}u(t)]+{1 \over \tau_0}e^{-t/ \tau_0}u(t)</math>

<math>=-{1 \over \tau_0^2}e^{-t/ \tau_0}u(t)+{1 \over \tau_0}e^{-t/ \tau_0}\delta(t)+{1 \over \tau_0}e^{-t/ \tau_0}u(t)</math>

<math>=({1 \over \tau_0}-{1 \over \tau_0^2})e^(-t/ \tau_0)u(t)+{1 \over \tau_0}\delta(t)</math>




----


== '''尋找系統單位脈衝響應的第二種方法''' ==


 ●已知系統的單位脈衝響應等於系統的單位步階響應的微分。

 ●當系統的微分方程式等號右邊沒有輸入訊號<math>x(t)</math>的微分時，可先求出系統的單位步階響應再將其微分而得系統的單位脈衝響應<math>h(t)</math>。




----


== '''範例3.24''' ==


 一LTI系統之微分方程式為<math>{d^2y(t) \over dt^2}+3{dy(t) \over dt}+2y(t)=x(t)</math>試求其單位脈衝響應<math>h(t)</math>。

 【解】先求單位步階響應<math>s(t)</math>

<math>{d^2s(t) \over dt^2}+3{ds(t) \over dt}+2s(t)=u(t)</math>

採用傳統工程數學解微分方程的方法 : 

 (1)對於<math>\mathbf t>0</math>

<math>\mathbf {d^2s(t) \over dt^2}+3{ds(t) \over dt^2}+2s(t)=1</math>

<math>\mathbf (D^2+3D+2)s(t)=1</math>


 (2)令<math>\mathbf \lambda^2+3\lambda +2=0</math>  <math>\mathbf \Rightarrow \lambda=-2,-1</math>

所以 齊性解<math>\mathbf s_h(t)=Ae^{-2t}+Be^{-t}</math>

        
 (3)利用未定係數法求特解<math>\mathbf s_p(t)=k</math>代入方程式

<math>{d^2s_p(t) \over dt^2}+3{ds_p(t) \over dt}+2s_p(t)=1</math>

<math>\mathbf 0+0+2k=1</math>

<math>\Rightarrow k={1 \over 2}</math>

故<math>s_p(t)={1 \over 2}</math>

 (4)由(2)(3)知

<math>\mathbf s(t)=s_h(t)=s_h(t)+s_p(t)</math>

<math>\mathbf =Ae^{-2t}+Be^{-t}+{1 \over 2}</math>


 (5)因為<math>\mathbf {d^2s(t) \over dt^2}+3{ds(t) \over dt}+2s(t)=u(t)</math>，故<math>\mathbf s(t)</math>及<math>\mathbf {ds(t) \over dt}</math>在<math>\mathbf t=0</math>處均 必須為連續。

 又假設系統的初始條件為0，即<math>s(0^-)=0</math>，<math>s\prime(0^-)=0</math>。
          
故<math>\mathbf s(0^+)=s(0^-)=0</math>

<math>s\prime(0^+)=s\prime(0^-)=0</math>

 將(5)的初始條件代入(4) 

<math>\Rightarrow \begin{cases} A+B+{1 \over 2}=0 \\ -2A-B=0 &  \end{cases}</math>

<math>\Rightarrow \begin{cases} A={1 \over 2} \\ B=-1 \end{cases}</math>

所以<math>s(t)=[{1 \over 2}e^{-2t}+{1 \over 2}u(t)]</math>

 (6)系統的單位脈衝響應 :

<math>h(t)={ds(t) \over dt}=[-e^{-2t}+e^{-t}u(t)]+[{1 \over 2}e^{-2t}+{1 \over 2}]\delta(t)</math> 

<math>\mathbf =[-e^{-2t}+e{-t}]u(t)</math>

----
== '''參考資料''' ==

●B. P. Lathi, Signal Processing and Linear Systems, Berkeley-Cambridge Press, 1998.

●G. E. Carlson, Signal and Linear System Analysis, 2nd ed., John Wiley & Sons, 1998.

●余兆棠、李志鵬，信號與系統， 滄海書局，2007。

●Edward W. Kamen and Bonnie S. Heck, Fundamentals of Signals and System Using the Web and Matlab, 2nd ed.,Prentice Hall International, 2000.

●Rodger E. Ziemer, William H. Tranter, D. Ronald Fannin, Signals & Systems:  Continuous and Discrete, 4th ed.,Prentice Hall International, 1998.

●Charls L. Phillips, John M. Parr, Eve A. Riskin, Signals, Systems, and Transforms, 3rd ed., Pearson Education, 2003.

●Rodger E. Ziemer and William H. Tranter, Principles of Communications, John Wiley & Sons, 2002.

●Simon Haykin, Communication Systems, 4th ed., John Wiley & Sons,  2001.

●John G. Proakis and M. Salehi, Communication Systems Engineering, 2nd ed., Prentice Hall International, 2002.

●Benoit Boulet, Fundamentals of Signals and Systems, Da Vinci Engineering Press, 2005.