方波分解為多次正弦波之和

1、目錄1 技術(shù)要求 12 基本原理 13 建立模型描述 24 模塊功能分析或源程序代碼 25 調(diào)試過程及結(jié)論 55.1 調(diào)試過程 55.2 結(jié)論 66 心得體會 87 參考文獻 9TCOS2(TMx)t)dt = J；對于任何一個周期為T的周期信號而將此函數(shù)在區(qū)間(t1 ， t1+T)表示為(t)，都可以求出它在上述各函數(shù)中的分量，從方波分解為多次正弦波之和的設(shè)計1技術(shù)要求已知某一周期性方波，用 matlab仿真軟件演示諧波合成情況，討論參數(shù)對分解和合 成波形的影響。2基本原理根據(jù)三角傅里葉級數(shù):sin2(n(nt)dt (1)J： + Tcos (nifdt)cos：；+rsin(7H

2、3;Jt)sin (ntnt)dt 0,m 尹 n(2)(3)J； + Tsjn (m(z)t)cos (najt)dt = 0n為任何整數(shù)/(0口(1=三 + axcos (wt) 4- aos (2u)t) 4- “ + ucos + + ALsin (wt) + i2sin (+ bnsin (?ui>)十= y + iw _ 1a)lcos十 垢更陽!吒iwM)(4)這就是函數(shù)f(t)在上述區(qū)間的三角傅里葉級數(shù)表達式方波函數(shù)表達式：/(0 = 10 <t<T/2T/2 < t < 7(6)(8)先把這個函數(shù)展開為三角級數(shù)，為此就要求出分量系數(shù)a和bo因此，

3、該非周期性方波在區(qū)間(0, T)可以表示為f (r) = .4 sin 3十)+ sin (3o)t)十 sin (5t)十n.I.J1(9)3建立模型描述4模塊功能分析或源程序代碼figure(1)ts=0.0001;t=0:ts:4*pi; %t的取值圍為0到4*pi，間隔0.0001取一個點 f=1/ts;N=length(t);y仁square(0.32*pi*t);%產(chǎn)生一個方波plot(t,y1 ) ;title('產(chǎn)生一個方波');pausehold onfigure(2)subplot 321y2=4/pi*sin(t) ;%頻率為 1 (f=1/2*pi)的正

4、弦基波 % plot(t,y2);title('正弦基波')subplot 322y3=4/pi*(sin(3*t)/3) ;%三次諧波plot(t,y3);title('三次諧波')subplot 323y4=4/pi*(sin(5*t)/5) ;%五次諧波plot(t,y4) ;title('五次諧波')subplot 324y5=4/pi*(sin(7*t)/7) ;%七次諧波plot(t,y5);title('七次諧波')subplot 325y6=4/pi*(sin(9*t)/9) ;%九次諧波plot(t,y6);ti

5、tle('九次諧波');pausefigure(3)subplot 221y7=4/pi*(sin(t)+sin(3*t)/3) ;%將前兩次諧波疊加plot(t,y7) ;title('前兩次諧波疊加')subplot 222y8=4/pi*(sin(t)+sin(3*t)/3+sin(5*t)/5);% 將前三次諧波疊加plot(t,y8);title('前三次諧波疊加')subplot 223;%將前五次諧波疊加y9=4/pi*(sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/t) plot

