matlab求傅里叶级数展开式_连续时间的傅里叶级数
生活随笔
收集整理的這篇文章主要介紹了
matlab求傅里叶级数展开式_连续时间的傅里叶级数
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
如果信號x(t)是周期信號,那么對于所有t,存在一個最小正數(shù)T,使得x(t)=x(t+T)其中T為這個周期信號的最小正周期。根據(jù)周期函數(shù)的周期性:x(t)=x(t+N*T)(N為整數(shù))
稱為這個信號的基波頻率
周期信號x(t)也可以用周期復(fù)指數(shù)信號表示則因?yàn)閤(t)是周期信號,那么所以那么周期信號T必須滿足倍數(shù)關(guān)系,也就是N*2π倍數(shù)關(guān)系(N為整數(shù)),而不能是非2π倍數(shù)關(guān)系。對于復(fù)指數(shù)信號的最小正周期為T,同時對應(yīng)的基波頻率w0,然而這個復(fù)指數(shù)信號可能含有其他諧波頻率的復(fù)指數(shù),而這些復(fù)指數(shù)是以T為基準(zhǔn)周期的,實(shí)際上比最小正周期T小。而諧波復(fù)指數(shù)信號的周期必須是最小正周期的T\K,(其中K為非零整數(shù)),復(fù)指數(shù)信號才可能含有諧波復(fù)指數(shù)信號。例如:這個周期函數(shù)最小正周期2π,但含有3次諧波的函數(shù),3次諧波周期為基波周期的2π\(zhòng)3,因此這個信號中不僅有基波信號,還有三次諧波信號。在一個最小正周期內(nèi),只有最小正周期是諧波周期的整數(shù)倍時,兩者波形的起點(diǎn)和終點(diǎn)才能重合。周期信號x(t)是以T為最小正周期的,與x(t)有關(guān)的周期諧波復(fù)指數(shù)信號集合可以表示如下:那么對于一個周期為T0的信號就可用周期復(fù)指數(shù)諧波信號疊加而成。上式的表達(dá)式就是傅里葉級數(shù),其中ak代表周期復(fù)指數(shù)諧波信號的振幅分量,又稱傅里葉系數(shù)或者頻譜系數(shù),w0代表信號的基波頻率,kw0代表諧波頻率,當(dāng)K=0時為常數(shù),當(dāng)K=1時為基波頻率,當(dāng)K=2時為2次諧波。以周期為T0的信號可以用復(fù)指數(shù)諧波信號疊加而成,那么應(yīng)該如何求解傅里葉系數(shù)呢?方程兩邊同時乘以-jnwt,兩邊同時積分。如果周期信號為實(shí)信號,傅里葉系數(shù)存在以下關(guān)系ak*=a-k(ak的共軛復(fù)數(shù)等于a-k)證明:如果x(t)是實(shí)周期函數(shù),那么有x(t)*=x(t)(共軛關(guān)系)令k=-k,目的是調(diào)換位置,使得復(fù)指數(shù)項(xiàng)從負(fù)到正。因?yàn)?比較兩式得:從而ak的共軛等于a-k。例題1:x(t)是以基波頻率為2π的周期信號,x(t)寫成以下形式。其中,a0=1,,a1=a-1=1/4,a2=a-2=1/2,a3=a-3=1/3,求其傅里葉展開式解:將(1-1)式中具有同一諧波頻率的諧波分量合在一起,重新組合得例題2:x(t)是sint,基波頻率為2π的實(shí)信號,求其傅里葉展開式。由歐拉公式得:則則k1=1/2j,k-1=-1/2j,ak=0(k≠±1)例題3,在電工學(xué)中,函數(shù)信號發(fā)生器可以產(chǎn)生方波,方波的幅值和相角可以調(diào)整,為了方便,設(shè)定原點(diǎn)為方波的對稱點(diǎn),幅值為1,最小正周期為2,波形以下展示。f(t)的周期為2,頻率為π,f(t)用傅里葉級數(shù)展開通過推導(dǎo)聯(lián)立兩個方程:用三角級數(shù)表示取k=17,自變量t為2π,Matlab代碼如下:>> clear all>> x=0:0.001:2*pi;>> y=4/pi*(sin(pi*x)+1/3*sin(3*pi*x)+1/5*sin(5*pi*x)+1/7*sin(7*pi*x)+1/9*sin(9*pi*x)+1/11*sin(11*pi*x)+1/13*sin(13*pi*x)+1/15*sin(15*pi*x)+1/17*sin(17*pi*x));>> plot(x,y,'b')輸出二維圖形如下:例題4,如圖所示是三角波周期函數(shù),最小正周期為2,求其傅里葉展開式。f(t)的周期為2,基波頻率為π,f(t)用傅里葉級數(shù)展開通過推導(dǎo)聯(lián)立兩個方程,用三角級數(shù)表示也可以用Matlab進(jìn)行傅里葉展開驗(yàn)證,取變量t為2π,k=11時,如下圖所示。代碼如下:>> clear all>> x=0:0.001:2*pi;>> y=0.5-4*(1/(pi*pi)*cos(pi*x)+1/(9*pi*pi)*cos(3*pi*x)+1/(25*pi*pi)*cos(5*pi*x)+1/(49*pi*pi)*cos(7*pi*x)+1/(81*pi*pi)*cos(9*pi*x)+1/(121*pi*pi)*cos(11*pi*x)+1/(169*pi*pi)*cos(13*pi*x));>> plot(x,y,'r')輸出二維圖形如下:參考內(nèi)容
[1]奧本海姆.信號與系統(tǒng)[M],第二版,北京:電子工業(yè)出版社,2013.
[2]王晶,翁國慶.張有兵.電力系統(tǒng)的MATLAB/SIMMULINK仿真與應(yīng)用[M],西安:西安電子科技大學(xué)出版社,2008.
總結(jié)
以上是生活随笔為你收集整理的matlab求傅里叶级数展开式_连续时间的傅里叶级数的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 计算机有必要报英语四级吗,我已工作了,现
- 下一篇: matlab如何实现降维,matlab怎