連續(xù)時間信號傅里葉級數(shù)分析及matlab實(shí)現(xiàn)

1、課程設(shè)計任務(wù)書學(xué)生姓名：專業(yè)班級：指導(dǎo)教師：工作單位：題目：連續(xù)時間信號傅里葉級數(shù)分析及MATLAB實(shí)現(xiàn)初始條件：MATLAB6.5要求完成的主要任務(wù)：深入研究連續(xù)時間信號傅里葉級數(shù)分析的理論知識，利用MATLAB強(qiáng)大的圖形處理功能，符號運(yùn)算功能以及數(shù)值計算功能，實(shí)現(xiàn)連續(xù)時間周期信號頻域分析的仿真波形。1.用MATLAB實(shí)現(xiàn)周期信號的傅里葉級數(shù)分解與綜合。2用MATLAB實(shí)現(xiàn)周期信號的單邊頻譜及雙邊頻譜。3用MATLAB實(shí)現(xiàn)典型周期信號的頻譜。4.撰寫MATLAB應(yīng)用實(shí)踐課程設(shè)計說明書。時間安排：學(xué)習(xí)MATLAB語言的概況第1天學(xué)習(xí)MATLAB語言的基本知識第2、3天學(xué)習(xí)MATLAB語言的應(yīng)用

2、環(huán)境，調(diào)試命令，繪圖能力第4、5天課程設(shè)計第6-9天答辯第10天指導(dǎo)教師簽名：年月日系主任（或責(zé)任教師）簽名:摘要I.ABSTRACT錯誤！未定義書簽。緒論I.1 MATLAB簡介II1.1 MATLAB語言功能II1.2 MATLAB語言特點(diǎn)II2連續(xù)時間周期信號的傅里葉級數(shù)JJJ2.1 連續(xù)時間周期信號的分解III2.1.1 三角形式的傅里葉級數(shù)III2.1.2 指數(shù)形式的傅里葉級數(shù)IV2.2 連續(xù)時間周期信號的傅里葉綜合V2.3 吉布斯現(xiàn)象V.3連續(xù)時間周期信號的頻譜分析VII3.1 單邊與雙邊頻譜關(guān)系VII3.2 以單邊幅度頻譜為例，研究脈沖寬度與頻譜的關(guān)系IX3.3 以單邊幅度頻譜為

3、例，研究脈沖周期與頻譜的關(guān)系X4典型周期脈沖的頻譜XI4.1 周期方波脈沖頻譜的MATLAB實(shí)現(xiàn)XI4.1.1 周期方波脈沖雙邊頻譜的MATLAB實(shí)現(xiàn)XII4.1.2 周期方波脈沖單邊頻譜的MATLAB實(shí)現(xiàn)XIV4.2 周期三角波脈沖頻譜的MATLAB實(shí)現(xiàn)XV4.2.1 周期三角波雙邊頻譜的MATLAB實(shí)現(xiàn)XVII4.2.2 周期三角波單邊頻譜的MATLAB實(shí)現(xiàn)XVIII5小結(jié)即心得體會XX致謝XXI參考文獻(xiàn)錯誤！未定義書簽。附錄XXIII摘要MATLAB目前已發(fā)展成為由MATLAB語言、MATLAB工作環(huán)境、MATLAB圖形處理系統(tǒng)、MATLAB數(shù)學(xué)函數(shù)庫和MATLAB應(yīng)用程序接口五大部分組

4、成的集數(shù)值計算、圖形處理、程序開發(fā)為一體的功能強(qiáng)大的系統(tǒng)。本次課程設(shè)計則在深入研究連續(xù)時間信號傅里葉級數(shù)分析理論知識的基礎(chǔ)上，利用MATLAB強(qiáng)大的圖形處理功能、符號運(yùn)算功能以及數(shù)值計算功能，通過MATLAB編程進(jìn)行圖形功能仿真，從而實(shí)現(xiàn)連續(xù)時間周期信號頻域分析的仿真波形，包括以下內(nèi)容：用MATLAB實(shí)現(xiàn)周期信號的傅里葉級數(shù)分解與綜合的波形；用MATLAB實(shí)現(xiàn)周期信號的單邊頻譜及雙邊頻譜的波形與分析；用MATLAB實(shí)現(xiàn)典型周期信號的頻譜的波形。關(guān)鍵詞：MATLAB;圖形處理；傅里葉級數(shù)；周期信號；頻譜緒論在科學(xué)技術(shù)飛速發(fā)展的今天，計算機(jī)正逐步將科技人員從繁重的計算工作中解脫出來。在進(jìn)行科學(xué)研究

