2022年信號與系統(tǒng)實(shí)驗(yàn)指導(dǎo)書

1、.精選文檔精選文檔.精選文檔?信號與系統(tǒng)?實(shí)驗(yàn)指導(dǎo)書重慶科技學(xué)院 電氣與信息工程學(xué)院2021年3月目 錄實(shí)驗(yàn)一 連續(xù)時間信號的MATLAB表示實(shí)驗(yàn)二 連續(xù)時間LTI系統(tǒng)的時域分析實(shí)驗(yàn)三 連續(xù)時間LTI系統(tǒng)的頻率特性及頻域分析實(shí)驗(yàn)四 連續(xù)時間LTI系統(tǒng)的零極點(diǎn)分析實(shí)驗(yàn)五 典型離散信號及其MATLAB實(shí)現(xiàn)實(shí)驗(yàn)一 連時間信號的MATLAB表示實(shí)驗(yàn)?zāi)康?掌握MATLAB語言的根本操作，學(xué)習(xí)根本的編程功能；2掌握MATLAB產(chǎn)生常用連續(xù)時間信號的編程方法；3觀察并熟悉常用連續(xù)時間信號的波形和特性。實(shí)驗(yàn)原理1. 連續(xù)信號MATLAB實(shí)現(xiàn)原理從嚴(yán)格意義上講，MATLAB數(shù)值計(jì)算的方法并不能處理連續(xù)時間信號。

2、然而，可用連續(xù)信號在等時間間隔點(diǎn)的取樣值來近似表示連續(xù)信號，即當(dāng)取樣時間間隔足夠小時，這些離散樣值能夠被MATLAB處理，并且能較好地近似表示連續(xù)信號。MATLAB提供了大量生成根本信號的函數(shù)。比方常用的指數(shù)信號、正余弦信號等都是MATLAB的內(nèi)部函數(shù)。為了表示連續(xù)時間信號，需定義某一時間或自變量的范圍和取樣時間間隔，然后調(diào)用該函數(shù)計(jì)算這些點(diǎn)的函數(shù)值，最后畫出其波形圖。2. 典型信號的MATLAB表示(1) 實(shí)指數(shù)信號：，式中K，a為實(shí)數(shù)。在MATLAB中可以利用exp()函數(shù)實(shí)現(xiàn)，其語句格式為：y=K*exp(a*t)。(2) 正弦信號：或，其中K為振幅，是角頻率；是初相位。在MATLAB中

3、可以利用sin()或cos()函數(shù)實(shí)現(xiàn)，其語句格式為K*sin(*t+phi)，K*cos(*t+phi)。(3) 單位階躍信號：單位階躍信號定義為，單位階躍信號在MATLAB中用“( t =0)產(chǎn)生。MATLAB表達(dá)式“y = ( t =0)的含義就是t 0時y =1，而當(dāng)t =0);保存為uCT.m后，就可以調(diào)用該函數(shù)。 自定義函數(shù)注意，在此定義的uCT函數(shù)是階躍信號數(shù)值表示方法，因此在數(shù)值計(jì)算中我們將調(diào)用uCT函數(shù)。而在MATLAB的MAPLE內(nèi)核中，將Heaviside函數(shù)定義為階躍信號符號表達(dá)式，在符號運(yùn)算過程中，假設(shè)要調(diào)用它必須用sym定義后，才能實(shí)現(xiàn)。(4) 復(fù)指數(shù)信號：，其中，

4、s=+j是復(fù)變量；、為實(shí)數(shù)。在MATLAB中可以利用exp()函數(shù)表示復(fù)指數(shù)信號，MATLAB中默認(rèn)變量i為虛部單位，求取實(shí)部用real()函數(shù)，求取虛部用imag()函數(shù)，求取模用abs()函數(shù)，求取相角用angle()函數(shù)。 (5) 抽樣信號:，在MATLAB中用與Sa(t)類似的sinc(t)函數(shù)表示，定義為。(6) 矩形脈沖信號：在MATLAB中可用rectpuls()函數(shù)產(chǎn)生，其語句格式為y=rectpuls(t,width)，該函數(shù)用于產(chǎn)生一個幅度為1、寬度為width，且以t=0為對稱軸的矩形脈沖信號，width的默認(rèn)值為1。周期性矩形波信號或方波在MATLAB中可用square

