周期信號的合成與分解實驗報告

1、武漢大學教學實驗報告電子信息學院 通信工程 專業(yè) 2017 年 9 月 14 日 實驗名稱 周期信號的合成與分解 指導教師 姓名 年級 學號 成績 一、 預習部分1. 實驗目的2. 實驗基本原理3. 主要儀器設備（含必要的元器件、工具）一、實驗目的 1在理論學習的基礎上，通過實驗深刻領會周期信號傅里葉級數分解的物理意義。 2理解實際應用中通常采用有限項級數來逼近無限項級數，此時方均誤差隨項數的增加而減小。 3觀察并初步了解 Gibbs 現象。 4深入理解周期信號的頻譜特點，比較不同周期信號頻譜的差異。 二、實驗基本原理 滿足 Dirichlet 條件的周期信號 f(t)可以分解成三角函數形式的

2、傅里葉級數，表達式為： 式中n為正整數；角頻率1由周期T1決定：。該式表明：任何滿足Dirichlet 條件的周期信號都可以分解成直流分量及許多正弦、余弦分量。這些正弦、余弦分量的頻率必定是基頻的整數倍。通常把頻率為的分量稱為基波，頻率為n的分量成為n次諧波。周期信號的頻譜只會出現在0，1,21，n1，等離散的頻率點上，這種頻譜稱為離散譜，是周期信號頻譜的主要特點。f(t)波形變化越劇烈，所包含的高頻分量的比重就越大；變化越平緩，所包含的低頻分量的比重就越大。 一般來說，將周期信號分解得到的三角函數形式的傅里葉級數的項數是無限的。也就是說，通常只有無窮項的傅里葉級數才能與原函數精確相等。但在實

3、際應用中，顯然無法取至無窮多項，而只能采用有限項級數來逼近無窮項級數。而且，所取項數越多，有限項級數就越逼近原函數，原函數與有限項級數間的方均誤差就越小，而且低次諧波分量的系數不會因為所取項數的增加而變化。當選取的傅里葉有限級數的項數越多，所合成的波形的峰起就越靠近 f(t)的不連續(xù)點。當所取得項數 N 很大時，該峰起值趨于一個常數，約等于總跳變值的 9%，這種現象稱為 Gibbs 現象。 三、需要掌握的 MATLAB 函數 結果的顯示會用到 plot 和 pause 函數，請參考 MATLAB 幫助。 二、 實驗操作部分1. 實驗數據、表格及數據處理2. 實驗操作過程（可用圖表示）3. 實驗

4、結論四、實驗內容 1.周期對稱方波信號的合成圖示方波既是一個奇對稱信號，又是一個奇諧信號。根據函數的對稱性與傅里葉系數的關系可知，它可以用無窮個奇次諧波分量的傅里葉級數來表示: 選取奇對稱周期方波的周期T=0.02s，幅度 E=6，請采用有限項級數替代無限項級數來逼近該函數。分別取前 1、10、50 和 200 項有限級數來近似，編寫程序并把結果顯示在一幅圖中，觀察它們逼近方波的過程。 MATLAB 程序如下： %奇對稱方波合成t=0:0.00001:0.1;sishu=12/pi;y=sishu*sin(100*pi*t);subplot(221)plot(t,y);axis(0,0.05,

5、-4,4);xlabel('time');ylabel('前1 項有限級數');y=0;for i=1:10y=y+sishu*(sin(2*i-1)*100*pi*t)/(2*i-1);endsubplot(222);plot(t,y);axis(0,0.05,-4,4);xlabel('time');ylabel('前10項有限級數');y=0;for i=1:50y=y+sishu*(sin(2*i-1)*100*pi*t)/(2*i-1);endsubplot(223);plot(t,y);axis(0,0.05,-4,4

6、);xlabel('time');ylabel('前50項有限級數');y=0;for i=1:200y=y+sishu*(sin(2*i-1)*100*pi*t)/(2*i-1);endsubplot(224);plot(t,y);axis(0,0.05,-4,4);xlabel('time');ylabel('前200 項有限級數');顯示結果如圖4-2所示：圖 4-2 奇對稱方波信號的合成2.觀察Gibbs現象分別取前 5、7、10和 20項有限級數來逼近奇對稱方波，觀察 Gibbs 現象。MATLAB程序如下：%觀察Gib

