連續(xù)時(shí)間信號(hào)傅里葉級(jí)數(shù)分析matlab課程設(shè)計(jì)

1、<p><b>　　課程設(shè)計(jì)任務(wù)書</b></p><p><b>　　題 目: </b></p><p>　　連續(xù)時(shí)間信號(hào)的傅利葉變換及MATLAB實(shí)現(xiàn) </p><p><b>　　初始條件：</b></p

2、><p>　　MATLAB軟件，微機(jī)</p><p>　　要求完成的主要任務(wù): </p><p>　　利用MATLAB強(qiáng)大的圖形處理功能，符號(hào)運(yùn)算功能和數(shù)值計(jì)算功能，實(shí)現(xiàn)連續(xù)時(shí)間非周期信號(hào)頻域分析的仿真波形；</p><p>　　1、用MATLAB實(shí)現(xiàn)典型非周期信號(hào)的頻域分析；</p><p>　　2、用MATLAB實(shí)現(xiàn)信

3、號(hào)的幅度調(diào)制；</p><p>　　3、用MATLAB實(shí)現(xiàn)信號(hào)傅立葉變換性質(zhì)的仿真波形；</p><p>　　4、寫出課程設(shè)計(jì)報(bào)告。</p><p><b>　　時(shí)間安排：</b></p><p>　　學(xué)習(xí)MATLAB語(yǔ)言的概況 第1天</p>

4、<p>　　學(xué)習(xí)MATLAB語(yǔ)言的基本知識(shí) 第2天</p><p>　　學(xué)習(xí)MATLAB語(yǔ)言的應(yīng)用環(huán)境，調(diào)試命令，繪圖能力 第3、4天</p><p>　　課程設(shè)計(jì) 第5-9天</p><p>　　答辯

5、 第10天</p><p>　　指導(dǎo)教師簽名： 年 月 日</p><p>　　系主任（或責(zé)任教師）簽名： 年 月 日</p><p><b>

6、;　　目 錄</b></p><p>　　摘要…………………………………………………………………………………III</p><p>　　ABSTRACT……………………………………………………………………………III</p><p>　　緒論…………………………………………………………………………………IV</p><p>　　

7、1傅里葉變換原理概述 ………………………………………………………………1</p><p>　　1.1 傅里葉變換及逆變換的MATLAB實(shí)現(xiàn) ………………………………………2</p><p>　　2 用MATLAB實(shí)現(xiàn)典型非周期信號(hào)的頻域分析 ………………………………………3</p><p>　　2.1 單邊指數(shù)信號(hào)時(shí)域波形圖、頻域圖 ……………………………

8、…………3</p><p>　　2.2 偶雙邊指數(shù)信號(hào)時(shí)域波形圖、頻域圖 ………………………………………4</p><p>　　2.3 奇雙邊指數(shù)信號(hào)時(shí)域波形圖、頻域圖 ………………………………………4</p><p>　　2.4 直流信號(hào)時(shí)域波形圖、頻域圖 …………………………………………5</p><p>　　2.

9、5 符號(hào)函數(shù)信號(hào)時(shí)域波形圖、頻域圖 ……………………………………………5</p><p>　　2.6 單位階躍信號(hào)時(shí)域波形圖、頻域圖 ……………………………………………6</p><p>　　2.7 單位沖激信號(hào)時(shí)域波形圖、頻域圖 ……………………………………………6</p><p>　　2.8 門函數(shù)信號(hào)時(shí)域波形圖、頻域圖 …………………………………

10、……………7</p><p>　　3 用MATLAB實(shí)現(xiàn)信號(hào)的幅度調(diào)制 ……………………………………………………8</p><p>　　3.1 實(shí)例1 …………………………………………………………………………8</p><p>　　3.2 實(shí)例2 ………………………………………………………………………10</p><p>　

11、　4 實(shí)現(xiàn)傅里葉變換性質(zhì)的波形仿真…………………………………………………11</p><p>　　4.1 尺度變換特性 ………………………………………………………………11</p><p>　　4.2 時(shí)移特性 ………………………………………………………………………14</p><p>　　4.3 頻移特性 ……………………………………………

12、…………………………16</p><p>　　4.4 時(shí)域卷積定理 ………………………………………………………………… 18 </p><p>　　4.5 對(duì)稱性質(zhì) …………………………………………………………………………20</p><p>　　4.6 微分特性 ………………………………………………………………………22</p>&l

13、t;p>　　心得體會(huì)…………………………………………………………………………………25</p><p>　　參考文獻(xiàn)…………………………………………………………………………………26</p><p>　　附錄………………………………………………………………………………………27</p><p><b>　　摘要</b></p>

14、<p>　　MATLAB和Mathematica、Maple并稱為三大數(shù)學(xué)軟件。MATLAB在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計(jì)算方面首屈一指。Simulink是MATLAB軟件的擴(kuò)展，它是實(shí)現(xiàn)動(dòng)態(tài)系統(tǒng)建模和仿真的一個(gè)軟件包。MATLAB具有強(qiáng)大的圖形處理功能、符號(hào)運(yùn)算功能和數(shù)值計(jì)算功能。其中系統(tǒng)的仿真（Simulink）工具箱是從底層開(kāi)發(fā)的一個(gè)完整的仿真環(huán)境和圖形界面。在這個(gè)環(huán)境中，用戶可以完成面向框圖系統(tǒng)仿真的全部過(guò)程，并且

15、更加直觀和準(zhǔn)確地達(dá)到仿真的目標(biāo)[1]。</p><p>　　本文主要介紹基于MATLAB的一階動(dòng)態(tài)電路特性分析。</p><p>　　關(guān)鍵字：MATLAB；仿真；圖形處理；一階動(dòng)態(tài)電路。</p><p><b>　　Abstract</b></p><p>　　MATLAB, and Mathematica, Maple

