基于MATLAB的諧波分析

目錄（1）Matlab6.5以上版本軟件； 1緒論 11公式分析及計算 21.1傅里葉變換的原理 21.2傅里葉變換的證明 31.3周期信號的分解……………..31.4方波的分解…………………...52建模與仿真 72.1建模 72.2仿真 83仿真結(jié)果分析 104小結(jié) 11參考文獻 13武漢理工大學《基礎技能強化訓練》課程設計說明書

緒論方波是一種非正弦曲線的波形，通常會于電子和訊號處理時出現(xiàn)。由于一般電子零件只有“高(1)”和“低(0)”兩個值，方波就自然產(chǎn)生，所以理想方波只有“高”和“低”這兩個值。電流的波形為矩形的電流即為方波電流。不論時間軸上下是不是對稱的，只要是矩形就可叫方波，必要時，可加“對稱”，“不對稱”加以說明。而在現(xiàn)實世界，方波只有有限的帶寬。因為方波可以快速從一個值轉(zhuǎn)至另一個(即0→1或1→0)，所以方波就用作\o"時鐘訊號"時鐘訊號來準確地觸發(fā)\o"同步電路"同步電路。但是如果用\o"頻率"頻率\o"定義域"定義域來表示方波，就會出然一連串的諧波。所以方波可用相應頻率的基波及其奇次諧波合成。在電路信號系統(tǒng)的分析中，隨著電路規(guī)模的加大，微分方程的階數(shù)以及聯(lián)立后所得的方程的個數(shù)也隨之加大，加上電器元件的多樣化，這些都給解題運算分析電路系統(tǒng)帶來了一定的困難。傳統(tǒng)的計算機編程語言，如FORTRAN、C語言等，雖然都可以幫助計算，但在處理高階微分方程和大規(guī)模的聯(lián)立方程組的問題時大量的時間和精力都花在矩陣處理和圖形的生成分析等繁瑣易錯的細節(jié)上。而MATLAB憑借其強大的矩陣運算能力、簡便的繪圖功能、可視化的仿真環(huán)境以及豐富的算法工具箱，已成為科研和工程技術(shù)人員的有力開發(fā)工具。利用MATLAB不僅可以簡單快速的求解電路方程，同時，MAYLAB提供的Simulink工具還可以直接建立電路模型，隨意改變模型的參數(shù)，并且還可以快速得到仿真擬結(jié)果，進一步省去了編程的步驟。MATLAB具有數(shù)值計算功能；圖形處理及可視化功能；可視化建模及動態(tài)仿真功能等等。它給用戶帶來的是最直觀，最簡潔的程序開發(fā)環(huán)境。它的語言簡潔緊湊，使用方便靈活，程序書寫形式自由，利用起豐富的庫函數(shù)避開繁雜的子程序編程任務，壓縮了一切不必要的編程工作。同時，它的運算符也很豐富。由于MATLAB是用C語言編寫的，MATLAB提供了和C語言幾乎一樣多的運算符，靈活使用MATLAB的運算符將使程序變得極為簡短。它的程序的可移植性，基本上不做修改就可以在各種型號的計算機和操作系統(tǒng)上運行。本文應用MATLAB來驗證定理：方波可用相應頻率的基波及其奇次諧波合成。1公式分析及計算1.1傅里葉變換的原理任何具有性質(zhì)周期為T的波函數(shù)都可以表示為三角函數(shù)所構(gòu)成的級數(shù)之和，即：（1）其中：t為時間，為角頻率。=（T為周期），第一項為直流分量。圖1方波所謂周期性函數(shù)的傅里葉變換（Fouriertransform）就是將周期性函數(shù)張凱成直流分量，基波和所有n次諧波的疊加。圖1所示的方波可以寫成函數(shù)形式：h(0≤t＜)-h(-≤t＜0)在這里，h為常數(shù)2。很明顯，此方波為奇函數(shù)，并且它沒有常數(shù)項，同時，它是一個周期為T的函數(shù)，所以我們可以用傅里葉級數(shù)來表示這個函數(shù)。（2）我們把它展開，可以得到：（3）1.2傅里葉變換的證明下面，我們要從數(shù)學角度來證明為什么公式（3）能成立。由于這是一個奇函數(shù),常數(shù)項可以用積分函數(shù)計算出來：所以其常數(shù)項不存在，即=0，下面開始計算與:由上式可知，當n=2，4，6，時，；當n=1，3，5，時，。然后我們將與都帶入公式（1），就可以得到公式（3）：如果我們?nèi)〉捻椩蕉?，就會越逼近原本的方波函?shù)。1.3周期信號的分解周期信號是定義在（-，）區(qū)間，每隔一定的時間T，按相同規(guī)律重復變化的信號，它可表示為f(t)=f(t+mT)式中m為任意整數(shù)。時間T稱為該信號的重復周期，簡稱為周期。需要指出的是，只有當周期信號滿足狄里赫利條件時，才能展開為傅里葉級數(shù)。狄里赫利條件是：1）函數(shù)在任意有限區(qū)間內(nèi)連續(xù)，或只有有限個第一類間斷點（當t從左或右趨于這個間斷點時，函數(shù)有有限的左極限和右極限）2）在一周期內(nèi)，函數(shù)有有限個極大值或極小值。設有周期信號f(t),它的周期是T，角頻率=2F=，它可分解為(2-1)上式中系數(shù)a,b稱為傅里葉系數(shù)。為簡便，式積分區(qū)間（t,t+T）取為（-，）或（0，T）?？紤]到正、余弦函數(shù)的正交條件，可得傅里葉系數(shù)a=(t)cos(nt)dt,n=0,1,2,……(2-2)b=(t)sin(nt)dt,n=0,1,2,……(2-3)式中T為函數(shù)f(t)的周期，=為角頻率，由上述兩式，傅里葉系數(shù)a和b都是n的函數(shù)，其中a是n的偶函數(shù)，即a=a；而b是n的奇函數(shù)，既有b=-b將式（2-1）中同頻率項合并，可寫成如下形式f(t)=……=(2-4)式中A=，n=1,2,……=-arctan()如將式（2-4）的形式化為（2-1）的形式，他們系數(shù)之間的關(guān)系為aa,n=1,2,……b,式（2-4）表明，任何滿足狄里赫利條件的周期函數(shù)可分解為直流和許多余弦（或正弦）分量。其中第一項是常數(shù)項，它是周期信號中所包涵的直流分量；式中第二項A稱為基波或一次諧波，它的角頻率與原周期信號相同，A是基波振幅，是基波初相角；式中第三項A稱為二次諧波，它的頻率是基波頻率的兩倍，A是二次諧波振幅，是其初相角。以此類推，還有三次、四次、……諧波。一般而言，A稱為n次諧波，A是n次諧波的振幅，是其初相角。式（2-4）表明，周期函數(shù)可以分解為各諧波分量。1.4方波的分解設方波信號f(t)的周期為T，寬度為，將其展開為傅里葉級數(shù)由式（2-2）和（2-3）可得=+考慮到，可得將它們代入到式（2-1），得到信號的傅里葉級數(shù)展開式為，n=1，3，5，……它只含一、三、五…奇次諧波分量。下圖中畫出了一個周期的方波組成情況，由圖可見，當它包含的諧波分量愈多時，波形就愈接近原來的方波信號（圖中虛線所示），其均方誤差愈小，還可以看出，頻率較低的諧波，其振幅較大，他們組成方波的主體，而頻率較高的高次諧波振幅較小，它們主要影響波形的細節(jié)，波形中所包含的高次諧波愈多，波形的邊緣愈陡峭。（a）基波（b）基波“+”三次諧波（c）基波+三次諧波+五次諧波（d）基波+三次諧波+五次諧波+七次諧波由圖中還可以看出，合成波形所包含的諧波，除間斷點附近外，它愈接近于原方波信號。在間斷點附近，隨著所含諧波次數(shù)的增高，合成波形的尖峰愈靠近間斷點，但尖峰幅度并未明顯減小?？梢宰C明，即使合成波形所含諧波次數(shù)n時，在間斷點處仍有約9%的偏差，這種現(xiàn)象稱為吉布斯現(xiàn)象。在傅里葉級數(shù)的項數(shù)取得很大時，間斷點處尖峰下的面積非常小以致趨近于零，因而在均方的意義上合成波形同原方波的真值之間沒有區(qū)別。2建模與仿真2.1建模上文中，我們證明了一個以原點為奇對稱中心的方波可以用奇次正弦波的疊加來逼近?？梢院喕癁椋海?）如果我們能驗證是方波，那么我們可以得出也是方波，只是的方波的幅值是幅值的倍。已知方波的寬度為，周期為，我們可以用MATLAB程序來檢驗這種逼近的程度與特征。程序如下：t=0:.01:2*pi;%設定一個時間數(shù)組,有101個點

