1、Matlab應用實踐課程設計課程設計任務書學生姓名： 專業班級：指導教師： 工作單位:題目：連續時間信號和系統時域分析及 MATLAB實現初始條件：MATLAB 6.5要求完成的主要任務：一、用MATLAB實現常用連續時間信號的時域波形（通過改變參數，分析其時 域特性）。1、單位階躍信號，2、單位沖激信號，3、正弦信號，4、實指數信號，5、虛指 數信號，6、復指數信號。二、用MATLAB實現信號的時域運算1、相加，2、相乘，3、數乘，4、微分，5、積分三、用MATLAB實現信號的時域變換（參數變化，分析波形變化）1、反轉，2、使移（超時，延時），3、展縮，4、倒相，5、綜合變化四、用MATLA

2、B實現信號簡單的時域分解1、信號的交直流分解，2、信號的奇偶分解五、用MATLAB實現連續時間系統的卷積積分的仿真波形給出幾個典型例子，對每個例子，要求畫出對應波形。六、用MATLAB實現連續時間系統的沖激響應、階躍響應的仿真波形。給出幾個典型例子，四種調用格式。七、利用MATLAB實現連續時間系統對正弦信號、實指數信號的零狀態響應的 仿真波形。給出幾個典型例子，要求可以改變激勵的參數，分析波形的變化時間安排：學習MATLAB語言的概況第1天學習MATLAB語言的基本知識第2、3天學習MATLAB語言的應用環境，調試命令，繪圖能力第4、5天課程設計第6-9天答辯第10天指導教師簽名：年月日系主

3、任（或責任教師）簽名：年月日摘要I.ABSTRACT II緒論11 MATLAB 簡介21.1 MATLAB語言功能21.2 MATLAB 語言特點 22常用連續時間信號的時域波形 32.1 單位階躍信號 32.2 單位沖激信號32.3 正弦信號42.4 實指數信號52.5 虛指數信號52.6 復指數信號63連續時間信號的時域運算 73.1 相加73.2 相乘73.3 數乘83.4 微分83.5 積分94 連續時間信號的時域變換 104.1 反轉104.2 時移104.3 展縮114.4 倒相114.5 綜合變化125連續時間信號簡單的時域分解 135.1 信號的交直流分解 135.2 信號的

4、奇偶分解146連續時間系統的卷積積分的仿真波形 157連續時間系統的沖激響應、階躍響應的仿真波形 167.1 IMPULSE ()函數177.2 STEP () 函數 198連續時間系統對正弦信號、實指數信號的零狀態響應的仿真波形 218.1 正弦信號的零狀態響應 218.2 實指數信號的零狀態響應 229小結即心得體會24致謝25參考文獻26附錄27Matlab應用實踐課程設計摘要MATLAB目前已發展成為由 MATLAB 語言、MATLAB 工作環境、MATLAB 圖形 處理系統、MATLAB數學函數庫和MATLAB應用程序接口五大部分組成的集數值計 算、圖形處理、程序開發為一體的功能強大

5、的系統。 本次課程設計則在深入研究連續時間 信號傅里葉級數分析理論知識的基礎上，利用MATLAB強大的圖形處理功能、符號運算功能以及數值計算功能，通過 MATLAB編程進行圖形功能仿真，從而實現連續時間周期 信號頻域分析的仿真波形，包括以下內容：用 MATLAB實現常用連續時間信號的時域波 形；用MATLAB實現信號的時域運算；用 MATLAB實現信號的時域變換；用 MATLAB 實現信號簡單的時域分解；用 MATLAB實現連續時間系統的卷積積分的仿真波形；用 MATLAB實現連續時間系統的沖激響應、階躍響應的仿真波形；用MATLAB實現連續時 問系統對正弦信號、實指數信號的零狀態響應的仿真波

6、形。關鍵詞：MATLAB ;圖形處理；連續時間信號和系統；時域2AbstractMATLAB 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, progra

7、m 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 function

8、al simulation MATLAB graphical programming in order to achieve continuous time periodic signal frequency domain analysis of the simulation waveforms, including the following: Time domain waveform of continuous time signal by MATLAB; time domain operation signal by MATLAB; realize the time domain sig

