訊號與系統/傅立葉轉換的定理

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

Template:NoteTA

傅立葉轉換定理(1) —線性(linearity)

也稱作重疊定理(superposition) 已知: χ(t)X(f)(Xω(ω)) y(t)Y(f)(Yω(ω))

則對於任意實數或複數常數a,b

aχ(t)+by(t)aX(f)+bX(f)(aXω(ω)+bYω(ω))

證明 aχ(t)+by(t)=[aχ(t)+by(t)]ej2πftdt

=aχ(t)ej2πftdt+by(t)ej2πftdt

=aX(f)+bY(f)


範例5.11

試求χ(t)=Bcos2πf0t的傅立葉轉換。 解χ(t)=Bcos2πf0t=B2ej2πf0t+B2ej2πf0t(尤拉公式)

依傅立葉轉換的線性定理知


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

χ(t)=B2ej2πf0t+B2ej2πf0t

=B2ej2πf0t+B2ej2πf0t

=B2(δ(ff0)+δ(f+f0))

=B2(δ(ωω02π)+δ(ω+ω02π))

                                 由前面範例知
                                     ej2πf0t
                                       =δ(ff0)
                                       ω=2πf
                                       ω0=2πf0

πBδ(ωω0)+(ω+ω0)


範例5.12

試求單位步階函數u(t)的傅立葉轉換。 【解】


                                 © Charls L. Phillips, John M. Parr, Eve A. Riskin, Signals, Systems, and Transforms, 3rd ed., Pearson Education, 2003.
          (1)由上圖知:u(t)=12[1+sgn(t)]
            (2)已知sgn(t)1jπf(=2jω)


(3)故u(t)=12[1+sgn(t)]
        =12[δ(f)+1jπf]
      1j2πf+12δ(f)(=1jω+πδ(ω))
                                   © G. E. Carlson, Signal and Linear System Analysis, 2nd ed., John Wiley & Sons, 1998.
                                                          

傅立葉轉換定理(2)-時間比例調整(time scaling)

已知:χ(t)X(f)(Xω(ω)) 則:χ(at)1aX(fa)(1aXω(ωa)) 信號在時域的時間參數t做等比例放大或縮小a倍,此程序在頻域的頻率參數f 縮小或放大1a倍,同時振幅大小也縮小或放大\frac{1}{\mid1\mid}倍。訊號在時間軸壓縮(a>1)則其頻譜會擴張;反之,信號在時間擴張 (a<1)則其頻譜會壓縮。


傅立葉轉換定理(2)-時間比例調整(time scaling)

【證明】(1) a>0 χ(at)=χ(at)ej2πft,dt

at=τ

=1aχ(τ)ej2πfaτ,dτ =1aX(fa)

(2) a<0 χ(at)=1aX(fa) χ(at)=1aX(fa) #


範例5.13

試繪出χ(t)=rect(ta)y(t)=χ(t2)=rect(t4)之頻譜。

【解】 χ(t)=rect(ta)X(f)=2sinc(2f) y(t)=χ(t2)Y(f)=2X(2f)=4sinc(4f)


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

傅立葉轉換原理(3.)—時間反轉(time reversal)

已知:χ(t)X(f)<Xω(ω)

χ(t)X(f)<Xω(ω) 明顯的,此定理為時間比例調整定理的特例。

訊號在時域的時間參數 t 反轉造成在頻域的頻率參數 f 也反轉。

χ(t)為實數,則χ(t)=X(f)=X*(f)


範例5.14

試求χ(t)=eat,a>0之傅立葉轉換。 【解】(1)χ(t)=eat=eatu(t)+eatu(t) =χ1(t)+χ2(t)

其中χ1(t)=eatu(t) χ2(t)=eatu(t)

明顯的,χ2(t)=ea(t)u(t) =eatu(t)=χ1(t)

χ1(t)=χ2(t) =X2(f)


(2)由前面範例可知: X2(f)=X2(t)=eatu(t)=1a+j2πf (3)X(f)=X(t)=X1(t)+X2(t)

                                           線性定理
 =X1(t)+X2(t)
 =X2(f)+X2(f)
 =1aj2πf+1a+j2πf
 =2aa2+(2πf)2=2aa2+(ω)2


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

傅立葉轉換定理(4) —乘上tn