16、, and known as the three major mathematical software. It is the application of technology in mathematics classes in numerical computing software, second to none. Simulink is an extension of MATLAB software, which is the

17、realization of dynamic system modeling and simulation of a package. MATLAB has a powerful graphics processing capabilities, symbolic computing and numerical computing functions. One system simulation (Simulink) toolbox f

18、rom the bottom of the developm</p><p>　　In this paper, MATLAB-based first-order characteristics of dynamic circuits.</p><p>　　Keywords: MATLAB；Simulation；Graphics；First Order Circuit。</p>

19、<p>　　緒論 </p><p>　　在科學(xué)技術(shù)飛速發(fā)展的今天，計(jì)算機(jī)正扮演著愈來(lái)愈重要的角色。在進(jìn)行科學(xué)研究與工程應(yīng)用的過(guò)程中，科技人員往往會(huì)遇到大量繁重的數(shù)學(xué)運(yùn)算和數(shù)值分析，傳統(tǒng)的高級(jí)語(yǔ)言Basic、Fortran 及C 語(yǔ)言等雖然能在一定程度上減輕計(jì)算量，但它們均用人員具有較強(qiáng)的編程能力和對(duì)算法有深入的研究。MATLAB 正是在這一應(yīng)用要求背景下產(chǎn)生的

20、數(shù)學(xué)類科技應(yīng)用軟件。MATLAB 是matrix 和laboratory 前三個(gè)字母的縮寫，意思是“矩陣實(shí)驗(yàn)室”，是Math Works 公司推出的數(shù)學(xué)類科技應(yīng)用軟件[2]。MATLAB 具有以下基本功能：（1）數(shù)值計(jì)算功能；（2）符號(hào)計(jì)算功能；（3）圖形處理及可視化功能；（3）可視化建模及動(dòng)態(tài)仿真功能。本文介紹了如何利用MATLAB強(qiáng)大的圖形處理功能、符號(hào)運(yùn)算功能以及數(shù)值計(jì)算功能，實(shí)現(xiàn)連續(xù)時(shí)間系統(tǒng)頻域分析。</p>&l

21、t;p>　　本次課程設(shè)計(jì)介紹了用MATLAB實(shí)現(xiàn)典型非周期信號(hào)的頻譜分析，用MATLAB實(shí)現(xiàn)信號(hào)的幅度調(diào)制以及用MATLAB實(shí)現(xiàn)信號(hào)傅里葉變換性質(zhì)的仿真波形。</p><p>　　1傅里葉變換原理概述</p><p>　　設(shè)有連續(xù)時(shí)間周期信號(hào)，它的周期為T，角頻率，且滿足狄里赫利條件，則該周期信號(hào)可以展開(kāi)成傅里葉級(jí)數(shù)，即可表示為一系列不同頻率的正弦或復(fù)指數(shù)信號(hào)之和。傅里葉級(jí)數(shù)有三角形

22、式和指數(shù)形式兩種[3]。 </p><p>　　1. 三角形式的傅里葉級(jí)數(shù)[2]： 式中系數(shù)，稱為傅里葉系數(shù)，可由下式求得： [</p><p>　　2. 指數(shù)形式的傅里葉級(jí)數(shù)[2]： </p><p>　　式中系數(shù)稱為傅里葉復(fù)系數(shù)，可由下式求得： </p><p>　　周期信號(hào)頻譜具有三個(gè)特點(diǎn)[1]：</p><p

23、>　?。?） 離散性，即譜線是離散的；</p><p>　?。?） 諧波性，即譜線只出現(xiàn)在基波頻率的整數(shù)倍上；</p><p>　?。?） 收斂性，即諧波的幅度隨諧波次數(shù)的增高而減小。</p><p>　　周期信號(hào)的傅里葉分解用Matlab進(jìn)行計(jì)算時(shí)，本質(zhì)上是對(duì)信號(hào)進(jìn)行數(shù)值積分運(yùn)算。在Matlab中有多種進(jìn)行數(shù)值積分運(yùn)算的方法，我們采用quadl函數(shù)，它有兩種

24、其調(diào)用形式。 (1) y＝quadl(‘func’, a, b)。 其中func是一個(gè)字符串，表示被積函數(shù)的.m文件名（函數(shù)名）；a、b分別表示定積分的下限和上限。 (2) y＝quadl(@myfun, a, b)。其中“@”符號(hào)表示取函數(shù)的句柄，myfun表示所定義函數(shù)的文件名。</p><p>　　1.1 傅里葉變換及逆變換的MATLAB實(shí)現(xiàn)</p><p>　　MATLAB 的

25、Symbolic Math Toolbox 提供了能直接求解傅里葉變換及逆變換的函數(shù)Fourier()及Fourier()[4]。</p><p>　　1.1 fourier 變換</p><p>　　1. (1) F=fourier(f);</p><p>　　(2) F=fourier(v);</p><p>　　(3) F=fouri

26、er(f,u,v);</p><p><b>　　說(shuō)明：</b></p><p>　　(1) F=fourier(f)是符號(hào)函數(shù)f 的Fourier 變換，缺省返回是關(guān)于ω的函數(shù)。如果 f=f(ω),則fourier 函數(shù)返回關(guān)于t 的函數(shù)。</p><p>　　(2)F=fourier(f,v)返回函數(shù)F 是關(guān)于符號(hào)對(duì)象v 的函數(shù)，而不是缺省

27、的ω</p><p>　　(3)F=fourier(f,u,v)對(duì)關(guān)于u 的函數(shù)f 進(jìn)行變換，返回函數(shù)F 是關(guān)于v 的函數(shù)。</p><p>　　1.2 fourier 逆變換</p><p>　　1. (1) f=ifourier(F);</p><p>　　(2) f=ifourier(F,u);</p><p>

28、　　(3) f=ifourier(F,v,u);</p><p><b>　　說(shuō)明：</b></p><p>　　(1) f=ifourier(F)中輸入?yún)⒘縁是傅里葉變換的符號(hào)表達(dá)式，缺省為符號(hào)變量w的函數(shù)，輸出參量f是F的傅里葉逆變換的符號(hào)表達(dá)式，缺省為符號(hào)變量x的函數(shù)。</p><p>　　(2)f=ifourier(F,u)中輸入?yún)⒘縁