5、與工程應(yīng)用中，往往需要大量的科學(xué)計算，一些科技人員曾經(jīng)嘗試使用傳統(tǒng)的高級語言Basic、Fortran及C語言編寫程序，以減輕工作量。但編制程序需要掌握高級語言的語法，還要對各種算法進(jìn)行了解，這對大多數(shù)科技人員來說是不大現(xiàn)實(shí)的，而且也是沒有沒有必要的。MATLAB正是在這一應(yīng)用要求背景下產(chǎn)生的數(shù)學(xué)類科技應(yīng)用軟件。它具有的頂尖的數(shù)值計算功能、強(qiáng)大的圖形可視化功能及簡潔易學(xué)的“科學(xué)便捷式”工作環(huán)境和編程語言，從根本上滿足了科技人員對工程數(shù)學(xué)計算的要求，并將科技人員從繁重的數(shù)學(xué)運(yùn)算中解放出來，因而越來越受到廣大科技工作者的普遍歡迎1。MATLAB是matrixftlaboratory前三個字母的縮寫

6、，意思是“矩陣實(shí)驗(yàn)室”，是MathWorks公司推出的數(shù)學(xué)類科技應(yīng)用軟件。其Dos版本(MATLAB1.0)發(fā)行于1984年，現(xiàn)已推出了Windows版本(MATLAB5.3)。經(jīng)過十多年的不斷發(fā)展與完善，MATLAB已發(fā)展成為由MATLAB語言、MATLAB工作環(huán)境、MATLAB圖形處理系統(tǒng)、MATLAB數(shù)學(xué)函數(shù)庫和MATLAB應(yīng)用程序接口五大部分組成的集數(shù)值計算、圖形處理、程序開發(fā)為一體的功能強(qiáng)大的系統(tǒng)。MATLAB由“主包”和三十多個擴(kuò)展功能和應(yīng)用學(xué)科性的工具箱(Toolboxs)組成。目前，MATLAB已經(jīng)成為國際上最流行的電子仿真計算機(jī)輔助設(shè)計的軟件工具，現(xiàn)在的MATLAB已經(jīng)不僅僅

7、是一個“矩陣實(shí)驗(yàn)室(MatrixLaboratory)",它已經(jīng)成為一種實(shí)用的、全新的計算機(jī)高級語言。正是由于MATLAB在數(shù)值計算及符號計算等方面的強(qiáng)大功能，使MATLAB一路領(lǐng)先，成為數(shù)學(xué)類科技應(yīng)用軟件中的佼佼者。目前，MATLAB已成為國際上公認(rèn)的最優(yōu)秀的科技應(yīng)用軟件。MATLAB的上述特點(diǎn)，使它深受工程技術(shù)人員及科技專家的歡迎，并很快成為應(yīng)用學(xué)科計算機(jī)輔助分析、設(shè)計、仿真、教學(xué)等領(lǐng)域不可缺少的基礎(chǔ)軟件。1 MATLAB簡介1.1 MATLAB語言功能MATLAB是一個高精度的科學(xué)計算語言，它將計算、可視化編程結(jié)合在一個容易使用的環(huán)境中，在這個環(huán)境中，用戶可以把提出的問題和解決

8、問題的辦法用熟悉的數(shù)學(xué)符號表示出來，它的典型使用包括：(1)數(shù)學(xué)和計算；(2)運(yùn)算法則；(3)建模、仿真；(4)數(shù)值分析、研究和可視化；(5)科學(xué)的工程圖形；(6)應(yīng)用程序開發(fā)，包括創(chuàng)建圖形用戶接口。1.2 MATLAB語言特點(diǎn)MATLAB是一個交互式系統(tǒng)，他的基本數(shù)據(jù)單元是數(shù)組，這個數(shù)組不要求固定的大小，因此可以讓用戶解決許多技術(shù)上的問題，特別是那些包含矩陣和矢量運(yùn)算的問題。MATLAB的指令表達(dá)與數(shù)學(xué)、工程中常用的習(xí)慣形式相似，與C、Fortran等高級語言相比，它的語法規(guī)則更簡單、表達(dá)更符合工程習(xí)慣，正因?yàn)槿绱耍藗冇肕ATLAB語言編寫程序就猶如在便箋上書寫公式和求解，因而MATLAB

