信號(hào)與系統(tǒng)課程設(shè)計(jì)報(bào)告---matlab在信號(hào)與系統(tǒng)中的使用

1、<p><b>　　《信號(hào)與系統(tǒng)》</b></p><p><b>　　——課程設(shè)計(jì)報(bào)告</b></p><p>　　報(bào)告題目： MATLAB在信號(hào)與系統(tǒng)中的使用</p><p>　　所在系部： 理學(xué)院 </p><p>　　所在專業(yè)：

2、 應(yīng)用物理 </p><p>　　所在班級(jí)： 物理092 </p><p>　　作者姓名： </p><p>　　作者學(xué)號(hào)： </p><p>　　指導(dǎo)教師： x

3、xxxx </p><p>　　完成時(shí)間： 2011.6 </p><p><b>　　目錄</b></p><p>　　內(nèi)容摘要 ———————————————————————— 3</p><p>　　MATLAB簡(jiǎn)介 ————————————————

4、——————— 3</p><p>　　第一章 MATLAB應(yīng)用功能簡(jiǎn)介</p><p>　　1.1 基本概念——————————————————————— 5</p><p>　　1.2 基本運(yùn)算——————————————————————— 6</p><p><b>　　1.3 可視化功能</b></p&

5、gt;<p>　　1.3.1 繪制二維圖形——————————————————— 7</p><p>　　1.3.2 繪制三維圖形——————————————————— 8</p><p><b>　　1.4 實(shí)例簡(jiǎn)介</b></p><p>　　1.4.1 應(yīng)用基礎(chǔ)————————————————————— 9</p&g

6、t;<p>　　1.4.2 繪制二維圖形——————————————————— 11</p><p>　　1.4.3 繪制三維圖形——————————————————— 12</p><p>　　第二章 信號(hào)與系統(tǒng)分析及MATLAB實(shí)現(xiàn)</p><p>　　2.1 信號(hào)的時(shí)域分析</p><p>　　2.1.1 信號(hào)的表示

7、及可視化———————————————— 13</p><p>　　2.1.2 信號(hào)的時(shí)域運(yùn)算、時(shí)域變換及MATLAB實(shí)現(xiàn)———————13</p><p>　　2.1.3 用MATLAB分析常用的時(shí)間信號(hào)特性——————————14</p><p>　　2.2 連續(xù)系統(tǒng)的是與分析及MATLAB實(shí)現(xiàn)</p><p>　　2.2.1 離散時(shí)間序

8、列卷積及MATLAB實(shí)現(xiàn)———————————15</p><p>　　2.2.2 連續(xù)時(shí)間信號(hào)卷積及MATLAB實(shí)現(xiàn)———————————15</p><p><b>　　2.3 實(shí)例簡(jiǎn)介</b></p><p>　　2.3.1 信號(hào)與系統(tǒng)分析及MATLAB實(shí)現(xiàn)————————————16</p><p>　　2.3.

9、2 連續(xù)系統(tǒng)的是與分析及MATLAB實(shí)現(xiàn)——————————21</p><p>　　第三章 課程設(shè)計(jì)總結(jié)</p><p>　　3.1課程設(shè)計(jì)體會(huì) ————————————————————— 26</p><p>　　3.2評(píng)語及成績(jī) —————————————————————— 27</p><p>　　參考文獻(xiàn) —————————————

10、———————————— 28</p><p><b>　　內(nèi)容摘要</b></p><p>　　通過本次課程設(shè)計(jì)，我首先學(xué)習(xí)了matlab7.0的應(yīng)用基礎(chǔ)。在此基礎(chǔ)之上又練習(xí)使用了運(yùn)用matlab7.0來分析信號(hào)與系統(tǒng)。</p><p>　　我學(xué)習(xí)并掌握了使用matlab 7.0的基本功能，包括數(shù)值計(jì)算，符號(hào)運(yùn)算，圖形控制及運(yùn)算方法?！缎盘?hào)與

11、系統(tǒng)》課程中學(xué)到的那些常用函數(shù)，加上對(duì)他們性質(zhì)和變換的理解，為運(yùn)用matlab 7.0進(jìn)行信號(hào)與系統(tǒng)分析打下了基礎(chǔ)。以matlab 7.0為工具，對(duì)信號(hào)與系統(tǒng)在時(shí)域，頻域及Z 域進(jìn)行了系統(tǒng)分析與計(jì)算機(jī)模擬。通過大量應(yīng)用實(shí)例介紹運(yùn)用matlab 7.0進(jìn)行信號(hào)與系統(tǒng)分析的具體方法。</p><p>　　關(guān)鍵詞： matlab 7.0、信號(hào)與系統(tǒng)</p><p><b>　　MATL

12、AB簡(jiǎn)介</b></p><p>　　MATLAB是matrix和laboratory前三個(gè)字母的縮寫，意思是：“矩陣實(shí)驗(yàn)室”，是Mathwork公司推出的數(shù)學(xué)類科技應(yīng)用軟件。用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計(jì)算的高級(jí)技術(shù)計(jì)算語言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。MATLAB和Mathematica、Maple并稱為三大數(shù)學(xué)軟件。它在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計(jì)算

13、方面首屈一指。MATLAB可以進(jìn)行矩陣運(yùn)算、繪制函數(shù)和數(shù)據(jù)、實(shí)現(xiàn)算法、創(chuàng)建用戶界面連接其他編程語言的程序等，主要應(yīng)用于工程計(jì)算、控制設(shè)計(jì)、信號(hào)處理與通訊、圖像處理、信號(hào)檢測(cè)、金融建模設(shè)計(jì)與分析等領(lǐng)域。具有以下基本功能：（1）數(shù)值計(jì)算功能；（2）符號(hào)計(jì)算功能；（3）可視化建模及動(dòng)態(tài)仿真功能。</p><p>　　MATLAB部分常用函數(shù)：</p><p>　　eps：浮點(diǎn)相對(duì)精度 </

14、p><p>　　exp：自然對(duì)數(shù)的底數(shù)e </p><p>　　i 或 j：基本虛數(shù)單位 </p><p>　　inf 或 Inf：無限大， 例如1/0 </p><p>　　nan或NaN：非數(shù)值（Not a number），例如0/0 </p><p>　　pi：圓周率 p（= 3.1415926...） </p