5、函數(shù)產(chǎn)生，其語句格式為square(t,DUTY)。該函數(shù)用于產(chǎn)生一個周期為、幅值為的周期性方波信號，其中，DUTY參數(shù)用來表示信號的占空比DUTY%，即在一個周期內(nèi)脈沖寬度正值局部與脈沖周期的比值。占空比默認(rèn)值為0.5。(7) 三角波脈沖信號：在MATLAB中可用tripuls()函數(shù)產(chǎn)生，其語句格式為y=tripuls(t,width,skew)，該函數(shù)用于產(chǎn)生一個幅度為1、寬度為width，且以t=0為中心左右各展開width/2大小、斜度為skew的三角波。width的默認(rèn)值為1，skew的取值范圍在-1+1之間。一般最大幅度1出現(xiàn)在t=(width/2)skew的橫坐標(biāo)位置，默認(rèn)時s

6、kew=0，此時產(chǎn)生對稱三角波。(8) 下面的函數(shù)也可以產(chǎn)生特定的信號y=fliplr(t)信號的翻轉(zhuǎn)；y=sawtooth(t,width)產(chǎn)生鋸齒波信號。實(shí)驗(yàn)內(nèi)容1. 實(shí)例分析與驗(yàn)證根據(jù)以上典型信號的MATLAB函數(shù)，分析與驗(yàn)證以下典型信號MATLAB程序，并實(shí)現(xiàn)各信號波形圖的顯示，連續(xù)信號的圖形顯示使用連續(xù)二維圖函數(shù)plot()。(1) 正弦信號：用MATLAB命令產(chǎn)生正弦信號，并會出時間0 t 3的波形圖。程序如下：K=2;w=2*pi;phi=pi/4;t=0:0.01:3;ft=K*sin(w*t+phi); plot(t,ft),grid on; axis(0,3,-2.2,2.

7、2)title(正弦信號)(2) 抽樣信號：用MATLAB中的sinc(t)函數(shù)命令產(chǎn)生抽樣信號Sa(t)，并會出時間為的波形圖。程序如下：t=-6*pi:pi/100:6*pi;ft=sinc(t/pi);plot(t,ft),grid on;axis(-20,20,-0.5,1.2)title(抽樣信號)(3) 矩形脈沖信號：用MATLAB中y=rectpuls(t,width)函數(shù)命令畫出以下矩形脈沖信號的波形圖程序如下：t=-0.5:0.01:3;t0=0.5;width=1;ft=2*rectpuls(t-t0,width);plot(t,ft),grid on;axis(-0.5,

8、3,-0.2,2.2)title(矩形脈沖信號)(4) 單位階躍信號：用MATLAB命令“y=(t=0)繪出-1 t 5單位階躍信號(t)。程序如下：t=-1:0.01:5;ft=(t=0);plot(t,ft),grid on;axis(-1,5,-0.5,1.5);title(單位階躍信號)也可通過定義階躍函數(shù)來實(shí)現(xiàn)，function f=uCT(t);f=(t=0);保存為uCT.m文件，上面實(shí)現(xiàn)階躍信號的程序變?yōu)閠=-1:0.01:5;ft=uCT(t);plot(t,ft),grid on;axis(-1,5,-0.5,1.5);title(單位階躍信號)2. 實(shí)踐編程(1) 復(fù)指數(shù)