9、被稱為“便箋式”的科學(xué)工程語言。MATLAB的最重要特征使他擁有解決特定應(yīng)用問題的程序組，也就是TOOLBOX（工具箱），如信號處理工具箱，控制系統(tǒng)工具箱、神經(jīng)網(wǎng)絡(luò)工具箱、模糊邏輯工具箱、通信工具箱和數(shù)據(jù)采集工具箱等許多專用工具箱，對大多數(shù)用戶來說，要想靈活、高效地運(yùn)用這些工具箱，通常都需要學(xué)習(xí)相應(yīng)的專業(yè)知識。止匕外，開放性也許是MATLA最重要和最受歡迎的特點(diǎn)之一。除內(nèi)部函數(shù)外，所有的MATLAB主要文件和各工具箱文件都是可讀的、可改的源文件，因?yàn)楣ぞ呦鋵?shí)際上是有一組復(fù)雜的MATLAB函數(shù)（M文件）組成，它擴(kuò)展了MATLAB的功能，用以解決待定的問題，因此用戶可以通過對源文件進(jìn)行修改和加入自

10、己編寫的文件去構(gòu)建新的專用工具箱。2連續(xù)時間周期信號的傅里葉級數(shù)頻域分析法即傅里葉分析法，它是變換域分析法的基石。其中，傅里葉級數(shù)是變換域分析法的理論基礎(chǔ)，傅里葉變換作為頻域分析法的重要數(shù)學(xué)工具，具有明確的物理意義,在不同的領(lǐng)域得到廣泛的應(yīng)用2.1 連續(xù)時間周期信號的分解以高等數(shù)學(xué)的知識，任何周期為T的周期函數(shù)f（t）,在滿足狄里赫利條件時，則該周期信號可以展開成傅里葉級數(shù)。傅里葉級數(shù)有三角形式和指數(shù)形式兩種2o2.1.1 三角形式的傅里葉級數(shù)三角形式的傅里葉級數(shù)為:f（t）=a0a1cos（'"）a2cosc,t）a3cos（"），3冶（',t）b2s（&

11、#39;什）0co（2-1）二a一二ancos（n't）bnsin（n'4）n=1,2,32n1n1式中系數(shù)an、bn稱為傅里葉系數(shù)，可由下式求得。121a0TftdtIr22an=-TftcosnL"dt22T氣bn=22Tftsinn4dtT，2二,rC=、，上八_，E&一一，、其中，T為基波頻率，nC為n次諧波頻率。如果將（2-1）式中同頻率的正弦和余弦分量合并，則三角形式的傅里葉級數(shù)可表示為：AftvAncosn"nn=1,2,3,2-32n4Ao=a0上式中An，a2bn2,n=1,2,2-4中n=-arctan包bna0=A0an=Aco

12、sn,n=1,2,bn=-AnSinQ,可以看出，傅里葉系數(shù)an和bn者B是n或（nQ）的函數(shù)，其中an和An是n或（nfi）的偶函數(shù)，即有a_n=4;而bn和Q是n或（nQ）的奇函數(shù)，即有b_n=-bn。2.1.2指數(shù)形式的傅里葉級數(shù)根據(jù)歐拉公式：cos（n-t：n）=lej（n'"n）通”.）22-5并考慮An和Q奇偶性可將（2-3）改寫為指數(shù)形式的傅里葉級數(shù)：ft）；=，,_Fnejn'n=0,一1,一2,一3，（2一6）n二即周期信號可分解為一系列不同頻率的虛指數(shù)信號之和，式中Fn稱為傅里葉復(fù)系數(shù)，可由下式求得：Fn2-72.2 連續(xù)時間周期信號的傅里葉綜合任

13、何滿足狄里赫里條件的周期信號，可以表示成式（21）或（26）的和式形式，（21）或（2-6W常稱為連續(xù)周期信號的傅里葉級數(shù)綜合公式。一般來說，傅里葉級數(shù)系數(shù)有無限個非零值，即任何具有有限個間斷點(diǎn)的周期信號都一定有一個無限項(xiàng)非零系數(shù)的傅里葉級數(shù)表示。但對數(shù)值計算來說，這是無法實(shí)現(xiàn)的。在實(shí)際的應(yīng)用中，但我們可以用有限項(xiàng)的傅里葉級數(shù)求和來逼近。為了比較有限項(xiàng)諧波的逼近情況，本次課設(shè)編寫了程序來繪制波形以給讀者一個直觀的感受。調(diào)用xiebo.m函數(shù)文件，即可繪出周期矩形波信號各次諧波的合成波形。如圖2.1所示?；?2次諧波+3次諧波1.5-基波+2次諧波+3次諧波+6次諧波1.5-10.5I1J0&