29、是傅里葉變換的符號(hào)表達(dá)式，缺省為符號(hào)變量w的函數(shù)，輸出參量f是F的傅里葉逆變換的符號(hào)表達(dá)式，為指定符號(hào)變量u的函數(shù)</p><p>　　(3)f=ifourier(F,v,u)中輸入?yún)⒘縁是傅里葉變換的符號(hào)表達(dá)式，為指定符號(hào)變量v的函數(shù)，輸出參量f是F的傅里葉逆變換的符號(hào)表達(dá)式，缺省為符號(hào)變量u的函數(shù)。</p><p>　　2 用MATLAB實(shí)現(xiàn)典型非周期信號(hào)的頻域分析</p>

30、<p>　　2.1單邊指數(shù)信號(hào)時(shí)域波形圖、頻域圖</p><p>　　的時(shí)域波形圖和頻譜圖如圖2.1.1 ：</p><p>　　2.2偶雙邊指數(shù)信號(hào)時(shí)域波形圖、頻域圖</p><p>　　偶雙邊指數(shù)信號(hào)時(shí)域波形圖、頻域圖如下圖圖2.2.1：</p><p>　　圖2.2.1偶雙邊指數(shù)信號(hào)</p><p>

31、;　　2.3奇雙邊指數(shù)信號(hào)時(shí)域波形圖、頻域圖</p><p>　　奇雙邊指數(shù)信號(hào)時(shí)域波形圖、頻域圖如下圖圖2.3.1：</p><p>　　圖2.3.1奇雙邊指數(shù)信號(hào)</p><p>　　2.4 直流信號(hào)時(shí)域波形圖、頻域圖</p><p>　　直流信號(hào)f（t）=A,不滿足絕對(duì)可積條件，但傅里葉變換卻存在?？梢园褑挝恢绷餍盘?hào)看做雙邊指數(shù)信號(hào)當(dāng)

32、a趨于0時(shí)的極限。直流信號(hào)時(shí)域波形圖、頻域圖如下圖2.4.1：</p><p>　　圖2.4.1直流信號(hào)</p><p>　　2.5 符號(hào)函數(shù)信號(hào)時(shí)域波形圖、頻域圖</p><p>　　符號(hào)函數(shù)信號(hào)時(shí)域波形圖、頻域圖如下圖2.5.1：</p><p>　　圖2.5.1符號(hào)函數(shù)信號(hào)</p><p>　　2.6 單位階

33、躍信號(hào)時(shí)域波形圖、頻域圖</p><p>　　單位階躍函數(shù)信號(hào)時(shí)域波形圖、頻域圖如下圖2.6.1：</p><p>　　圖2.6.1單位階躍函數(shù)信號(hào)</p><p>　　2.7 單位沖激信號(hào)時(shí)域波形圖、頻域圖</p><p>　　單位沖激函數(shù)信號(hào)時(shí)域波形圖、頻域圖如下圖2.7.1：</p><p>　　圖2.7.1單

34、位沖激函數(shù)信號(hào)</p><p>　　2.8 門函數(shù)信號(hào)時(shí)域波形圖、頻域圖</p><p>　　門函數(shù)信號(hào)時(shí)域波形圖、頻域圖如下圖2.8.1：</p><p>　　圖2.8.1門函數(shù)信號(hào)</p><p>　　3用MATLAB實(shí)現(xiàn)信號(hào)的幅度調(diào)制</p><p>　　設(shè)信號(hào)f (t) 的頻譜為F( jw) ,現(xiàn)將f (t)

35、 乘以載波信號(hào)cos (w0t) ，得到高頻的已調(diào)信號(hào)y(t ) ，即：y(t ) = f (t) cos (w0t)</p><p>　　從頻域上看，已調(diào)制信號(hào)y(t ) 的頻譜為原調(diào)制信號(hào)f (t) 的頻譜搬移到0 ±w 處，幅度降為原F( jw) 的1/2，即</p><p>　　上式即為調(diào)制定理，也是傅里葉變換性質(zhì)中“頻移特性”的一種特別情形。</p>&l

36、t;p>　　MATLAB 提供了專門的函數(shù)modulate（）用于實(shí)現(xiàn)信號(hào)的調(diào)制。調(diào)用格式為：</p><p>　　y=modulate(x,Fc,Fs,'method')[y,t]=modulate(x,Fc,Fs)</p><p>　　其中，x 為被調(diào)信號(hào)，F(xiàn)c 為載波頻率，F(xiàn)s 為信號(hào)x 的采樣頻率，method 為所采用的調(diào)制方式，若采用幅度調(diào)制、雙邊帶調(diào)制、

37、抑制載波調(diào)制，則＇method＇為＇am＇或amdsd-sc＇。其執(zhí)行算法為y=x*cos(2*pi*Fc*t)</p><p>　　其中y 為已調(diào)制信號(hào)，t 為函數(shù)計(jì)算時(shí)間間隔向量。</p><p>　　涉及到一個(gè)函數(shù)，暫時(shí)不容易理解，因此查閱工具書，特在此說(shuō)明：</p><p>　　MATLAB 的“信號(hào)處理工具箱函數(shù)”中的估計(jì)信號(hào)的功率譜密度函數(shù)psd(),其

38、格式是：[Pxx，f]=psd(x，Nfft，F(xiàn)s，window，noverlap，dflag)</p><p>　　其中，x 是被調(diào)制信號(hào)（即本例中的f (t) ），Nfft 指定快速付氏變換FFT 的長(zhǎng)度，F(xiàn)s</p><p>　　為對(duì)信號(hào)x 的采樣頻率。后面三個(gè)參數(shù)的意義涉及到信號(hào)處理的更深的知識(shí)，在此暫不介紹。</p><p>　　3.1 信號(hào)調(diào)制實(shí)例1&