9、信號：用MATLAB命令畫出0 t 3復(fù)指數(shù)信號的實(shí)部、虛部、模及相角隨時間變化的曲線，并觀察其時域特性。程序如下：t=0:0.01:3;k=2;a=-1.5;b=10;ft=k*exp(a+i*b)*t);subplot(2,2,1);plot(t,real(ft);title(實(shí)部);axis(0,3,-2,2);grid on;subplot(2,2,2);plot(t,imag(ft);title(虛部);axis(0,3,-2,2);grid on;subplot(2,2,3);plot(t,abs(ft);title(模);axis(0,3,0,2);grid on;subplot

10、(2,2,4);plot(t,angle(ft);title(相角);axis(0,3,-4,4);grid on;(2) 用MATLAB中y=square(t,DUTY)產(chǎn)生頻率為10Hz、占空比為30%的周期方波信號。程序如下：t=0:0.001:3;y=square(2*pi*10*t,30);plot(t,y);grid on;axis(0,0.3,-1.2,1.2);title(周期方波信號)(3) 用上面定義的uCT階躍函數(shù)實(shí)現(xiàn)幅度為1、寬度為1的門函數(shù)。程序如下：t=-1:0.01:1;ft=uCT(t+0.5)-uCT(t-0.5);plot(t,ft);grid on;axi

11、s(-1,1,-0.2,1.2)title(門函數(shù))實(shí)驗(yàn)分析觀察實(shí)驗(yàn)結(jié)果，掌握、分析典型的連續(xù)時間信號波形特點(diǎn)。實(shí)驗(yàn)總結(jié)總結(jié)實(shí)驗(yàn)認(rèn)識、過程、效果、問題、收獲、體會、意見和建議。實(shí)驗(yàn)二 連續(xù)時間LTI系統(tǒng)的時域分析實(shí)驗(yàn)?zāi)康?運(yùn)用MATLAB符號求解連續(xù)系統(tǒng)的零輸入響應(yīng)和零狀態(tài)響應(yīng)； 2運(yùn)用MATLAB數(shù)值求解連續(xù)系統(tǒng)的零狀態(tài)響應(yīng)；3運(yùn)用MATLAB求解連續(xù)系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)；4運(yùn)用MATLAB卷積積分法求解系統(tǒng)的零狀態(tài)響應(yīng)。實(shí)驗(yàn)原理1. 連續(xù)時間系統(tǒng)零輸入響應(yīng)和零狀態(tài)響應(yīng)的符號求解LTI連續(xù)系統(tǒng)可用線性常系數(shù)微分方程描述，即MATLAB符號工具箱提供了dsolve函數(shù)，可實(shí)現(xiàn)常系數(shù)微分方程

12、的符號求解，其調(diào)用格式為dsolve(eq1,eq2,cond1,cond2,v)其中，參數(shù)eq1，eq2、表示各微分方程，它與MATLAB符號表達(dá)式的輸入根本相同，微分或?qū)?shù)的輸入是用Dy、D2y、D3y、來表示y的一階導(dǎo)數(shù)、二階導(dǎo)數(shù)、三階導(dǎo)數(shù)、；參數(shù)cond1、cond2、表示各初始條件或起始條件；參數(shù)v表示自變量，默認(rèn)為是變量t。可利用dsolve函數(shù)來求解系統(tǒng)微分方程的零輸入響應(yīng)和零狀態(tài)響應(yīng)，進(jìn)而求出完全響應(yīng)。2. 連續(xù)時間系統(tǒng)零狀態(tài)響應(yīng)的數(shù)值求解實(shí)際工程中用的較多求解微分方程的方法是數(shù)值求解。零輸入響應(yīng)的求解可通過函數(shù)initial來實(shí)現(xiàn)，initial函數(shù)中的參量必須是狀態(tài)變量所描

13、述的系統(tǒng)模型。對于零狀態(tài)響應(yīng)，MATLAB控制系統(tǒng)工具箱提供了對LTI系統(tǒng)的零狀態(tài)響應(yīng)進(jìn)行數(shù)值仿真的函數(shù)lsim，該函數(shù)可求解零初始條件下微分方程的數(shù)值解，語句格式為y=lsim(sys,f,t)其中，t表示計(jì)算系統(tǒng)響應(yīng)的時間抽樣點(diǎn)向量；f是系統(tǒng)的輸入信號向量；sys表示LTI系統(tǒng)模型，用來表示微分方程、差分方程或狀態(tài)方程。在求解微分方程時，sys是由MATLAB的tf函數(shù)根據(jù)微分方程系數(shù)生成的系統(tǒng)函數(shù)對象，其語句格式為sys=tf(b,a)其中，b和a分別為微分方程右端和左端的系數(shù)向量。注意，如果微分方程的左端或右端表達(dá)式中有缺項(xiàng)，那么其向量a或b中的對應(yīng)元素應(yīng)為0，不能省略不寫，否那么會出