6、(t,y9);title（'前五次諧波疊加'）%重新定義y,把各次波形數(shù)據(jù)存為一個三維數(shù)組y=zeros(10,max(size(t);x=zeros(size(t)for k=1:2:15x=x+4/pi*(si n( k*t)/k);y(k+1)/2,:)=x;end%將各波形疊合繪出pause,figure(4),plot(t,y(1:9,:),title('各波形疊合')5調(diào)試過程及結(jié)論5.1調(diào)試過程在調(diào)試過程中，編寫方波時由于當(dāng)時沒有注意方波掃描的精確度，將ts的數(shù)值選取的過于大了而導(dǎo)致方波的波形出現(xiàn)了大的失真，經(jīng)過反復(fù)幾次的調(diào)試后終于得出了一個比較

7、理想的方波圖形。由于在設(shè)計之初是要將方波與后來要形成的多次正弦波疊加后的波形放 到一起，所以對于方波周期的選取以及方波振幅的確定也嘗試了好多次。最終才將這兩個 數(shù)值確定下來。本事課程設(shè)計的目的是要證明方波能夠分解為多次正弦波之和，所以程序 的設(shè)計目的在于體現(xiàn)分解的過程以及最后n次諧波疊加后所形成的波形是否為方波。根據(jù)上述公式9得出的方波表達式，令® =1后分別畫出y=4/pi*si n(t)的正弦基波、 y=4/pi*(sin(3*t)/3)的三次諧波、y=4/pi*(sin(5*t)/5)的五次諧波、y=4/pi*(sin(7*t)/7)的七次諧 波以及y=4/pi*(sin(9*

8、t)/9)的九次諧波，為了能夠充分的顯示風(fēng)波分解為多次正弦波之和的 過程，又分別作出了前兩、三、五次諧波疊加的波形。最后為了能夠完美的證明方波分解 為多次正弦波之和，利用了一個for循環(huán)語句畫出了經(jīng)過前n次諧波疊加后的波形(其中 令n=15)，并且與之前程序最開始畫出的方波波形相互疊加。由于經(jīng)過三角級數(shù)傅里葉分解得到的方波表達式是一個非周期性函數(shù)，故只有當(dāng)無限次諧波疊加后才能形成正真的方波，所以本次的實驗是存在一定的誤差的，本想在程序中作出其誤差的圖形表示以及計算出不同次數(shù)的方波疊加形成的方波所造成的誤差，但由于知識以及能力有限，經(jīng)過多次嘗試無果后放棄了計算誤差這部分項目。 造成了一定的遺憾5

9、.2結(jié)論經(jīng)過此次的實驗，充分的證明了方波能夠分解為多次正弦波之和，疊加次數(shù)越多誤差越小越接近方波真實的波形。所得波形如下所示：產(chǎn)主一個方波圖5-1方波波形20150.500510150wnn51015220-2105IQ0510150 20 50 2圖5-2各次諧波波形-0 5-0.2三次諧液510七次諧液正弦基波51015方次諧液51015九;欠諧液前兩歡諧液軽加前三次諧波鬆加前五炭諧波蠡擔(dān)圖5-3疊加后的波形各液冊養(yǎng)合1 iIIIiniini246a101214圖5-4前n次疊加后的波形6心得體會通過此次的課程設(shè)計可以知道利用Matlab所得出的結(jié)果與利用理論知識通過計算得出的結(jié)果是一樣的

10、，方波，正弦波，以及正弦波疊加后的波形都能夠利用程序編寫出來的， 并且通過利用Matlab所得出的波形能夠很清晰的表達出各個波形的物理意義。Matlab是一款非常強大非常實用的軟件，利用Matlab我們不僅能夠證明許多數(shù)學(xué)上的問題，而且許多我們通過理論知識所無法計算的問題都能夠解決，比如說我們這次的課程設(shè)計，通過 人工計算我們能夠得出方波的公式，但是這個公式卻是一個非周期性的，我們無法準確的 算出方波的波形。但是利用 Matlab卻是能夠做到。通過利用 Matlab進行仿真，僅僅用一 個循環(huán)語句就能夠讓我們的波形無限的接近方波的真實波形，而且我們還能利用Matlab計算出不同次正弦波疊加后所產(chǎn)

11、生的誤差并且通過圖形模式表示出來以及我們還可以通過建立一個三維坐標(biāo)模型將方波分解為多次正弦波之和的趨勢表示出來，使我們更加直觀 的了解方波是如何分解為多次正弦波之和的。這次的課程設(shè)計讓我真的很難忘， 經(jīng)過兩個多星期的努力， 終于順利完成了課程設(shè)計。 開始做課程設(shè)計不知道從何入手，困難很多，經(jīng)過查閱資料，和同學(xué)討論，終于了解了許 多。在做課程設(shè)計的過程中，我學(xué)會了很多，最主要的是我對 Matlab 的運用更加熟練了， 同時對信號這門課的知識又弄懂了不少，增加了對信號課程的學(xué)習(xí)興趣。課程設(shè)計是每個 大學(xué)生必須面臨的一項綜合素質(zhì)的考驗，如果說在我們的學(xué)習(xí)階段是一個知識的積累過 程，那么現(xiàn)在的課程設(shè)計就是對過去所學(xué)的知識的綜合應(yīng)用，是對理論進行深化和重新認 識的實踐活動。在這期間，我們有艱辛的付出，當(dāng)然也有豐收的喜悅。首先，學(xué)習(xí)能力和 解決問題的信心都得到了提高。通過這次課程設(shè)計，我不僅對理論有了更深一步的認識， 還培養(yǎng)了自學(xué)能力和解決問題的能力，更重要的是，培養(yǎng)了克服困難的勇氣和信心。7 參考文