39、lt;/p><p>　　例1： f(t)=sin(100πt) f=400Hz，繪出原信號(hào)f(t)以及調(diào)制信號(hào)y(t)=f(t)coswt的實(shí)域波形圖、頻鋪圖以及功率譜。</p><p><b>　　程序如下：</b></p><p><b>　　Fm=50;</b></p><p>　　Fc=400;

40、%載波頻率</p><p>　　Fs=1000;%信號(hào)x的抽樣頻率</p><p><b>　　N=1000;</b></p><p>　　k=0:N-2;%采樣點(diǎn)</p><p>　　t=k/Fs;%采樣時(shí)間</p><p>　　x=sin(2*pi*Fm*t);%被調(diào)制信號(hào)</p>

41、<p>　　subplot(221);</p><p>　　plot(t,x);%畫出被調(diào)制信號(hào)的波形</p><p>　　xlabel('t(s)');</p><p>　　ylabel('x');</p><p>　　title('被調(diào)制信號(hào)的波形');</p>

42、<p>　　axis([0 0.1 -1 1]);%坐標(biāo)系范圍 t取值范圍不能大，因?yàn)椴蓸宇l率很高，不便于觀察</p><p>　　Nfft=1024;</p><p>　　window=hamming(512);</p><p>　　noverlap=256;</p><p>　　dflag='none';<

43、/p><p>　　[Pxx,f]=psd(x,Nfft,Fs,window,noverlap,dflag);%估算被調(diào)信號(hào)的功率譜密度Nfft是快速傅里葉變換的長(zhǎng)度</p><p>　　subplot(222);</p><p>　　plot(f,Pxx);%畫出被調(diào)信號(hào)的功率譜密度的波形</p><p>　　ylabel('功率譜(x

44、)');</p><p>　　xlabel('f(hz)');</p><p>　　%axis([0 600 0 100]);%坐標(biāo)系的范圍</p><p>　　title('被調(diào)信號(hào)的功率譜密');</p><p><b>　　grid</b></p><p&

45、gt;　　y=modulate(x,Fc,Fs,'am');%得到調(diào)制信號(hào)</p><p>　　subplot(223);</p><p>　　plot(t,y);%會(huì)出調(diào)制信號(hào)的波形</p><p>　　xlabel('t(s)');</p><p>　　ylabel('y');</p&

46、gt;<p>　　title('已調(diào)信號(hào)');</p><p>　　axis([0 0.1 -1 1]);%坐標(biāo)系的范圍t取值范圍不能大，因?yàn)椴蓸宇l率很高，不便于觀察</p><p>　　[Pxx,f]=psd(y,Nfft,Fs,window,noverlap,dflag);%估算被調(diào)信號(hào)的功率譜密度Nfft是快速傅里葉變換的長(zhǎng)度</p>&l

47、t;p>　　subplot(224);</p><p>　　plot(f,Pxx);%畫出被調(diào)信號(hào)的功率譜密度的波形</p><p>　　ylabel('功率譜(y)');</p><p>　　xlabel('f(hz)');</p><p>　　%axis([0 600 0 100]);%坐標(biāo)系的范圍

48、</p><p>　　title('已調(diào)信號(hào)功率譜');</p><p><b>　　grid</b></p><p>　　圖3.1.1調(diào)制信號(hào)與被調(diào)信號(hào)</p><p>　　3.2 信號(hào)調(diào)制實(shí)例2</p><p>　　例2：設(shè) ,繪出原信號(hào)f(t)以及調(diào)制信號(hào)y(t)=f(t)

49、coswt的實(shí)域波形圖、頻譜圖以及功率譜。</p><p><b>　　解：</b></p><p><b>　　n=0.005;</b></p><p>　　t=-1.5:n:1.5;</p><p>　　f=Heaviside(t+1)-Heaviside(t-1);</p>&l

50、t;p>　　ft=f.*cos(10*pi*t);%FT為已調(diào)信號(hào)，要滿足矩陣相乘規(guī)則，點(diǎn)乘，.w</p><p>　　subplot(221);</p><p>　　plot(t,f);%畫出被調(diào)制信號(hào)波形</p><p>　　xlabel('t');</p><p>　　ylabel('f(t)');

51、</p><p>　　title('被調(diào)制信號(hào)波形');</p><p>　　subplot(222);</p><p>　　plot(t,ft);%畫出已調(diào)制信號(hào)波形</p><p>　　xlabel('t');</p><p>　　ylabel('ft(t)');&l

52、t;/p><p>　　title('已調(diào)制信號(hào)波形');</p><p><b>　　w1=40;</b></p><p><b>　　N=1000;</b></p><p><b>　　k=-N:N;</b></p><p><b&g

53、t;　　w=w1*k/N;</b></p><p>　　Fw=f*exp(-j*t'*w)*n;%得到被調(diào)制信號(hào)頻譜</p><p>　　Ftw=ft*exp(-j*t'*w)*n;%得到已調(diào)制信號(hào)頻譜</p><p>　　Fwr=real(Fw);%熱∈擋?</p><p>　　Ftwr=real(Ftw);&l

54、t;/p><p>　　subplot(223);</p><p>　　plot(w,Fwr);%畫出被調(diào)制信號(hào)頻譜</p><p>　　xlabel('w');</p><p>　　ylabel('F(jw)');</p><p>　　title('被調(diào)制信號(hào)頻譜');<

55、;/p><p>　　subplot(224);</p><p>　　plot(w,Ftwr);%畫出已調(diào)制信號(hào)頻譜</p><p>　　xlabel('w');</p><p>　　ylabel('Ft(jw)');</p><p>　　title('已調(diào)信號(hào)頻譜');<