y=sin(t);plot(t,y),figure(gcf),pause%頻率為w=1(f=1/2π)的正弦基波

y=sin(t)+sin(3*t)/3;plot(t,y),pause%疊加三次諧波

y=sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/9;

plot(t,y)%用1,3,5,7,9次諧波疊加%為了繪制三維曲面,要把各次波形數(shù)據(jù)存為一個三維數(shù)組,因此必須重新定義y,重編程。y=zeros(10,max(size(t)));x=zeros(size(t));

fork=1:2:19

x=x+sin(k*t)/k;y((k+1)/2,:)=x;

end

pause,figure(1),plot(t,y(1:9,:)),grid%將各波形迭合繪出

line([0,pi+0.5],[pi/4,pi/4])%加上方波幅度線及標注

text(pi+0.5,pi/4,'pi/4')

halft=ceil(length(t)/2);

pause,figure(2),

mesh(t(1:halft),[1:10],y(:,1:halft)),%只用正半周波形

pause,clc2.2仿真仿真結(jié)果如圖2至圖6所示，其中所有的橫坐標為時間t，縱坐標為幅值y：圖2頻率為的正弦基波圖3疊加三次諧波圖4用1,3,5,7,9次諧波疊加圖5諧波合成的二維曲線圖6諧波合成的三維曲面從圖6可以看出，所取的階次越高，合成結(jié)果就越接近與方波。至于總是消除不了邊緣上的尖峰，這個就是吉布斯效應。3仿真結(jié)果分析根據(jù)周期信號的傅里葉展開式可知，任何方波周期性信號，只要滿足公式（1）條件，都可以分解為一個直流分量和由基波及奇次諧波分量的疊加。