9、nal by MATLAB; MATLAB time domain signal simple decomposition; simulation waveform of convolution integrals of continuous time system with MATLAB; MATLAB impact simulation waveform impulse response, the step response of the continuous time system to achieve zero state; simulation waveform of sine si

10、gnal, the real exponential signal responseof continuous time system with MATLAB.Keywords: MATLAB; image processing; continuous time signals and systems; time domainMatlab應用實踐課程設計緒論在科學技術飛速發展的今天，計算機正逐步將科技人員從繁重的計算工作中解脫出 來。在進行科學研究與工程應用中，往往需要大量的科學計算，一些科技人員曾經嘗試使 用傳統的高級語言Basic、Fortran及C語言編寫程序，以減輕工作量。但編制程序需

11、要掌 握高級語言的語法，還要對各種算法進行了解，這對大多數科技人員來說是不大現實的， 而且也是沒有沒有必要的。MATLAB正是在這一應用要求背景下產生的數學類科技應用 軟件。它具有的頂尖的數值計算功能、強大的圖形可視化功能及簡潔易學的“科學便捷式” 工作環境和編程語言，從根本上滿足了科技人員對工程數學計算的要求，并將科技人員從 繁重的數學運算中解放出來，因而越來越受到廣大科技工作者的普遍歡迎1。MATLAB 是matrix和laboratory前三個字母的縮寫，意思是“矩陣實驗室”，是 MathWorks公司推出的數學類科技應用軟件。其Dos版本(MATLAB 1.0 )發行于1984年， 現

12、已推出了 Windows版本(MATLAB 5.3 )。經過十多年的不斷發展與完善，MATLAB 已 發展成為由MATLAB語言、MATLAB工作環境、MATLAB 圖形處理系統、MATLAB數 學函數庫和MATLAB應用程序接口五大部分組成的集數值計算、圖形處理、程序開發為 一體的功能強大的系統。MATLAB由“主包”和三十多個擴展功能和應用學科性的工具 箱(Toolboxs)組成。目前，MATLAB已經成為國際上最流行的電子仿真計算機輔助設計的軟件工具，現在 的MATLAB已經不僅僅是一個“矩陣實驗室(Matrix Laboratory) ”，它已經成為一種實 用的、全新的計算機高級語言。

13、正是由于MATLAB在數值計算及符號計算等方面的強大功能，使 MATLAB 一路領 先，成為數學類科技應用軟件中的佼佼者。目前， MATLAB已成為國際上公認的最優秀 的科技應用軟件。MATLAB的上述特點，使它深受工程技術人員及科技專家的歡迎，并 很快成為應用學科計算機輔助分析、設計、仿真、教學等領域不可缺少的基礎軟件。61 MATLAB 簡介1.1 MATLAB語言功能MATLAB是一個高精度的科學計算語言，它將計算、可視化編程結合在一個容易使用 的環境中，在這個環境中，用戶可以把提出的問題和解決問題的辦法用熟悉的數學符號表 示出來，它的典型使用包括：(1)數學和計算；(2)運算法則；(3

14、)建模、仿真；(4)數值分析、研究和可視化；(5)科學的工程圖形；(6)應用程序開發，包括創建圖形用戶接口。1.2 MATLAB語言特點MATLAB 是一個交互式系統，他的基本數據單元是數組，這個數組不要求固定的大 小，因此可以讓用戶解決許多技術上的問題，特別是那些包含矩陣和矢量運算的問題。MATLAB的指令表達與數學、工程中常用的習慣形式相似，與 C、Fortran、等高級語言相 比，它的語法規則更簡單、表達更符合工程習慣，正因為如此，人們用 MATLAB語言編寫 程序就猶如在便箋上書寫公式和求解，因而 MATLAB被稱為“便箋式”的科學工程語言。MATLAB的最重要特征使他擁有解決特定應用