14、錯。3. 連續(xù)時間系統(tǒng)沖激響應(yīng)和階躍響應(yīng)的求解在MATLAB中，對于連續(xù)LTI系統(tǒng)的沖激響應(yīng)和階躍響應(yīng)的數(shù)值解，可分別用控制系統(tǒng)工具箱提供的函數(shù)impulse和step來求解。其語句格式為y=impulse(sys,t)y=step(sys,t)其中，t表示計(jì)算系統(tǒng)響應(yīng)的時間抽樣點(diǎn)向量，sys表示LTI系統(tǒng)模型。4. 利用卷積積分法求系統(tǒng)的零狀態(tài)響應(yīng)由卷積積分公式可以得到，LTI系統(tǒng)對于任意輸入信號的零狀態(tài)響應(yīng)，可由系統(tǒng)的單位沖激響應(yīng)與輸入信號的卷積積分得到。預(yù)先建立連續(xù)時間信號卷積運(yùn)算的函數(shù)ctsconv.m，其MATLAB源程序?yàn)閒unctionf,t=ctsconv(f1,f2,t1,t

15、2,dt) % f1和f2分別是參與計(jì)算的兩個時間信號f=conv(f1,f2); % t1,t2是兩信號時間變量，dt為時間采樣間隔f=f*dt;ts=min(t1)+min(t2);te=max(t1)+max(t2);t=ts:dt:te;subplot(221)plot(t1,f1);grid on;axis(min(t1),max(t1),min(f1)-abs(min(f1)*0.2),max(f1)+abs(max(f1)*0.2)title(f1(t);xlabel(t)subplot(222)plot(t2,f2);grid on;axis(min(t2),max(t2),m

16、in(f2)-abs(min(f2)*0.2),max(f2)+abs(max(f2)*0.2)title(f2(t);xlabel(t)subplot(212)plot(t,f);grid on;axis(min(t),max(t),min(f)-abs(min(f)*0.2),max(f)+abs(max(f)*0.2)title(f(t)= f1(t)* f2(t);xlabel(t)在實(shí)際計(jì)算中，調(diào)入即可應(yīng)用。實(shí)驗(yàn)內(nèi)容1. 實(shí)例分析與驗(yàn)證(1) 連續(xù)時間系統(tǒng)零輸入響應(yīng)和零狀態(tài)響應(yīng)的符號求解試用MATLAB命令求解微分方程，當(dāng)輸入，起始條件為、時系統(tǒng)的零輸入響應(yīng)、零狀態(tài)響應(yīng)及完全響應(yīng)。源

17、程序?yàn)椋篹q=D2y+3*Dy+2*y=0;cond=y(0)=1,Dy(0)=2;yzi=dsolve(eq,cond); yzi=simplify(yzi)運(yùn)行結(jié)果：yzi = -3*exp(-2*t)+4*exp(-t)eq1=D2y+3*Dy+2*y=Dx+3*x;eq2=x=exp(-3*t)*Heaviside(t);cond=y(-0.001)=0,Dy(-0.001)=0;yzs=dsolve(eq1,eq2,cond);yzs=simplify(yzs.y)運(yùn)行結(jié)果：yzs = heaviside(t)*(-exp(-2*t)+exp(-t) yt=simplify(yzi+

18、yzs)運(yùn)行結(jié)果： yt = -3*exp(-2*t)+4*exp(-t)-exp(-2*t)*heaviside(t)+exp(-t)*heaviside(t) subplot(311)ezplot(yzi,0,8);grid on;xlabel(t)title(零輸入響應(yīng))subplot(312)ezplot(yzs,0,8);grid onxlabel(t)title(零狀態(tài)響應(yīng))subplot(313)ezplot(yt,0,8);grid onxlabel(t)title(完全響應(yīng))(2) 連續(xù)時間系統(tǒng)零狀態(tài)響應(yīng)的數(shù)值求解試用MATLAB數(shù)值求解微分方程，當(dāng)輸入時系統(tǒng)的零狀態(tài)響應(yīng)。源

19、程序?yàn)椋簍s=0;te=8;dt=0.01;sys=tf(1,3,1,3,2);t=ts:dt:te;f=exp(-3*t).*uCT(t);y=lsim(sys,f,t);plot(t,y),grid on;axis(0 8 -0.02 0.27)xlabel(t(s),ylabel(y(t)title(零狀態(tài)響應(yīng))(3) 連續(xù)時間系統(tǒng)沖激響應(yīng)和階躍響應(yīng)的求解某LTI系統(tǒng)的微分方程，試用MATLAB命令繪出范圍內(nèi)系統(tǒng)的沖激響應(yīng)h(t)和階躍響應(yīng)s(t)。源程序?yàn)椋簍=0:0.001:4;sys=tf(1,16,1,2,32);h=impulse(sys,t);s=step(sys,t);su

20、bplot(211);plot(t,h),grid onxlabel(t(s),ylabel(h(t)title(沖激響應(yīng))subplot(212);plot(t,s),grid onxlabel(t(s),ylabel(s(t)title(階躍響應(yīng))(4) 利用卷積積分法求系統(tǒng)的零狀態(tài)響應(yīng)某LTI系統(tǒng)的微分方程，其中，。試用MATLAB卷積積分方法繪出系統(tǒng)零狀態(tài)響應(yīng)y(t)的波形圖。程序如下：dt=0.01;t1=0:dt:4;f1=exp(-2*t1);t2=t1;sys=tf(1,16,1,2,32);f2=impulse(sys,t2);t,f=ctsconv(f1,f2,t1,t2,

21、dt);2. 實(shí)踐編程(1) 系統(tǒng)的微分方程和鼓勵信號為，試用MATLAB命令繪出系統(tǒng)零狀態(tài)響應(yīng)的時域仿真波形圖。源程序?yàn)椋簍s=0;te=8;dt=0.01;sys=tf(1,3,1,4,4);t=ts:dt:te;f=exp(-t).*uCT(t);y=lsim(sys,f,t);plot(t,y),grid on;xlabel(t(s),ylabel(y(t)title(零狀態(tài)響應(yīng))(2) 系統(tǒng)的微分方程為，試用MATLAB命令求系統(tǒng)沖激響應(yīng)和階躍響應(yīng)的數(shù)值解，并繪出沖激響應(yīng)和階躍響應(yīng)的時域仿真波形圖。源程序?yàn)椋簍=0:0.001:4;sys=tf(1,1,3,2);h=impulse(

22、sys,t);s=step(sys,t);subplot(211)plot(t,h),grid onxlabel(t(s),ylabel(h(t)title(沖激響應(yīng))subplot(212)plot(t,s),grid onxlabel(t(s),ylabel(s(t)title(階躍響應(yīng))實(shí)驗(yàn)分析觀察實(shí)驗(yàn)結(jié)果，掌握、分析連續(xù)時間系統(tǒng)零狀態(tài)響應(yīng)和零輸入響應(yīng)的求解方法。實(shí)驗(yàn)總結(jié)總結(jié)實(shí)驗(yàn)認(rèn)識、過程、效果、問題、收獲、體會、意見和建議。實(shí)驗(yàn)三 連續(xù)時間LTI系統(tǒng)的頻率特性及頻域分析實(shí)驗(yàn)?zāi)康?運(yùn)用MATLAB分析連續(xù)系統(tǒng)的頻率特性；2運(yùn)用MATLAB進(jìn)行連續(xù)系統(tǒng)的頻域分析。實(shí)驗(yàn)原理1. 連續(xù)時間LT

23、I系統(tǒng)的頻率特性一個連續(xù)時間LTI系統(tǒng)的數(shù)學(xué)模型通常用常系數(shù)線性微分方程來描述，即對上式兩邊取傅里葉變換，并根據(jù)傅里葉變換的時域微分特性，得到系統(tǒng)的頻率響應(yīng)為MATLAB信號處理工具箱提供的freqs函數(shù)可直接計(jì)算系統(tǒng)的頻率響應(yīng)的數(shù)值解，其語句格式為H=freqs(b,a,w)其中，b和a分別表示H(w)的分子和分母多項(xiàng)式的系數(shù)向量；w為系統(tǒng)頻率響應(yīng)的頻率范圍，其一般形式為w1：p: w2，w1為頻率起始值，w2為頻率終止值，p為頻率取樣間隔。H返回w所定義的頻率點(diǎn)上系統(tǒng)響應(yīng)頻率響應(yīng)的樣值。注意，H返回的樣值可能為包含實(shí)部和虛部的復(fù)數(shù)。因此，如果想得到系統(tǒng)的幅頻特性或相頻特性，還需利用abs和

24、angle函數(shù)來分別求得。2. 連續(xù)時間LTI系統(tǒng)的頻域分析連續(xù)LTI系統(tǒng)的頻域分析法，也稱為傅里葉變換分析法。該方法是基于信號頻譜分析的概念，討論信號作用于線性系統(tǒng)時在頻域中求解響應(yīng)的方法。傅里葉分析法的關(guān)鍵是求系統(tǒng)的頻率響應(yīng)。傅里葉分析法主要用來分析系統(tǒng)的頻率響應(yīng)特性，或分析輸出信號的頻譜，也可用來求解正弦信號作用下的穩(wěn)態(tài)響應(yīng)。對于周期信號鼓勵而言，可首先將周期信號進(jìn)行傅里葉級數(shù)展開，然后求系統(tǒng)在各傅里葉級數(shù)分解的頻率分量作用下系統(tǒng)的穩(wěn)態(tài)響應(yīng)分量，再由系統(tǒng)的線性性質(zhì)將這些穩(wěn)態(tài)響應(yīng)分量疊加，從而得出系統(tǒng)總的響應(yīng)。該方法的理論根底是基于正弦信號作用下系統(tǒng)的正弦穩(wěn)態(tài)響應(yīng)。對于正弦鼓勵信號，當(dāng)經(jīng)過

25、系統(tǒng)，其穩(wěn)態(tài)響應(yīng)為實(shí)驗(yàn)內(nèi)容1. 實(shí)例分析與驗(yàn)證(1) 一個連續(xù)時間LTI系統(tǒng)的微分方程為求系統(tǒng)的頻率響應(yīng)，并用MATLAB繪出其幅頻特性和相頻特性圖。解：對微分方程取傅里葉變換，得因此，頻率響應(yīng)為MATLAB源程序：w=-3*pi:0.01:3*pi;b=13,7;a=1,10,8,5;H=freqs(b,a,w);subplot(211)plot(w,abs(H),grid onxlabel(omega(rad/s),ylabel(|H(omega)|)title(H(w)的頻率特性)subplot(212)plot(w,angle(H),grid onxlabel(omega(rad/s)

26、,ylabel(phi(omega)|)title(H(w)的相頻特性)(2) 如圖為RC低通濾波器電路，在輸入端參加矩形脈沖u1(t)。利用傅里葉分析法求輸出端電壓u2(t)。Ru2(t)u1(t)tt(2)120kRu(C10Fu) 解：RC低通濾波器的頻率響應(yīng)為 其中鼓勵信號的傅里葉變換為 因此，響應(yīng)的傅里葉變換為MATLAB源程序：w=-6*pi:0.01:6*pi;b=5;a=1,5;H1=freqs(b,a,w);plot(w,abs(H1),grid onxlabel(omega(rad/s),ylabel(|H(omega)|)title(RC低通濾波電路的幅頻特性)u1t=s

27、ym(Heaviside(t)-Heaviside(t-1);U1w=simplify(fourier(u1t);figuresubplot(221),ezplot(u1t,-0.2,2),grid ontitle(矩形脈沖信號)xlabel(Time(sec),ylabel(u1(t)subplot(222),ezplot(abs(U1w),-6*pi 6*pi),grid ontitle(矩形脈沖的頻譜)xlabel(omega(rad/s),ylabel(X(omega)U2w=sym(5*(1-exp(-i*w)/(5*i*w-w2);u2t=simplify(ifourier(U2w

28、);subplot(223),ezplot(u2t,-0.2 2),grid ontitle(響應(yīng)的時域波形)xlabel(Time(sec),ylabel(u2(t) subplot(224),ezplot(abs(U2w),-6*pi 6*pi),grid ontitle(響應(yīng)的頻譜)xlabel(omega(rad/s),ylabel(U2(omega)2. 實(shí)踐編程設(shè)系統(tǒng)的頻率響應(yīng)為 QUOTE ，假設(shè)外加鼓勵，用MATLAB求其穩(wěn)態(tài)響應(yīng)。解：MATLAB源程序：t=0:0.1:20;w1=1;w2=10; H1=1/(-w12+j*3*w1+2);H2=1/(-w22+j*3*w2+2);f=5*cos(w1*t)+2*cos(w2*t);y=5*abs(H1)*cos(w1*t+angle(H1)+2*abs(H2)*cos(w2*t+angle(H2);subplot(211);plot(t,f);grid onylabel(f(t),xlabel(Time(sec)title(輸入信號的波形)subplot(212);plot(t,