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

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

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

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

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

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

6、三角波。(8) 下面的函數(shù)也可以產(chǎn)生特定的信號(hào)y=fliplr(t)信號(hào)的翻轉(zhuǎn)；y=sawtooth(t,width)產(chǎn)生鋸齒波信號(hào)。實(shí)驗(yàn)內(nèi)容1. 實(shí)例分析與驗(yàn)證根據(jù)以上典型信號(hào)的MATLAB函數(shù)，分析與驗(yàn)證下列典型信號(hào)MATLAB程序，并實(shí)現(xiàn)各信號(hào)波形圖的顯示，連續(xù)信號(hào)的圖形顯示使用連續(xù)二維圖函數(shù)plot()。(1) 正弦信號(hào)：用MATLAB命令產(chǎn)生正弦信號(hào)，并會(huì)出時(shí)間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.2)title(正弦信號(hào)

7、)(2) 抽樣信號(hào)：用MATLAB中的sinc(t)函數(shù)命令產(chǎn)生抽樣信號(hào)Sa(t)，并會(huì)出時(shí)間為的波形圖。程序如下：t=-6*pi:pi/100:6*pi;ft=sinc(t/pi);plot(t,ft),grid on;axis(-20,20,-0.5,1.2)title(抽樣信號(hào))(3) 矩形脈沖信號(hào)：用MATLAB中y=rectpuls(t,width)函數(shù)命令畫(huà)出下列矩形脈沖信號(hào)的波形圖程序如下：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,3,-0.2,2.2)t

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

9、令畫(huà)出0 t 3復(fù)指數(shù)信號(hào)的實(shí)部、虛部、模及相角隨時(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(2,2,4);plo

10、t(t,angle(ft);title(相角);axis(0,3,-4,4);grid on;(2) 用MATLAB中y=square(t,DUTY)產(chǎn)生頻率為10Hz、占空比為30%的周期方波信號(hào)。程序如下：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(周期方波信號(hào))(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;axis(-1,1,-0.2

11、,1.2)title(門函數(shù))實(shí)驗(yàn)分析觀察實(shí)驗(yàn)結(jié)果，掌握、分析典型的連續(xù)時(shí)間信號(hào)波形特點(diǎn)。實(shí)驗(yàn)總結(jié)總結(jié)實(shí)驗(yàn)認(rèn)識(shí)、過(guò)程、效果、問(wèn)題、收獲、體會(huì)、意見(jiàn)和建議。實(shí)驗(yàn)二 連續(xù)時(shí)間LTI系統(tǒng)的時(shí)域分析實(shí)驗(yàn)?zāi)康?運(yùn)用MATLAB符號(hào)求解連續(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ù)時(shí)間系統(tǒng)零輸入響應(yīng)和零狀態(tài)響應(yīng)的符號(hào)求解LTI連續(xù)系統(tǒng)可用線性常系數(shù)微分方程描述，即MATLAB符號(hào)工具箱提供了dsolve函數(shù)，可實(shí)現(xiàn)常系數(shù)微分方程的符號(hào)求解，其調(diào)用格式

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

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

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

15、分別是參與計(jì)算的兩個(gè)時(shí)間信號(hào)f=conv(f1,f2); % t1,t2是兩信號(hào)時(shí)間變量，dt為時(shí)間采樣間隔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),min(f2)-abs(mi

16、n(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ù)時(shí)間系統(tǒng)零輸入響應(yīng)和零狀態(tài)響應(yīng)的符號(hào)求解試用MATLAB命令求解微分方程，當(dāng)輸入，起始條件為、時(shí)系統(tǒng)的零輸入響應(yīng)、零狀態(tài)響應(yīng)及完全響應(yīng)。源程序?yàn)椋篹q=D2y+3*

17、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+yzs)運(yùn)行結(jié)果： yt

18、= -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ù)時(shí)間系統(tǒng)零狀態(tài)響應(yīng)的數(shù)值求解試用MATLAB數(shù)值求解微分方程，當(dāng)輸入時(shí)系統(tǒng)的零狀態(tài)響應(yīng)。源程序?yàn)椋簍s=0;te=8

19、;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ù)時(shí)間系統(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);subplot(211);

20、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,dt);2. 實(shí)踐

21、編程(1) 已知系統(tǒng)的微分方程和激勵(lì)信號(hào)為，試用MATLAB命令繪出系統(tǒng)零狀態(tài)響應(yīng)的時(shí)域仿真波形圖。源程序?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)的時(shí)域仿真波形圖。源程序?yàn)椋簍=0:0.001:4;sys=tf(1,1,3,2);h=impulse(sys,t

22、);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ù)時(shí)間系統(tǒng)零狀態(tài)響應(yīng)和零輸入響應(yīng)的求解方法。實(shí)驗(yàn)總結(jié)總結(jié)實(shí)驗(yàn)認(rèn)識(shí)、過(guò)程、效果、問(wèn)題、收獲、體會(huì)、意見(jiàn)和建議。實(shí)驗(yàn)三 連續(xù)時(shí)間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ù)時(shí)間LTI系統(tǒng)的頻

23、率特性一個(gè)連續(xù)時(shí)間LTI系統(tǒng)的數(shù)學(xué)模型通常用常系數(shù)線性微分方程來(lái)描述，即對(duì)上式兩邊取傅里葉變換，并根據(jù)傅里葉變換的時(shí)域微分特性，得到系統(tǒng)的頻率響應(yīng)為MATLAB信號(hào)處理工具箱提供的freqs函數(shù)可直接計(jì)算系統(tǒng)的頻率響應(yīng)的數(shù)值解，其語(yǔ)句格式為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和angle

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

25、態(tài)響應(yīng)為實(shí)驗(yàn)內(nèi)容1. 實(shí)例分析與驗(yàn)證(1) 已知一個(gè)連續(xù)時(shí)間LTI系統(tǒng)的微分方程為求系統(tǒng)的頻率響應(yīng)，并用MATLAB繪出其幅頻特性和相頻特性圖。解：對(duì)微分方程取傅里葉變換，得因此，頻率響應(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),yl

26、abel(phi(omega)|)title(H(w)的相頻特性)(2) 如圖為RC低通濾波器電路，在輸入端加入矩形脈沖u1(t)。利用傅里葉分析法求輸出端電壓u2(t)。Ru2(t)u1(t)tt(2)120kRu(C10Fu) 解：RC低通濾波器的頻率響應(yīng)為 其中激勵(lì)信號(hào)的傅里葉變換為 因此，響應(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=sym(

27、Heaviside(t)-Heaviside(t-1);U1w=simplify(fourier(u1t);figuresubplot(221),ezplot(u1t,-0.2,2),grid ontitle(矩形脈沖信號(hào))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);s

28、ubplot(223),ezplot(u2t,-0.2 2),grid ontitle(響應(yīng)的時(shí)域波形)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)為 ，若外加激勵(lì)，用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(輸入信號(hào)的波形)subplot(212);plot(t,y);