通過數(shù)學工具MATLAB中編程可得，在正弦波（圖2）的基礎上實現(xiàn)三次諧波的疊加出現(xiàn)波形如圖3所示。依次對正弦波信號進行奇次諧波疊加過程（圖4），最終結(jié)果我們看到了圖5中的圖像，波形穩(wěn)定在了的周期性方波。圖六為圖五的的三維圖像分布圖。通過圖像的特性我們證明了：當奇次諧波的次數(shù)越高時，方波可由正弦波形產(chǎn)生諧波構(gòu)成。同理可知，一定條件下，周期性方波信號可由相應頻率的基波和奇次諧波合成。

4小結(jié)通過這次基礎強化訓練讓我認識了MATLAB這個有著強大功能的軟件。經(jīng)過這一段時間的學習，我對MATLAB有了較為系統(tǒng)的了解。它的應用領域相當廣泛：微積分、矩陣代數(shù)、應用數(shù)學、物理、力學、信號與系統(tǒng)、電子線路、電機學、機械振動、自動控制和通信技術(shù)等。它只需幾筆簡單的程序，就可以完成繁瑣的計算。它的擴展性強，在學好其基礎部分之后，還有幾十種工具箱可以用于各類科研需要，這可以縮短學習和實踐工作的距離。拿到老師給的題目后，我針對題目進行了分析。方波可用相應頻率的基波及其奇次諧波合成，并用MATLAB實現(xiàn)仿真。仔細思考了課題，發(fā)現(xiàn)相應頻率的基波及其奇次諧波的合成，就是傅里葉的展開式。顯然，這要求我要深入研究傅里葉級數(shù)展開分析的理論知識，利用MATLAB強大的圖形處理功能，符號運算功能以及數(shù)值計算功能，實現(xiàn)相應頻率的基波及其奇次諧波合成方波。在這次訓練中，我用MATLAB實現(xiàn)方波信號的傅里葉級數(shù)分解與綜合，繪出包含不同諧波次數(shù)的合成波形，觀察合成波形與原矩形波形之間的關(guān)系及吉布斯現(xiàn)象，證明課題的正確性。于是我使用MATLAB，輸入相關(guān)的程序，并且運行它，檢查它的正確性，并且修正。在程序的調(diào)試過程中遇到很多的問題。有的程序我根本就不知道它錯在哪里，也不知道怎么去改正。只有不斷的翻書，不斷的修改程序，不斷的調(diào)試。終于繪制出了所需要的圖像。在寫相關(guān)源程序的時候，我還收索了大量的網(wǎng)站，在網(wǎng)上收索了很多關(guān)于MATLAB的資料。在這個過程中我發(fā)現(xiàn)網(wǎng)上有很多有用的知識。以后應該多注意，充分合理的利用網(wǎng)絡，通過網(wǎng)絡來學習東西。這次基礎強化訓練使我明白了在知識的領域里我還有很多很多的不足，并且再一次的深深的體會到理論和實踐之間還有很到的差別。在以后的學習中應該多多的注意實踐知識的訓練和積累。在以后的學習生活中要不斷的開拓自己的動手能力，不斷的訓練自己的動手能力。這次基礎強化訓練讓我學習了了數(shù)字信號系統(tǒng)處理里的相關(guān)知識，也復習了高數(shù)中的相關(guān)知識，并且我還對mathtype數(shù)學公式編輯器有了一定的了解，并且會用它編輯公式。21世紀人類的知識正以指數(shù)規(guī)律飛速增長，我們已經(jīng)借助計算機輔助設計和制造，設計業(yè)和制造業(yè)已經(jīng)大大的提高了效率，創(chuàng)造了空前的財富。但是在我們的教和學的過程中，我們還在做繁瑣重復的計算這種機械勞動，如何從其中解放出來？就成了提高學習效率的關(guān)鍵環(huán)節(jié)?，F(xiàn)在我們雖然在大學中學習計算機課程，也只是為了以后的就業(yè)需要，很少對學生在學校學習有直接幫助。目前我們還處在“計算器水平”。正因為這樣，我們沒有足夠時間用于概念的思考，知識的擴充和思維的創(chuàng)新。就像學習信號與系統(tǒng)一樣，