56、;/p><p>　　圖3.2.1 原信號(hào)f (t) 、調(diào)制信號(hào)ft( t) 的波形及其頻譜F( jw) 、Ft( jw)</p><p>　　4用MATLAB實(shí)現(xiàn)信號(hào)傅立葉變換性質(zhì)的仿真波形</p><p>　　4.1 傅里葉變換的尺度變換特性</p><p>　　若f (t) « F( jw) ，則傅里葉變換的尺度變換特性為[5]

57、：</p><p>　　例1: 設(shè) ，即門寬為τ=2 的門信號(hào)，用MATLAB 求</p><p>　　的頻譜Y ( jw) ，并與f (t) 的頻譜F( jw) 進(jìn)行比較。 </p><p><b>　　%尺度變換</b></p><p>　　n=0.02;%采樣間隔</p><p>　　t=-

58、2:n:2;%采樣范圍</p><p>　　f=Heaviside(t+1)-Heaviside(t-1);%脈寬為2的門信號(hào)</p><p>　　h=Heaviside(2*t+1)-Heaviside(2*t-1);%脈寬為1的門信號(hào)</p><p>　　w1=5*2*pi;</p><p><b>　　N=500;</b

59、></p><p><b>　　k=-N:N;</b></p><p><b>　　w=k*w1/N;</b></p><p>　　F=f*exp(-j*t'*w)*n;%求出Fw</p><p>　　H=h*exp(-j*t'*w)*n;%求出Hw</p>&l

60、t;p>　　subplot(221);</p><p>　　plot(t,f);%畫出脈寬為2的門信號(hào)的時(shí)域波形</p><p>　　xlabel('t');</p><p>　　ylabel('f(t)');</p><p>　　title('脈寬為2的門信號(hào)的時(shí)域波形');</p

61、><p>　　axis([-2.5 2.5 0 1.1]);</p><p>　　subplot(222);</p><p>　　plot(t,h);%畫出脈寬為1的門信號(hào)的時(shí)域波</p><p>　　xlabel('t');</p><p>　　ylabel('h(t)');</p&

62、gt;<p>　　title('脈寬為1的門信號(hào)的時(shí)域波形');</p><p>　　axis([-2.5 2.5 0 1.1]);</p><p>　　subplot(223);</p><p>　　plot(w,F);%畫出脈寬為2的門信號(hào)的頻域波</p><p>　　xlabel('w');

63、</p><p>　　ylabel('F(w)');</p><p>　　title('脈寬為2的門信號(hào)的頻域波形');</p><p>　　axis([-5*pi 5*pi -0.5 2.1]);</p><p>　　subplot(224);</p><p>　　plot(w,H);

64、%畫出脈寬為1的門信號(hào)的頻域波</p><p>　　xlabel('w');</p><p>　　ylabel('H(w)');</p><p>　　title('脈寬為1的門信號(hào)的頻域波形');</p><p>　　axis([-5*pi 5*pi -0.5 2.1]);</p>

65、<p>　　圖4.1.1 傅里葉變換的尺度變換特性</p><p>　　由圖4.1.1，y(t ) 信號(hào)相當(dāng)于原信號(hào)f (t) 在時(shí)域上壓縮一倍，即y(t ) = f (2t) ，a = 2 ，按式，Y ( jw) 的頻域?qū)挾葢?yīng)是F( jw) 的兩倍，而幅度下降為F( jw) 的一半。</p><p>　　4.2 傅里葉變換的時(shí)移變換特性</p><p&g

66、t;　　若f (t) « F( jw) ，則傅里葉變換的時(shí)移特性為：</p><p>　　例2： 設(shè),試用MATLAB 繪出f (t-t0),f (t+t0) 及其頻譜（幅度譜及相位譜）。</p><p>　　t0=0.2;%時(shí)移大小</p><p>　　n=0.02;%采樣間隔</p><p>　　t=-5:n:5;%采樣范圍&l

67、t;/p><p>　　f1=1/2*exp(-2*t).*Heaviside(t);%定義函數(shù)f1</p><p>　　f2=1/2*exp(-2*(t-t0)).*Heaviside(t-t0);%定義函數(shù)f2，時(shí)域右移t0</p><p>　　f3=1/2*exp(-2*(t+t0)).*Heaviside(t+t0);%定義函數(shù)f1，時(shí)域左移t0</p>

68、;<p>　　subplot(311);</p><p>　　plot(t,f1);%畫出f1,f2,f3的時(shí)域波形</p><p>　　xlabel('t');</p><p>　　ylabel('f(t)');</p><p><b>　　hold on</b></

69、p><p>　　plot(t,f2,'-.');</p><p>　　plot(t,f3,':');</p><p>　　axis([-6 6 0 0.7]);</p><p>　　legend('f1(t)','f2(t)','f3(t)');標(biāo)注f1f2f3<

70、;/p><p><b>　　N=300;</b></p><p>　　w1=5*pi*2;%頻譜寬度</p><p><b>　　k=-N:N;</b></p><p><b>　　w=k*w1/N;</b></p><p>　　F1=n*f1*exp(-j

71、*t'*w);%f1的傅里葉變換</p><p>　　F1f=abs(F1);%f1的幅度頻譜</p><p>　　F1a=angle(F1);%f1的相位頻譜</p><p>　　F2=n*f2*exp(-j*t'*w);%f2的傅里葉變換</p><p>　　F2f=abs(F2);%f2的幅度頻譜</p>

