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

1、武漢理工大學(xué)專(zhuān)業(yè)課程設(shè)計(jì)3（信號(hào)與線性系統(tǒng)）課程設(shè)計(jì)說(shuō)明書(shū)目錄1 技術(shù)要求12 基本原理13 建立模型描述24 模塊功能分析或源程序代碼25 調(diào)試過(guò)程及結(jié)論45.1調(diào)試過(guò)程45.2 結(jié)論56 心得體會(huì)77參考文獻(xiàn)8方波分解為多次正弦波之和的設(shè)計(jì)1 技術(shù)要求已知某一周期性方波，用matlab仿真軟件演示諧波合成情況，討論參數(shù)對(duì)分解和合成波形的影響。2 基本原理根據(jù)三角傅里葉級(jí)數(shù)： t1t1+Tcos2ntdt=t1t1+Tsin2ntdt=T2 （1） t1t1+Tcosmtcosntdt=t1t1+Tsinmtsinntdt=0,mn （2） t1t1+Tsinmtcosntdt=0 m、n為

2、任何整數(shù) （3）對(duì)于任何一個(gè)周期為T(mén)的周期信號(hào)f（t），都可以求出它在上述各函數(shù)中的分量，從而將此函數(shù)在區(qū)間（t1，t1+T）內(nèi)表示為ft=a02+a1cost+a1cos2t+ancosnt+b1sint+b2sin2t+ +bnsinn+=a02+n=1ancosnt+bnsin(nt) （4）這就是函數(shù)f(t)在上述區(qū)間內(nèi)的三角傅里葉級(jí)數(shù)表達(dá)式。方波函數(shù)表達(dá)式：ft=1 0<t<T/2 ft T/2<t<T （5）先把這個(gè)函數(shù)展開(kāi)為三角級(jí)數(shù)，為此就要求出分量系數(shù)a和b。 a0=2T0Tf(t)dt=2T0T2dt-T2Tdt=0 （6） an=2T0Tf(t)co

3、s(nt)dt=2T0T2cosntdt-T2Tcosnt=0 （7） bn=2T0Tf(t)sin(nt)dt=2T0T2sinntdt-T2Tsinntdt=4n 當(dāng)n為奇數(shù)0 當(dāng)n為偶數(shù) （8）因此，該非周期性方波在區(qū)間（0，T）內(nèi)可以表示為 ft=4sint+13sin3t+15sin5t+ （9）3 建立模型描述分別作出各次諧波的波形根據(jù)公式（9）令=1利用Matlab畫(huà)出方波波形與方波比較將各波形圖疊加作圖分別作出不同次諧波疊加后的波形得出結(jié)論：方波可以分解為多次正弦波之和4 模塊功能分析或源程序代碼figure(1)ts=0.0001;t=0:ts:4*pi;%t的取值范圍為0到

4、4*pi，間隔0.0001取一個(gè)點(diǎn)f=1/ts;N=length(t);y1=square(0.32*pi*t);%產(chǎn)生一個(gè)方波plot(t,y1);title('產(chǎn)生一個(gè)方波');pausehold onfigure(2)subplot 321y2=4/pi*sin(t);%頻率為1（f=1/2*pi）的正弦基波%plot(t,y2);title('正弦基波')subplot 322y3=4/pi*(sin(3*t)/3);%三次諧波plot(t,y3);title('三次諧波')subplot 323y4=4/pi*(sin(5*t)/5);

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);title('九次諧波');pausefigure(3)subplot 221y7=4/pi*(sin(t)+sin(3*t)/3);%將前兩次諧波疊加plot(t,y7);title('前兩次諧波疊加')subplot 222y8=4/pi*(sin(t)+

6、sin(3*t)/3+sin(5*t)/5);%將前三次諧波疊加plot(t,y8);title('前三次諧波疊加')subplot 223y9=4/pi*(sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/t);%將前五次諧波疊加plot(t,y9);title('前五次諧波疊加')%重新定義y，把各次波形數(shù)據(jù)存為一個(gè)三維數(shù)組y=zeros(10,max(size(t);x=zeros(size(t);for k=1:2:15 x=x+4/pi*(sin(k*t)/k);y(k+1)/2,:)=x;end%將各波

7、形疊合繪出pause,figure(4),plot(t,y(1:9,:),title('各波形疊合')5 調(diào)試過(guò)程及結(jié)論5.1調(diào)試過(guò)程在調(diào)試過(guò)程中，編寫(xiě)方波時(shí)由于當(dāng)時(shí)沒(méi)有注意方波掃描的精確度，將ts的數(shù)值選取的過(guò)于大了而導(dǎo)致方波的波形出現(xiàn)了大的失真，經(jīng)過(guò)反復(fù)幾次的調(diào)試后終于得出了一個(gè)比較理想的方波圖形。由于在設(shè)計(jì)之初是要將方波與后來(lái)要形成的多次正弦波疊加后的波形放到一起，所以對(duì)于方波周期的選取以及方波振幅的確定也嘗試了好多次。最終才將這兩個(gè)數(shù)值確定下來(lái)。本事課程設(shè)計(jì)的目的是要證明方波能夠分解為多次正弦波之和，所以程序的設(shè)計(jì)目的在于體現(xiàn)分解的過(guò)程以及最后n次諧波疊加后所形成的波形

8、是否為方波。根據(jù)上述公式9得出的方波表達(dá)式，令=1后分別畫(huà)出y=4/pi*sin(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*t)/9)的九次諧波，為了能夠充分的顯示風(fēng)波分解為多次正弦波之和的過(guò)程，又分別作出了前兩、三、五次諧波疊加的波形。最后為了能夠完美的證明方波分解為多次正弦波之和，利用了一個(gè)for循環(huán)語(yǔ)句畫(huà)出了經(jīng)過(guò)前n次諧波疊加后的波形（其中令n=15），并且與之前程序最開(kāi)始畫(huà)出的方波波形相互疊加。由于經(jīng)過(guò)三角級(jí)數(shù)傅里葉分解得到的方波表達(dá)

9、式是一個(gè)非周期性函數(shù)，故只有當(dāng)無(wú)限次諧波疊加后才能形成正真的方波，所以本次的實(shí)驗(yàn)是存在一定的誤差的，本想在程序中作出其誤差的圖形表示以及計(jì)算出不同次數(shù)的方波疊加形成的方波所造成的誤差，但由于知識(shí)以及能力有限，經(jīng)過(guò)多次嘗試無(wú)果后放棄了計(jì)算誤差這部分項(xiàng)目。造成了一定的遺憾。5.2 結(jié)論經(jīng)過(guò)此次的實(shí)驗(yàn)，充分的證明了方波能夠分解為多次正弦波之和，疊加次數(shù)越多誤差越小越接近方波真實(shí)的波形。所得波形如下所示：圖5-1 方波波形圖5-2 各次諧波波形圖5-3 疊加后的波形圖5-4 前n次疊加后的波形6 心得體會(huì)通過(guò)此次的課程設(shè)計(jì)可以知道利用Matlab所得出的結(jié)果與利用理論知識(shí)通過(guò)計(jì)算得出的結(jié)果是一樣的，方

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

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