15、><p>　　realmax：系統(tǒng)所能表示的最大數(shù)值 </p><p>　　realmin：系統(tǒng)所能表示的最小數(shù)值 </p><p>　　nargin: 函數(shù)的輸入引數(shù)個(gè)數(shù) </p><p>　　nargout: 函數(shù)的輸出引數(shù)個(gè)數(shù) </p><p>　　lasterr：存放最新的錯(cuò)誤信息 </p><

16、p>　　lastwarn：存放最新的警告信息 </p><p>　　sqrt(x)：開平方 </p><p>　　real(z)：復(fù)數(shù)z的實(shí)部 </p><p>　　imag(z)：復(fù)數(shù)z的虛部 </p><p>　　round(x)：四舍五入至最近整數(shù) </p><p>　　fix(x)：無論正負(fù)，舍去小數(shù)

17、至最近整數(shù) </p><p>　　sign(x)：符號(hào)函數(shù) (Signum function)。 </p><p>　　rem(x,y)：求x除以y的余數(shù) </p><p>　　gcd(x,y)：整數(shù)x和y的最大公因數(shù) </p><p>　　lcm(x,y)：整數(shù)x和y的最小公倍數(shù) </p><p>　　exp(x

18、) ：自然指數(shù) </p><p>　　pow2(x)：2的指數(shù) </p><p>　　log(x)：以e為底的對(duì)數(shù)，即自然對(duì)數(shù)或 </p><p>　　log2(x)：以2為底的對(duì)數(shù) </p><p>　　log10(x)：以10為底的對(duì)數(shù) </p><p>　　sin(x)：正弦函數(shù) </p><

19、p>　　cos(x)：余弦函數(shù) </p><p>　　tan(x)：正切函數(shù) </p><p>　　asin(x)：反正弦函數(shù) </p><p>　　acos(x)：反余弦函數(shù) </p><p>　　atan(x)：反正切函數(shù) </p><p>　　min(x): 向量x的元素的最小值 </p>

20、<p>　　max(x): 向量x的元素的最大值 </p><p>　　MATLAB基本繪圖函數(shù)： </p><p>　　plot: x軸和y軸均為線性刻度（Linear scale） </p><p>　　loglog: x軸和y軸均為對(duì)數(shù)刻度（Logarithmic scale） </p><p>　　fplot 較精確的函

21、數(shù)圖形 </p><p><b>　　stem 針狀圖 </b></p><p><b>　　fill 實(shí)心圖 </b></p><p>　　第一章 MATLAB應(yīng)用功能簡(jiǎn)介</p><p>　　1、MATLAB 的數(shù)值計(jì)算功能</p><p><b>　　1.1

22、 基本概念</b></p><p><b>　　變量 </b></p><p>　　和其它高級(jí)語言一樣，MATLAB 也是使用變量來保存信息。變量由變量名表示，變量 的命名應(yīng)遵循如下規(guī)則：</p><p>　?。ǎ保┳兞棵仨氁宰帜搁_頭；</p><p>　?。ǎ玻┳兞棵梢杂勺帜浮?shù)字和下劃線混合組成；&l

23、t;/p><p>　?。ǎ常┳兞棵麉^(qū)分字母大小寫；</p><p>　?。ǎ矗┳兞棵淖址L(zhǎng)度不應(yīng)超過 31 個(gè)。</p><p><b>　　函數(shù) </b></p><p><b>　?。ǎ保┍磉_(dá)式</b></p><p>　　（２）變量=表達(dá)式 表達(dá)式由變量名、常數(shù)、函數(shù)和運(yùn)

24、算符構(gòu)成。</p><p>　　語句結(jié)尾若不加“；”，則表示在語句執(zhí)行后，在將計(jì)算結(jié)果存入內(nèi)存的 同時(shí)，還將運(yùn)算結(jié)果顯示出來</p><p><b>　　矩陣 </b></p><p><b>　?。ǎ保┲苯虞斎敕?；</b></p><p>　?。ǎ玻├?MATLAB 內(nèi)部函數(shù)創(chuàng)建；</p&

25、gt;<p>　　（３）從外部數(shù)據(jù)文件（*.mat）裝載并創(chuàng)建矩陣。</p><p><b>　　向量</b></p><p>　?。?）利用冒號(hào) “ ：” 運(yùn)算生成向量 </p><p>　?。?）利用函數(shù) linspace ()生 成向量 </p><p>　?、賚inspace(m,n)</p&

26、gt;<p>　?、趌inspace(m,n,s)</p><p><b>　　命令</b></p><p>　　Clear 清除當(dāng)前工作空間全部變量</p><p>　　Clear a b c 清除當(dāng)前工作空間變量a、b、c</p><p>　　who 或 who

27、s 了解已賦值變量的有關(guān)信息</p><p>　　clc 擦除工作窗中所顯示的全部?jī)?nèi)容</p><p>　　clf 擦除當(dāng)前工作窗中的圖像</p><p><b>　　1.2 基本運(yùn)算</b></p><p><b>　　矩陣加減與數(shù)組加

28、減</b></p><p>　　（１）若參與運(yùn)算的兩矩陣（數(shù)組）的維數(shù)相同，則加減運(yùn)算的結(jié)果是將兩矩陣的對(duì)應(yīng)元 素進(jìn)行加減 </p><p>　　（２）若參與運(yùn)算的兩矩陣的之一為標(biāo)量（1×1 的矩陣），則加減運(yùn)算的結(jié)果是將矩陣（數(shù) 組）的每一元素與該標(biāo)量逐一相加減</p><p><b>　　矩陣乘與數(shù)組乘</b><

29、/p><p>　　矩陣乘與數(shù)組乘有著較大差別，運(yùn)算結(jié)果也完全不同。矩陣乘的運(yùn)算符為“*”，運(yùn)算是按矩陣的乘法規(guī)則進(jìn)行的，即參與乘運(yùn)算的兩矩陣的內(nèi)維必須相同。設(shè) A、 B 為參與乘 運(yùn)算的兩矩陣，C 為 A 和 B 矩陣乘的結(jié)果，則它們必須滿足關(guān)系 Cm×n=Am×k Bk×n。因此，參與運(yùn)算的兩矩陣的順序不能任意調(diào)換，因?yàn)?A*B 和 B*A 計(jì)算結(jié)果很可能是完全不一樣 的。數(shù)組乘的運(yùn)算