72、<p>　　F2a=angle(F2);%f2的相位頻譜</p><p>　　F3=n*f3*exp(-j*t'*w);%f3的傅里葉變換</p><p>　　F3f=abs(F3);%f3的幅度頻譜</p><p>　　F3a=angle(F3);%f3的相位頻譜</p><p>　　subplot(312);</p

73、><p>　　plot(w,F1f);</p><p>　　xlabel('w');</p><p>　　ylabel('F(w)');</p><p><b>　　hold on</b></p><p>　　plot(w,F2f,'-.');</

74、p><p>　　plot(w,F3f,':');%畫出f1f2f3的幅度譜</p><p>　　axis([-6 6 0 0.7]);</p><p>　　legend('F1f(w)','F2f(w)','F3f(w)');</p><p>　　subplot(313);</

75、p><p>　　plot(w,F1a*180/pi);</p><p>　　xlabel('w');</p><p>　　ylabel('P(度)');</p><p><b>　　hold on</b></p><p>　　plot(w,F2a*180/pi,'

76、;-.');</p><p>　　plot(w,F3a*180/pi,':');%畫出f1f2f3的相位</p><p>　　axis([-6 6 -200 200]);</p><p>　　圖4.2 .1f (t), f (t-t0),f (t+t0)及其幅頻特性與相頻特性</p><p>　　從圖4.2.1可以看出

77、，信號(hào)時(shí)移后其幅度頻譜并沒(méi)有改變，只是相位頻譜發(fā)生了改變，增加0.2w或減少0.2w。</p><p>　　4.3 傅里葉變換的頻移變換特性</p><p>　　若f (t)的傅里葉變換為 F( jw) ，則傅里葉變換的頻移特性為：</p><p>　　例4:設(shè)f (t) = e (t +1) - e (t - 1) ，試用MATLAB 繪出</p>

78、<p>　　的頻譜F1( jw) 及F 2 ( jw) ， 并與f (t) 的頻譜F( jw) 進(jìn)行比較。</p><p><b>　　w0=30;</b></p><p><b>　　n=0.02;</b></p><p><b>　　t=-1:n:1;</b></p>&

79、lt;p>　　g=Heaviside(t+1)-Heaviside(t-1);</p><p>　　gt=1/2*g.*(exp(j*w0*t)+exp(-j*w0*t));</p><p>　　subplot(221);</p><p>　　plot(t,g);</p><p>　　xlabel('t');</p

80、><p>　　ylabel('g');</p><p>　　title('脈寬為2的門函數(shù)的時(shí)域波形','FontSize',15);</p><p>　　axis([-1.5 1.5 0 1.2]);</p><p>　　subplot(222);</p><p>　　pl

81、ot(t,gt);</p><p>　　xlabel('t');</p><p>　　ylabel('gt');</p><p>　　title('調(diào)制信號(hào)的時(shí)域波形','FontSize',15);</p><p>　　axis([-1.5 1.5 0 1.2]);</p

82、><p><b>　　N=500;</b></p><p><b>　　w1=2*w0;</b></p><p><b>　　k=-N:N;</b></p><p><b>　　w=k*w1/N;</b></p><p>　　gw=n*

83、g*exp(j*t'*w);</p><p>　　gtw=n*gt*exp(j*t'*w);</p><p>　　gwf=abs(gw);</p><p>　　subplot(223);</p><p>　　plot(w,gwf);</p><p>　　xlabel('w');</

84、p><p>　　ylabel('gwf');</p><p>　　title('門函數(shù)的幅度譜','FontSize',15);</p><p>　　axis([-60 60 -0.5 2.5]);</p><p>　　gtwf=abs(gtw);</p><p>　　sub

85、plot(224);</p><p>　　plot(w,gtwf);</p><p>　　xlabel('w');</p><p>　　ylabel('gtwf');</p><p>　　title('調(diào)制信號(hào)的幅度譜','FontSize',15);</p>&l

86、t;p>　　axis([-60 60 -0.5 2.5]); </p><p>　　圖4.4.1 傅里葉變換的頻移特性</p><p>　　由圖4.4.1可見(jiàn)，對(duì)比的結(jié)果可知 g(jw) 及) gt ( jw) 是將F( jw) 分別搬移到</p><p>　　w 0= -20及w0 = 20 處的頻譜。</p><p>　　4.

87、4 傅里葉變換的時(shí)域卷積定理</p><p>　　變換的時(shí)域卷積定理如下：</p><p>　　若信號(hào)f1( t) , f2( t) 的傅里葉變換分別為, F1( jw) F2( jw) ，則： f1( t) ??f2( t) ?F1( jw) ??F2 (jw) </p><p>　　例5：設(shè)f (t) ?u(t ?1) ??u (t ?1), y(

88、t ) ??f (t) * f (t )</p><p>　　試用MATLAB 給出f(t)、y(t)、F(jω)、F(jω)?F(jω)及Y（jω)的圖形，驗(yàn)證式（9-13）的時(shí)域卷積定理。</p><p><b>　　n=0.005;</b></p><p><b>　　t=-2:n:2;</b></p>

89、<p>　　f=Heaviside(t+1)-Heaviside(t-1);%脈寬為2的門函數(shù)</p><p>　　subplot(221);</p><p>　　plot(t,f);</p><p>　　xlabel('t');</p><p>　　ylabel('f(t)');</p>

90、;<p>　　axis([-2 2 0 1.1]);</p><p>　　title('脈寬為2的門函數(shù)','FontSize',15);</p><p>　　y=conv(f,f)*n;%對(duì)f f求積</p><p>　　dt=-4:n:4;%y的取值范圍為f取值范圍的2倍</p><p>　　