14、#39;j.-0.5c-50510.5Ij10'。-J'：丁7-0.5c-505圖2.1周期矩形脈沖信號的合成由圖2.1可見，當(dāng)它所包含的諧波分量越多時，合成波形愈接近于原來的矩形波脈沖（由圖2.1還可以看到，合成波形所包含的諧波分量愈多時，除間斷點(diǎn)附近外，它越接近于原矩形波脈沖。在間斷點(diǎn)附近，隨著所含諧波次數(shù)的增加，合成波形的尖峰愈接近間斷點(diǎn)，但尖峰幅度并未明顯減少?？梢宰C明，即使合成波形所含諧波次數(shù)nr®時，在間斷點(diǎn)處仍有約9%的偏差，這種現(xiàn)象稱為吉布斯(Gibbs)現(xiàn)象。在傅里葉級數(shù)的項(xiàng)數(shù)取得很大時，間斷點(diǎn)處尖峰下的面積非常小以致趨近于零，因而在均方的意義上合成

15、波形同原波形的真值之間沒有區(qū)別4。2.3 吉布斯現(xiàn)象上一節(jié)中我們提到了吉布斯現(xiàn)象，本節(jié)我們將作重點(diǎn)來討論。我們知道滿足狄里赫利條件的周期函數(shù)表示成的傅立葉級數(shù)都收斂。狄里赫利條件如下：1 .在任何周期內(nèi)，x(t)必須絕對可積；2 .在任一有限區(qū)間中，x(t)只能取有限個最大值或最小值；3 .在任何有限區(qū)間上，x(t)只能有有限個第一類間斷點(diǎn)。所謂的吉布斯現(xiàn)象就是：在x(t)的不可導(dǎo)點(diǎn)上，如果我們只取x(t)等式右邊的無窮級數(shù)中的有限項(xiàng)作和X(t),那么X(t)在這些點(diǎn)上會有起伏1。具體現(xiàn)象如下圖所示，以下分別為諧波次數(shù)為N=50,N=10QN=500合成波的情況。諧波次數(shù)N=50時的合成波形諧

16、波次數(shù)N=100時的合成波形1!T11-iT¥ITL1-«F11rb0.50-1.5-1-0.500.511.5-2TLLkJL1_t1rrrr諧波次數(shù)N=500時的合成波形10.50-1.5-1-0.500.511.5-2圖2.2不同時N值時的合成波從上面的圖像中可以看出，當(dāng)N=500的時候，合成波與原來的方波擬合得非常好，但是在不可導(dǎo)的點(diǎn)上，即為x=-1.5,x=-0.5,x=0.5,x=1.5這樣的點(diǎn)的時候，合成波會有較大的波動，這就是非常明顯的吉布斯現(xiàn)象。3連續(xù)時間周期信號的頻譜分析3.1 單邊與雙邊頻譜關(guān)系如前所述，周期信號可以分解成一系列正弦（余弦）信號或虛指數(shù)

17、信號之和，為了直觀地表示出信號所含各分量的振幅An或|Fn|,隨頻率的變化情況，通常以角頻率為橫坐標(biāo)，以各次諧波的振幅人或虛指數(shù)函數(shù)|Fn|的幅度為縱坐標(biāo)，畫出如圖3.1和3.2所示的各諧波的振幅An或|匕|與角頻率的關(guān)系圖，稱為周期信號的幅度（振幅）頻譜，簡稱幅度譜。圖中每條豎線代表該頻率分量的幅度，稱為譜線。各譜線頂點(diǎn)連線的曲線（如圖中原點(diǎn)所示）稱為頻譜包絡(luò)線，它反映了各諧波分量幅度隨頻率變化的情況。圖3.1中幅度譜為單邊幅度譜（用An繪制的頻譜）。圖3.2中幅度譜為雙邊幅度譜（用|Fn|繪制的頻譜）。類似地，也可畫出各諧波初相角Q與角頻率的關(guān)系圖，如圖3.1和3.2中各諧波初相角中n與角

