連續(xù)時(shí)間信號(hào)和系統(tǒng)時(shí)域分析及MATLAB實(shí)現(xiàn)

1、；. MATLAB 課程設(shè)計(jì)任務(wù)書課程設(shè)計(jì)任務(wù)書 姓名：王姓名：王* 學(xué)號(hào)：學(xué)號(hào)：2010*010 題目題目: 連續(xù)時(shí)間信號(hào)和系統(tǒng)時(shí)域分析及 MATLAB 實(shí)現(xiàn) 初始條件：初始條件： MATLAB 7.5.0 ，Windows XP 系統(tǒng) 實(shí)驗(yàn)任務(wù)：實(shí)驗(yàn)任務(wù)： 一、用 MATLAB 實(shí)現(xiàn)常用連續(xù)時(shí)間信號(hào)的時(shí)域波形（通過改變參數(shù)，分析其時(shí) 域特性） 。 1、單位階躍信號(hào)，2、單位沖激信號(hào)，3、正弦信號(hào)，4、實(shí)指數(shù)信號(hào)，5、虛指 數(shù)信號(hào)，6、復(fù)指數(shù)信號(hào)。 二、用 MATLAB 實(shí)現(xiàn)信號(hào)的時(shí)域運(yùn)算 1、相加 ，2、相乘 ，3、數(shù)乘，4、微分，5、積分 三、用 MATLAB 實(shí)現(xiàn)信號(hào)的時(shí)域變換（參數(shù)變

2、化，分析波形變化） 1、反轉(zhuǎn)，2、使移（超時(shí)，延時(shí)） ，3、展縮，4、倒相，5、綜合變化 四、用 MATLAB 實(shí)現(xiàn)信號(hào)簡(jiǎn)單的時(shí)域分解 1、信號(hào)的交直流分解，2、信號(hào)的奇偶分解 五、用 MATLAB 實(shí)現(xiàn)連續(xù)時(shí)間系統(tǒng)的卷積積分的仿真波形 給出幾個(gè)典型例子，對(duì)每個(gè)例子，要求畫出對(duì)應(yīng)波形。 六、用 MATLAB 實(shí)現(xiàn)連續(xù)時(shí)間系統(tǒng)的沖激響應(yīng)、階躍響應(yīng)的仿真波形。 給出幾個(gè)典型例子，四種調(diào)用格式。 七、利用 MATLAB 實(shí)現(xiàn)連續(xù)時(shí)間系統(tǒng)對(duì)正弦信號(hào)、實(shí)指數(shù)信號(hào)的零狀態(tài)響應(yīng)的 仿真波形。 . ；. 目目 錄錄 1 MATLAB 簡(jiǎn)介.1 1.1 MATLAB 設(shè)計(jì)目的.1 1.2 MATLAB 語言特點(diǎn)

3、.1 2 常用連續(xù)時(shí)間信號(hào)的時(shí)域波形.1 2.1 單位階躍信號(hào).1 2.2 單位沖激信號(hào).2 2.3 正弦信號(hào).3 2.4 實(shí)指數(shù)信號(hào).4 2.5 虛指數(shù)信號(hào).5 2.6 復(fù)指數(shù)信號(hào).6 3 連續(xù)時(shí)間信號(hào)的時(shí)域運(yùn)算.7 3.1 相加.7 3.2 相乘.8 3.3 數(shù)乘.9 3.4 微分.10 3.5 積分.11 4.1 反轉(zhuǎn).12 4.2 時(shí)移.13 4.3 展縮.14 4.4 倒相.15 . ；. 4.5 綜合變化.16 5 連續(xù)時(shí)間信號(hào)簡(jiǎn)單的時(shí)域分解.17 5.1 信號(hào)的交直流分解.17 5.2 信號(hào)的奇偶分解.18 6 連續(xù)時(shí)間系統(tǒng)的卷積積分的仿真波形.20 7 連續(xù)時(shí)間系統(tǒng)的沖激響應(yīng)、