91、subplot(222);</p><p>　　plot(dt,y);</p><p>　　xlabel('t');</p><p>　　ylabel('f(t)');</p><p>　　axis([-4 4 0 2.3]);</p><p>　　title('y=f(t)*f

92、(t)','FontSize',15);</p><p><b>　　N=500;</b></p><p>　　w1=2*pi*5;</p><p><b>　　k=-N:N;</b></p><p><b>　　w=k*w1/N;</b></p&

93、gt;<p>　　F=n*f*exp(j*t'*w);</p><p><b>　　F=abs(F);</b></p><p><b>　　F1=F.*F;</b></p><p>　　subplot(223);</p><p>　　plot(w,F);</p>

94、<p>　　xlabel('w');</p><p>　　ylabel('F(w)的幅度頻譜');</p><p>　　title('F(w)的幅度頻譜','FontSize',15);</p><p>　　subplot(224);</p><p>　　plot(w,

95、F1);%畫出F1(w)的幅度頻</p><p>　　xlabel('w');</p><p>　　ylabel('F1(w)的幅度頻譜');</p><p>　　title('F1(w)的幅度頻譜','FontSize',15);</p><p>　　圖4.5.1 時(shí)域卷積定理

96、</p><p>　　由圖4.5.1可見(jiàn)，F(xiàn)（jω）與F(jω)??F(jω)的圖形一致，而y(t)的波形正是我們熟知的</p><p>　　t)*f(t)的波形，F(xiàn)（jω）也是熟知的y(t)的付氏變換，從而驗(yàn)證時(shí)域卷積定理。</p><p>　　4.5 傅里葉變換的對(duì)稱性</p><p>　　例5 設(shè)f (t ) = Sa(t ), 已知

97、信號(hào)f (t) 的傅里葉變換為： 用MATLAB 求f2( t) ??pg（t）的傅里葉變換F1 (jw) , 并驗(yàn)證對(duì)稱性。</p><p>　　解：MATLAB程序?yàn)椋?lt;/p><p><b>　　n=0.1;</b></p><p>　　t=-20:n:20;</p><p>　　f1=sinc(t/pi);%抽

98、樣信號(hào)</p><p>　　f2=pi*(Heaviside(t+1)-Heaviside(t-1));%門函數(shù)</p><p>　　subplot(221);</p><p>　　plot(t,f1);%畫出抽樣信號(hào)的時(shí)域波形</p><p>　　xlabel('t');</p><p>　　ylab

99、el('f1');</p><p>　　title('抽樣信號(hào)的時(shí)域波形','FontSize',15);</p><p>　　axis([-21 21 -1 1.5]);</p><p>　　subplot(222);</p><p>　　plot(t,f2);%畫出門函數(shù)的時(shí)域波形</

100、p><p>　　xlabel('t');</p><p>　　ylabel('f2');</p><p>　　title('門信號(hào)的時(shí)域波形','FontSize',15);</p><p>　　axis([-2 2 -3 5]);</p><p>　　N=5

101、00;%采樣數(shù)</p><p><b>　　k=-N:N;</b></p><p><b>　　w1=20*pi;</b></p><p><b>　　w=k*w1/N;</b></p><p>　　F1=f1*exp(-j*t'*w)*n;%抽樣信號(hào)的傅里葉變換<

102、;/p><p>　　F2=f2*exp(-j*t'*w)*n;%門函數(shù)的傅里葉變換</p><p>　　subplot(223);</p><p>　　plot(w,F1);%畫出抽樣信號(hào)的頻譜</p><p>　　xlabel('w');</p><p>　　ylabel('F1(w)&#

103、39;);</p><p>　　title('抽樣信號(hào)的頻域','FontSize',15);</p><p>　　axis([-2 2 -3 7]);</p><p>　　subplot(224);</p><p>　　plot(w,F2);%畫出門函數(shù)的頻譜</p><p>　　xl

104、abel('w');</p><p>　　ylabel('F2(w)');</p><p>　　title('門函數(shù)信號(hào)的頻域','FontSize',15);</p><p>　　axis([-20 20 -3 7]);</p><p>　　圖4.5.1 傅里葉變換對(duì)稱性實(shí)例&

105、lt;/p><p>　　由圖4.5.1可見(jiàn)，f (t ) ??Sa(t ) 的傅里葉變換為</p><p><b>　　的傅里葉變換為</b></p><p>　　考慮到Sa(w) 是w 的偶函數(shù)，因此我</p><p>　　們有： F(t)=2*pi*f (jw) ，即驗(yàn)證了傅里葉變換的對(duì)稱性。</p>&l

106、t;p>　　4.6 傅里葉變換的時(shí)域微分性質(zhì)</p><p>　　傅里葉變換的時(shí)域微分特性為：</p><p>　　若f (t)的傅里葉變換為F(jw)，則:</p><p>　　例6: 已知f (t) 的波形如圖9.13 所示，試用MATLAB求f (t) 及df (t)/ dt的傅里葉變換，</p><p>　　F(jw) 及F

107、'(jw)，并驗(yàn)證時(shí)域微分特性。</p><p>　　圖4.7 f(t)的波形</p><p>　　解：在MATLAB中，有專門的三角波形生成函數(shù)sawtooth()，其格式為：</p><p>　　f = sawtooth(t, width)</p><p>　　其中width(0<width≤1的標(biāo)量)用于確定最大值的位置，

108、即當(dāng)t從0到2π´ width變化時(shí)，</p><p>　　f從-1上升到+1，然后當(dāng)t從2π´ width至4π時(shí)f(t)又線性地從+1下降到-1，周而復(fù)始。當(dāng)</p><p>　　width=0.5時(shí)，可產(chǎn)生一對(duì)稱的標(biāo)準(zhǔn)三角波。利用此三角波與一門信號(hào)g2π(t)相乘，再進(jìn)</p><p>　　行必要的幅度調(diào)整（乘系數(shù)2/π ），并時(shí)移（左移π

109、）可得到f(t)：</p><p>　　又設(shè)f 1 ( t) = df( t)/dt，其波形為：</p><p>　　f1( t) 可用階躍函數(shù)Heaviside()生成：</p><p><b>　　即驗(yàn)證：</b></p><p><b>　　n=0.02;</b></p><

110、;p><b>　　t=-6:n:6;</b></p><p>　　h1=Heaviside(t+pi)-Heaviside(t-pi);</p><p>　　f1=Heaviside(t+pi)-2*Heaviside(t)+Heaviside(t-pi);%三角波信號(hào)的一介倒數(shù)</p><p>　　f=pi/2*(sawtooth(t+

111、pi,0.5)+1).*h1;%三角波函數(shù)表達(dá)式</p><p>　　N=500;%采樣點(diǎn)數(shù)</p><p><b>　　k=-N:N;</b></p><p>　　w1=10*pi;%采樣頻譜寬度</p><p><b>　　w=k*w1/N;</b></p><p>　　

112、F=n*f*exp(-j*t'*w);%三角波函數(shù)的傅里葉變換</p><p>　　F1=n*f1*exp(-j*t'*w);%三角波一介函數(shù)的傅里葉變換</p><p>　　F2=F1./(j*w);%三角波一介函數(shù)的傅里葉變換除以jw</p><p>　　subplot(221);</p><p>　　plot(t,f);

113、%畫出三角波函數(shù)的時(shí)域波形</p><p>　　set(gca,'box','off');</p><p>　　xlabel('t');</p><p>　　ylabel('f(t)');</p><p>　　title('三角波時(shí)域波形','FontSi

114、ze',15);</p><p>　　subplot(222);</p><p>　　plot(t,f1);%畫出三角波一介導(dǎo)數(shù)的時(shí)域波形</p><p>　　set(gca,'box','off');</p><p>　　xlabel('t');</p><p>

115、;　　ylabel('f1(t)');</p><p>　　title('三角波一介導(dǎo)數(shù)時(shí)域波形f1(t)','FontSize',15);</p><p>　　subplot(223);</p><p>　　plot(w,F);%畫出三角波的頻譜</p><p>　　set(gca,'

116、box','off');</p><p>　　xlabel('w');</p><p>　　ylabel('F(w)');</p><p>　　axis([-40 40 0 10]);</p><p>　　title('三角波的頻譜','FontSize'

117、,15);</p><p>　　subplot(224);</p><p>　　plot(w,F2);%畫出三角波一介導(dǎo)數(shù)除jw的頻譜</p><p>　　set(gca,'box','off');</p><p>　　xlabel('w');</p><p>　　yla

118、bel('F1(w)');</p><p>　　axis([-40 40 0 10]);</p><p>　　title('F1(w)./jw','FontSize',15);</p><p>　　圖4.6.1 時(shí)域微分特性</p><p>　　程序運(yùn)行結(jié)果如圖4.6.1所示。結(jié)果表明，F(xiàn)( j

119、w) 與一致</p><p>　　從而驗(yàn)證了時(shí)域微分特性。</p><p><b>　　心得體會(huì)</b></p><p>　　第一次接觸matlab是大二上學(xué)期的實(shí)驗(yàn)課，但是只是簡(jiǎn)單接觸，只知道它畫圖很強(qiáng)大，也沒(méi)怎么認(rèn)真學(xué)，因?yàn)槲抑鲃?dòng)去接觸的未知的東西我才會(huì)全身心的投入學(xué)習(xí)。后來(lái)也下了一些教程，看了一點(diǎn)，但沒(méi)有堅(jiān)持，因?yàn)樯婕暗念I(lǐng)域太多，找不到方