18、頻率的關(guān)系圖，稱為相位頻譜，簡稱相位譜。圖3.1中相位譜為單邊相位譜。圖3.2中相位譜為雙邊相位譜。如果Fn為實(shí)數(shù)，那么可用Fn的正負(fù)來表示Q為0或n也可把幅度譜和相位譜畫在一張圖上。由圖可見，周期信號的譜線只出現(xiàn)在頻率為0,。,2夏,.等原周期信號頻率的整數(shù)倍的離散頻率上，即周期信號的頻譜是離散譜。圖3.1周期信號的單邊幅度譜和相位譜0.10.050050100150-150-100-50圖3.2周期信號的雙邊幅度譜和相位譜由此可見周期信號頻譜具有三個特點(diǎn)：(1)離散性，即譜線是離散的；(2)諧波性，即譜線只出現(xiàn)在基波頻率的整數(shù)倍上；(3)收斂性，即諧波的幅度隨諧波次數(shù)的增高而減小30單邊頻

19、譜和雙邊頻譜的區(qū)別就是求值的范圍不同，單邊頻譜求的是頻率大于0的情況,而雙邊頻譜求的是所有頻率的情況，即包括頻率小于0的情況，這個區(qū)別在上面的兩張圖中可以非常明顯地看出來。3.2 以單邊幅度頻譜為例，研究脈沖寬度與頻譜的關(guān)系首先令方波首期T=5。改變脈沖寬度，就是在圖3.3中T值不變的情況下，改變的p值的大小，同時p必須小于T。在MATLAEB件里可以比較方便地改變這個值。xsqual=(x)1/2.*(x=-1/2)+1.*(x>-1/2&x<1/2)+1/2.*(x=-1/2);這個語句是控制p值的，現(xiàn)在的參數(shù)是1/2。下面是比較三種不同p值的矩形脈沖單邊頻譜圖像。&l

20、t;>*0"容易看出，在T不變的情況下，減小T®,可以使頻譜變得更密集，增大P值則可以使頻譜變得稀疏，因此，需要在不同的情況下選擇不同的T®,才能是系統(tǒng)變得更加符合實(shí)際需要。由于周期T相同，因而相鄰譜線的間隔相同；脈沖寬度窄，其頻譜包絡(luò)線第一個零點(diǎn)的頻率愈高，即信號帶寬愈寬，頻帶內(nèi)所含的分量愈多。可見，信號的頻帶寬度與脈沖寬度P成反比2。3.3 以單邊幅度頻譜為例，研究脈沖周期與頻譜的關(guān)系上面是改變P值來觀察頻譜的變化情況，現(xiàn)在來改變T值以達(dá)到改變頻譜的目的。在MATLA賀碼中a=-5;b=5;T0=b-a;這幾句代碼是用來控制方波的T值的，b-a就是方波的

21、周期T,在上面的討論中使用的參數(shù)是a=-5,b=5,現(xiàn)在將p的參數(shù)，即這句xsqual=(x)1/2.*(x=-1/2)+1.*(x>-1/2&x<1/2)+1/2.*(x=-1/2)固定為1/2,然后分別將a,b值變?yōu)椋篴=-4,b=4和a=-6,b=6,來研究方波周期對其頻譜的影響。通過觀察以上三個圖像中第一個零點(diǎn)的位置，不難看出：當(dāng)方波的周期越大，頻譜就越密集，周期越小，頻譜就越稀疏，其實(shí)這點(diǎn)也不難理解。因?yàn)閜值不變，改變T值就等于改變了T=ar中比例系數(shù)a的大小。由于周期脈沖信號的時域?qū)挾炔蛔儯@時頻譜包絡(luò)線的零點(diǎn)所在位置不變，而當(dāng)周期增長時，相鄰譜線的間隔減少，頻

22、譜變密。如果周期無限增長（這時就成為非周期信號），那么，相鄰譜線的間隔將趨近于零，周期信號的離散頻譜就過渡到非周期信號的連續(xù)頻譜。隨著周期的增長，各諧波分量的幅度也相應(yīng)減少。脈沖周期T愈大，譜線間隔愈小，頻譜越稠密；反之，則越稀疏。4典型周期脈沖的頻譜4.1周期方波脈沖頻譜的MATLAB實(shí)現(xiàn)周期方波脈沖信號如圖4.1所示，其幅度為1,脈沖寬度占空比：duty=0.5,周期T=510.80.60.40.20- 0.2- 0.4- 0.6- 0.8-1周期萬波脈沖-10-8-6-4-20246810圖4.1周期方波脈沖編寫fangbo.m函數(shù)文件，源程序文件見附錄程序四。調(diào)用函數(shù)fangbo.m,

