訊號與系統/指數傅立葉級數

来自testwiki
跳转到导航 跳转到搜索

Template:NoteTA

指數傅立葉級數

假設X(t)為一週期訊號,基本週期為0,則X(t)可表示成複指數訊號的和: χ(𝐭)=...+X2ej2ω0t+X1ejω0t+𝐗𝟎+X1ejω0t+X2ej2ω0t+...

=n=Xnejnω0t


其中 ω0=2π𝑻0


指數傅立葉級數的係數

三角傅立葉級數

     χ(𝐭)=a0+n=1(ancosnω0t+bnsinnω0t)


尤拉公式 :

    sinnω0t=ejnω0tejnω0t2j


    cosnω0t=ejnω0t+ejnω0t2


x(t)=a0+n=1(anejnω0t+ejnω0t2+bnejnω0tejnω0t2j)

  =a0+n=112(anjbn)ejnω0t+n=112(an+jbn)ejnω0t
  =n=112(an+jbn)ejnω0t+a0+n=112(anjbn)ejnω0t


X(t)=...+X2ej2ω0t+X1ejω0t+X0+X1ejω0t+X2ej2ω0t+...

  =n=Xnejnω0t




其中 Xn={12(an+jbn)n<0a0n=012(anjbn)n>0


以上為指數傅立葉級數的係數Xn與三角傅立葉級數的係數anbn 之關係式


三角傅立葉級數第二式

    X(t)=C0+n=1Cncos(nω0t+θn)


尤拉公式


    Cncos(nω0t+θn)=Cn2[ej(nω0t+θn)+ej(nω0t+θn)]
               =(Cn2ejθn)ejnω0t+(Cn2ejθn)ejnω0t
                              
   X(t)=C0+n=1(Cn2ejθn)ejnω0t+n=1(Cn2ejθn)ejnω0t   
   =n=Xnejnω0t        




        其中     Xn={Cn2ejθnn<0C0n=0Cn2ejθnn>0

以上為指數傅立葉級數的係數Xn與三角傅立葉級數第二式的係數Cnθn之關係式


由上說明可知,三角傅立葉級數與指數傅立葉級數是完全等效的。



指數傅立葉級數係數的直接算法

假設一週期訊號X(t)的基本週期T0,此一訊號可表示成指數傅立葉級數 :

                X(t)=n=Xnejnω0t


其中 ω0=2π𝑻0



直接求係數 Xk:

   等號兩邊同乘上ejkω0t再對時間 t 積分,積一個基本週期T0的時間範圍 : 
   T0X(t)ejkω0tdt=T0(n=Xnejnω0t)ejkω0tdt
   =n=XnT0ej(nk)ω0tdt={0nkT0n=k
  =XkT0 


Xk=1T0T0X(t)ejkω0tdt


範例4.12

下圖為一方波週期訊號X(t)之時域波形,其週期為T0(ω0=2ΠT0),求此訊號之指數傅立葉級數。

...............圖


一方波週期訊號之時域波形 ©余兆棠、李志鵬,信號與系統, 滄海書局,2007。



【解】

      X0=1T0T02T02X(t)dt=1T0T04T041dt=12
     Xn=1T0T02T02X(t)ejnω0tdt=1T0T04T04ejnω0tdt


    =1j2πn[ejnπ2ejnπ2]=1nπ[ejnπ2ejnπ2j2]


    =sin(nπ2)nπ={0,n=2K0(1)k1(2k+1)π,n=2K+1



整理並改寫上述係數為:

X0=12; X2k=0,k0; X2k+1=(1)k1(2K+1)π


最後將此方波週期訊號表示成以下之指數傅立葉級數式展開式。

X(t)= n=Xnejnω0t=12+n=(1)k(2K+1)πej(2k+1)ω0t


訊號X(t)的雙邊振幅頻譜和相位頻譜


............圖


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

範例4.13

請將訊號表示成指數傅立葉級數展開式

Xc(t)=3ej(2000πt+π6)+4ej(4000πt+π3)+ej(6000πt+π6)

【解】

Xc(t)=n=ejnω0

          =3ej(2000πt+π6)+4ej(4000πt+π3)+ej(6000πt+π6)

訊號由3項複指數函數組合而成,其頻率分別為1000、2000與3000Hz,先找出3項頻率的最大公因數為1000,因此訊號是一週期訊號且其基本頻率為1000 Hz(ω0=2000πrad/s),仔細觀察可發現其表示式已是指數傅立葉級數之型式,因此只要改寫成為



Xc(t) = 3ej(2000πt+π6)+4ej(4000πt+π3)+ej(6000πt+π6)

= 3ejπ/6ej2000πt+4ejπ/3ej4000πt+ejπ/6ej6000πt

= 3ejπ/6ej2πf0t+4ejπ/3ej2π2f0t+ejπ/6ej2π3f0t,(f0=1000)

= 3ejπ/6ejω0t+4ejπ/3ej2ω0t+ejπ/6ej3ω0t(ω0=2000π)

X1=3ejπ/6;X2=4ejπ/3;X3=ejπ/6;Xn=0,n1,2,3


注意:Xc(t)為一複數週期訊號,故不滿足於後面將討論的指數傅立葉級數的對稱性(即Xn=Xn*) 。


範例4.14

經過半波整流後的週期弦波訊號可表示如下 : X(t)={Asinω0t0tT0/20T0/2tT0

...........................圖

此訊號仍為週期訊號(X(t)=X(t+T0)),試求其指數傅立葉級數

【解】 週期T0ω0=2πT0

           指數傅立葉級數 :
              X(t)=Xnejnω0t

Xn=1T00T0/2Asinω0tejnω0tdt

=A2jT00T0/2(ejω0tejω0t)ejnω0tdt

=A2jT0[0T0/2ejω0(1n)tdt-0T0/2ejω0(1+n)tdt]

=A4π[ej(1n)π11n+ej(1+n)π11+n] ,n±1


說明:

   ej(1±n)π=cos(1±n)π+jsin(1±n)π
   =(1)n


Xn={0nodd,n±1Aπ11n2neven?n=±1

X1=A2jT00T0/2(ejω0tejω0t)ejω0tdt


 =A2jT00T0/2(1ej2ω0t)dt
 =A4j=jA4

同法可得 :

X1=A4j=jA4

範例4.15

週期脈衝序列(periodic impulse train) δT0(t)的定義為

δT0(t)=m=δ(tmT0)試求其傅立葉級數


..............................圖


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



【解】 指數傅立葉級數


         δT0(t)=n=Dnejnω0tω0=2πT0


其中

Dn=1T0T0δT0(t)ejnω0tdt 取積分範圍由T02T02


=1T0T0/2T0/2δ(t)ejnω0tdt δT0(t)=δ(t)

=1T0T0/2T0/2δ(t)dt在此一範圍內

=1T0


         δT0(t)=1T0n=ejnω0t
                





       δT0(t)的雙邊振幅頻譜如下 :


.....................圖

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






    利用尤拉公式ejnω0t+ejnω0t=2cosnω0t

δT0(t)=1T0n=ejnω0t=1T0[1+2(cosω0t+cos2ω0t+...)]

=1T0[1+n=12cos(nω0t)]


δT0(t)的單邊頻譜如下 :


..................................圖


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



常見訊號的傅立葉級數

................................................

................................................

.....................圖.....................

................................................


指數傅立葉級數係數的對稱性

指數傅立葉級數

          X(t)=n=Xnejnω0t


     其係數 : 
         Xn=1T0T0x(t)ejnω0tdt


Xn=1T0T0x(t)ej(n)ω0tdt

   ={[1T0T0x(t)ej(n)ω0tdt]*}*      *表示共軛複數
   ={1T0T0x*(t)ej(+n)ω0tdt}*x(t) 為實數函數
   ={1T0T0x(t)ejnω0tdt}*
  =Xn*




     Xn 一般而言是複數,用極座標表示 :
            Xn=|Xn|ejXn    Xn表示Xn的相角


根據上述條件 :

             Xn=Xn*
        |Xn|=|Xn*|=|Xn|取共軛複數絕對值等
        Xn=Xn*=Xn共軛複數的相角差一個負號


結論 :當x(t) 為實數週期訊號的條件下,其指數傅立葉級數的係數Xn

             有底下特性:  
           (1)Xn的絕對值對 n 而言為偶對稱。
           (2)Xn的相位角對 n 而言為奇對稱。

考慮X(t)為實數偶函數 :

           指數傅立葉級數的係數Xn
    Xn=1T0T0x(t)ejnω0tdt
    =1T0T0/2T0/2x(t)[cosnω0tjsinnω0t]dt
    =1T0[T0/2T0/2x(t)cosnω0tdt-jT0/2T0/2x(t)sinnω0tdt



X(t)是偶函數,cosnω0t也是偶函數 X(t)是偶函數,sinnω0t是奇函數


X(t)cosnω0t仍為偶函數 X(t)sinnω0t為奇函數


         T0/2T0/2x(t)sinnω0tdt=0
=2T00T0/2x(t)cosnω0tdt




x(t)cosnω0t均為實數,故積分結果仍為實數。

        結論 : 當x(t)為實數偶函數時,其指數傅立葉級數的係數Xn 為實數 。

同法可證明,當x(t) 為實數奇函數時,x(t) 的指數傅立葉級數的係數Xn 為純虛數或 0。




x(t)為半波奇對稱(half-wave odd symmetrical ,即 :

         x(tT02)=x(t)Xn=0n=0,±2,±4,...

證明請同學自行練習


級數 係數 對稱性

1.三角正弦-餘弦


..........................圖




2.複數指數


..........................圖



傅立葉級數的變換(1)

振幅轉換 :

            已知  X(t)=n=Xnejnω0t
                           另一訊號 y(t)=αx(t)+β

                            則y(t)的傅立葉級數為
                                 y(t)=n=Y0ejnω0t
                                  =α[n=Xnejnω0t]+β
                                   =(αX0+β)+n=n0αXnejnω0t




{Y0=αX0+βYn=αXnn0



範例4.16

試求週期訊號y(t)的傅立葉級數 :


.........................圖


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

【解】仔細觀察上圖可以發現,此一圖形與標準的鋸齒波X(t)(如下圖)存

      在一關係 :
               y(t)=4Ax(t)+1=αx(t)+β


............................圖


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

經查表可知,標準鋸齒波x(t)的傅立葉級數為 :

          x(t)=A2+n=n0A2πnejπ2ejnω0t



y(t)=n=Ynejnω0t


     其中  Y0=αX0+β
 
                (4A)A2+1=1
          Yn=αXn=(4A)A2πnejπ/2
          2πnejπ/2n0

傅立葉級數的變換(2)

時間變換

   (a) y(t) 為 x(t) 延遲t0 時間,即y(t)=x(tt0)
      y(t)=x(tt0)=n=Xnejnω0(tt0)
                  =n=[Xnejnω0t0]ejnω0t
                  =n=Ynejnω0t
   Yn=Xnejnω0t0




傅立葉級數的變換

(b) y(t) 為 x(t) 左右翻轉,即y(t)=x(-t)

y(t)=x(-t)=n=Xnejnω0(t)

         =n=Xnejnω0t

         =n=Xn*ejnω0t
          =n=Ynejnω0t


Yn=Xn*


範例4.17

同範例4.16,改用時間變換搭配振幅變換求解。



【解】首先定義z(t)=x(-t)


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




z(t)=x(-t)=n=Xnejnω0(t)

         =n=Xnejnω0t

         =n=Xn*ejnω0t
          =n=Znejnω0t


Zn=Xn*={A2n=0A2πnejπ/2n0


觀察y(t)z(t)可得

    y(t)=4Az(t)3

{Y0=4AZ03=4AA23=1Yn=4AZn=4AA2πnejπ/2=2πnejπ/2n0



大師的風範

同學辛苦了,休息一下並瞻仰一下大師風範!


.........................圖



Jean Baptiste Joseph Fourier J. Willard Gibbs