已知:χ(t)X(f)(Xω(ω))

tnχ(t)(1j2π)ndnX(f)dfn(jndXωn(ω)dωn) 【證明】證明n=1的情形:

                (1)依傅立葉轉換的公式:

X(f)=χ(t)ej2πftdt



(2)等號兩邊對f微分

dX(f)df=ddfχ(t)ej2πftdt

=χ(t)[ddfej2πft]dt

=χ(t)j2πtej2πft]dt

=(j2π)[tχ(t)]ej2πft]dt

=(j2π)tχ(t)

tχ(t)1j2πdX(f)df


範例5.15

試求下圖x(t)的傅立葉轉換。


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

【解】

      (1)χ(t)=(t) rect(f2)
     (2)已知{rect(f2)}=2sinc(2f)  (2sa(ω))
     (3)故 {χ(t)}= {t rect(f2)}
       =1j2πddf2sinc(2f)       ((j)ddω2sa(ω)))

=j2πfcos(2πf)sin(2πf)2π2f2 (j2ωcosωsinωω2)



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

傅立葉轉換定理(5) —共軛複數

已知:χ(t)X(f)(Xω(ω))

χ*(t)X*(f)(Xω*(ω))

                      明顯的,當x(t)為實數時,χ*(t)=χ(t)X*(f)=X(f)X(f)=X*(f0)

【證明】X(f)=χ(t)ej2πftdt

     X(f)=χ(t)ej2πftdt
    X*(f)=[χ(t)ej2πftdt]*
 =χ*(t)ej2πft]dt
 = {χ*(t)}

傅立葉轉換定理(6) —對偶性(duality)

已知:χ(t)X(f) (Xω(ω))

X(t)χ(f)

Xω(t)2πx(ω)

【證明】】χ(t)=χ(f)ej2πftdf

     χ(t)=X(t)ej2πftdf
                                            變數t與f互換
    χ*(f)=X(t)ej2πftdt
   = {X(t)}

範例5.16

試求的傅立葉轉換。 χb(t)=2Wsinc(2Wt) 【解】:(1)已知χa(t)=rect(tτ) Xa(f)=τsinc(τf)

(2)依據對偶性知:Xa(t)=τsinc(τt)χa(f)=rect(fτ)=rect(fτ)rect為偶函數

(3)將上式τ用2W代可得 χb(t)=Xa(t)τ=2W <math>Xb(f)=χa(f)τ=2W =2Wsinc(2Wt) =rect(f2W)


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

傅立葉轉換定理(7) —時移(time shift)

已知:

X(T)x(f) (Xω(ω))

x(tt0)X(f)ej2πft0(Xω(ω)ejωt0)

訊號在時間軸上平移(訊號超前或延遲)在頻域的效果相當於在原訊號的相位頻譜加上一個線性變化量2πft0,此變化量稱為傅立葉轉換X(f)的線性相位平移(linear phase shift) 。


{x(tt0)} =x(tt0)ej2πftdt =(τ=tt0,t=τ+t0dτ=dt)x(τ)ej2πf(τ+t0)dτ =ej2πft0x(τ)ej2πfτdτ =ej2πft0X(f)


範例5.17

已知 {δ(t)} = 1,試求δ(tt0) 的傅立葉轉換。

【解】

{δ(tt0)} =ej2πft0 {δ(t)} =ej2πft0*1 =ej2πft0


範例5.18

重複範例5.4,試求xb(t)=rect(t+12)rect(t12) 的傅立葉轉換。

【解】:

(1)已知rect(t)sinc(f)

(2)根據時移定理知:

rect(t+12)sinc(f)ej2πf(12) =sinc(f)ejπf

rect(t12)sinc(f)ej2πf(12) =sinc(f)ejπf

(3)故

{xb(t)} = {rect(t+12)} - {rect(t12)} =sinc(f)[ejπfejπf] =sinc(f)2jsin(πf) =sinc(f)2jπfsin(πf)πf =j2πfsinc2(f)


時移定理補充說明

由時移定理知,對於t0時間的延遲將會造成2πft0的相移,此一相移量與頻率 f 成正比。也就是說,針對t0時間的延遲,訊號的高頻成分會有較大的相位移,而低頻部分則相移較小。


傅立葉轉換定理(8) —頻移(frequency shift)

已知:

x(t)X(f)(Xω(ω)

x(t)ej2πf0tX(ff0)

x(t)ejω0tXω(ωω0))

訊號在時域乘上ㄧ複指數訊號ej2πf0t,在頻域的效果相當於訊號的頻譜在頻率軸上平移f0

頻移定理與時移定理互為對偶定理。

【證明】

{x(t)ej2πf0t} =[x(t)ej2πf0t]ej2πftdt ==[x(t)ej2π(ff0)t]dt =X(ff0)


範例5.19—調變原理(modulation)

已知x(t)X(f) ,試求(a)y1(t)=x(t)cos2πf0t(b)y2(t)=x(t)sin2πf0t 的傅立葉轉換。

【解】 (a)

(1)y1(t)=x(t)cos2πf0t =x(t)ej2πf0t+ej2πf0t2 =12[x(t)j2πf0t+x(t)j2πf0t]

(2) {y1(t)} =12[ {x(t)ej2πf0t}+ {x(t)ej2πf0t}] =12 {X(ff0)+X(f+f0)}



(b)

(1)y2(t)x(t)sin2πf0t =x(t)ej2πf0tej2πf0t2j =12j[x(t)j2πf0tx(t)j2πf0t]

(2) {y2(t)} =12j {X(ff0)X(f+f0)}


範例5.20

試求x1(t) =rect(t2)cos20πt的傅立葉轉換。

【解】 已知 rect(t2)cos20πt2sinc(2f)

X1(f) = {x1(t)} =sinc[2(f10)]+sinc[2(f+10)]


傅立葉轉換定理(9) —旋積定理(convolution)

已知:

x(t)X(f)(Xω(ω)

y(t)Y(f)(Yω(ω)

x(t)*y(t) =x(λ)y(tλ)dλX(f)Y(f)(Xω(ω)Yω(ω))

兩個訊號在時域做旋積運算相當於此二訊號在頻域相乘。

【證明】︰

(1)令z(t)=x(t)*y(t)=x(λ)y(tλ)dλZ(f) = {z(t)} =[x(λ)y(tλ)dλ]ej2πftdt交換積分順序Z(f)=x(λ)[y(tλ)ej2πftdt]dλ


【證明】

(2)根據時移定理:y(tλ)ej2πftdt=Y(f)ej2πfλ

(3)故Z(f)=x(λ)Y(f)ej2πfλdλ=[x(λ)ej2πfλdλ]Y(f)=X(f)Y(f)


範例5.21

試求方波x(t)=rect(tτ)自己做旋積運算後的傅立葉轉換。

【解】

(1)已知rect(tτ)τsinc(τf)

(2)令y(t)=x(t)*x(t)=rect(tτ)*rect(tτ)

根據旋積定理知:Y(f) = {y(t)} = {x(t)*x(t)} = X(f)X(f) = τ2sinc2(τf)


範例5.22—三角波的傅立葉轉換

定義:

x(t) = Λ(tτ){1|t|/τ,|t|<τ0,other

試求其傅立葉轉換。

【解】

(1)根據旋積運算公式及三角波的定義知:rect(tτ)*rect(tτ)=τΛ(tτ)

(2)又由範例5.21可知::rect(tτ)*rect(tτ)τ2sinc2(τf)

(3)故τΛ(tτ)τ2sinc2(τf)Λ(tτ)τsinc2(τf)


系統的轉換函數(transfer function)

一線性非時變系統對輸入訊號 x(t) 的響應可表示為

其中h(t) 為系統的單位脈衝響應

y(t)=x(t)*h(t)取傅立葉轉換可得:Y(f)=X(f)H(f)

其中H(f) = {h(t)}

H(f)稱為系統的轉換函數(transfer function) ;或稱為系統的頻率響應(frequency response) 。


傅立葉轉換定理(10) —乘積定理(multiplication)

已知:

x(t)X(f)(Xω(ω))

y(t)Y(f)(Yω(ω))

x(t)y(t)X(f)*Y(f)(12πXω(ω)*Yω(ω))

明顯的,乘積定理與旋積定理互為對偶定理。

【證明】︰

{x(t)y(t)} = x(t)y(t)ej2πftdt = [X(λ)ej2πλtdt]y(t)ej2πftdt = X(λ)[y(t)ej2π(fλ)tdt]dλ = X(λ)Y(fλ)dλ = X(f)*Y(f)


範例5.23

試求餘弦脈波函數(cosinusoidal pulse)

x(t)=Arect(tτ)cos2πf0t

的傅立葉轉換。

【解】

(1)已知Arect(tτ)Aτsinc(τf) cos(2πf0t)12 {δ(ff0)+δ(f+f0)}

(2)依據乘積定理可知:

X(f) = {x(t)} = {Arect(fractτ)} * {cos(2πf0t)} = Aτsinc(τf) * 12 {δ(ff0)+δ(f+f0)}(f(x)*δ(x±x0)=f(x±x0) = Aτ2 {sinc[τ(ff0)]+sinc[τ(f+f0)]}



傅立葉轉換定理(11) —時域微分(time-domain differentiation)

已知:

x(t)X(f)(Xω(ω))

若 x(t) 可微分,則

dnx(t)dtn(j2πf)nX(f)((jω)nXω(ω))

在時域對 x(t) 作微分相當於在頻域乘上(j2πf) 。由於(j2πf)的大小與頻率 f 成正比,故頻率越高的成分將會乘上越大的倍數。也就是說,微分的動作會對高頻有放大的效果。

時域微分定理與定理(4)乘上tn定理互為對偶定理。

【證明】

(1)已知x(t)=X(f)ej2πftdf

(2)等號兩邊分別對 t 微分:dx(t)dt=ddt[X(f)ej2πftdf] = X(f)ddt(ej2πft)df = X(f)(j2πf)ej2πftdf = [(j2πf)X(f)]ej2πftdf

【證明】

(3)故dx(t)dt(j2πf)X(f)

(4)重複上述步驟可得dnx(t)dtn(j2πf)nX(f)


傅立葉轉換定理(12) —時域積分(time-domain integration)

已知:x(t)X(f)(Xω(ω))

tx(λ)dλ1j2πfX(f)+X(0)δ(0)

(1jωXω(ω)+Xω(0)δ(ω))

在時域對 x(t) 作積分相當於在頻域除以(j2πf),故積分會衰減訊號的高頻部份。

公式中,X(0)=x(t) ,dt。 若X(0)=0,則公式可簡化為:tx(λ)dλ1j2πfX(f)(1jωXω(ω))

【證明】

(1)根據旋積運算的定義:x(t)*u(t)=x(λ)u(tλ) ,dλ=tx(λ)dλ

(2) {tx(λ)dλ } = {x(t)*u(t)}

【證明】︰

(3)根據旋積定理: {tx(λ)dλ } = {x(t)*u(t)} = {x(t)} {u(t)} = X(f)[1j2πf+12δ(f)] = 1j2πfX(f)+12X(f)δ(f) = 1j2πfX(f)+12X(0)δ(f)


範例5.24

試求下圖 x(t) 的傅立葉轉換。



【解法1】 (1)令y(t)=dx(t)dt,z(t)=dy(t)dtx(t)=ty(λ)dλ,y(t)=tz(τ)dτ


(2)由上圖知:z(t)=K[δ(t+b)δ(t+a)δ(ta)+δ(tb)]

其中(K=dx(A)ba)

Z(f) = {z(t)} = K[ej2πfbej2πfaej2πfa+ej2πfb] = 2K[cos2πfbcos2πfa]


(3)根據時域積分定理:

y(t)=tz(τ)dτ

Y(f)=1j2πfZ(f)+12Z(0)δ(f)=2K[cos2πfbcos2πfaj2πf]+0=2K[cos2πfbcos2πfaj2πf]


(4)因x(t)=ty(λ)dλ,再一次利用時域積分定理:

X(f)=1j2πfY(f)+12Y(0)δ(f)=2K[cos2πfbcos2πfaj2πf]+0(Y(0)=y(t)dt=0,cos2θ=12sin2θ) = 2K(12sin2πfb)(12sin2πfa)(2πf)2 = Ksin2πfbsin2πfa(πf)2 = K[ b2sinc2(fb)a2sinc2(fa) ]


【解法2】

(1)由圖可知:X(t)=BΛ(tb)(BA)Λ(ta)

其中B=Ab(ba)

(2)故X(f) = {x(t)} = B {Λ(tb)} - (BA) {Λ(ta} = Bbsinc2(bf)(BA)asinc2(af)