23、即可繪出方波脈沖的雙邊頻譜，其中周期T和占空比duty可變，修改程序即可得到單邊頻譜。將在下一小節(jié)中給出不同參數(shù)時的頻譜圖。4.1,1周期方波脈沖雙邊頻譜的MATLAB實(shí)現(xiàn)圖4.2周期方波脈沖的雙邊頻譜a連續(xù)時間函數(shù)周期方波脈沖的雙邊幅度譜0.40.20L-20口.I!/_L-15-51I1L_一r_L，r、lJI。，LlUi二5_LII-10101520圖4.3周期方波脈沖的雙邊頻譜b圖4.4周期方波脈沖的雙邊頻譜c由上面三個圖可以看出，當(dāng)T一定時占空比越大頻譜主瓣的寬度越大，當(dāng)占空比一定時周期越小頻譜的主瓣寬度越大。周期方波信號頻譜與周期矩形脈沖信號具有相同的規(guī)律，這里不再贅述。4.1.2

24、周期方波脈沖單邊頻譜的MATLAB實(shí)現(xiàn)圖4.5周期方波脈沖的單邊頻譜a周期萬波脈沖的單邊幅度譜0.40.2002468101214161820圖4.6周期方波脈沖的單邊頻譜b連續(xù)時間函數(shù)周期方波脈沖的雙邊幅度譜圖4.7周期方波脈沖的單邊頻譜c單邊頻譜就是雙邊頻譜正半軸部分，其具有的規(guī)律也與雙邊頻譜相同4.2周期三角波脈沖頻譜的MATLAB實(shí)現(xiàn)周期三角波脈沖如圖4.8所示，周期T=5,其幅度為1圖4.8周期三角波脈沖MATLAB內(nèi)置有產(chǎn)生三角波的函數(shù)sawtooth(t),其調(diào)用格式為:x=sawtooth(t,width);根據(jù)width值的不同產(chǎn)生不同形狀的三角波，參數(shù)width是01之間的

25、標(biāo)量，指定在一個周期之間最大值的位置，width是該位置的橫坐標(biāo)和周期的比值.因而,當(dāng)width=0.5時產(chǎn)生標(biāo)準(zhǔn)的對稱三角波，當(dāng)width不等于0.5時(可缺省)產(chǎn)生鋸齒波。在附錄的源程序五中，只需要給出不同的T值和width值，就會得到三角波的不同的雙邊和單邊頻譜圖。4.2.1周期三角波雙邊頻譜的MATLAB實(shí)現(xiàn)T=5,脈寬為0.5的周期三角波脈沖-30-20-100102030圖4.9周期三角波脈沖的雙邊頻譜a連續(xù)時間函數(shù)周期三角脈沖的雙邊幅度譜0.40.20一trII11I1IFlII-20-15-10-5T忻05101520圖4.11周期三角波脈沖的雙邊頻譜c由上面三個圖可以看出，當(dāng)

26、三角波為脈寬width=0.5的對稱三角波時頻譜圖在0點(diǎn)的幅值為零。主瓣寬度與周期和脈沖寬度的關(guān)系與方波時的規(guī)律基本一致，這里不再贅述。4.2.2周期三角波單邊頻譜的MATLAB實(shí)現(xiàn)T=10,脈寬為0.5的周期三角波脈沖連續(xù)時間函數(shù)周期三角脈沖的單邊幅度譜024681012141618200.20.10圖4.12周期三角波脈沖的單邊頻譜aT=5,脈寬為0.5的周期三角波脈沖-30-20-100102030連續(xù)時間函數(shù)周期三角脈沖的單邊幅度譜1 ITCIEC0.2-0.1-71-中a.II¥/Tlf1平守仃G中匯rir.簾002468101214161820圖4.13周期三角波脈沖的單

27、邊頻譜bT=10,脈寬為1的周期三角波脈沖0.5-0.510連續(xù)時間函數(shù)周期三角脈沖的雙邊幅度譜0.4JI1牛19聞'下。IirI。531r.IhLL0.20024681012141618205小結(jié)即心得體會本次課程設(shè)計至此已經(jīng)接近尾聲，一周的時間雖然很短暫，但在這一個星期的設(shè)計過程中收獲頗多。設(shè)計的核心內(nèi)容就是利用MATLAB強(qiáng)大的圖形處理功能，符號運(yùn)算功能以及數(shù)值計算功能，實(shí)現(xiàn)連續(xù)時間周期信號頻域分析的仿真波形。整個設(shè)計過程中首先對所學(xué)的信號與系統(tǒng)與數(shù)字信號處理有了更深的了解，比如傅立葉級數(shù)、信號頻譜等；其次,實(shí)現(xiàn)過程是通過MATLAB軟件完成的，MATLAB的圖形功能強(qiáng)大，具有良