120、向，也不像單片機(jī)Fpga得出的實(shí)驗(yàn)現(xiàn)象那么直觀。但是這次課設(shè)，可謂給我一個(gè)方向。Matlab語(yǔ)言，初看起來(lái)不那么順眼，但是一個(gè)小小的plot就能畫出那么形象生動(dòng)的sinx波形，這時(shí)我才明白，課本上的那些復(fù)雜的波形都可以通過(guò)matlab畫出，剛開(kāi)始寫程序的確很痛苦，因?yàn)閙atlab的界面不是那么友好，字體，標(biāo)記不如其它編譯軟件那么順手，后來(lái)更改了字體和編輯工具，感覺(jué)對(duì)matlab的興趣有如對(duì)Fpga一般。初學(xué)都是很痛苦的，需要很大的耐性和

121、良好的身體素質(zhì)，在電腦面前一呆好幾個(gè)小時(shí)那是很正常的。但是呆的越久，寫的越多，領(lǐng)悟也就越多。特別是那些在課本上晦澀難懂的傅里葉變換，那些當(dāng)初怎么想也想不到的圖形就這么在電腦前出現(xiàn)了。對(duì)那些變換性質(zhì)有個(gè)更直觀的理解，我想在學(xué)習(xí)信號(hào)與系統(tǒng)時(shí)，就應(yīng)該好好同步學(xué)matlab。Matlab有很多應(yīng)用的地方，數(shù)學(xué)，信號(hào)，我想現(xiàn)在我找到了學(xué)matlab</p><p><b>　　參考文獻(xiàn)</b><

122、/p><p>　　[1]MATLAB及在電子信息課程中的應(yīng)用，陳懷琛主編，電子工業(yè)出版社，2003.07（第二版）</p><p>　　[2]信號(hào)與系統(tǒng)分析及MATLAB實(shí)現(xiàn)，梁虹編，電子工業(yè)出版社，2002.02</p><p><b>　　參考書</b></p><p>　　[3]MATLAB6.X信號(hào)處理，鄒鯤等編，清

123、華大學(xué)出版，2002.05（第一版）</p><p>　　[4]精通MATLAB6.5版，張志涌編，北京航空航天大學(xué)出版社，2003.03（第一版）</p><p>　　[5]MATLAB M語(yǔ)言高級(jí)編程，陳永春編，清華大學(xué)出版社，2004.01（第一版）</p><p>　　[6]MATLAB程序設(shè)計(jì)，阮沈勇、王永利等編，電子工業(yè)出版社，2004.01（第一版）&

124、lt;/p><p><b>　　附錄</b></p><p>　　MATLAB 常用函數(shù)表</p><p>　　本科生課程設(shè)計(jì)成績(jī)?cè)u(píng)定表</p><p><b>　　指導(dǎo)教師簽字：</b></p><p><b>　　年 月 日</b></p&