30、符為“.*”，運(yùn)算符中的點(diǎn)號(hào)千萬不能遺漏，也不能隨意加空格符。參與 數(shù)組乘運(yùn)算的兩數(shù)組的大小必須相等（即為同維數(shù)組）。</p><p>　　數(shù)組乘的結(jié)果是將兩同維數(shù)組（矩 陣）的對(duì)應(yīng)元素逐一相乘，因此，A.*B 和 B.*A 計(jì)算結(jié)果是完全相同的.由于，矩陣運(yùn)算和數(shù)組運(yùn)算的差異，能進(jìn)行數(shù)組乘運(yùn)算的兩矩陣，不一定能進(jìn)行矩陣乘運(yùn)算。</p><p><b>　　矩陣除與數(shù)組除<

31、/b></p><p>　　矩陣除與數(shù)組除也有著較大差別，運(yùn)算結(jié)果也完全不同。 矩陣除分為矩陣右除和矩陣左除兩種情況。矩陣右除的運(yùn)算符為“/”，設(shè) A、B 為兩矩陣，則“A/B”是指方程 X * B = A 的解矩陣 X。顯然，矩陣右除運(yùn)算對(duì)參與運(yùn)算的兩矩陣的維數(shù)是有一定要求的，即矩陣 A 和 B 的列數(shù)必須相等。矩陣右除允許參與右除運(yùn)算的矩陣 B 為標(biāo)量，這時(shí)矩陣右除運(yùn)算的結(jié)果是將矩陣 A 的，每一元素逐

32、一與該標(biāo)量進(jìn)行除法運(yùn)算矩陣左除的運(yùn)算符為“\”，設(shè) A、B 為兩矩陣，則“ A\ B”是指方程 B * X = A 的解矩陣 X。顯然，矩陣左除運(yùn)算對(duì)參與運(yùn)算的兩矩陣的維數(shù)是也有一定的要求，即矩陣 A 和 B的行數(shù)必須相等。</p><p>　　數(shù)組右除的運(yùn)算符為“. /”，左除的運(yùn)算符為“. \”。數(shù)組右除和左除的運(yùn)算結(jié)果是完全 等效的。設(shè) A、B 為兩同維矩陣 ，則“A. /B”的運(yùn)算結(jié)果是將矩陣 A 的每一

33、個(gè)元素與矩陣 B 的對(duì)應(yīng)元素相除。需要就意的是，參與數(shù)組運(yùn)算的兩矩陣（數(shù)組）的大小必須相等。</p><p><b>　　1.3可視化功能</b></p><p>　　1.3.1 繪制二維圖形</p><p>　　1、繪制簡(jiǎn)單的二維曲線</p><p>　　MATLAB 中最常用的繪圖函數(shù)為 plot。根據(jù) plot 函

34、數(shù)參數(shù)的不同，可以在平面上繪制 不同的曲線。plot 函數(shù)是將各個(gè)數(shù)據(jù)點(diǎn)通過連折線的方式來繪制二維圖形的，若對(duì)曲線細(xì) 分的話，曲線可以看成是由直線連接而成的。plot 命令的格式有以下幾種：</p><p>　　（1） plot(y) 當(dāng) y 為一向量時(shí)，以 y 的序號(hào)作為 X 軸，按向量 y 的值繪制曲線。</p><p>　?。?） plot(x,y) x,y 均為向量時(shí)，以 x

35、向量作為 X 軸，向量 y 作為 Y 軸繪制曲線。</p><p>　?。?） plot(x,y1,‘option’,x,y2,‘option’,…)</p><p>　?。?） plot(x1,y1,‘option’,x2,y2,‘option‘,…)</p><p>　　2、離散序列圖的繪制</p><p>　?。?）stem(y)&

36、lt;/p><p>　?。?）stem(x,y,‘option’)</p><p>　　（3）stem(x,y,‘filled’)</p><p><b>　　3、標(biāo)識(shí)坐標(biāo)軸名稱</b></p><p>　　通過 xlable（‘string’）和 ylable（‘string’）命令給 X 軸和 Y 軸加上標(biāo)注。title

37、（‘string’）命令給圖形加上標(biāo)題。使用 grid on 或 grid off 命令在所畫出的圖形中添加或去掉網(wǎng)絡(luò)線。使用坐標(biāo)軸確定文字位置的 text 命令；使用鼠標(biāo)確定文字位置的 gtext命令。</p><p>　　(1) text(x,y,‘string’,‘option’)</p><p>　　(2) gtext(‘string’)</p><p>&

38、lt;b>　　4、調(diào)整坐標(biāo)軸狀態(tài)</b></p><p>　　axis(‘string’)這個(gè)命令將坐標(biāo)軸的狀態(tài)調(diào)整為字符串‘string’的指定的狀態(tài)。</p><p>　　1.3.2 繪制三維圖形</p><p>　　MATLAB 提供了許多實(shí)現(xiàn)三維數(shù)據(jù)可視化的函數(shù)，可以在三維空間中繪制曲線或曲面。</p><p>　　

39、1、三維折線及曲線的基本繪圖命令</p><p>　?。?）plot3(x1,y1,z1,’option1‘,x2,y2,z2,‘option2’,…)</p><p>　?。?）axis(xmin xmax ymin ymax zmin zmax)</p><p>　　2、繪制三維網(wǎng)格曲面</p><p>　　(1) mesh(X,Y,Z

40、)命令繪制網(wǎng)格曲面</p><p>　　參數(shù) X，Y，Z 都是矩陣，X 矩陣的行向量相同，Y 矩陣的列向量相同。</p><p>　　(2) mesh(x,y,Z)命令繪制網(wǎng)格曲面</p><p>　　參數(shù) x 和 y 分別是長(zhǎng)度為 n 和 m 的向量，Z 是 m*n 矩陣。</p><p>　　(3) mesh(Z) 命令繪制網(wǎng)格曲面<

41、;/p><p>　　若 Z 是 m*n 矩陣，則柵格數(shù)據(jù)點(diǎn)的取法是：x=1:n，y=1:m</p><p>　　3、圖形繪制中的視角定義</p><p>　　view(2) 設(shè)定 [AZ ， EL]=[0 ， 90]， 即觀測(cè) x-y 平面的顯示效果。 view(3) 設(shè)定[AZ ，EL]=[-37.5，30]，即缺省視點(diǎn)位置。</p><p

42、>　　1.3.3 圖形的控制和設(shè)置</p><p>　　1、創(chuàng)建或打開圖形窗口</p><p>　?。?）figure 命令 每調(diào)用一次就打開一個(gè)新的圖形窗口。</p><p>　?。?）figure（n）命令</p><p><b>　　2、圖形重疊</b></p><p>　　hold