28、好的人機(jī)界面，此次設(shè)計過程中熟練了MATLAB的編程，掌握了很多函數(shù)的作用及使用方法；最后,通過此次課程設(shè)計，我對設(shè)計所用到的軟件MATLAB有了更加深刻地了解，MATLAB不管在數(shù)值計算方面的功能很強(qiáng)大，而且其圖形仿真功能更能滿足各個領(lǐng)域的需要，因此我們以后更要經(jīng)常運(yùn)用MATLAB軟件，使其成為自己不可或缺的工具。在寫相關(guān)源程序的時候，我還收索了大量的網(wǎng)站，在網(wǎng)上收索了很多關(guān)于MATLAB勺資料。在這個過程中我發(fā)現(xiàn)網(wǎng)上有很多有用的知識。以后應(yīng)該多注意，充分合理的利用網(wǎng)絡(luò)，通過網(wǎng)絡(luò)來學(xué)習(xí)東西。在收集資料的階段我復(fù)習(xí)了數(shù)字信號系統(tǒng)處理里的相關(guān)知識。對以前的理論知識有了更進(jìn)一步的認(rèn)識和理解。通過這

29、次課程設(shè)計我還對mathtype數(shù)學(xué)公式編輯器有了一定的了解，并且會用它編輯公式。對word也有了進(jìn)一步的掌握。雖然我順利完成了課程設(shè)計的要求，但是我感覺到我對MATLAB的理解我掌握還停留在比較淺的層次。要想真正掌握它還需要繼續(xù)努力學(xué)習(xí)它。這次課程設(shè)計也使我明白了在知識的領(lǐng)域里我還有很多很多的不足，并且再一次的深深的體會到理論和實(shí)踐之間還有很到的差別。在以后的學(xué)習(xí)中應(yīng)該多多的注意實(shí)踐知識的訓(xùn)練和積累。在以后的學(xué)習(xí)生活中要不斷的開拓自己的動手能力，不斷的訓(xùn)練自己的動手能力。這次課程設(shè)計讓我深深的明白了自己以后該做什么，該怎么去做。致謝感謝學(xué)校給我們這次MATLAB課程設(shè)計的機(jī)會，不僅讓我們更加

30、學(xué)會了MATLAB的強(qiáng)大圖形處理方法，掌握了MATLAB的編程技術(shù)，而且也鍛煉了我們的動手能力。通過這次課設(shè)讓我明白了理論聯(lián)系實(shí)踐的重要性，書本上的理論知識學(xué)了不少，我們必須得應(yīng)用到實(shí)踐當(dāng)中，做到學(xué)以致用，這樣我們才能有不斷的創(chuàng)新。這次課程設(shè)計也感謝指導(dǎo)老師在設(shè)計過程中的輔導(dǎo)以及同學(xué)們的幫助。沒有他們的幫助我不會那么快克服那些困難，也不會這么快學(xué)到這么多的知識。附錄注意：由于大部分程序都很相似，這里只給出五個主要的源程序。對于那些只需要該參數(shù)就可以實(shí)現(xiàn)的程序也給只出主要程序。源程序一：(連續(xù)周期信號的分解與綜合一一諧波分析)functionA_sym,B_sym=CTFS1symstnkxT=

31、5;tao=0.2*T;a=0.5;ifnargin<4;Nf=6;endifnargin<5;Nn=32;endx=time_fun_x(t);A0=2*int(x,t,-a,T-a)/T;%求傅里葉級數(shù)展開式的系數(shù)。As=int(2*x*cos(2*pi*n*t/T)/T,t,-a,T-a);Bs=int(2*x*sin(2*pi*n*t/T)/T,t,-a,T-a);A_sym(1)=vpa(A0,Nn);fork=1:Nf%Nf為最高次諧波的次數(shù)A_sym(k+1)=vpa(subs(As,n,k),Nn);B_sym(k+1)=vpa(subs(Bs,n,k),Nn);e