4、階躍響應(yīng)的仿真波形.22 7.1 IMPULSE（）函數(shù) .23 7.2 STEP（）函數(shù).27 8 連續(xù)時(shí)間系統(tǒng)對(duì)正弦信號(hào)、實(shí)指數(shù)信號(hào)的零狀態(tài)響應(yīng)的仿真 波形.29 8.1 正弦信號(hào)的零狀態(tài)響應(yīng).30 8.2 實(shí)指數(shù)信號(hào)的零狀態(tài)響應(yīng).31 9 小結(jié).32 . ；. 1 MATLAB 簡(jiǎn)介簡(jiǎn)介 1.1 MATLAB 設(shè)計(jì)目的設(shè)計(jì)目的 深入研究連續(xù)時(shí)間信號(hào)和系統(tǒng)時(shí)域分析的理論知識(shí)。利用MATLAB強(qiáng)大的圖形處理 功能、符號(hào)運(yùn)算功能以及數(shù)值計(jì)算功能，實(shí)現(xiàn)連續(xù)時(shí)間信號(hào)和系統(tǒng)時(shí)域分析的仿真波形 1.2 MATLAB 語言特點(diǎn)語言特點(diǎn) MATLAB和Mathematica、Maple并稱為三大數(shù)學(xué)軟件。

5、它在數(shù)學(xué)類科技應(yīng)用軟件中 在數(shù)值計(jì)算方面首屈一指。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)域。MATLAB的最重要特征 使他擁有解決特定應(yīng)用問題的程序組，也就是TOOLBOX(工具箱)，如信號(hào)處理工具箱， 控制系統(tǒng)工具箱、神經(jīng)網(wǎng)絡(luò)工具箱、模糊邏輯工具箱、通信工具箱和數(shù)據(jù)采集工具箱等 許多專用工具箱，對(duì)大多數(shù)用戶來說，要想靈活、高效地運(yùn)用這些工具箱，通常都需要 學(xué)習(xí)相應(yīng)的專業(yè)知識(shí)。 2 常用連續(xù)時(shí)間信號(hào)的時(shí)域波形常用連續(xù)時(shí)間信號(hào)的時(shí)域波形 連續(xù)信號(hào)又

6、稱為模擬信號(hào)，其信號(hào)存在于整個(gè)時(shí)間范圍內(nèi)，包括單位階躍信號(hào)，單 位沖激信號(hào)，正弦信號(hào)，實(shí)指數(shù)信號(hào)，虛指數(shù)信號(hào)，復(fù)指數(shù)信號(hào)。 2.1 單位階躍信號(hào)單位階躍信號(hào) 單位階躍信號(hào)的定義如下： () =0,( 0) ? 單位階躍信號(hào)的MATLAB程序: t=-0.5:0.01:5; t0=1.0; q=stepfun(t,t0); plot(t,q); . ；. axis equal 其信號(hào)圖如下： 2.2 單位沖激信號(hào)單位沖激信號(hào) MATLAB實(shí)現(xiàn)程序如下: t=-5:0.01:5; a=(t=0); plot(t,a); 信號(hào)圖如下： . ；. 2.3 正弦信號(hào)正弦信號(hào) 正弦信號(hào)其 MATLAB 實(shí)

7、現(xiàn)程序如下： t=-1:0.0001:1; A=6; f=5; b=1; u=A*sin(2*pi*f*t+b); plot(t,u) axis(-1 1 -6.5 6.5) 其信號(hào)圖如下： . ；. 2.4 實(shí)指數(shù)信號(hào) 實(shí)指數(shù)信號(hào)可由下面的表達(dá)式來表示： () = MATLAB 實(shí)現(xiàn)程序如下： t=0:0.002:3; A=3; a=0.5; b=A*exp(a*t); plot(t,b) axis(-0.2 3.1 -0.2 14) 其信號(hào)圖如下： . ；. 2.5 虛指數(shù)信號(hào) 虛指數(shù)信號(hào)可由下面的表達(dá)式來表示： () = A=2，的虛指數(shù)信號(hào)的 MATLAB 實(shí)現(xiàn)程序如下 = /4()