43、on 命令：保留當(dāng)前窗口的圖形；hold off 命令：解除 hold on 命令。</p><p>　　利用圖形窗口分割函數(shù) subplot（），也可以在同一個(gè)圖形窗口中繪制多幅圖形。</p><p>　　當(dāng)我們調(diào)用 plot 命令繪制二維曲線時(shí)，MATLAB 的執(zhí)行過程大致如下：</p><p>　　（1） ‘使用 figure 命令，在屏幕（root）對(duì)象上生

44、成一個(gè)圖形窗口（Figure 對(duì)象）；</p><p>　　（2）使用 axis 命令，在圖形窗口內(nèi)生成一個(gè)繪圖區(qū)域（Axes 對(duì)象）；</p><p>　?。?）最后用 line 命令在 Axes 指定的區(qū)域內(nèi)繪制線條（Line 對(duì)象）</p><p><b>　　3、屬性</b></p><p>　?。?）獲取屬性P

45、ropertyValue=get(handle,’ProperName’)</p><p>　?。?）屬性的設(shè)置set(handle,‘ProtertyName1’,ProtertyValue1,’PropertyName2’,PropertyValue2,…)</p><p>　　通過 set()函數(shù)設(shè)置坐標(biāo)軸對(duì)象和字符對(duì)象的屬性</p><p><b>

46、;　　1.4 實(shí)例簡(jiǎn)介</b></p><p>　　1.4.1 應(yīng)用基礎(chǔ)</p><p><b>　　變量</b></p><p>　　a='happy new year' b=365</p><p><b>　　函數(shù)</b></p><p>

47、<b>　　a=rand(4)</b></p><p>　　b=rand(2,5)</p><p><b>　　運(yùn)行結(jié)果為：</b></p><p><b>　　a =</b></p><p>　　0.9501 0.89130.82140.9218</p>

48、<p>　　0.2311 0.76210.44470.7382</p><p>　　0.6068 0.45650.61540.1763</p><p>　　0.4860 0.01850.79190.4057 </p><p><b>　　b=</b></p><p>　　0.3028

49、 0.15090.3784 0.8537 0.4966</p><p>　　0.5417 0.69790.8600 0.5936 0.8998</p><p><b>　　矩陣</b></p><p>　　a=[1 2 3 ;4 5 6;7 8 9]</p><p><b>　　運(yùn)行結(jié)果

50、為</b></p><p><b>　　a=</b></p><p><b>　　1 2 3 </b></p><p><b>　　4 5 6</b></p><p><b>　　7 8 9</b></p><p

51、><b>　　向量</b></p><p><b>　　a=1:10</b></p><p><b>　　運(yùn)行結(jié)果為</b></p><p><b>　　a =</b></p><p>　　1 2 3 4 5 6 7 8 9 10</p>

52、<p><b>　　矩陣加減與數(shù)組加減</b></p><p><b>　　A=[1 1 1</b></p><p><b>　　2 2 2</b></p><p><b>　　3 3 3]; </b></p><p><b>　　

53、B=A;</b></p><p><b>　　A+B</b></p><p><b>　　ans =</b></p><p><b>　　2 2</b></p><p><b>　　4 4</b></p><p>&

54、lt;b>　　6 6 6</b></p><p><b>　　A-B</b></p><p><b>　　ans =</b></p><p><b>　　0 0 0</b></p><p><b>　　0 0 0</b><

55、;/p><p><b>　　0 0 0</b></p><p><b>　　矩陣乘與數(shù)組乘</b></p><p>　　a=[1 1 1;2 2 2;3 3 3];</p><p><b>　　b=a; </b></p><p><b>　　a

56、*b </b></p><p><b>　　ans =</b></p><p><b>　　6 6</b></p><p>　　12 12 12</p><p>　　18 18 18</p><p>　　f=ones(1,3); g=ones(3,1);

57、</p><p><b>　　f*g</b></p><p><b>　　ans =</b></p><p><b>　　g*f</b></p><p><b>　　ans =</b></p><p><b>　　3<

58、;/b></p><p><b>　　1 1 1</b></p><p><b>　　1 1 1</b></p><p><b>　　1 1 1</b></p><p><b>　　矩陣除與數(shù)組除</b></p>

59、<p>　　A=[2 2 3 3 4 4</p><p>　　1 1 2 2 3 3</p><p>　　4 4 5 5 6 6]; </p><p>　　B=[1 1 3 3 2 2</p><p>　　1 1 1 1 1 1</p><p>　　2 2 5 5 3 3]; </p><

60、;p><b>　　A./B</b></p><p><b>　　ans =</b></p><p>　　2 2 1 1 2 2</p><p>　　1 1 2 2 3 3</p><p>　　2 2

61、 1 1 2 2</p><p><b>　　B./A</b></p><p><b>　　ans =</b></p><p>　　0.5000 0.5000 1.0000 1.0000 0.5000 0.5000 </p><p>　　1

62、.0000 1.0000 0.5000 0.5000 0.3333 0.3333</p><p>　　0.5000 0.5000 1.0000 1.0000 0.5000 0.5000</p><p>　　1.4.2 繪制二維圖形</p><p>　　1、連續(xù)序列圖的繪制</p><

63、p>　　function [f]=aa(t)</p><p>　　f=sym('sin(pi/4*t)')</p><p>　　ezplot(f,[-16,16])</p><p><b>　　圖形如下：</b></p><p><b>　　n=0:15;</b></p&

64、gt;<p>　　2、離散序列圖的繪制</p><p>　　% steam plot</p><p>　　x=0:0.1:4; y=sin(x.^2).*exp(-x); stem(x,y)</p><p><b>　　圖形如右：</b></p><p>　　1.4.3繪制三維圖形</p>&

65、lt;p><b>　　1、</b></p><p>　　x=-8:0.5:8; </p><p><b>　　y=x; </b></p><p>　　X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps; </p><p>　　Z=sin(R)./R; &

66、lt;/p><p>　　mesh(X,Y,Z) </p><p><b>　　grid on </b></p><p>　　axis([-10 10 -10 10 -1 1])</p><p><b>　　2、視角變換</b></p><p>　　[X,Y] = meshgrid

67、([-2:.25:2]); </p><p>　　Z = X.*exp(-X.^2 -Y.^2);</p><p>　　surf(X,Y,Z); view([0 0]); </p><p>　　view([180 0]); </p><p>　　第二