32、ndifnargout=0c=A_sym;disp(c)d=B_sym;disp(d)t=-8*a:0.01:T-a;f1=2*(0.2/2+0.1871.*cos(2*pi*1*t/5)+0.*sin(2*pi*1*t/5);f2=2*(0.1514.*cos(2*pi*2*t/5)+0.*sin(2*pi*2*t/5);f3=2*(0.1009.*cos(2*pi*3*t/5)+0.*sin(2*pi*3*t/5);f4=2*(0.0468.*cos(2*pi*4*t/5)+0.*sin(2*pi*4*t/5);f5=2*(-0.0312.*cos(2*pi*6*t/5)+0.*sin(2

33、*pi*6*t/5);f6=f1+f2;f7=f6+f3;f8=f7+f4+f5;subplot(2,2,1)plot(t,f1),holdontitle('基波')subplot(2,2,2)plot(t,f6),holdontitle('基波+2次諧波')subplot(2,2,3)plot(t,f7),holdontitle('基波+2次諧波+3次諧波')subplot(2,2,4)plot(t,f8),holdontitle('基波+2次諧波+3次諧波+6次諧波')endfunctiony=time_fun_ea=0.5;

34、T=5;h=1;tao=0.2叮；t=-8*a:0.01:T-a;e1=1/2+1/2.*sign(t+tao/2);e2=1/2+1/2.*sign(t-tao/2);y=h.*(e1-e2);%定義矩形脈沖函數(shù)（繪圖用）functionx=time_fun_x(t)h=1;x1=sym('Heaviside(t+0.5)*h;x=x1-sym('Heaviside(t-0.5)')*h;%定義矩形脈沖函數(shù)（積分用）源程序二：(吉布斯效應(yīng))%N為輸入要達(dá)到的最高次諧波的次數(shù)t=-2:0.001:2;N=input('N=');c0=0.5;fN=c0*

35、ones(1,length(t);forn=1:2:NfN=fN+cos(pi*n*t)*sinc(n/2);endfigureplot(t,fN)axis(-22-0.21.2)源程序三：(矩形脈沖的頻譜分析)functionA_sym,B_sym=CTFS2symstnyifnargin<3;Nf=input('pleasInput所需展開的最高諧波次數(shù)：Nf=');endT=input('pleasInput信號的周期T=');ifnargin<5;Nn=32;endy=fun_in(t);A0=2*int(y,t,0,T)/T;As=int(

36、2*y*cos(2*pi*n*t/T)/T,t,0,T);Bs=int(2*y*sin(2*pi*n*t/T)/T,t,0,T);A_sym(1)=double(vpa(A0,Nn);fork=1:NfA_sym(k+1)=double(vpa(subs(As,n,k),Nn);B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn);endifnargout=0S1=fliplr(A_sym)S1(1,k+1)=A_sym(1)S2=fliplr(1/2*S1)S3=fliplr(1/2*B_sym)S3(1,k+1)=0S4=fliplr(S3)S5=S2-i*S4;S

37、6=fliplr(S5);N=Nf*2*pi/T;k2=-N:2*pi/T:N;S7=S6,S5(2:end);x=fun_mcsubplot(3,1,2)stem(k2,abs(S7);axis(-150,150,0,0.12)subplot(3,1,3)stem(k2,angle(S7);axis(-150,150,-4,4)endfunctiony=fun_in(t)symsaa1T=input('pleasInput信號的周期T=');M=input('周期與脈沖寬度之比M=');A=1;tao=T/M;a=tao/2;y1=sym('Heavi

38、side(t+a1)*A;y=y1-sym('Heaviside(t-a1)*A;y=subs(y,a1,a);y=simple(y);functionx=fun_mcT=5;tao=T/5;n=4;t=-n*T:0.01:n*T;x=rectpuls(t,1);fori=1:n;x=x+rectpuls(t-i*T,1)+rectpuls(t+i*T,1);endsubplot(3,1,1)plot(t,x)holdonaxis(-20,20,0,1.2)源程序四：(典型周期脈沖的頻譜分析一一方波)functionA_sym,B_sym=fangbosymstnkyT=10;ifna

39、rgin<4;Nf=input('pleasInput所需展開的最高諧波次數(shù)Nf=');endifnargin<5;Nn=32;endy=fangbo_1;A0=2*int(y,t,0,T)/T;As=int(2*y*cos(2*pi*n*t/T)/T,t,0,T);Bs=int(2*y*sin(2*pi*n*t/T)/T,t,0,T);A_sym(1)=double(vpa(A0,Nn);fork=1:NfA_sym(k+1)=double(vpa(subs(As,n,k),Nn);B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn);endifnargout=0S1=fliplr(A_sym)S1(1,k+1)=A_sym(1)S2=fliplr(1/2*S1)S3=fliplr(1/2*B_sym)S3(1,k+1)=0S4=fliplr(S3)S5=S2-i*S4;S6=fliplr(S5);N=Nf*2*pi/T;k2=-N:2*pi/T:N;S7=S6,S5(2:end);subplot(2,1,1)x=fangbo_2T=5;t=-2*T:0.01:2*T;plot(t,x)title(&