matlab課程設(shè)計(jì)（含完整程序）.doc

matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 課程設(shè)計(jì)任務(wù)書(shū)課程設(shè)計(jì)任務(wù)書(shū) 學(xué)生姓名學(xué)生姓名: 王偉 _ 專(zhuān)業(yè)班級(jí)專(zhuān)業(yè)班級(jí): 電子科學(xué)與技術(shù) 0703 班 指導(dǎo)教師指導(dǎo)教師: 鐘毅 _ 工作單位工作單位: 信息工程學(xué)院 題題 目目: 連續(xù)時(shí)間信號(hào)傅里葉級(jí)數(shù)分析及 matlab 實(shí)現(xiàn) 初始條件：初始條件： matlab 6.5 要求完成的主要任務(wù)：要求完成的主要任務(wù)： 深入研究連續(xù)時(shí)間信號(hào)傅里葉級(jí)數(shù)分析的理論知識(shí)，利用 matlab 強(qiáng)大的圖 形處理功能，符號(hào)運(yùn)算功能以及數(shù)值計(jì)算功能，實(shí)現(xiàn)連續(xù)時(shí)間周期信號(hào)頻域分析 的仿真波形。 1.用 matlab 實(shí)現(xiàn)周期信號(hào)的傅里葉級(jí)數(shù)分解與綜合。 2.用 matlab 實(shí)現(xiàn)周期信號(hào)的單邊頻譜及雙邊頻譜。 3.用 matlab 實(shí)現(xiàn)典型周期信號(hào)的頻譜。 4.撰寫(xiě)matlab 應(yīng)用實(shí)踐課程設(shè)計(jì)說(shuō)明書(shū)。 時(shí)間安排：時(shí)間安排： 學(xué)習(xí) matlab 語(yǔ)言的概況 第 1 天 學(xué)習(xí) matlab 語(yǔ)言的基本知識(shí) 第 2、3 天 學(xué)習(xí) matlab 語(yǔ)言的應(yīng)用環(huán)境，調(diào)試命令，繪圖能力 第 4、5 天 課程設(shè)計(jì) 第 6-9 天 答辯 第 10 天 指導(dǎo)教師簽名：指導(dǎo)教師簽名： 年年 月月 日日 系主任（或責(zé)任教師）簽名：系主任（或責(zé)任教師）簽名： 年年 月月 日日 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 目 錄 摘要i abstractii 緒論1 1 matlab 內(nèi)容簡(jiǎn)介.2 1.1 matlab 語(yǔ)言功能 2 1.2 matlab 語(yǔ)言特點(diǎn) 2 2 連續(xù)時(shí)間周期信號(hào)的傅里葉級(jí)數(shù)ctfs3 2.1 連續(xù)時(shí)間周期信號(hào)的分解 .3 2.1.1 三角形式的傅里葉級(jí)數(shù).3 2.1.2 指數(shù)形式的傅里葉級(jí)數(shù).4 2.2 連續(xù)時(shí)間周期信號(hào)的傅里葉綜合 .4 3 連續(xù)時(shí)間周期信號(hào)的頻譜分析6 4 周期信號(hào)的傅里葉級(jí)數(shù)分解與綜合的仿真波形8 4.1 實(shí)現(xiàn)流程 .8 4.2 matlab 算法提示 8 4.3 程序運(yùn)行結(jié)果 .9 5 連續(xù)時(shí)間周期信號(hào)的頻譜分析的仿真波形.10 5.1 實(shí)現(xiàn)流程10 5.2 matlab 算法提示 11 5.3 程序運(yùn)行結(jié)果及分析 .12 5.3.1 程序運(yùn)行結(jié)果.12 5.3.2 脈沖寬度與頻譜的關(guān)系13 5.3.3 脈沖周期與頻譜的關(guān)系15 6 典型周期脈沖的頻譜18 6.1 周期方波脈沖頻譜的 matlab 實(shí)現(xiàn).18 6.2 周期三角波脈沖頻譜的 matlab 實(shí)現(xiàn)20 結(jié)束語(yǔ)23 致謝24 參考文獻(xiàn)25 附錄26 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) i 摘要 matlab 目前已發(fā)展成為由 matlab 語(yǔ)言、matlab 工作環(huán)境、matlab 圖形 處理系統(tǒng)、matlab 數(shù)學(xué)函數(shù)庫(kù)和 matlab 應(yīng)用程序接口五大部分組成的集數(shù)值計(jì)算、 圖形處理、程序開(kāi)發(fā)為一體的功能強(qiáng)大的系統(tǒng)。本次課程設(shè)計(jì)則在深入研究連續(xù)時(shí)間信 號(hào)傅里葉級(jí)數(shù)分析理論知識(shí)的基礎(chǔ)上，利用 matlab 強(qiáng)大的圖形處理功能、符號(hào)運(yùn)算功 能以及數(shù)值計(jì)算功能，通過(guò) matlab 編程進(jìn)行圖形功能仿真，從而實(shí)現(xiàn)連續(xù)時(shí)間周期信 號(hào)頻域分析的仿真波形，包括以下內(nèi)容：用 matlab 實(shí)現(xiàn)周期信號(hào)的傅里葉級(jí)數(shù)分解與 綜合的波形；用 matlab 實(shí)現(xiàn)周期信號(hào)的單邊頻譜及雙邊頻譜的波形與分析；用 matlab 實(shí)現(xiàn)典型周期信號(hào)的頻譜的波形。 關(guān)鍵詞：matlab；圖形處理；傅里葉級(jí)數(shù)；周期信號(hào)；頻譜 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) ii abstract matlab now evolved into matlab language, matlab working environment, matlab graphics processing systems, matlab math library and the matlab application program interface has five major components of the set of numerical computation, graphics processing, program development as one powerful system. the curriculum design, in-depth study fourier series analysis of continuous-time signal on the basis of theoretical knowledge, using matlab a powerful graphics processing capabilities, symbolic computing and numerical computing capabilities, through the functional simulation matlab graphical programming in order to achieve continuous time periodic signal frequency domain analysis of the simulation waveforms, including the following: realization of periodic signals using matlab fourier series decomposition and integration of the waveform; periodic signals using matlab to achieve unilateral and bilateral spectrum waveform and spectrum analysis; using matlab to achieve a typical cycle of the signal wave spectrum. keywords: matlab; graphics processing; fourier series; periodic signal; spectrum matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 1 緒論 在科學(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ì)算量，但它們均要求應(yīng)用人 員具有較強(qiáng)的編程能力和對(duì)算法有深入的研究。另外，在運(yùn)用這些高級(jí)語(yǔ)言進(jìn)行計(jì)算結(jié) 果的可視化分析及圖形處理方面，對(duì)非計(jì)算機(jī)專(zhuān)業(yè)的普通用戶(hù)來(lái)說(shuō)，仍存在著一定的難 度。matlab 正是在這一應(yīng)用要求背景下產(chǎn)生的數(shù)學(xué)類(lèi)科技應(yīng)用軟件。它具有的頂尖的 數(shù)值計(jì)算功能、強(qiáng)大的圖形可視化功能及簡(jiǎn)潔易學(xué)的“科學(xué)便捷式”工作環(huán)境和編程語(yǔ) 言，從根本上滿(mǎn)足了科技人員對(duì)工程數(shù)學(xué)計(jì)算的要求，并將科技人員從繁重的數(shù)學(xué)運(yùn)算 中解放出來(lái)，因而越來(lái)越受到廣大科技工作者的普遍歡迎1。 matlab 是matrix 和laboratory 前三個(gè)字母的縮寫(xiě)，意思是“矩陣實(shí)驗(yàn)室”，是 mathworks 公司推出的數(shù)學(xué)類(lèi)科技應(yīng)用軟件。其dos 版本（matlab 1.0）發(fā)行于1984 年，現(xiàn)已推出了windows 版本（matlab 5.3）。經(jīng)過(guò)十多年的不斷發(fā)展與完善， matlab 已發(fā)展成為由matlab 語(yǔ)言、matlab 工作環(huán)境、matlab 圖形處理系統(tǒng)、 matlab 數(shù)學(xué)函數(shù)庫(kù)和matlab 應(yīng)用程序接口五大部分組成的集數(shù)值計(jì)算、圖形處理、 程序開(kāi)發(fā)為一體的功能強(qiáng)大的系統(tǒng)。matlab 由“主包”和三十多個(gè)擴(kuò)展功能和應(yīng)用學(xué) 科性的工具箱（toolboxs）組成。 matlab 語(yǔ)言是以矩陣計(jì)算為基礎(chǔ)的程序設(shè)計(jì)語(yǔ)言，語(yǔ)法規(guī)則簡(jiǎn)單易學(xué)，用戶(hù)不用 花太多時(shí)間即可掌握其編程技巧。其指令格式與教科書(shū)中的數(shù)學(xué)表達(dá)式非常相近，用 matlab 編寫(xiě)程序尤如在便箋上列寫(xiě)公式和求解，因而被稱(chēng)為“便箋式”的編程語(yǔ)言。另 外，matlab 還具有功能豐富和完備的數(shù)學(xué)函數(shù)庫(kù)及工具箱，大量繁雜的數(shù)學(xué)運(yùn)算和分 析可通過(guò)調(diào)用matlab 函數(shù)直接求解，大大提高了編程效率，其程序編譯和執(zhí)行速度遠(yuǎn) 遠(yuǎn)超過(guò)了傳統(tǒng)的c 和fortran 語(yǔ)言，因而用matlab 編寫(xiě)程序，往往可以達(dá)到事半功倍 的效果。在圖形處理方面，matlab 可以給數(shù)據(jù)以二維、三維乃至四維的直觀(guān)表現(xiàn)，并 在圖形色彩、視角、品性等方面具有較強(qiáng)的渲染和控制能力，使科技人員對(duì)大量原始數(shù) 據(jù)的分析變得輕松和得心應(yīng)手。 正是由于 matlab 在數(shù)值計(jì)算及符號(hào)計(jì)算等方面的強(qiáng)大功能，使matlab一路領(lǐng) 先，成為數(shù)學(xué)類(lèi)科技應(yīng)用軟件中的佼佼者。目前，matlab 已成為國(guó)際上公認(rèn)的最優(yōu)秀 的科技應(yīng)用軟件。matlab 的上述特點(diǎn)，使它深受工程技術(shù)人員及科技專(zhuān)家的歡迎，并 很快成為應(yīng)用學(xué)科計(jì)算機(jī)輔助分析、設(shè)計(jì)、仿真、教學(xué)等領(lǐng)域不可缺少的基礎(chǔ)軟件。 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 2 1 matlab 內(nèi)容簡(jiǎn)介 1.1 matlab 語(yǔ)言功能 matlab功能豐富，可擴(kuò)展性強(qiáng)。matlab軟件包括基本部分和專(zhuān)業(yè)擴(kuò)展兩大部分 的功能?；静糠职ǎ壕仃嚨倪\(yùn)算和各種變換；代數(shù)和超越方程的求解；數(shù)據(jù)處理和 傅立葉變換；數(shù)值部分等等，可以充分滿(mǎn)足大學(xué)理工科本科的計(jì)算需要。擴(kuò)展部分稱(chēng)為 工具箱。它實(shí)際上是用matlab的基本語(yǔ)句辯稱(chēng)的各種子程序集，用于解決某一方面的 專(zhuān)門(mén)問(wèn)題，或?qū)崿F(xiàn)某一類(lèi)的新算法。 matlab 具有以下基本功能： （1）數(shù)值計(jì)算功能； （2）符號(hào)計(jì)算功能； （3）圖形處理及可視化功能； （3）可視化建模及動(dòng)態(tài)仿真功能。 1.2 matlab 語(yǔ)言特點(diǎn) matlab 給用戶(hù)帶來(lái)的是最直觀(guān)、最簡(jiǎn)潔的程序開(kāi)發(fā)環(huán)境。它具有以下特點(diǎn)： （1）語(yǔ)言簡(jiǎn)潔緊湊，使用方便靈活，庫(kù)函數(shù)極其豐富。matlab 程序書(shū)寫(xiě)形式自 由，利用起豐富的庫(kù)函數(shù)避開(kāi)繁雜的子程序編程任務(wù)，壓縮了一切不必要的編程工作。 由于庫(kù)函數(shù)都由本領(lǐng)域的專(zhuān)家編寫(xiě)，用戶(hù)不必?fù)?dān)心函數(shù)的可靠性。 （2）運(yùn)算符豐富。由于matlab 是用c 語(yǔ)言編寫(xiě)的，matlab 提供了和c語(yǔ)言幾 乎一樣多的運(yùn)算符，靈活使用matlab 的運(yùn)算符將使程序變得極為簡(jiǎn)短。 （3）matlab 既具有結(jié)構(gòu)化的控制語(yǔ)句（如for 循環(huán)，while 循環(huán)，break 語(yǔ)句和if 語(yǔ)句），又有面向?qū)ο缶幊痰奶匦浴?（4）程序限制不嚴(yán)格，程序設(shè)計(jì)自由度大。例如，在matlab 里，用戶(hù)無(wú)需對(duì)矩 陣預(yù)定義就可使用。 （5）程序的可移植性很好，基本上不做修改就可以在各種型號(hào)的計(jì)算機(jī)和操作系統(tǒng) 上運(yùn)行。 （6）matlab 的圖形功能強(qiáng)大。在fortran 和c 語(yǔ)言里，繪圖都很不容易，但 在matlab 里，數(shù)據(jù)的可視化非常簡(jiǎn)單。matlab 還具有較強(qiáng)的編輯圖形界面的能力。 （7）功能強(qiáng)大的工具箱是matlab 的另一特色。matlab 包含兩個(gè)部分：核心部 分和各種可選的工具箱。核心部分中有數(shù)百個(gè)核心內(nèi)部函數(shù)。其工具箱又分為兩類(lèi)：功 能性工具箱和學(xué)科性工具箱。功能性工具箱主要用來(lái)擴(kuò)充其符號(hào)計(jì)算功能，圖示建模仿 真功能，文字處理功能以及與硬件實(shí)時(shí)交互功能，而學(xué)科性工具箱是專(zhuān)業(yè)性比較強(qiáng)的， 如control，toolbox，signl proceessing， toolbox，commumnication toolbox 等。 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 3 2 連續(xù)時(shí)間周期信號(hào)的傅里葉級(jí)數(shù)ctfs 周期信號(hào)是定義在區(qū)間，按一定時(shí)間間隔(周期)不斷重復(fù)的信號(hào)。它可表-+，t 示為，式中為任意整數(shù)，為周期，周期的倒數(shù)稱(chēng)為該信號(hào)的頻率2。 f t =f t+mtmt 2.1 連續(xù)時(shí)間周期信號(hào)的分解 設(shè)有周期信號(hào) ，它的周期為t ，角頻率，且滿(mǎn)足狄里赫里條件，則 f t 2 =2 f= t 該周期信號(hào)可以展開(kāi)成傅里葉級(jí)數(shù)，即可表示為一系列不同頻率的正弦或復(fù)指數(shù)信號(hào)之 和。傅里葉級(jí)數(shù)有三角形式和指數(shù)形式兩種2。 2.1.1 三角形式的傅里葉級(jí)數(shù) 三角形式的傅里葉級(jí)數(shù)為： 0 12 2 tcoscos 2 a fatat 12 sinsin 2btbt 0 11 cossin 2 nn nn a an tbn t 1,2,3,n 2 1 式中系數(shù)、稱(chēng)為傅里葉系數(shù)，可由下式求得。 n a n b 2 2 2 2 2 2 0 1 2 cos 2 sin t t t t t t n n af t dt t af tn t dt t bf tn t dt t 22 如果將式中同頻率的正弦和余弦分量合并，則三角形式的傅里葉級(jí)數(shù)可表示為：2 1 0 1 tcos 2 nn n a fan t 1,2,3,n 23 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 4 上式中 00 22, 1,2, arctan nnn n n n aa aabn a b 24 00 cos,1,2, sin nnn nnn aa aan ba 由式可見(jiàn)；傅里葉系數(shù)和都是或的函數(shù)，其中是或的偶22 n a n bnn n ann 函數(shù)，即有；而是或的奇函數(shù)，即有。 n a n a n bnn n b n b 由式可見(jiàn)，是或的偶函數(shù)，即有；而是或的奇函24 n ann n a n a n nn 數(shù)，即有。 n n 式表明，任何滿(mǎn)足狄里赫里條件的周期信號(hào)可分解為一系列不同頻率的余弦23 （或正弦）分量的疊加。其中第一項(xiàng)是常數(shù)項(xiàng)，它是周期信號(hào)中所包含的直流分量； 0 2a 第二項(xiàng)稱(chēng)為基波或一次諧波，它的角頻率與原周期信號(hào)相同，是基波振 11 cosat 1 a 幅，是基波初相角；第三項(xiàng)稱(chēng)為二次諧波，它的頻率是基波頻率的二 1 22 cos 2at 倍， 是二次諧波振幅，是其初相角；依此類(lèi)推，還有三次、四次、 等諧波。一 2 a 2 般而言，稱(chēng)為次諧波，是次諧波振幅，是其初相角。式cos nn an t n n an n 表明，周期信號(hào)可以分解為各次諧波分量的疊加。23 2.1.2 指數(shù)形式的傅里葉級(jí)數(shù) 指數(shù)形式的傅里葉級(jí)數(shù)表達(dá)式為： ,0, 1, 2, 3, jn t n n f tf en 25 即周期信號(hào)可分解為一系列不同頻率的虛指數(shù)信號(hào)之和，式中稱(chēng)為傅里葉復(fù)系數(shù)， n f 可由下式求得： 2 2 1 t t jn t n ff t edt t 26 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 5 2.2 連續(xù)時(shí)間周期信號(hào)的傅里葉綜合 任何滿(mǎn)足狄里赫里條件的周期信號(hào)，可以表示成式或的和式形式，2 125 或式常稱(chēng)為ctfs綜合公式。2 125 一般來(lái)說(shuō)，傅里葉級(jí)數(shù)系數(shù)有無(wú)限個(gè)非零值，即任何具有有限個(gè)間斷點(diǎn)的周期信號(hào) 都一定有一個(gè)無(wú)限項(xiàng)非零系數(shù)的傅里葉級(jí)數(shù)表示。但對(duì)數(shù)值計(jì)算來(lái)說(shuō)，這是無(wú)法實(shí)現(xiàn)的。 在實(shí)際的應(yīng)用中，但我們可以用有限項(xiàng)的傅里葉級(jí)數(shù)求和來(lái)逼近。即對(duì)有限項(xiàng)和 n jn t n nn f tf e 0 11 cossin 2 nn nn nn a an tbn t 27 當(dāng)值取得較大時(shí)，上式就是原周期信號(hào)的一個(gè)很好的近似。式常稱(chēng)作n f t27 的截?cái)喔道锶~級(jí)數(shù)表示。 f t matlab的符號(hào)積分函數(shù)可以幫助我們求出連續(xù)時(shí)間周期信號(hào)的截?cái)喔道锶~級(jí) int 數(shù)及傅里葉表示。 求積指令的具體使用格式如下：int 給出符號(hào)表達(dá)式對(duì)指定變量的(不帶積分常數(shù)的)不定積分。intint,fff 給出符號(hào)表達(dá)式對(duì)指定變量的定積分。intint, , ,ffa bf matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 6 3 連續(xù)時(shí)間周期信號(hào)的頻譜分析 如前所述，周期信號(hào)可以分解成一系列正弦（余弦）信號(hào)或虛指數(shù)信號(hào)之和，即 n jn t n nn f tf e 0 11 cossin 2 nn nn nn a an tbn t 3 1 其中， 11 22 n j nnnn fa eajb 32 或 幅度和相位 22 11 | 22 arctan nnnn n n n faab b a 為了直觀(guān)地表示出信號(hào)所含各分量的振幅或，隨頻率的變化情況，通常以角 n a| n f 頻率為橫坐標(biāo)，以各次諧波的振幅或虛指數(shù)函數(shù)的幅度為縱坐標(biāo)，畫(huà)出如圖3.1和3.2 n a| n f 所示的各諧波的振幅或與角頻率的關(guān)系圖，稱(chēng)為周期信號(hào)的幅度（振幅）頻譜， n a| n f 簡(jiǎn)稱(chēng)幅度譜。圖中每條豎線(xiàn)代表該頻率分量的幅度，稱(chēng)為譜線(xiàn)。各譜線(xiàn)頂點(diǎn)連線(xiàn)的曲線(xiàn) （如圖中原點(diǎn)所示）稱(chēng)為頻譜包絡(luò)線(xiàn)，它反映了各諧波分量幅度隨頻率變化的情況。圖 3.1中幅度譜為單邊幅度譜（用繪制的頻譜）。圖3.2中幅度譜為雙邊幅度譜（用繪 n a| n f 制的頻譜）。 類(lèi)似地，也可畫(huà)出各諧波初相角與角頻率的關(guān)系圖，如圖3.1和3.2中各諧波初相角 n 與角頻率的關(guān)系圖，稱(chēng)為相位頻譜，簡(jiǎn)稱(chēng)相位譜。圖3.1中相位譜為單邊相位譜。圖3.2 n 中相位譜為雙邊相位譜。如果為實(shí)數(shù)，那么可用的正負(fù)來(lái)表示為0或也可把幅度 n f n f n 譜和相位譜畫(huà)在一張圖上。 由圖可見(jiàn)，周期信號(hào)的譜線(xiàn)只出現(xiàn)在頻率為等原周期信號(hào)頻率的整數(shù)倍的0,2 ,. 離散頻率上，即周期信號(hào)的頻譜是離散譜。 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 7 圖3.1 周期信號(hào)的幅度譜和相位譜（單邊） 圖3.2 周期信號(hào)的幅度譜和相位譜（雙邊） 由此可見(jiàn)周期信號(hào)頻譜具有三個(gè)特點(diǎn)： （1）離散性，即譜線(xiàn)是離散的； （2）諧波性，即譜線(xiàn)只出現(xiàn)在基波頻率的整數(shù)倍上； （3）收斂性，即諧波的幅度隨諧波次數(shù)的增高而減小3。 matlab的符號(hào)積分函數(shù)、一維數(shù)組的尋訪(fǎng)概念和可以幫助我們求出 int fliplr 連續(xù)時(shí)間周期信號(hào)的頻譜（傅里葉級(jí)數(shù)分析法）。 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 8 4 周期信號(hào)的傅里葉級(jí)數(shù)分解與綜合的仿真波形 matlab強(qiáng)大的符號(hào)運(yùn)算功能為進(jìn)行周期信號(hào)的分析提供了強(qiáng)有力的工具。以周期 矩形脈沖信號(hào)為例，采用三角形式傅里葉級(jí)數(shù)分解與綜合形式，利用matlab來(lái)實(shí)現(xiàn)周 期信號(hào)的分解與綜合過(guò)程。 周期矩形脈沖信號(hào)如圖4.1所示。其幅度為1，脈沖寬度為1，周期。tao 5*ttao 用式求出傅里葉級(jí)數(shù)分解系數(shù)，運(yùn)用matlab的符號(hào)運(yùn)算功能，用式實(shí)現(xiàn)信223 1 號(hào)的綜合，諧波的階數(shù)nf =6。 圖4.1 周期矩形脈沖信號(hào) 4.1 實(shí)現(xiàn)流程 利用matlab實(shí)現(xiàn)上述分析過(guò)程的流程如下： （1）編寫(xiě)子函數(shù)x=time_fun_x(t)，用符號(hào)表達(dá)式表示出周期信號(hào)在第一個(gè)周期內(nèi)的 符號(hào)表達(dá)式，并賦值給返回符號(hào)變量x。 （2）編寫(xiě)子函數(shù) y=time_fun_e(t)，求出該周期信號(hào)在繪圖區(qū)間內(nèi)的信號(hào)樣值，并賦 值給返回變量y。 （3）編寫(xiě)求解信號(hào)傅里葉系數(shù)及繪制合成波形圖的通用函數(shù)ctfs1.m，該函數(shù)的流 程如下： 調(diào)用函數(shù)time_fun_x(t)，獲取周期信號(hào)的符號(hào)表達(dá)式。 求出信號(hào)的傅里葉系數(shù)。 求出各次諧波。 繪制各次諧波疊加波形圖。 調(diào)用函數(shù)time_fun_e(t)，繪制原信號(hào)波形圖。 4.2 matlab 算法提示 (1) 采用符號(hào)積分 4求 內(nèi)時(shí)間函數(shù)的三角級(jí)數(shù)展開(kāi)系數(shù)：，int0,t 0 a 0 a n aas ，即計(jì)算式。 n bbs22 (2) 用循環(huán)語(yǔ)句求出三角級(jí)數(shù)展開(kāi)系數(shù)，的數(shù)值：分別為，.for end n a n b_asym matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 9 。_bsym (3) 用語(yǔ)句輸出三角級(jí)數(shù)展開(kāi)系數(shù)，。 disp_asym_bsym (4) 用傅里葉三角級(jí)數(shù)展開(kāi)式合成(綜合)連續(xù)時(shí)間信號(hào)。3 1 (5) 化簡(jiǎn)表達(dá)式，據(jù)函數(shù)的奇偶性，可知，若為奇函數(shù)，則。若為偶 f t n a0 f t 函數(shù)，則，以此化簡(jiǎn)三角級(jí)數(shù)展開(kāi)式。0 n b 源程序見(jiàn)附錄中程序一。 4.3 程序運(yùn)行結(jié)果 調(diào)用ctfs1.m函數(shù)文件，即可繪出周期矩形波信號(hào)各次諧波的合成波形。指令如下： 在matlab命令窗口鍵入ctfs1，并回車(chē)，即可繪出周期矩形波信號(hào)各次諧波的合成波 形，如圖4.2所示。 圖4.2 周期矩形脈沖信號(hào)的合成 由圖4.2可見(jiàn)，當(dāng)它所包含的諧波分量越多時(shí)，合成波形愈接近于原來(lái)的矩形波脈沖 （圖中虛線(xiàn)）。由圖4.2還可以看到，合成波形所包含的諧波分量愈多時(shí)，除間斷點(diǎn)附近 外，它越接近于原矩形波脈沖。在間斷點(diǎn)附近，隨著所含諧波次數(shù)的增加，合成波形的 尖峰愈接近間斷點(diǎn)，但尖峰幅度并未明顯減少。可以證明，即使合成波形所含諧波次數(shù) 時(shí)，在間斷點(diǎn)處仍有約9%的偏差，這種現(xiàn)象稱(chēng)為吉布斯（gibbs）現(xiàn)象。在傅里葉n 級(jí)數(shù)的項(xiàng)數(shù)取得很大時(shí)，間斷點(diǎn)處尖峰下的面積非常小以致趨近于零，因而在均方的意 義上合成波形同原波形的真值之間沒(méi)有區(qū)別4。 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 10 5 連續(xù)時(shí)間周期信號(hào)的頻譜分析的仿真波形 由周期信號(hào)頻譜分析只要求出了周期信號(hào)傅里葉級(jí)數(shù)（或）及，我們就可根 n a n f n 據(jù)（或）及隨頻率的變化關(guān)系畫(huà)出信號(hào)的振幅和相位頻譜。以圖5.1 n a n f n n 2 t 所示的周期矩形脈沖信號(hào)為例，利用matlab繪制周期信號(hào)的頻譜圖，并對(duì)周期信號(hào)的 頻譜特性進(jìn)行分析。 圖5.1 周期矩形脈沖信號(hào) 由于繪制頻譜的前提是必須先求出周期信號(hào)的傅里葉系數(shù)，因此只需對(duì)求周期信號(hào) 傅里葉級(jí)數(shù)的函數(shù)ctfs1.m進(jìn)行適當(dāng)修改，即可編寫(xiě)出繪制周期信號(hào)頻譜的通用函數(shù)。由 于周期信號(hào)的頻譜是離散的，故在繪制頻譜時(shí)，采用stem命令而不是plot命令3。 5.1 實(shí)現(xiàn)流程 采用三角形式傅里葉級(jí)數(shù)分解形式，用式求出傅里葉級(jí)數(shù)分解系數(shù)和，22 n a n b 再用式求出傅里葉復(fù)指數(shù)系數(shù)，并畫(huà)出的振幅和相位頻譜。諧波的階數(shù)nf 可32 n f n f 任意指定為60。 利用matlab實(shí)現(xiàn)上述分析過(guò)程的流程如下： （1）編寫(xiě)子函數(shù)y=time_fun_s(t)，用符號(hào)表達(dá)式表示出周期信號(hào)在第一個(gè)周期內(nèi)的 符號(hào)表達(dá)式，并賦值給返回符號(hào)變量y。 （2）編寫(xiě)子函數(shù)x=time_fun_e，求出該周期信號(hào)在繪圖區(qū)間內(nèi)的信號(hào)樣值，并賦值 給返回變量x。 （3）編寫(xiě)求解信號(hào)傅里葉復(fù)指數(shù)系數(shù)及繪制頻譜圖的通用函數(shù)，該函數(shù)的流程如 n f 下： 調(diào)用函數(shù)time_fun_s(t)，獲取周期信號(hào)的符號(hào)表達(dá)式。 求出信號(hào)的三角級(jí)數(shù)形式的傅里葉級(jí)數(shù)展開(kāi)系數(shù)和。 n a n b matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 11 求出信號(hào)的復(fù)指數(shù)形式的傅里葉級(jí)數(shù)展開(kāi)系數(shù)。 n f 繪制的振幅頻譜圖和相位頻譜圖。 n f 調(diào)用函數(shù)time_fun_e，繪制信號(hào)波形圖。 5.2 matlab 算法提示 （1）采用符號(hào)積分求內(nèi)時(shí)間函數(shù)的三角級(jí)數(shù)展開(kāi)系數(shù)：，int0,t 0 a 0 a n aas 。即計(jì)算式。 n bbs22 （2）用循環(huán)語(yǔ)句求出三角級(jí)數(shù)展開(kāi)系數(shù)，的數(shù)值：并賦值給變量.for end n a n b ，如圖5.2所示。1sym k _1bsym k 圖5.2 三角級(jí)數(shù)展開(kāi)系數(shù)，頻譜（n=k+1：nf） n a n b （3）從三角級(jí)數(shù)展開(kāi)系數(shù)，得到復(fù)指數(shù)展開(kāi)系數(shù)為了得到復(fù)指數(shù)展開(kāi)系數(shù) n a n b n f 必須先求出三角形式的傅里葉級(jí)數(shù)展開(kāi)系數(shù)和，如圖5.2所示，再據(jù)式求出 n f n a n b32 。但要注意、和的自變量取值情況，即： n f n a n b n f 從上式的自變量取值情況及圖5.2可見(jiàn)，三角級(jí)數(shù)展開(kāi)系數(shù)，的變量的取值范 n a n bn 圍為，而指數(shù)形式展開(kāi)系數(shù)的變量的取值范圍為,為0,1,2,3,.,nn n fn0, 1, 2,.,nn matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 12 了從和得到，需要用到matlab的反折函數(shù)fliplr來(lái)實(shí)現(xiàn)頻譜的反折。 n a n b n f 雙邊、單邊（幅度，相位）頻譜分析源程序見(jiàn)附錄程序二。 5.3 程序運(yùn)行結(jié)果及分析 5.3.1 程序運(yùn)行結(jié)果 調(diào)用ctfs2.m函數(shù)文件，即可繪出周期矩形波信號(hào)的頻譜。指令如下：在matlab 命令窗口鍵入ctfs2，并回車(chē)，命令窗口將出現(xiàn)：pleas input 所需展開(kāi)的最高諧波次數(shù) nf=？，輸入nf=60 然后命令窗口將出現(xiàn)：pleas input 信號(hào)的周期t=？，輸入周期t=5， 然后命令窗口又出現(xiàn)：周期與脈沖寬度之比m=？，輸入m =5，即可繪出周期矩形波信號(hào) 的雙邊頻譜，如圖5.3所示。修改繪制周期矩形脈沖雙邊頻譜的函數(shù)文件ctfs2.m，將其 改為繪制周期矩形脈沖單邊頻譜的函數(shù)文件ctfs21.m，則可繪出單邊頻譜如圖5.4所示。 圖5.3 周期t=5，脈寬tao=1 的矩形脈沖及雙邊頻譜圖 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 13 圖5.4 周期t=5，脈寬tao=1 的矩形脈沖及單邊頻譜圖 由圖5.3和5.4可見(jiàn)，周期矩形脈沖的頻譜具有一般周期信號(hào)頻譜的共同特點(diǎn)，即它們 的頻譜都是離散的。僅含有的各分量，其相鄰兩譜線(xiàn)的間隔是，同時(shí)單n 2t 邊譜僅為雙邊譜的一半波形。 5.3.2 脈沖寬度與頻譜的關(guān)系 重復(fù)調(diào)用周期矩形脈沖單邊頻譜的函數(shù)文件，當(dāng)窗口出現(xiàn)輸入信號(hào)諧波次數(shù)nf，周 期t，周期與脈沖寬度之比m時(shí)，分別將信號(hào)設(shè)置為表5-1所示三種情況，則可繪出表5-1 所示各情況的信號(hào)波形及頻譜圖，如圖5.5所示。 表5-1 周期矩形脈沖脈寬取值參數(shù)表 周期t555 脈沖寬度tao 被展開(kāi)函數(shù)的時(shí)間 區(qū)間的左端a=tao/2 0.6250.31250.15625 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 14 a 周期t=5，脈沖寬度tao=t/4=1.25，a=tao/2=0.625 的矩形脈沖頻譜 b 周期t=5，脈沖寬度tao=t/8=0.625, a=tao/2=0.3125 的矩形脈沖頻譜 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 15 c 周期t=5，脈沖寬度tao=t/16=0.3125, a=tao/2=0.15625 的矩形脈沖頻譜 圖5.5 脈沖寬度與頻譜的關(guān)系 由圖5.5可見(jiàn)，由于周期t相同，因而相鄰譜線(xiàn)的間隔相同；脈沖寬度tao愈窄，其頻 譜包絡(luò)線(xiàn)第一個(gè)零點(diǎn)的頻率愈高，即信號(hào)帶寬愈寬，頻帶內(nèi)所含的分量愈多?？梢?jiàn)，信 號(hào)的頻帶寬度與脈沖寬度tao成反比2。 5.3.3 脈沖周期與頻譜的關(guān)系 重復(fù)調(diào)用繪制周期矩形脈沖單邊頻譜的函數(shù)文件，當(dāng)窗口出現(xiàn)輸入信號(hào)諧波次數(shù) nf，周期t，周期與脈沖寬度之比m 時(shí)，分別將信號(hào)設(shè)置為表5-2所示各種情況（即信號(hào) 時(shí)域?qū)挾缺３植蛔儯瑒t可繪出表5-2所示各情況下的信號(hào)波形及頻譜圖如圖5.6所示。 表5-2 周期矩形脈沖周期取值參數(shù)表 脈沖寬度tao111 周期t 被展開(kāi)函數(shù)的時(shí)間 區(qū)間的左端a=tao/2 0.50.50.5 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 16 a 脈沖寬度tao=1，周期t=4*tao =4的矩形脈沖頻譜 b 脈沖寬度tao=1，周期t=8*tao =8 的矩形脈沖頻譜 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 17 c 脈沖寬度tao=1，周期t=16*tao =16 的矩形脈沖頻譜 圖5.6 周期與頻譜的關(guān)系 由圖5.6可見(jiàn)，由于周期脈沖信號(hào)的時(shí)域?qū)挾炔蛔?，這時(shí)頻譜包絡(luò)線(xiàn)的零點(diǎn)所在位置 不變，而當(dāng)周期增長(zhǎng)時(shí)，相鄰譜線(xiàn)的間隔減少，頻譜變密。如果周期無(wú)限增長(zhǎng)（這時(shí)就 成為非周期信號(hào)），那么，相鄰譜線(xiàn)的間隔將趨近于零，周期信號(hào)的離散頻譜就過(guò)渡到 非周期信號(hào)的連續(xù)頻譜。隨著周期的增長(zhǎng)，各諧波分量的幅度也相應(yīng)減少。脈沖周期t 愈長(zhǎng)，譜線(xiàn)間隔愈小，頻譜越稠密；反之，則越稀疏2。 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 18 6 典型周期脈沖的頻譜 6.1 周期方波脈沖頻譜的 matlab 實(shí)現(xiàn) 周期方波脈沖信號(hào)如圖6.1所示，其幅度為1，脈沖寬度占空比：duty=1/2，周期 t=5。 圖6-1 周期方波脈沖 編寫(xiě)ctfs3.m函數(shù)文件，源程序文件見(jiàn)附錄程序三。調(diào)用函數(shù)ctfs3.m，即可繪出方 波脈沖的雙邊頻譜如圖6.2所示，其中周期t和占空比duty可變，修改程序即可得到單邊頻 譜如圖6.3所示。 a 周期為t=5，占空比duty=50的方波脈沖雙邊頻譜 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 19 b 周期為t=10，占空比duty=80的方波脈沖雙邊頻譜 圖6.2 方波脈沖雙邊幅度譜 a 周期為t=5，占空比duty=50的方波脈沖單邊頻譜 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 20 b 周期為t=10，占空比duty=80的方波脈沖單邊頻譜 圖6.3 方波脈沖單邊幅度譜 由圖可以看出，周期方波信號(hào)頻譜與周期矩形脈沖信號(hào)具有相同的規(guī)律，由于方波 的周期與脈寬比，因此頻譜的第一個(gè)過(guò)零點(diǎn)內(nèi)只有兩根譜線(xiàn)。2 t 6.2 周期三角波脈沖頻譜的 matlab 實(shí)現(xiàn) 周期三角波脈沖如圖6.4所示，周期t=5，其幅度為1。 圖6.4 周期三角波脈沖 matlab 內(nèi)置有產(chǎn)生三角波的函數(shù)sawtooth(t)，其調(diào)用格式為:x= sawtooth(t,width): 根據(jù)width值的不同產(chǎn)生不同形狀的三角波,參數(shù)width 是01 之間的標(biāo)量,指定在一個(gè)周 期之間最大值的位置,width是該位置的橫坐標(biāo)和周期的比值.因而,當(dāng)width=0.5 時(shí)產(chǎn)生標(biāo)準(zhǔn) 的對(duì)稱(chēng)三角波，當(dāng)width時(shí)（可缺?。┊a(chǎn)生鋸齒波。 編寫(xiě)ctfs4.m函數(shù)文件，源程序文件見(jiàn)附錄程序四。調(diào)用函數(shù)ctfs4.m，即可繪出方 波脈沖的雙邊頻譜如圖6.5所示，其中周期t和width可變，修改程序即可得到單邊頻譜如 圖6.6所示。 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 21 a 周期t=5，width=0.5的三角脈沖雙邊頻譜 b 周期t=10，width=1的三角脈沖雙邊頻譜 圖6.5 周期三角波雙邊幅度頻譜 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 22 a 周期t=5，width=0.5的三角脈沖單邊頻譜 b 周期t=10，width=1的三角脈沖單邊頻譜 圖6.6周期三角波單邊幅度頻譜 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 23 結(jié)束語(yǔ) 本次課程設(shè)計(jì)至此已經(jīng)接近尾聲，一周的時(shí)間雖然很短暫，但在這一個(gè)星期的設(shè)計(jì) 過(guò)程中收獲頗多。設(shè)計(jì)的核心內(nèi)容就是利用 matlab 強(qiáng)大的圖形處理功能，符號(hào)運(yùn)算功 能以及數(shù)值計(jì)算功能，實(shí)現(xiàn)連續(xù)時(shí)間周期信號(hào)頻域分析的仿真波形。整個(gè)設(shè)計(jì)過(guò)程中首 先對(duì)所學(xué)的信號(hào)與系統(tǒng)與數(shù)字信號(hào)處理有了更深的了解，比如傅立葉級(jí)數(shù)、信號(hào)頻譜等； 其次，實(shí)現(xiàn)過(guò)程是通過(guò) matlab 軟件完成的，matlab 的圖形功能強(qiáng)大，具有良好的 人機(jī)界面5，此次設(shè)計(jì)過(guò)程中熟練了 matlab 的編程，掌握了很多函數(shù)的作用及使用方 法；最后，通過(guò)此次課程設(shè)計(jì)，我對(duì)設(shè)計(jì)所用到的軟件 matlab 有了更加深刻地了解， matlab 不管在數(shù)值計(jì)算方面的功能很強(qiáng)大，而且其圖形仿真功能更能滿(mǎn)足各個(gè)領(lǐng)域的 需要，因此我們以后更要經(jīng)常運(yùn)用 matlab 軟件，使其成為自己不可或缺的工具。 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 24 致謝 感謝學(xué)校給我們這次matlab課程設(shè)計(jì)的機(jī)會(huì)，不僅讓我們更加學(xué)會(huì)了matlab的 強(qiáng)大圖形處理方法，掌握了matlab的編程技術(shù)，而且也鍛煉了我們的動(dòng)手能力。通過(guò) 這次課設(shè)讓我明白了理論聯(lián)系實(shí)踐的重要性，書(shū)本上的理論知識(shí)學(xué)了不少，我們必須得 應(yīng)用到實(shí)踐當(dāng)中，做到學(xué)以致用，這樣我們才能有不斷的創(chuàng)新。這次課程設(shè)計(jì)也感謝指 導(dǎo)老師在設(shè)計(jì)過(guò)程中的輔導(dǎo)以及同學(xué)們的幫助。 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 25 參考文獻(xiàn) 1 陳懷琛，吳大正，高西全.matlab 及在電子信息課程中的應(yīng)用z. 北京：電子工業(yè) 出版社，2005 2 劉泉，江雪梅.信號(hào)與系統(tǒng)z. 北京：高等教育出版社，2006 3 劉泉，闕大順，郭志強(qiáng).數(shù)字信號(hào)處理原理與實(shí)現(xiàn)z. 北京：電子工業(yè)出版社，2009 4 梁虹. 信號(hào)與系統(tǒng)分析及 matlab 實(shí)現(xiàn)z. 北京：電子工業(yè)出版社，2002 5 羅建軍. matlab 教程z. 北京：電子工業(yè)出版社，2005 matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 26 附錄 源程序一 周期信號(hào)的傅里葉級(jí)數(shù)分解與綜合的源程序如下： functiona_sym,b_sym=ctfs1 syms t n k x t=5;tao=0.2*t;a=0.5; if nargin4;nf=6;end if nargin5;nn=32;end x=time_fun_x(t); a0=2*int(x,t,-a,t-a)/t; 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)=double(vpa(a0,nn); for k=1:nf a_sym(k+1)=double(vpa(subs(as,n,k),nn); b_sym(k+1)=double(vpa(subs(bs,n,k),nn); end if nargout=0 c=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*pi*6*t/5); f6=f1+f2; f7=f6+f3; f8=f7+f4+f5; subplot(2,2,1) plot(t,f1),hold on y=time_fun_e plot(t,y,:) title(周期矩形波的形成基波) subplot(2,2,2) plot(t,f6),hold on y=time_fun_e plot(t,y,:) title(周期矩形波的形成基波+2次諧波) subplot(2,2,3) plot(t,f7),hold on y=time_fun_e matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 27 plot(t,y,:) title(周期矩形波的形成基波+2次諧波+3次諧波) subplot(2,2,4) plot(t,f8),hold on y=time_fun_e plot(t,y,:) title(周期矩形波的形成基波+2次諧波+3次諧波+6次諧波) end function x=time_fun_x(t) h=1; x1=sym(heaviside(t+0.5)*h; x=x1-sym(heaviside(t-0.5)*h; function y=time_fun_e a=0.5;t=5;h=1;tao=0.2*t;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); 源程序二 雙邊、單邊（幅度，相位）頻譜分析源程序如下： 雙邊頻譜分析源程序： function a_sym,b_sym=ctfs2 syms t n y if nargin3;nf=input(pleas input 所需展開(kāi)的最高諧波次數(shù): nf=);end t=input(pleas input 信號(hào)的周期=); if nargin5;nn=32;end y=time_fun_s(t); 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); for k=1:nf a_sym(k+1)=double(vpa(subs(as,n,k),nn); b_sym(k+1)=double(vpa(subs(bs,n,k),nn);end if nargout=0 s1=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)=0 s4=fliplr(s3) s5=s2-i*s4; s6=fliplr(s5); matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 28 n=nf*2*pi/t; k2=-n:2*pi/t:n; s7=s6,s5(2:end); x=time_fun_e subplot(3,1,2) stem(k2,abs(s7); title(連續(xù)時(shí)間函數(shù)周期矩形脈沖的雙邊幅度譜) axis(-80,80,0,0.12) line(-80,80,0,0) line(0,0,0,0.12) subplot(3,1,3) stem(k2,angle(s7); title(連續(xù)時(shí)間函數(shù)周期矩形脈沖的雙邊相位譜) axis(-80,80,-4,4) line(-80,80,0,0) line(0,0,-4,4) end function y=time_fun_s(t) syms a a1 t=input(pleas input 信號(hào)的周期t=); m=input(周期與脈沖寬度之比m=); a=1;tao=t/m;a=tao/2; y1=sym(heaviside(t+a1)*a; y=y1-sym(heaviside(t-a1)*a; y=subs(y,a1,a); y=simple(y); function x=time_fun_e t=5;t=-2*t:0.01:2*t;tao=t/5; x=rectpuls(t,1); subplot(3,1,1) plot(t,x) hold on x=rectpuls(t-5,1); plot(t,x) hold on x=rectpuls(t+5,1); plot(t,x) title(周期為t=5,脈寬tao=1的矩形脈沖) axis(-10,10,0,1.2) 單邊頻譜分析源程序： function a_sym,b_sym=ctfs21 syms t n k y matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 29 if nargin4;nf=input(pleas input 所需展開(kāi)的最高諧波次數(shù): nf=);end t=input(pleas input 信號(hào)的周期t=); if nargin5;nn=32;end y=time_fun_s(t); 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); for k=1:nf a_sym(k+1)=double(vpa(subs(as,n,k),nn); b_sym(k+1)=double(vpa(subs(bs,n,k),nn);end if nargout=0 s1=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)=0 s4=fliplr(s3) s5=s2-i*s4; n=nf*2*pi/t; k2=0:2*pi/t:n; x=time_fun_e subplot(3,1,2) stem(k2,abs(s5); title(連續(xù)時(shí)間函數(shù)周期矩形脈沖的單邊幅度譜) axis(0,80,0,0.12) line(0,80,0,0) line(0,0,0,0.12) subplot(3,1,3) stem(k2,angle(s5); title(連續(xù)時(shí)間函數(shù)周期矩形脈沖的單邊相位譜) axis(0,80,-5,5) line(0,80,0,0) line(0,0,-5,5) end function y=time_fun_s(t) syms a a1 t=input(pleas input 信號(hào)的周期t=); m=input(周期與脈沖寬度之比m=); a=1;tao=t/m;a=tao/2; y1=sym(heaviside(t+a1)*a; y=y1-sym(heaviside(t-a1)*a; y=subs(y,a1,a); y=simple(y); matlab 應(yīng)用實(shí)踐課程設(shè)計(jì) 30 function x=time_fun_e t=5;t=-2*t:0.01:2*t;tao=t/5; x=rectpuls(t,1); subplot(3,1,1) plot(t,x) hold on x=rectpuls(t-5,1); plot(t,x) hold on x=rectpuls(t+5,1); plot(t,x) title(周期為t=5,脈