68、章 信號(hào)與系統(tǒng)分析及MATLAB實(shí)現(xiàn)</p><p>　　2.1 信號(hào)的時(shí)域分析</p><p>　　2.1.1 信號(hào)的表式及可視化 </p><p><b>　?。?）連續(xù)時(shí)間信號(hào)</b></p><p><b>　　向量表示法：</b></p><p>　　對(duì)于連續(xù)

69、時(shí)間信號(hào) f (t) ， 我們可以用兩個(gè)行向量 f 和 t 來表示 ， 其中向量 t 是形如 t=t1：p : t2 的 MATLAB 命令定義的時(shí)間范圍向量， t1 為信號(hào)起始時(shí)間， t 2 為終止時(shí)間，p 為時(shí)間間隔。向量 f 為連續(xù)信號(hào) f (t) 在向量 t 所定義的時(shí)間點(diǎn)上的樣值。例如對(duì)于連續(xù)信號(hào) f (t)=Sa (t)=sin( t)/t</p><p><b>　　符號(hào)運(yùn)算表示法：&l

70、t;/b></p><p>　　用ezplot（）命令繪制出信號(hào)波形</p><p><b>　?。?）離散時(shí)間信號(hào)</b></p><p>　　一般說來，離散時(shí)間信號(hào)用 f (k ) 表示，其中變量 k 為整數(shù)，代表離散的采樣時(shí)間點(diǎn)。f (k ) 可表示為：</p><p>　　f (k )= {……f (2 )

71、, f (1), f (0), f (1), f (2)……} </p><p>　　在用MATLAB表示離散序列并將其可視化時(shí)，注意：</p><p>　　與連續(xù)時(shí)間信號(hào)不同，離散時(shí)間信號(hào)無法用符號(hào)運(yùn)算來表示；</p><p>　　第二，由于在MATLAB中，矩陣的元素個(gè)數(shù)是有限的，因此，MATLAB無法表示無序列；</p><p>　　第

72、三，在繪制離散信號(hào)波形時(shí)，要使用專門繪制離散數(shù)據(jù)的stem命令，而不是plot命令。</p><p>　　2.1.2 信號(hào)的時(shí)域運(yùn)算、時(shí)域變換及MATLAB實(shí)現(xiàn)</p><p>　?。?）連續(xù)信號(hào)的時(shí)域運(yùn)算、時(shí)域變換</p><p><b>　　1． 相加</b></p><p>　　連續(xù)信號(hào)的相加，是指兩信號(hào)的對(duì)應(yīng)時(shí)刻

73、值相加，即</p><p>　　f (t )=f1 (t)+f 2 (t ) </p><p>　　2. 相乘 連續(xù)信號(hào)的相乘，是指兩信號(hào)的對(duì)應(yīng)時(shí)刻值相乘，</p><p>　　即f (t)=f1 (t )*f2 (t)</p><p><b>　　3. 移位</b></p><p>　　連續(xù)信

74、號(hào)的移位也稱平移。對(duì)于連續(xù)信號(hào) f (t) ，若有常數(shù) t0</p><p>　　>0 ，延時(shí)信號(hào) f (t-t0 )是將原信號(hào)沿正 t 軸方向平移時(shí)間 t0 ，而 f (t+t0 ) 是將原信號(hào)沿負(fù) t 軸方向移動(dòng)時(shí)間 t0 。</p><p><b>　　4. 反折</b></p><p>　　連續(xù)信號(hào)的反折，是指將信號(hào)以縱坐標(biāo)為軸反

75、折，即將信號(hào) f (t) 中的自變量 t 換為-t 。</p><p><b>　　5. 尺度變換</b></p><p>　　連續(xù)信號(hào)的尺度變換，是指將信號(hào)的橫坐標(biāo)進(jìn)行展寬或壓縮變換，即將信號(hào)f (t) 中的自變量 t 換為 at ，當(dāng) a>1 時(shí)，信號(hào)f (at ) 以原點(diǎn)為基準(zhǔn)，沿橫軸壓縮到原來的 1/ a ；當(dāng)0<a<1 時(shí)，信號(hào) f (a

76、t) 將沿橫軸展寬至原來的 1/ a 倍。</p><p><b>　　6．倒相</b></p><p>　　連續(xù)信號(hào)的倒相，是指將信號(hào) f (t)以橫軸為對(duì)稱軸對(duì)折得到f (t)</p><p>　?。?）離散信號(hào)的時(shí)域運(yùn)算、時(shí)域變換</p><p>　　對(duì)于離散序列來說，序列相加、相乘是將兩序列對(duì)應(yīng)時(shí)間序號(hào)的值逐項(xiàng)相

77、加或相乘， 平移、反折、及倒相變換與連續(xù)信號(hào)的定義完全相同，這里就不再累述。但需要注意，與 連續(xù)信號(hào)不同的是，在MATLAB中，離散序列的時(shí)域運(yùn)算和變換不能用符號(hào)運(yùn)算來實(shí)現(xiàn)，而 必須用向量表示的方法，即在MATLAB中離散序列的相加、相乘需表示成兩個(gè)向量的相加、 相乘，因而參加運(yùn)算的兩序列向量必須具有相同的維數(shù)。</p><p>　　2.1.3 用MATLAB分析常用的時(shí)間信號(hào)特性</p><

78、p><b>　　1、連續(xù)時(shí)間信號(hào)</b></p><p>　　正弦信號(hào)定義為：f (t )=A cos(wt+φ )</p><p>　　實(shí)指數(shù)信號(hào)： f (t)=Ce^(at) ，其中 c 和 a 為實(shí)常數(shù)</p><p>　　虛指數(shù)信號(hào)：f (t )=Ae^(jwt), A 為 常數(shù), w 為虛指數(shù)信號(hào)的角頻率。</p>

79、<p>　　復(fù)指數(shù)信號(hào)：f (t )=Ae^(st) ，其中 s=σ+jw 為復(fù)常數(shù)。</p><p><b>　　2、離散時(shí)間信號(hào)</b></p><p>　　離散時(shí)間正弦序列的表達(dá)式為：f (k )= A cos(kw+j ) ：式中 k 為無量綱的整數(shù)， w 和 j 以弧度為單位， w 稱為離散正弦序列的數(shù)字角頻率，j 為初相位。</p>