15、問題的程序組， 也就是TOOLBOX(工具 箱)，如信號處理工具箱，控制系統工具箱、神經網絡工具箱、模糊邏輯工具箱、通信工具 箱和數據采集工具箱等許多專用工具箱，對大多數用戶來說，要想靈活、高效地運用這些 工具箱，通常都需要學習相應的專業知識。止匕外，開放性也許是MATLA最重要和最受歡迎的特點之一。除內部函數外，所有的 MATLAB主要文件和各工具箱文件都是可讀的、可改的源文件，因為工具箱實際上是有一組復雜的MATLAB函數(M文件)組成，它擴展了 MATLAB的功能，用以解決待定的問題， 因此用戶可以通過對源文件進行修改和加入自己編寫的文件去構建新的專用工具箱。2常用連續時間信號的時域波形

16、連續信號又稱為模擬信號，具信號存在于整個時間范圍內，包括單位階躍信號，單 位沖激信號，正弦信號，實指數信號，虛指數信號，復指數信號。2.1 單位階躍信號單位階躍信號的定義如下:Q (t < 0)L(t > 0)單位階躍信號的Matlab實現程序見附錄，其信號圖如下:圖2.1單位階躍信號2.2 單位沖激信號在連續時間系統中，單位沖激是一種重要的信號。任何一種模擬信號都能通過沖激給 予近似，通過系統對沖激輸入的響應可以求的所有其他輸入信號的響應。單位沖激信號8 (t)也稱為狄拉克(Dirac)分布，定義如下：/ K.1 5(0 = 1(-06(火。=1工,0第一個條件表明8 在所有t不

17、為0是取值為0,第二個條件是沖擊下的面積為1,因 此8 信號具有單位面積的特性。特別需要指出的是，s ce在t=o點的值5（。）是沒有定義的，& 并不等于無窮。 沖激信號6 （e可以近似地用一個位于原點、幅度為 A、持續時間為1/A的脈沖來表示，這 里A是一個很大的正值。t=i/A=i/50,單位脈沖5 0）的Matlab實現程序見附錄，其信號圖如下:圖2.2單位沖激信號2.3 正弦信號正弦信號和余弦信號二者僅在相位上相差九/2 ,經常統稱為正弦信號，一般寫作加）-A cos（2jrft +或f（t） = A cos（wt+中）cos幅度A=3,頻率f=5,相移 = 1的正弦信號其Ma

18、tlab實現程序見附錄，其信號圖如下:圖2.3正弦信號2.4 實指數信號實指數信號可由下面的表達式來表示：f(t) = AeQt式中e是自然數2.718，a和A是實數。若a>0,信號將隨時間而增長，若a<0,信 號將隨時間而衰減，若a=0,信號不隨時間而變化，成為直流信號。常數 A表示指數信號 在t=0點的初始值。A=3, a=0.5的實指數信號，什)=3/元的Matlab實現程序見附錄，其信號圖如下：圖2.4實指數信號2.5 虛指數信號虛指數信號可由下面的表達式來表示：f(t) = A/.2.6 復指數信號復指數信號可由下面的表達式來表示：f(t) = AeCtt+itd)tA=

19、1 , a=-l , & = 10的復指數信號網分=a/一通的Matlab實現程序見附錄，其信 號圖如下：圖2.6復指數信號Matlab應用實踐課程設計在信號的傳輸和處理過程中往往需要進行信號的運算，它包括信號的相加、相乘、數乘、微分、積分。3.1 相加要實現兩信號的相加,即f (t) =fl (t) +f2 (t)fl (t)為單位階躍信號，f2 (t)為正弦信號，兩信號相加的實現程序加附錄，其信號 圖如下：圖3.1相加3.2 相乘要實現兩信號的相乘,即f (t) =fl (t) *f2 (t)fl (t)為單位階躍信號，f2 (t)為正弦信號，兩信號相乘的實現程序加附錄，其信號 圖

20、如下：3.3 數乘要實現信號的數乘，即f (t) =A*(t)A=2, fi (t)為單位階躍信號，信號數乘的實現程序加附錄，其信號圖如下:圖3.3數乘3.4 微分微分即求信號的導數。對函數f (t) =t2求一階微分的實現程序見附錄，其信號圖如下:圖3.4微分93.5積分Matlab應用實踐課程設計4連續時間信號的時域變換4.1 反轉信號的反轉就是將信號的波形以某軸為對稱軸翻轉180?,將信號f (t)中的自變量t替換成-t即可得到其反轉信號。信號f (t) =t的反轉實現程序見附錄，其信號圖如下：4.2 時移實現連續時間信號的時移即f (t-to)或者f (t+t0),常數t0>0。

21、正弦信號的時移實現程序見附錄，其信號圖如下：4.3 展縮信號的展縮即將信號f (t)中的自變量t替換為at, aw0 正弦信號的展縮實現程序見附錄，其信號圖如下：圖4.3展縮4.4 倒相連續信號的倒相是指將信號f(t)以橫軸為對稱軸對折得到-f(t) 正弦信號的展縮實現程序見附錄，其信號圖如下：圖4.4倒相11Matlab應用實踐課程設計4.5 綜合變化將f(t尸sin(t)/t通過反褶、移位、尺度變換由f的波形得到f(-2t+3)的波形。該變化的 實現程序見附錄，其信號圖如下：36圖4.5綜合變化5連續時間信號簡單的時域分解5.1 信號的交直流分解信號的交直流分解即將信號分解成直流分量和交流

22、分量兩部分之和，其中直流分量定義為fD (t) =f/ (t) /t交流分量定義為fA (t) =f (t) -fD (t)例如對函數f (t) =sin (t) +2進行交直流分解。MATLAB命令見附錄，分解波形圖如圖5.1所示圖5.1信號的交直流分解5.2 信號的奇偶分解信號的奇偶分解即將信號分解成偶分量和奇分量兩部分之和，偶分量定義為fe (t) =fe (-t)奇分量定義為fo (t) =- fo (- t)則任意信號f (t)可寫成f (t) = -+-/(£)-£i£i上式第一部分是偶分量，第二部分是奇分量，即倔=加。+/(-明2泌人一切2例如對函數

23、f (t) =sin (t-0.1) +t進行交直流分解MATLAB命令見附錄，分解波形圖如圖5.2所示圖5.2奇偶分解6連續時間系統的卷積積分的仿真波形卷積積分在信號與線形系統分析中具有非常重要的意義，是信號與系統分析的基本方法之一。連續時間信號fl(t)和f2(t)的卷積積分(簡稱為卷積)f定義為：f(t尸 fl(t)* f2(t)= -1(t)f2(t- ).：.由此可得到兩個與卷積相關的重要結論，即是：(D f(t)= fi(t)* fi(t),即連續信號可分解為一系列幅度由f (t)決定的沖激信號5。)及其平移信號之和；(2)線形時不變連續系統，設其輸入信號為f，單位響應為h (t

24、),其零狀態響應為 y，則有：y (t ) = f (t) ?h (t)。用MATLAB實現連續信號f i(t)與f2(t)卷積的過程如下：(1)將連續信號f l(t)與f2(t)以時間問隔？進行取樣，得到離散序列f l(k?)和f2(k?)；(2)構造與fi(k?)和f2(k?)相對應的時間向量ki和k2 ；(3)調用conv()函數計算卷積積分f (t)的近似向量f (n?)；(4)構造f (n?)對應的時間向量k0卷積實現程序見附錄。實現程序如下:p=0.1；k1=0:p:2;f1=0.5*k1;k2=k1;f2=f1;f,k=sconv(f1,f2,k1,k2,p)例二:圖6.2例二實

25、現程序如下:p=0.1;k1=0:p:2;f1=rectpuls(k1-1,length(k1);k2=k1;f2=f1;f,k=sconv(f1,f2,k1,k2,p)7連續時間系統的沖激響應、階躍響應的仿真波形對于連續時間系統，求解系統的沖激響應 h(t)和階躍響應g(t)對我們進行連續系統的 分析具有非常重要的意義。MATLAB為用戶提供了專門用于求連續系統沖激響應和階躍 響應并繪制其時域波形的函數impulse ()和step ()。在調用impulse ()和step ()函數時，我們需要用向量來對連續時間系統進行分析。 設描述連續系統的微分方程為：Jiy(i)(t尸x(j)(t)則

26、我們可用向量A和B來表示該系統，即：A=A n,An-i,Ai,AoB=Bn,Bn-i,Bi,Bo注意，向量A和B的元素一定要以微分方程中時間求導的降幕次序來排列，且缺項 J要用0來補齊。例如，對微分方程 式。+知仕)+2y年)=fQt) + f(t),則表示該系統的 對應向量應為 A=1 3 2 , B=1 0 107.1 impulse ()函數函數impulse ()將繪出由向量 a和b表示的連續系統在指定時間范圍內的沖激響 應h(t)的時域波形圖，并能求出指定時間范圍內沖激響應的數值解。impulse ()函數有如下四種調用格式：(1) impulse(b,a):該調用格式以默認方式繪

27、出向量A和B定義的連續系統的沖激響應的時域波形。例如描述連續系統的微分方程為DIIIy(t)+5y(t)+6M = 3而+2f(t)運行如下MATLAB 命令：a=1 5 6;b=3 2;impulse(b,a);則繪出系統的沖激響應波形，如圖 7.1.1所示。圖7.1.1連續系統的沖激響應1(2) impulse(b,a,t):繪出系統在0t時間范圍內沖激響應的時域波形。對上例，若 運行命令impulse(b,a,10),則繪出系統在010秒范圍內沖激響應的時域波形，如圖7.1.2 所示.Tlime < s-e- c )圖7.1.2連續系統的沖激響應2(3) impulse(b,a,h

28、p:t2):繪出在t1t2時間范圍內，且以時間間隔 p均勻取樣的沖 激響應波形。對上例，若運行命令 impulse(b,a,1:0.1:2),則繪出12秒內，每隔0.1秒 取樣的沖激響應的時域波形，如圖 7.1.3所示。(4) y=impulse(b,a,ti：p:t2):不繪出波形，而是求出系統沖激響應的數值解。對上例， 若運行命令y=impulse(b,a,0:022),則運行結果為：y =3.00001.16040.3110-0.0477-0.1726-0.1928-0.1716-0.1383-0.1054-0.0777-0.0559圖7.1.3連續系統的沖激響應37.2 step ()

29、函數step ()函數可繪出連續系統的階躍響應 g(t)在指定時間范圍的時域波形并能求出其 數值解，和impulse ()函數一樣也有四種調用格式。(1) step(b,a):該調用格式以默認方式繪出向量 A和B定義的連續系統的階躍響應 的時域波形。例如描述連續系統的微分方程為|II州+5而通二3而十2陽運行如下MATLAB 命令：a=1 5 6;b=3 2;step(b,a);則繪出系統的階躍響應波形，如圖 7.2.1所示圖7.2.1連續系統的階躍響應 1(2) step(b,a,t):繪出系統在0t時間范圍內階躍響應的時域波形。對上例，若運 行命令step(b,a,10),則繪出系統在01

30、0秒范圍內階躍響應的時域波形，如圖7.2.2所示.圖7.2.2連續系統的階躍響應 2(3) step(b,a,1：p:t2):繪出在tit2時間范圍內，且以時間間隔p均勻取樣的階躍響 應波形。對上例，若運行命令 step(b,a,1:0.1:2)則繪出12秒內，每隔0.1秒取樣的階 躍響應的時域波形，如圖7.2.3所示。圖7.2.3連續系統的階躍響應 3(4) y=step(b,a,l：p：t2):不繪出波形，而是求出系統階躍響應的數值解。對上例，若 運行命令y=step(b,a,0:022),則運行結果為：y= 00.3930.5290.5500.5250.4880.4510.4200.39

31、60.3770.3648連續時間系統對正弦信號、實指數信號的零狀態響應的仿真波形MATLAB中的函數lsim ()能對微分方程描述的LTI連續時間系統的響應進行仿真。 該函數能繪制連續時間系統在指定的任意時間范圍內系統響應的時域波形圖，還能求出連續時間系統在指定的任意時間范圍內系統響應的數值解，函數 lsim ()的調用格式如下：lsim(b,a,x,t)在該調用格式中，a和b是由描述系統的微分方程系統決定的表示該系統的兩個行向 量。x和t則是表示輸入信號的行向量，其中t為表示輸入信號時間范圍的向量，x則是 輸入信號在向量t定義的時間點上的抽樣值。該調用格式將繪出向量b和a所定義的連續 系統在

32、輸入量為向量x和t所定義的信號時，系統的零狀態響應的時域仿真波形，且時間范圍與輸入信號相同。8.1正弦信號的零狀態響應描述某連續時間系統的微分方程為111r(t)+2r(t)+r(t)=e(t)+2e(t)當輸入信號為e(t)=sin 2腐£(。時，該系統的零狀態響應r (t)MATLAB命令如下：clc;a=1,2,1;b=1,2;p=0.5;t=0:p:5;x=sin(2*pi*t);lsim(b,a,x,t);hold on;p=0.2;t=0:p:5;x=sin(2*pi*t);lsim(b,a,x,t);p=0.01;t=0:p:5;x=sin(2*pi*t);lsim(b

33、,a,x,t);hold off;LIp_w-dLLhFLnear Simulation results圖8.1正弦信號的零狀態響應8.2實指數信號的零狀態響應描述某連續時間系統的微分方程為而+2r(t)+r(t)=“+2崩當輸入信號為二廣/時，該系統的零狀態響應r (t)MATLAB命令如下：clc;a=1,2,1;b=1,2;p=0.5;t=0:p:5;x=exp(-2*t);lsim(b,a,x,t);hold on;p=0.3;t=0:p:5;x=exp(-2*t);lsim(b,a,x,t);p=0.01;t=0:p:5;x=exp(-2*t);lsim(b,a,x,t);hold

34、off;Linear Sknuiatdon ResultsTime E等=圖8.2實指數信號的零狀態響應圖8.1、8.2中藍線、綠線、紅線分別代表 p=0.5、p=0.3、p=0.01o顯然可以看出，函 數lsim ()對系統響應進行仿真的效果取決于向量t的時間間隔的密集程度。圖8.1、8.2繪出了上述系統在不同抽樣時間問隔時函數lsim ()仿真的情況，可見抽樣時問問隔越小仿真效果越好。9小結即心得體會本次課程設計至此已經接近尾聲，一周的時間雖然很短暫，但在這一個星期的設計過 程中收獲頗多。設計的核心內容就是利用 MATLAB強大的圖形處理功能，符號運算功能 以及數值計算功能，實現連續時間周

35、期信號頻域分析的仿真波形。整個設計過程中首先對 所學的信號與系統與數字信號處理有了更深的了解，比如傅立葉級數、信號頻譜等；其次,實現過程是通過MATLAB軟件完成的，MATLAB的圖形功能強大，具有良好的人機界 面，此次設計過程中熟練了 MATLAB的編程，掌握了很多函數的作用及使用方法； 最后, 通過此次課程設計，我對設計所用到的軟件MATLAB有了更加深刻地了解，MATLAB不管在數值計算方面的功能很強大，而且其圖形仿真功能更能滿足各個領域的需要，因此我們以后更要經常運用 MATLAB軟件，使其成為自己不可或缺的工具。在寫相關源程序的時候，我還收索了大量的網站，在網上收索了很多關于MATL

36、AB的資料。在這個過程中我發現網上有很多有用的知識。以后應該多注意，充分合理的利用 網絡，通過網絡來學習東西。在收集資料的階段我復習了數字信號系統處理里的相關知識。 對以前的理論知識有了更進一步的認識和理解。通過這次課程設計我還對mathtype數學公式編輯器有了一定的了解，并且會用它編輯公式。對word也有了進一步的掌握。雖然我順利完成了課程設計的要求，但是我感覺到我對MATLAB的理解我掌握還停留在比較淺的層次。要想真正掌握它還需要繼續努力學習它。這次課程設計也使我明白了在知識的領域里我還有很多很多的不足，并且再一次的深深的體會到理論和實踐之間還有很到的差別。 在以后的學習中應該多多的注意

37、實踐知識的 訓練和積累。在以后的學習生活中要不斷的開拓自己的動手能力，不斷的訓練自己的動手能力。這次課程設計讓我深深的明白了自己以后該做什么，該怎么去做。致謝感謝學校給我們這次MATLAB課程設計的機會，不僅讓我們更加學會了 MATLAB的強 大圖形處理方法，掌握了 MATLAB的編程技術，而且也鍛煉了我們的動手能力。通過這次 課設讓我明白了理論聯系實踐的重要性，書本上的理論知識學了不少，我們必須得應用到 實踐當中，做到學以致用，這樣我們才能有不斷的創新。這次課程設計也感謝指導老師在 設計過程中的輔導以及同學們的幫助。沒有他們的幫助我不會那么快克服那些困難，也不 會這么快學到這么多的知識。參考

38、文獻1陳懷琛，吳大正，高西全.MATLAB及在電子信息課程中的應用Z.北京：電子工業出版社，20052劉泉，江雪梅.信號與系統Z.北京：高等教育出版社，20063劉泉，闕大順，郭志強.數字信號處理原理與實現Z.北京：電子工業出版社，20094梁虹.信號與系統分析及MATLAB實現Z.北京：電子工業出版社，20025羅建軍.MATLAB 教程Z.北京：電子工業出版社，20056施陽.MATLAB語言工具箱-ToolBox實用指南.西安：西北工業大學出版社，1999口 鄧微.MATLAB函數速查手冊.北京：人民郵電出版社，2010附錄階躍信號clc;t=-0.5:0.001:1;t0=0;u=st

39、epfun(t,t0);plot(t,u);axis(-0.5 1-0.2 1.2)沖激信號clc;t = -3:0.01:3;y = (t=0);plot(t,y);正弦信號clc;t=-0.5:0.001:1;A=3;f=5;fai=1;u=A*sin(2*pi*f*t+fai); plot(t,u)axis(-0.5 1-3.2 3.2)實指數信號clc;t=0:0.002:3;A=3;a=0.5;u=A*exp(a*t);plot(t,u)axis(-0.2 3.1 -0.2 14)虛指數信號clc;t=0:0.001:15;a=2;w=pi/4;z=a*exp(i*w*t);subp

40、lot(2,2,1),plot(t,real(z),axis(0, 15,-2.5,2.5),title('實部') subplot(2,2,3),plot(t,imag(z),axis(0,15,-2.5,2.5),title('虛部') subplot(2,2,2),plot(t,abs(z),axis(0,15,1.5,2.5),title('模') subplot(2,2,4),plot(t,angle(z),axis(0,15,-4,4),title('相角')復指數信號clc;t=0:0.01:3;a=-1;A=1 ；

41、b=10;z=A*exp(a+i*b)*t);subplot(2,2,1),plot(t,real(z),title('實部')subplot(2,2,3),plot(t,imag(z),title('虛部')subplot(2,2,2),plot(t,abs(z),title('模')subplot(2,2,4),plot(t,angle(z),title('相角')相加clc;t=-0.5:0.0001:2.5;t0=1;u=stepfun(t,t0);y=sin(2*pi*t);f=y+u;plot(t,f)axis(-0.

42、5 2.5 -1.5 2.5)相乘clc;t=0:0.0001:3;t0=1;u=stepfun(t,t0);y=sin(2*pi*t);f=u.*y;plot(t,f);axis(0 3 -1.5 1.5);數乘clc;t=0:0.0001:3;a=2;t0=1;u=stepfun(t,t0);f=a*u;plot(t,f);axis(0 3 0 2.5);微分clc;t=-1:0.02:1;g=t.*t；d=diff(g);subplot(211);plot(t,g,'-');subplot(212);plot(d,'-');積分clc;t=-1:0.2:1

43、;syms t;f=t*t;g=int(f);subplot(211);ezplot(f);subplot(212);ezplot(g);反轉clc;t=-1:0.2:1;f=t;g=fliplr(f);h=flipud(f); subplot(311);plot(t,f);axis(-1 1-1 1);title('原函數')； subplot(312);Plot(t,g)；axis(-1 1-1 1);title('左右反轉);subplot(313);plot(t,h);axis(-1 1 -1 1); title('上下反轉');時移clc;t=0:0.0001:2;y=sin(2*pi*t);y1=sin(2*pi*(t-0.2);plot(t,y,'-',t,y1);axis(0 2 -1.5 1.5);展縮clc;t=0:0.0001:2;a=2;y=sin(2*pi*t);y1=subs(y,t,a*t);subplot(211);ezplot(y);subplot(212);ezplot(yl);倒相clc;t=0:0.0001:2;y=sin(2*pi*