7、bs現象t=0:0.00001:0.1;y=0;for i=1:5y=y+sishu*(sin(2*i-1)*100*pi*t)/(2*i-1);endsubplot(221);plot(t,y);axis(0,0.05,-4,4);xlabel('time');ylabel('前5項有限級數');g=(max(y)-3)/6;legend(sprintf('Gibbs:%f',g);y=0;for i=1:7y=y+sishu*(sin(2*i-1)*100*pi*t)/(2*i-1);endsubplot(222);plot(t,y);axi

8、s(0,0.05,-4,4);xlabel('time');ylabel('前7項有限級數');g=(max(y)-3)/6;legend(sprintf('Gibbs:%f',g);y=0;for i=1:10y=y+sishu*(sin(2*i-1)*100*pi*t)/(2*i-1);endsubplot(223);plot(t,y);axis(0,0.05,-4,4);xlabel('time');ylabel('前10項有限級數');g=(max(y)-3)/6;legend(sprintf('G

9、ibbs:%f',g);y=0;for i=1:20y=y+sishu*(sin(2*i-1)*100*pi*t)/(2*i-1);endsubplot(224);plot(t,y);axis(0,0.05,-4,4);xlabel('time');ylabel('前20項有限級數');g=(max(y)-3)/6;legend(sprintf('Gibbs:%f',g);顯示結果如圖4-3所示：圖4-3 Gibbs現象3.周期對稱三角信號的合成設計采用有限項級數逼近偶對稱周期三角信號的實驗，編制程序并顯示結果。4周期信號的頻譜 

10、;分析奇對稱方波信號與偶對稱三角信號的頻譜，編制程序并顯示結果，深入 討論周期信號的頻譜特點和兩信號頻譜的差異。  五、實驗要求 1. 輸入實驗內容1中提供的奇對稱方波信號合成的MATLAB程序，生成M文件，編譯并運行，觀察合成結果。  2. 輸入實驗內容2中提供的有限項級數逼近方波信號的MATLAB程序生成M文件，編譯并運行，觀察Gibbs現象。  3. 自行編制完整的MATLAB程序，完成實驗內容3中偶對稱三角信號的合成。在實驗報告中給出程序和顯示結果。  該信號的傅里葉級數表示為: 選取偶

11、對稱周期三角信號T=0.02s，幅度E=6，采用有限項級數替代無限項級數來逼近該函數。分別取前1、5、10和100項有限級數來近似。 MATLAB程序如下： %偶對稱周期三角波 t=0:0.001:0.1;sishu=24/pi2;y=3+sishu*cos(100*pi*t);subplot(221)plot(t,y);axis(0,0.05,-1,7);xlabel('time');ylabel('前1 項有限級數');y=0;for i=1:5y=y+sishu*(sin(i*pi/2)2*cos(i*100*pi*t)/i2

12、;endy=y+3;subplot(222);plot(t,y);axis(0,0.05,-1,7);xlabel('time');ylabel('前5項有限級數');y=0;for i=1:10y=y+sishu*(sin(i*pi/2)2*cos(i*100*pi*t)/i2;endy=y+3;subplot(223);plot(t,y);axis(0,0.05,-1,7);xlabel('time');ylabel('前10項有限級數');y=0;for i=1:100y=y+sishu*(sin(i*pi/2)2*cos(

13、i*100*pi*t)/i2;endy=y+3;subplot(224);plot(t,y);axis(0,0.05,-1,7);xlabel('time');ylabel('前100 項有限級數');顯示結果如圖4-4所示：圖4-4 偶對稱三角波信號的合成4. 自行編制完整的MATLAB程序，完成實驗內容4中奇對稱方波信號和偶對稱三角波信號的頻譜分析。在實驗報告中給出程序和顯示結果，討論周期信號的頻譜特點和兩信號頻譜的差異。 為了把奇對稱方波信號和偶對稱三角波信號的頻譜做一個對比，修改圖4-2中t的步長，MATLAB程序如下：t=0:0.001:0.1;y=0

14、;for i=1:100y=y+sishu*(sin(2*i-1)*100*pi*t)/(2*i-1);endsubplot(224);plot(t,y);axis(0,0.05,-4,4);subplot(211);plot(t,y);xlabel('time');ylabel('奇對稱周期方波信號');N=100;X=fft(y,N);f=1/0.1*(-N/2:(N/2-1);subplot(212);stem(f,abs(fftshift(X);xlabel('Frequency(Hz)'); ylabel('magnitude&

15、#39;)結果顯示如圖4-5所示：圖4-5 奇對稱方波信號及其頻譜圖%接圖4-4程序：subplot(211);plot(t,y);xlabel('time');ylabel('偶對稱周期三角波信號');N=100;X=fft(y,N);f=1/0.1*(-N/2:(N/2-1);subplot(212);stem(f,abs(fftshift(X);xlabel('Frequency(Hz)');ylabel('magnitude');顯示結果如圖4-6所示：圖 4-6 偶對稱三角波信號及其頻譜圖三、 實驗效果分析（包括儀器設備

16、等使用效果）六、實驗結果分析1.由圖4-2和4-4可觀察發(fā)現，采用傅里葉有限項級數替代無限項級數來逼近這兩種函數時，隨著有限級數的增加，所得到的波形越來越接近原函數波形。2.圖4-3展現了Gibbs現象，即當選取的傅里葉有限級數的項數越多，所合成的波形的峰起就越靠近 f(t)的不連續(xù)點。當所取得項數 N 很大時，該峰起值趨于一個常數，約等于總跳變值的 9%。3.由圖4-5和圖4-6，可總結出周期信號的頻譜具有如下特點： (1)離散性。周期信號的頻譜是由不連續(xù)的譜線組成，每條譜線代表一個諧波分量。 (2)諧波性。頻譜中每條譜線只出現在基波頻率的整數倍上。 (3)收斂性。各頻率分量的譜線高度表示各次諧波分量的幅值或相位角。 兩信號頻譜的差異: 由以上周期性方波和三角波信號的頻譜分析可知，周期性三角波信號的各次諧波幅值衰減比周期性方波的頻譜衰減快得多，這說明三角波的頻率結構中低頻成分較多，而方波的高頻成分比較多。4.誤差分析：1）圖形曲線不連續(xù)是因為matlab中作圖時是取的有限的點，無法做到連續(xù)連線，故畫出的圖形曲線會出現間斷或轉折等情況。 2）所作出的圖形不是完全標準的方波或三角波是因為我們是用有限項傅里葉級數去逼近的，無法到達用無窮項去逼近作圖的效果。七、思考題 1. 利用有限項的指數形式的傅里葉級數重復奇對稱方