80、<p>　　離散時(shí)間實(shí)指數(shù)序列的一般形式為:f (k)=ca ^k ，其中 c 和 a 為實(shí)常數(shù)。</p><p>　　離散時(shí)間虛指數(shù)序列的一般形式為:f (k )=e^(jwk)其中 w 為角頻率</p><p>　　復(fù)指數(shù)序列的一般形式為：f (k )=r^(k)*e^(jwk)</p><p>　　2.2 連續(xù)系統(tǒng)的是與分析及MATLAB實(shí)現(xiàn)&

81、lt;/p><p>　　2.2.1 離散時(shí)間序列卷積及MATLAB實(shí)現(xiàn)</p><p>　　MATLAB 的 conv()函數(shù)可以幫助我們快速求出兩個(gè)離散序列的卷積和。conv 函數(shù)的 調(diào)用格式為：f = conv (f1, f2 )其中 f1 為包含序列 f1(k ) 的非零樣值點(diǎn)的行向量，f2 為包含序列 f 2 (k ) 的非零樣值點(diǎn)的行向量，向量 f 則返回序列 f (k )=f1 (

82、k )*f 2 (k ) 的所有非零樣值點(diǎn)行向量</p><p>　　2.2.2 連續(xù)時(shí)間信號(hào)卷積及MATLAB實(shí)現(xiàn)</p><p>　　兩個(gè)與卷積相關(guān)的重要結(jié)論：</p><p>　?。?）f (t ) = f (t) *d (t )，即連續(xù)信號(hào)可分解為一系列幅度由f (t )決定的沖激信號(hào)d (t)及其平移信號(hào)之和。</p><p>　　

83、（2）線性時(shí)不變連續(xù)系統(tǒng)，設(shè)其輸入信號(hào)為f (t)，單位響應(yīng)為h(t)，其零狀態(tài)響應(yīng)為y(t )，則有： y (t) = f (t ) * h(t)</p><p>　　用 MATLAB 實(shí)現(xiàn)連續(xù)信號(hào) f1 (t) 與 f 2 (t ) 卷積的過程如下：</p><p>　?。?）將連續(xù)信號(hào) f1 (t) 與 f 2 (t ) 以時(shí)間隔Δ進(jìn)行取樣，得到離散序列 f1 (kΔ) 和 f2 (

84、kΔ) </p><p>　?。?）構(gòu)造與 f1 (k) 和 f 2 (k) 相對(duì)應(yīng)的時(shí)間向量f1 (kΔ)和 f2 (kΔ)（注意，此時(shí)時(shí)間序號(hào)向量 k1、k2 的元素不再是整數(shù)，而是取樣時(shí)間隔Δ的整數(shù)倍的時(shí)間間隔點(diǎn)）；</p><p>　?。?）調(diào)用 conv()函數(shù)計(jì)算卷積積分 f (t) 的近似向量 f (n ) ； </p><p>　?。?）構(gòu)造 f

85、(n) 對(duì)應(yīng)的時(shí)間向量 k。</p><p>　　LTI 系統(tǒng)當(dāng)輸入為沖激信號(hào)d (t)時(shí)產(chǎn)生的零狀態(tài)響應(yīng)稱為系統(tǒng)的沖激響應(yīng)，用h(t)表示。若輸入為單位階躍信號(hào)e (t)時(shí)系統(tǒng)產(chǎn)生的零狀態(tài)響應(yīng)則稱為系統(tǒng)的階躍響應(yīng)，記為g (t)，如圖所示：對(duì) LTI 連續(xù)系統(tǒng)，設(shè)其輸入信號(hào)為f (t)，沖激響應(yīng)為h(t)，零狀態(tài)響應(yīng)為y(t )，則有：y (t) = f (t ) * h(t)即h(t)包含了連續(xù)系統(tǒng)的固有特性

86、，與系統(tǒng)的輸入無關(guān)。</p><p>　　MATLAB 為用戶提供了專門用于求連續(xù)系統(tǒng)沖激響應(yīng)及階躍響應(yīng)，并繪制其時(shí)域波形 的函數(shù) impulse 和 step。在調(diào)用函數(shù) impulse()和 step()時(shí)，我們需要用向量來對(duì)連續(xù)系統(tǒng)進(jìn)行表示。</p><p>　　Impulse()函數(shù)有如下幾種調(diào)</p><p>　?。?）impulse( b, a )<

87、;/p><p>　　（2）impulse( b, a, t )</p><p>　?。?）impulse( b, a, t1:p:t2 )</p><p>　?。?）impulse( b, a, t1:p:t2 )</p><p>　　函數(shù) step()將繪出由向量 a 和 b 表示的連續(xù)系統(tǒng)的階躍響應(yīng) g (t) 在指定時(shí)間范圍內(nèi)的 波形圖，并

88、能求出其數(shù)值解。和 impulse()函數(shù)一樣，step()函數(shù)也有如下四種調(diào)用格式：</p><p>　?。?）step(b,a)</p><p>　?。?）step(b,a,t)</p><p>　?。?）step(b,a,t1:p:t2)</p><p>　?。?）y=step(b,a,t1:p:t2)</p><p

89、>　　函數(shù) impz( )能繪出向量 a 和 b 定義的離散系統(tǒng)在指定時(shí)間范圍內(nèi)單位響應(yīng)的時(shí)域波形， 并能求出系統(tǒng)單位響應(yīng)在指定時(shí)間范圍內(nèi)的數(shù)值解。Impz()函數(shù)有如下幾種調(diào)用格式：</p><p>　?。?）impz(b,a)</p><p>　　（2）impz(b,a,n)</p><p>　?。?）impz(b,a,n1:n2)</p>

90、<p>　?。?）impz(b,a,n1:n2)</p><p>　　MATLAB 的函數(shù) lsim( )能對(duì)上述微分方程描述的 LTI 連續(xù)系統(tǒng)的響應(yīng)進(jìn)行仿真。lsim( ) 函數(shù)能繪制連續(xù)系統(tǒng)在指定的任意時(shí)間范圍內(nèi)系統(tǒng)響應(yīng)的時(shí)域波形圖，還能求出連續(xù)系統(tǒng) 在指定的任意時(shí)間范圍內(nèi)系統(tǒng)響應(yīng)的數(shù)值解。lsim()函數(shù)有如下兩種調(diào)用格式：</p><p>　?。?）lsim(b,a,x

91、,t)</p><p>　　（2）lsim(b,a,x,t)</p><p>　　MATLAB 為用戶提供了求 LTI 離散系統(tǒng)響應(yīng)的專用函數(shù) filter()。該函數(shù)能求出由差分 方程描述的離散系統(tǒng)在指定時(shí)間范圍內(nèi)的輸入序列時(shí)所產(chǎn)生的響應(yīng)序列的數(shù)值解。</p><p>　　filter()函數(shù)的調(diào)用格式：filter(b,a,x)</p><p