8、= t=0:0.001:20; a=2; w=pi/4; b=a*exp(i*w*t); subplot(221),plot(t,real(b),axis(0, 20,-4,4),title(實(shí)部) subplot(222),plot(t,imag(b),axis(0,20,-4,4),title(虛部) subplot(223),plot(t,abs(b),axis(0,20,1,4),title(模) subplot(224),plot(t,angle(b),axis(0,20,-4,4),title(相角) . ；. 其信號(hào)圖如下： 2.6 復(fù)指數(shù)信號(hào) 復(fù)指數(shù)信號(hào)可由下面的表達(dá)式來表示：

9、 () = ( + ) MATLAB 實(shí)現(xiàn)程序如下： t=0:0.01:4; a=-1; A=1; b=12; c=A*exp(a+i*b)*t); . ；. subplot(221),plot(t,real(c),title(實(shí)部) subplot(222),plot(t,abs(c),title(模) subplot(223),plot(t,imag(c),title(虛部) subplot(224),plot(t,angle(c),title(相角) 其信號(hào)圖如下： 3 連續(xù)時(shí)間連續(xù)時(shí)間信號(hào)的時(shí)域運(yùn)算信號(hào)的時(shí)域運(yùn)算 3.1 相加相加 要實(shí)現(xiàn)兩信號(hào)的相加，即 f（t）=f1（t）+f2（t

10、） . ；. f1（t）為單位階躍信號(hào)，f2（t）為正弦信號(hào)，兩信號(hào)相加的實(shí)現(xiàn)程序如下， t=-6:0.0001:10; t0=2; a=stepfun(t,t0); b=sin(2*pi*t); f=b+a; plot(t,f) axis(-6 10 -3 3) 其信號(hào)圖如下： 3.2 相乘相乘 要實(shí)現(xiàn)兩信號(hào)的相乘，即 f（t）=f1（t）*f2（t） f1（t）為單位階躍信號(hào)，f2（t）為正弦信號(hào)，兩信號(hào)相乘的實(shí)現(xiàn)程序如下： t=0:0.0001:5; . ；. t0=0.5; a=stepfun(t,t0); b=sin(2*pi*t); f=a.*b; plot(t,f); axis(

11、0 5 -2 2); 其信號(hào)圖如下： 3.3 數(shù)乘數(shù)乘 要實(shí)現(xiàn)信號(hào)的數(shù)乘，即 f（t）=A*f1（t） A=3，f1（t）為單位階躍信號(hào)，信號(hào)數(shù)乘的實(shí)現(xiàn)程序如下： t=0:0.0001:5; . ；. a=3; t0=1; b=stepfun(t,t0); f=a*b; plot(t,f); axis(-2 5 0 5); 其信號(hào)圖如下： 3.4 微分微分 微分即求信號(hào)的導(dǎo)數(shù)。 對(duì)函數(shù) f（t）=t2求一階微分的實(shí)現(xiàn)程序如下： t=-40:0.002:40; a=t.*t; . ；. d=diff(a); subplot(211); plot(t,a,-); subplot(212); plo

12、t(d,-); 其信號(hào)圖如下： 3.5 積分積分 對(duì) f（t）=t2函數(shù)的一次積分的實(shí)現(xiàn)程序如下： t=-2:0.1:2; syms t; f=t*t; a=int(f); . ；. subplot(211); ezplot(f); subplot(212); ezplot(a); 其信號(hào)圖如下： 4.1 反轉(zhuǎn)反轉(zhuǎn) 信號(hào)的反轉(zhuǎn)就是將信號(hào)的波形以某軸為對(duì)稱軸翻轉(zhuǎn)180 信號(hào)f（t）=t的反轉(zhuǎn)MATLAB程序如下： t=-10:1:10; f=t; a=fliplr(f); h=flipud(f); . ；. subplot(311); plot(t,f); axis(-2 2 -2 2); ti

13、tle(原); subplot(312); plot(t,a); axis(-5 5 -5 5); title(上下); subplot(313); plot(t,h); axis(-5 5 -5 5); title(左右); 其信號(hào)圖如下： . ；. 4.2 時(shí)移時(shí)移 實(shí)現(xiàn)連續(xù)時(shí)間信號(hào)的時(shí)移即f（t-t0）或者f（t+t0），常數(shù)t00。 余弦信號(hào)的時(shí)移實(shí)現(xiàn)程序如下： t=0:0.0001:5; y=cos(2*pi*t); y1=cos(2*pi*(t-0.2); plot(t,y,-,t,y1); axis(0 5 -1.5 1.5); 其信號(hào)圖如下： 4.3 展縮展縮 信號(hào)的展縮即將信

14、號(hào)f（t）中的自變量t替換為at，a0。 . ；. 正切信號(hào)的展縮實(shí)現(xiàn)程序如下： t=0:0.001:2; a=1; y=sin(2*pi*t); y1=subs(y,t,a*t); subplot(211); ezplot(y); subplot(212); ezplot(y1); 其信號(hào)圖如下： 4.4 倒相倒相 連續(xù)信號(hào)的倒相是指將信號(hào) f(t)以橫軸為對(duì)稱軸對(duì)折得到-f(t)。 正弦信號(hào)的展縮實(shí)現(xiàn)程序如下： t=0:0.0001:4; y=cos(2*pi*t); y1=-y; subplot(211); plot(t,y); . ；. axis(0 4 -1.5 1.5); subp

15、lot(212); plot(t,y1); axis(0 4 -1.5 1.5); 其信號(hào)圖如下： 4.5 綜合變化綜合變化 將 f(t)=cos(t)/t 通過反褶、移位、尺度變換由 f(t)的波形得到 f(-2t+3)的波形。該變化的 實(shí)現(xiàn)程序如下 syms t; f=sym(cos(t)/t); f1=subs(f,t,t+2); f2=subs(f1,t,1*t); . ；. f3=subs(f2,t,-t); subplot(221);ezplot(f,-10,10); subplot(222);ezplot(f1,-10,10); subplot(223);ezplot(f2,-1

16、0,10); subplot(224);ezplot(f3,-10,10); 其信號(hào)圖如下： 5 連續(xù)時(shí)間信號(hào)簡(jiǎn)單的時(shí)域分解連續(xù)時(shí)間信號(hào)簡(jiǎn)單的時(shí)域分解 5.1 信號(hào)的交直流分解信號(hào)的交直流分解 信號(hào)的交直流分解即將信號(hào)分解成直流分量和交流分量?jī)刹糠种?，其中直流分?定義為 fD（t）=/t f（t） . ；. 交流分量定義為 fA（t）=f（t）-fD（t） 例如對(duì)函數(shù) f（t）=cos（t）+2 進(jìn)行交直流分解。 程序如下： t=-5:0.002:5; f=cos(2*pi*t)+2; a=mean(f); b=f-a; subplot(311); plot(t,f); axis(-2 2

17、 0.5 5); subplot(312); plot(t,a); axis(-2 2 -2 5); subplot(313); plot(t,b); axis(-2 2 -1.5 1.5); 圖像如下： . ；. 5.2 信號(hào)的奇偶分解信號(hào)的奇偶分解 對(duì)函數(shù) f（t）=cos（t-5）+t 進(jìn)行交直流分解。 MATLAB 程序如下： syms t; f=sym(cos(t- 5)+t ); f1=subs(f,t,-t) a=1/2*(f+f1); b=1/2*(f-f1); subplot(311); ezplot(f,-8,8); . ；. subplot(312); ezplot(a,

18、-8,8); subplot(313); ezplot(b,-8,8); 波形圖如下圖： 6 連續(xù)時(shí)間系統(tǒng)的卷積積分的仿真波形連續(xù)時(shí)間系統(tǒng)的卷積積分的仿真波形 卷積積分在信號(hào)與線形系統(tǒng)分析中具有非常重要的意義，是信號(hào)與系統(tǒng)分析 的基本方法之一。 連續(xù)時(shí)間信號(hào) f1(t)和 f2(t)的卷積積分（簡(jiǎn)稱為卷積）f(t)定義為： f(t)= f1(t)* f2(t)= 1(t)f2(t- ) 由此可得到兩個(gè)與卷積相關(guān)的重要結(jié)論，即是： （1） f(t)= f1(t)* (t)，即連續(xù)信號(hào)可分解為一系列幅度由 f (t) 決定的沖激 . ；. 信號(hào) (t) 及其平移信號(hào)之和； （2）線形時(shí)不變連續(xù)系統(tǒng)

19、，設(shè)其輸入信號(hào)為 f (t) ，單位響應(yīng)為 h (t ) ，其零狀態(tài)響應(yīng)為 y (t)，則有：y (t ) = f (t)h (t)。 用 MATLAB 實(shí)現(xiàn)連續(xù)信號(hào) f 1(t)與 f2(t)卷積的過程如下： （1）將連續(xù)信號(hào) f 1(t)與 f2(t)以時(shí)間間隔進(jìn)行取樣，得到離散序列 f 1(k)和 f2(k)； （2）構(gòu)造與 f 1(k)和 f2(k)相對(duì)應(yīng)的時(shí)間向量 k1 和 k2 ； （3）調(diào)用 conv()函數(shù)計(jì)算卷積積分 f (t) 的近似向量 f (n)； （4）構(gòu)造 f (n)對(duì)應(yīng)的時(shí)間向量 k。 卷積實(shí)現(xiàn)程序如下： function f,k=sconv(f1,f2,k1,k

20、2,p) f=conv(f1,f2); %計(jì)算序列 f1 與 f2 的卷積和 f f=f*p; k0=k1(1)+k2(1); %計(jì)算序列 f 非零樣值的起點(diǎn)位置 k3=length(f1)+length(f2)-2; %計(jì)算卷積和 f 的非零樣值的寬度 k=k0:p:k3*p; %確定卷積和 f 非零樣值的時(shí)間向量 subplot(2,2,1) plot(k1,f1) %在子圖 1 繪 f1(t)時(shí)域波形圖 title(f1(t) xlabel(t) ylabel(f1(t) subplot(2,2,2) plot(k2,f2) %在子圖 2 繪 f2(t)時(shí)波形圖 title(f2(t)

21、xlabel(t) ylabel(f2(t) subplot(2,2,3) plot(k,f); %畫卷積 f(t)的時(shí)域波形 h=get(gca,position); . ；. h(3)=2.5*h(3); set(gca,position,h) %將第三個(gè)子圖的橫坐標(biāo)范圍擴(kuò)為原來的 2.5 倍 title(f(t)=f1(t)*f2(t) xlabel(t) ylabel(f(t) 例一： 實(shí)現(xiàn)程序如下： p=0.1; k1=0:p:2; f1=0.5*k1; k2=k1; f2=f1; f,k=doc sconv(f1,f2,k1,k2,p) 例二： . ；. 實(shí)現(xiàn)程序如下： p=0.1

22、; k1=0:p:2; f1=rectpuls(k1-1,length(k1); k2=k1; f2=f1; f,k=sconv(f1,f2,k1,k2,p) 7 連續(xù)時(shí)間系統(tǒng)的沖激響應(yīng)、階躍響應(yīng)的仿真波形連續(xù)時(shí)間系統(tǒng)的沖激響應(yīng)、階躍響應(yīng)的仿真波形 對(duì)于連續(xù)時(shí)間系統(tǒng)，求解系統(tǒng)的沖激響應(yīng) h(t)和階躍響應(yīng) g(t)對(duì)我們進(jìn)行連續(xù)系統(tǒng)的 分析具有非常重要的意義。MATLAB 為用戶提供了專門用于求連續(xù)系統(tǒng)沖激響應(yīng)和階躍 響應(yīng)并繪制其時(shí)域波形的函數(shù) impulse（）和 step（）。 在調(diào)用 impulse（）和 step（）函數(shù)時(shí)，我們需要用向量來對(duì)連續(xù)時(shí)間系統(tǒng)進(jìn)行分析。 設(shè)描述連續(xù)系統(tǒng)的微分

23、方程為： iy(i)(t)=jx(j)(t) = 0 = 0 則我們可用向量 A 和 B 來表示該系統(tǒng)，即： A=AN,AN-1,A1,A0 B=BN,BN-1,B1,B0 注意，向量 A 和 B 的元素一定要以微分方程中時(shí)間求導(dǎo)的降冪次序來排列，且缺項(xiàng) 要用 0 來補(bǔ)齊。例如，對(duì)微分方程，則表示該系統(tǒng)的對(duì) () + 3() + 2() = () + () 應(yīng)向量應(yīng)為 A=1 3 2，B=1 0 1。 . ；. 7.1 impulse（）函數(shù)（）函數(shù) 函數(shù) impulse（）將繪出由向量 a 和 b 表示的連續(xù)系統(tǒng)在指定時(shí)間范圍內(nèi)的沖激響 應(yīng) h(t)的時(shí)域波形圖，并能求出指定時(shí)間范圍內(nèi)沖激響

24、應(yīng)的數(shù)值解。 impulse（）函數(shù)有如下四種調(diào)用格式： （1）impulse(b,a)：該調(diào)用格式以默認(rèn)方式繪出向量 A 和 B 定義的連續(xù)系統(tǒng)的沖激 響應(yīng)的時(shí)域波形 。例如描述連續(xù)系統(tǒng)的微分方程為 () + 5() + 6() = 3() + 2() 運(yùn)行如下 MATLAB 命令： a=1 5 6; b=3 2; impulse(b,a); 則繪出系統(tǒng)的沖激響應(yīng)波形，如圖所示： . ；. （2）impulse(b,a,t)：運(yùn)行命令 impulse(b,a,6)，則繪出系統(tǒng)在 06 秒范圍內(nèi)沖激響 應(yīng)的時(shí)域波形，如圖所示 （3）impulse(b,a,t1:p:t2)：若運(yùn)行命令 impu

25、lse(b,a,1:0.1:3)，則繪出 13 秒內(nèi)，每隔 0.1 秒取樣的沖激響應(yīng)的時(shí)域波形，如圖所示： . ；. （4）y=impulse(b,a,t1:p:t2)：若運(yùn)行命令 y=impulse(b,a,0:0.2:3)，則運(yùn)行結(jié)果為： . ；. 7.2 step（）函數(shù)（）函數(shù) step（）函數(shù)可繪出連續(xù)系統(tǒng)的階躍響應(yīng) g(t)在指定時(shí)間范圍的時(shí)域波形并能求出其 數(shù)值解，和 impulse（）函數(shù)一樣也有四種調(diào)用格式。 () + 5() + 6() = 3() + 2() 運(yùn)行如下 MATLAB 命令： a=1 5 6; b=3 2; step(b,a); 則繪出系統(tǒng)的階躍響應(yīng)波形，如

26、圖所示。 . ；. （2）step(b,a,t)：運(yùn)行命令 step(b,a,6)，繪出在 06 秒范圍內(nèi)階躍響應(yīng)的時(shí)域波形 . ；. （3）step(b,a,t1:p:t2)：運(yùn)行命令 step(b,a,1:0.1:3)，繪出 13 秒內(nèi)，每隔 0.1 秒取樣 的階躍響應(yīng)的時(shí)域波形 . ；. （4）y=step(b,a,t1:p:t2)：運(yùn)行命令 y=step(b,a,0:0.2:3)，則運(yùn)行結(jié)果為： . ；. 8 連續(xù)時(shí)間系統(tǒng)對(duì)正弦信號(hào)、實(shí)指數(shù)信號(hào)的零狀態(tài)響應(yīng)的仿真連續(xù)時(shí)間系統(tǒng)對(duì)正弦信號(hào)、實(shí)指數(shù)信號(hào)的零狀態(tài)響應(yīng)的仿真 波形波形 MATLAB 中的函數(shù) lsim（）能對(duì)微分方程描述的 LTI

27、 連續(xù)時(shí)間系統(tǒng)的響應(yīng)進(jìn)行仿真。 該函數(shù)能繪制連續(xù)時(shí)間系統(tǒng)在指定的任意時(shí)間范圍內(nèi)系統(tǒng)響應(yīng)的時(shí)域波形圖，還能求出 連續(xù)時(shí)間系統(tǒng)在指定的任意時(shí)間范圍內(nèi)系統(tǒng)響應(yīng)的數(shù)值解，函數(shù) lsim（）的調(diào)用格式如下： lsim(b,a,x,t) 在該調(diào)用格式中，a 和 b 是由描述系統(tǒng)的微分方程系統(tǒng)決定的表示該系統(tǒng)的兩個(gè)行向 量。x 和 t 則是表示輸入信號(hào)的行向量，其中 t 為表示輸入信號(hào)時(shí)間范圍的向量，x 則是 輸入信號(hào)在向量 t 定義的時(shí)間點(diǎn)上的抽樣值。該調(diào)用格式將繪出向量 b 和 a 所定義的連續(xù) 系統(tǒng)在輸入量為向量 x 和 t 所定義的信號(hào)時(shí)，系統(tǒng)的零狀態(tài)響應(yīng)的時(shí)域仿真波形，且時(shí)間 范圍與輸入信號(hào)相同。

28、 8.1 正弦信號(hào)的零狀態(tài)響應(yīng)正弦信號(hào)的零狀態(tài)響應(yīng) MATLAB 命令如下： a=1,2,1; b=1,2; p=0.5; t=0:p:5; x=sin(3*pi*t); lsim(b,a,x,t); hold on; p=1.0; t=0:p:8; x=sin(3*pi*t); lsim(b,a,x,t); p=0.01; t=0:p:5; . ；. x=sin(3*pi*t); lsim(b,a,x,t); hold off; 8.2 實(shí)指數(shù)信號(hào)的零狀態(tài)響應(yīng)實(shí)指數(shù)信號(hào)的零狀態(tài)響應(yīng) MATLAB 命令如下： a=1,3,1; b=1,3; p=0.5; t=0:p:5; x=exp(-2*t); lsim(b,a,x,t); . ；. hold on; p=0.2; t=0:p:5; x=e