92、><b>　　2.3 實(shí)例簡(jiǎn)介</b></p><p>　　2.3.1 信號(hào)的時(shí)域分析及MATLAB實(shí)現(xiàn)</p><p>　　例1:設(shè)信號(hào) f(t) =(1+t/2)*[ε(t+2)-ε(t-2)],用MATLAB求f(t+2)，f(t-2), f(-t), f(2t),- f(t)并繪出其時(shí)域波形。 </p><p>　　解：根據(jù)前面的

93、介紹，我們可用符號(hào)運(yùn)算來實(shí)現(xiàn)上述過程，MATLAB命令如下： </p><p><b>　　syms t </b></p><p>　　f=sym('(t/2+1)*(heaviside(t+2)-heaviside(t-2))') </p><p>　　subplot(2,3,1),ezplot(f,[-3,3]) </

94、p><p>　　y1=subs(f,t,t+2) </p><p>　　subplot(2,3,2),ezplot(y1,[-5,1]) </p><p>　　y2=subs(f,t,t-2)</p><p>　　subplot(2,3,3),ezplot(y2,[-1,5]) </p><p>　　y3=subs(f,t

95、,-t) </p><p>　　subplot(2,3,4),ezplot(y3,[-3,3]) </p><p>　　y4=subs(f,t,2*t) </p><p>　　subplot(2,3,5),ezplot(y4,[-2,2]) </p><p><b>　　y5=-f </b></p><

96、;p>　　subplot(2,3,6),ezplot(y5,[-3,3]) </p><p>　　命令執(zhí)行后得到f,y1,y2,y3,y4,y5的符號(hào)表達(dá)式如下： </p><p>　　f = (t/2+1)*(heaviside(t+2)-heaviside(t-2)); </p><p>　　y1 = (1/2*t+2)*(heaviside(t+4)-

97、heaviside(t)); </p><p>　　y2 = 1/2*t*(heaviside(t)-heaviside(t-4)) ; </p><p>　　y3 = (-1/2*t+1)*(heaviside(-t+2)-heaviside(-t-2)) ;</p><p>　　y4 = (t+1)*(heaviside(2*t+2)-heaviside(

98、2*t-2)) ;</p><p>　　y5 = -(1/2*t+1)*(heaviside(t+2)-heaviside(t-2));</p><p><b>　　圖形如下：</b></p><p>　　例 2：已知如圖所示信號(hào) f1（t）及信號(hào)f2（t）=sin(2兀t ) </p><p>　　用 MATLAB繪

99、出滿足下列要求的信號(hào)波形。</p><p>　　(1)f3（t）=f1（-t）+f1（t） (2) f4（t）=-[f1（-t）+ f1（t）]</p><p>　　(3)f5（t）=f2（t）* f3（t） (4) f6（t）=f1（t）* f2（t）</p><p>　　解：信號(hào)f1（t）可表示為f1（t）=(-t

100、+4)[ε(t)-ε(t-4)]，這樣我們即可用MATLAB 的符號(hào)運(yùn)算功能來實(shí)現(xiàn)題目要求的時(shí)域運(yùn)算： </p><p><b>　　syms t </b></p><p>　　f1=sym('(-1*t+4)*(heaviside(t)-heaviside(t-4))') </p>

101、<p>　　subplot(2,3,1),ezplot(f1) </p><p>　　f2=sym('sin(2*pi*t)') </p><p>　　subplot(2,3,4),ezplot(f2,[-4,4]) </p><p>　　y1=subs(f1,t,-t) </p><p><b>　　f3

102、=f1+y1 </b></p><p>　　subplot(2,3,2),ezplot(f3) </p><p><b>　　f4=-f3 </b></p><p>　　subplot(2,3,3),ezplot(f4) </p><p><b>　　f5=f2*f3 </b></

103、p><p>　　subplot(2,3,5),ezplot(f5) </p><p><b>　　f6=f1*f2 </b></p><p>　　subplot(2,3,6),ezplot(f6)</p><p>　　命令執(zhí)行后得到f1(t) ,f2(t) ,f3(t) ,f4(t), f5(t), f6(t)的符號(hào)表達(dá)式如

104、下：f1 =(-1*t+4)*(heaviside(t)-heaviside(t-4))</p><p>　　f2 =sin(2*pi*t)</p><p>　　y1 =(t+4)*(heaviside(-t)-heaviside(-t-4)) f3=(-t+4)*(heaviside(t)-heaviside(t-4))+(t+4)*(heaviside(-t)-heaviside(-t

105、-4)) </p><p>　　f4 = -(-t+4)*(heaviside(t)-heaviside(t-4))-(t+4)*(heaviside(-t)-heaviside(-t-4))</p><p>　　f5 = sin(2*pi*t)*((-t+4)*(heaviside(t)-heaviside(t-4))+(t+4)*(heaviside(-t)-heaviside(-t-

106、4))) </p><p>　　f6 = (-t+4)*(heaviside(t)-heaviside(t-4))*sin(2*pi*t);圖形如下：</p><p>　　例3：有兩離散序列，f 1 (k) ={ 2, 1,0,1, 2,f 2 },f2((k ) ={ 1,1, 1 }用MATLAB繪出它們的波形及f 1 (k )+f 2 (k) 的波形。 <

107、;/p><p>　　解： MATLAB命令為： </p><p><b>　　f1=-2:2； </b></p><p><b>　　k1=-2:2； </b></p><p>　　f2=[1 1 1]； </p><p><b>　　k2=-1:1； </b&g

108、t;</p><p>　　stem(k1,f1),axis(-3,3,-2.5,2.5) </p><p>　　stem(k2,f2),axis(-3,3,-2.5,2.5) </p><p>　　[f,k]=lsxj (f1,f2,k1,k2) </p><p>　　程序運(yùn)行結(jié)果如下，繪制的波形如圖所示： </p><p

109、>　　f= -2 0 1 2 2 </p><p>　　k= -2 –1 0 1 2 </p><p>　　例4：畫出復(fù)指數(shù)信號(hào)f (t)= e^(-t)* e^(10it)的實(shí)部、虛部、模及相角隨時(shí)間變化的曲線， 并觀察其時(shí)域特性。 </p><p>　　解：對(duì)應(yīng)的MA

110、TLAB程序如下：</p><p>　　t=0:0.01:3 </p><p>　　a=-1;b=10;</p><p>　　z=exp((a+i*b)*t); </p><p>　　subplot(2,2,1),plot(t,real(z)),title('實(shí)部'); subplot(2,2,3),plot(t,ima

111、g(z)),title('虛部'); subplot(2,2,2),plot(t,abs(z)),title('模'); subplot(2,2,4),plot(t,angle(z)),title('相角'); </p><p>　　例5：畫出復(fù)指數(shù)信號(hào)f (t)= r^(k)* e^(jwk)的實(shí)部、虛部、模及相角隨時(shí)間變化的曲線， 并觀察其時(shí)域特性&

112、lt;/p><p><b>　　解：創(chuàng)建M文件：</b></p><p>　　function dfzsu(n1,n2,r,w) </p><p>　　%n1:繪制波形的虛指數(shù)序列的起始時(shí)間序號(hào) </p><p>　　%n2：繪制波形的虛指數(shù)序列的終止時(shí)間序號(hào) </p><p>　　%w：虛指數(shù)序列的

113、角頻率 </p><p>　　%r: 指數(shù)序列的底數(shù) </p><p><b>　　k=n1:n2; </b></p><p>　　f=(r*exp(i*w)).^k; </p><p>　　Xr=real(f); </p><p>　　Xi=imag(f); </p><p&

114、gt;　　Xa=abs(f); </p><p>　　Xn=angle(f); </p><p>　　subplot(2,2,1), stem(k,Xr,'filled'),title('實(shí)部'); </p><p>　　subplot(2,2,3), stem(k,Xi,'filled'),title('虛部

115、'); </p><p>　　subplot(2,2,2), stem(k,Xa,'filled'),title('模'); </p><p>　　subplot(2,2,4),stem(k,Xn,'filled'),title（’相角’）;</p><p>　　輸入dfzsu(0,20,0.9,pi/4) ，

116、結(jié)果如圖：</p><p>　　2.3.2 離散時(shí)間序列卷積及MATLAB實(shí)現(xiàn)</p><p>　　例6：已知某LTI 離散系統(tǒng)，其單位響應(yīng)h(k) = e(k) -e(k - 4)，求該系統(tǒng)在激勵(lì)為f (k) = e(k) -e(k -3)時(shí)的零狀態(tài)響應(yīng)y(k)，并繪出其時(shí)域波形圖。</p><p>　　解：由LTI 變離散系統(tǒng)的分析可得</p>&

117、lt;p>　　y(k) = h(k) * f (k)</p><p>　　因此，我們可調(diào)用dconv()函數(shù)來解決此問題，相應(yīng)的MATLAB 命令為：</p><p>　　f1=ones(1,4);</p><p><b>　　k1=0:3;</b></p><p>　　f2=ones(1,3);</p>

118、;<p><b>　　k2=0:2;</b></p><p>　　[f,k]=dconv(f1,f2,k1,k2)</p><p><b>　　運(yùn)行結(jié)果為：</b></p><p><b>　　f =</b></p><p>　　1 2 3 3 2 1</p

119、><p><b>　　k =</b></p><p>　　0 1 2 3 4 5 時(shí)域波形圖</p><p>　　例7：已知兩連續(xù)時(shí)間信號(hào)如下圖所示，試用MATLAB 求f（t）??f1（t） ??f2（ t），并繪出f (t)的時(shí)域波形圖。</p>

120、;<p>　　解: 我們可以調(diào)用前述的函數(shù)sconv( )來解決此問題，即首先設(shè)定取樣時(shí)間間隔p，并對(duì)連續(xù)信號(hào)f1（t）和f2（t）的非零值區(qū)間以時(shí)間間隔p 進(jìn)行抽樣，產(chǎn)生離散序列f1和f2，然后構(gòu)造離散序列f1 和f2 所對(duì)應(yīng)的時(shí)間向量k1 和k2，最后再調(diào)用sconv 函數(shù)即可求出f1（t） ??f2（ t）的數(shù)值近似，并繪出其時(shí)域波形圖。實(shí)現(xiàn)上述過程的 MATLAB 命令如下：</p><p>

121、;<b>　　p=0.01;</b></p><p><b>　　k1=0:p:2;</b></p><p>　　f1=0.5*k1;</p><p><b>　　k2=k1;</b></p><p><b>　　f2=f1;</b></p>

122、<p>　　[f,k]=sconv(f1,f2,k1,k2,p)</p><p>　　上述命令繪制的波形如圖所示。圖中分別給出了取樣時(shí)間間隔p=0.5 和p=0.01時(shí)的處理效果?？梢?，當(dāng)取樣時(shí)間p 足夠小時(shí)，函數(shù)sconv( )的計(jì)算結(jié)果就是連續(xù)時(shí)間卷積f（t）=f1（t）*f2（ t）</p><p><b>　　的較好的數(shù)值近似</b></p&

123、gt;<p>　　例8：已知描述某連續(xù)系統(tǒng)的微分方程為：2 y??(t ) ??y?(t) ??8 y(t) ??f (t )</p><p>　　試用 MATLAB 繪出該系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)的波形。</p><p>　　解：直接調(diào)用函數(shù)impulse()和step()即可解決此問題，對(duì)應(yīng)的MATLAB 命令如下：</p><p><b&g

124、t;　　b=[1];</b></p><p>　　a=[2 1 8];</p><p>　　subplot(1,2,1)</p><p>　　impulse(b,a)</p><p>　　subplot(1,2,2)</p><p><b>　　step(b,a）</b></p&

125、gt;<p>　　上述命令繪制的系統(tǒng)沖激響應(yīng)和階躍響應(yīng)的波形如圖所示。</p><p>　　例9：已知描述某離散系統(tǒng)的差分方程如下：</p><p>　　2 y(k) ??2y(k ??1) ??y(k ??2) ??f (k ) ??3 f (k ?1) ??2 f (k ??2)</p><p>　　試用 MATLAB 繪出該系統(tǒng)0~50 時(shí)間范圍

126、內(nèi)單位響應(yīng)的波形。</p><p>　　解：首先用向量a 和b 表示出該離散系統(tǒng)，然后再調(diào)用impz()函數(shù)即可解決此問題。實(shí)現(xiàn)上述過程的MATLAB 命令如下：</p><p>　　a=[2 -2 1]</p><p>　　b=[1 3 2];</p><p><b>　　impz(b,a)</b></p>