信號實驗報告

1、信號與系統(tǒng)實驗報告課程名稱：信號與系統(tǒng)實驗實驗時間：10 至 11 月任課教師：王群實驗地點：4-427實驗教師：何冰松, 范哲意實驗類型： 原理驗證 綜合設(shè)計 創(chuàng)新學(xué)生尹迪學(xué)號/班級：1120131225/05111361組號：學(xué)院：信息與電子學(xué)院同組搭檔：專業(yè)：信息工程成績：信號與系統(tǒng)實驗報告實驗名稱：信號的時域描述與運算課程名稱：信號與系統(tǒng)實驗實驗時間：10 月 15 日任課教師：王群實驗地點：4-427實驗教師：何冰松, 范哲意實驗類型： 原理驗證 綜合設(shè)計 創(chuàng)新學(xué)生尹迪學(xué)號/班級：1120131225/05111361組號：學(xué)院：信息與電子學(xué)院同組搭檔：專業(yè)：信息工程成績：一1.2.

2、3.二1.2.實驗?zāi)康恼莆招盘柕谋磉_(dá)及其可視化方法掌握信號基本時域運算的實現(xiàn)方法。利用實驗原理利用利用分析常用信號，加深對信號時域特性的理解表達(dá)連續(xù)時間信號進(jìn)行連續(xù)時間信號的時域運算（相加，相乘，微分， 尺度變換）表達(dá)離散時間信號進(jìn)行離散時間信號的時域運算,移位利用利用反轉(zhuǎn)3.4.三1.相應(yīng)程序的代碼及實驗結(jié)果1.1t=0:0.01:10;x=(1-exp(-0.5*t).*heaviside(t); plot(t,x)1.2t=0:0.01:10;x=cos(pi*t).*(heaviside(t)-heaviside(t-2); plot(t,x)1.3t=-2:0.01:10;x=(ab

3、s(t)/2).*cos(pi*t).*(heaviside(t+2)-heaviside(t-2); plot(t,x)1.4t=0:0.01:10;x=exp(-t).*sin(2*pi*t).*(heaviside(t)-heaviside(t-3);plot(t,x)2.2.1n=0:10;x=heaviside(n-3); stem(n,x,'filled');xlabel('n');title('x');2.2n=-10:10;x=(-0.5).(n).*heaviside(n); stem(n,x,'filled')

4、;xlabel('n');title('x');2.3n=-10:10;x=n.*(heaviside(n)-heaviside(n-5); stem(n,x,'filled');xlabel('n');title('x');2.4n=-10:10;x=sin(n.*pi/2).*heaviside(n); stem(n,x,'filled');xlabel('n');title('x');3.t=0:0.01:6;x=square(3*t); plot(t,x)4.

5、4.1t=0:0.01:4;x1=-t+4;x2=sin(2*pi*t); x3=x1+x2; plot(x3)4.2t=0:0.01:4;x1=-t+4;x2=sin(2*pi*t); x3=(x1).*(x2);plot(x3)4.3syms t;x1=(-t+4)*(heaviside(t)-heaviside(t-4); x2=subs(x1,t,-t);x3=(x1)+(x2);ezplot(x3,-4,4)4.4syms t;x1=(-t+4)*(heaviside(t)-heaviside(t-4); x2=sin(2*pi*t);x3=(x1)+(x2);x31=subs(x3

6、,t,t-1); x6=x2*x31; ezplot(x6,-4,10)5.5.1n=-3:4;x=0 1 2 3 3 3 3 0;stem(n,x,'filled');xlabel('n');title('x');5.2n=-3:4;x=0 1 2 3 3 3 3 0;stem(-n,x,'filled'); xlabel('n');title('x');5. 3n=-3:4;x=0 1 2 3 3 3 3 0;stem(n-2,x,'filled'); xlabel('n

7、');title('x');5.4n=-3:4;x=0 1 2 3 3 3 3 0;stem(n+2,x,'filled');xlabel('n');title('x');6.6.1t=0:0.1:20;x=1+cos(pi/4*t-pi/3)+2*cos(pi/2*t-pi/4)+cos(2*pi*t); plot(t,x)周 期 函 數(shù)6.2t=0:0.01:20;x=sin(t)+2*sin(pi*t); plot(t,x)非 周 期 函 數(shù)6.3t=0:10;x=2+3*sin(2*t*pi/3-pi/8); st

8、em(t,x,'filled');xlabel('t');title('x');周 期 函 數(shù)6.4t=0:10;x=cos(t*pi/6)+sin(t*pi/3)+cos(t*pi/2); stem(t,x,'filled');xlabel('t');title('x');周 期 函 數(shù)四試驗中所遇到的問題及解決方法實驗中首次運用進(jìn)行實際編程，遇到了不熟練，操作不明確等問1.題，經(jīng)過的努力摸索和查閱書籍，逐漸克服了此類問題。在計算過程中，發(fā)現(xiàn) U（n）函數(shù)在零點的取值為 0.5，與實際情況不符，2

9、.通過與交流得知，在中，零點的值等于前值和后值之和的一半，雖然與書本不同，但是不影響實驗結(jié)果。在進(jìn)行反轉(zhuǎn)，移位等操作時，出現(xiàn)定義域不明確的問題，通過仔細(xì)閱讀書3.籍改正了心得體會的錯誤。五1.2.是一款非常強(qiáng)大的軟件的時候，由于其所保有的函數(shù)庫非常大，基本涵蓋了我們在使用所有的應(yīng)用，所以要勤于查閱，多使用 HELP 等幫助性功能，有助于我們學(xué)習(xí)。信號與系統(tǒng)實驗報告實驗名稱：LTI 系統(tǒng)的時域分析課程名稱：信號與系統(tǒng)實驗實驗時間：10 月 22 日任課教師：王群實驗地點：4-427實驗教師：何冰松, 范哲意實驗類型： 原理驗證 綜合設(shè)計 創(chuàng)新學(xué)生尹迪學(xué)號/班級：1120131225/051113

10、61組號：學(xué)院：信息與電子學(xué)院同組搭檔：專業(yè)：信息工程成績：一1.2.3.4.實驗?zāi)康恼莆绽脤ο到y(tǒng)進(jìn)行時域分析的方法掌握連續(xù)時間系統(tǒng)零狀態(tài)響應(yīng)，沖擊響應(yīng)和階躍響應(yīng)的求解方法掌握求解離散時間系統(tǒng)響應(yīng)，抽樣響應(yīng)的方法加深對卷積 和計算的方法實驗原理和卷積和的理解。掌握利用計算機(jī)進(jìn)行卷積和卷積二1.2.3.三利用利用實現(xiàn)連續(xù)時間系統(tǒng)的時域分析實現(xiàn)離散時間系統(tǒng)的時域分析卷積和與卷積相應(yīng)程序的代碼及實驗結(jié)果第一題：1.1 b=1;a=1 sqrt(2) 1;sys=tf(b,a);t=0:0.01:10;impulse(sys)Impulse Response0.60.50.40.30.20.10-0

11、.1012345678910Time (seconds)1.2b=1 0 0;a=1 sqrt(2) 1;sys=tf(b,a);t=0:0.01:10;impulse(sys)AmplitudeImpulse Response0.20-0.2-0.4-0.6-0.8-1-1.2-1.4-1.60123456789Time (seconds)1.3b=1 0;a=1 1 1;sys=tf(b,a);t=0:0.01:10;impulse(sys)Impulse Response10.80.60.40.20-0.2-0.40246Time (seconds)810121.4a=1 1 1;b=2

12、 0 1;sys=tf(b,a);t=0:0.01:10;impulse(sys)AmplitudeAmplitudeImpulse Response0.50-0.5-1-1.5-20246Time (seconds)81012第二題:2.1（階躍響應(yīng)） a=1 1 6;b=1;sys=tf(b,a); step(sys)Step Response0.350.30.250.20.150.10.05002468101214Time (seconds)(沖激響應(yīng)) a=1 1 6;b=1;sys=tf(b,a); impulse(sys)AmplitudeAmplitudeImpulse Resp

13、onse0.50.40.30.20.10-0.1-0.202468101214Time (seconds)2.2沖激響應(yīng)的波形在在負(fù)無窮到正無窮的求和不是無窮， 所以系統(tǒng)是穩(wěn)定的。2.3a=1 1 6;b=1;sys=tf(b,a);t=0:0.01:10;x=exp(-t);lsim(sys,x,t)Linear Simulation Results1.210.80.60.40.20-0.2012345678910Time (seconds)第三題3.1a=1 3 2;b=1;impz(b,a,0:10)AmplitudeAmplitudeImpulse Response2500200015

14、0010005000-500-1000-1500012345n (samples)678910不穩(wěn)定系統(tǒng)3.2a=1 -0.5 0.8;b=1 -3;impz(b,a,0:10)Impulse Response2 521 510 50-0 5-1-1 5-2-2 5012345n (samples)678910穩(wěn)定系統(tǒng)第四題抽樣響應(yīng)：a=1 1 0.25;b=1;impz(b,a,0:10)AmplitudeAmplitudeImpulse Response10 80 60.40 20-0 2-0.4-0 6-0 8-1012345n (samples)678910階躍響應(yīng)：n=0:10;a=

15、1 1 0.25;b=1;x=heaviside(n); y=filter(b,a,x);stem(n,y,'filled');xlabel('n');title('y');y0 50.450.40 350 30 250 20.150.10 050012345n678910第五題n=-3:4;x1=1,2,1,1; x2=1 1 1 1 1;x=conv(x1,x2); stem(n,x,'filled');xlabel('n');title('x')Amplitude0 50-3第六題n=0:40

16、;n1=0:80;h=sin(0.5*n);x=sin(0.2*n);y=conv(x,h); stem(n1,y,'filled');xlabel('n1');title('y');0-1-2-3-40第七題t1=-1:0.001:1;t2=-2:0.001:2;t3=-3:0.001:3; x1=2+t1-t1; x2=1+t2-t2;y=conv(x1,x2)./2000; plot(t3,y);2.521.510.50-3-2-10123四1.試驗中所遇到的問題及解決方法實驗中運用進(jìn)行實際編程，遇到了不熟練，操作不明確等問題，經(jīng)過的努力摸

17、索和查閱書籍，逐漸克服了此類問題。2.在進(jìn)行卷積運算的時候，會出現(xiàn)卷積結(jié)果定義域不明確的問題，需要自己注意。心得體會是一款非常強(qiáng)大的軟件五1.2.在使用的時候，由于其所保有的函數(shù)庫非常大，基本涵蓋了我們所有的應(yīng)用，所以要勤于查閱，多使用 HELP 等幫助性功能，有助于我們學(xué)習(xí)。3.在卷積部分，對于卷積結(jié)果定義域的確定是非常重要的。信號與系統(tǒng)實驗報告實驗名稱：信號的頻域分析課程名稱：信號與系統(tǒng)實驗實驗時間：10 月 23 日任課教師：王群實驗地點：4-427實驗教師：何冰松, 范哲意實驗類型： 原理驗證 綜合設(shè)計 創(chuàng)新學(xué)生尹迪學(xué)號/班級：1120131225/05111361組號：學(xué)院：信息與電

18、子學(xué)院同組搭檔：專業(yè)：信息工程成績：一實驗?zāi)康?. 深入理解信號頻譜的概念，掌握信號的頻域分析方法。2. 觀察典型周期信號和非周期信號的頻譜，掌握其頻譜特性。實驗原理與方法二1.2.3.4.利用利用利用利用分析連續(xù)周期信號的頻譜 分析連續(xù)非周期信號的頻譜分析離散周期時間信號的頻域情況 分析離散非周期時間信號的頻域情況三相應(yīng)程序代碼及實驗結(jié)果1.1.2t=-10:0.01:10;N=input('N='); x=zeros(size(t); x=0.5;for n=1:1:N;x=x+(2/(pi*n).*sin(n*pi/2).*cos(0.5*pi*n*t);end plot

19、(t,x);當(dāng)T=4;A=1;寬度為2時當(dāng)T=6;A=1;寬度為4，N=10的時候1.210.80.60.40.20-0.2-10-8-6-4-20246810當(dāng) T=6;A=1;寬度為 4，N=100 的時候1.210.80.60.40.20-0.2-10-8-6-4-20246810當(dāng) T=6;A=1;寬度為 4，N=1000 的時候1.210.80.60.40.20-0.2-10-8-6-4-20246810當(dāng) T=6;A=1;寬度為 4，N=1000 的時候1.210.80.60.40.20-0.2-10-8-6-4-202468101.3N = 10;T = input('T

20、=');tt = input('tt=');A = 1;w = 2*pi/T; n1 = -N:-1;c1 = (2*A*l/T)*sin(n1*w*l/2)/w/l./n1*2; c0 = E*l/T;n2 = 1:N;c2 = (2*A*l/T)*sin(n1*w*l/2)/w/l./n1*2; cn = c1 c0 c2;n = -N:N;subplot(211); stem(n,abs(cn),'filled'); xlabel('omega/omega_0'); title('Magnitude of ck');

21、 subplot(212); stem(n,angle(cn),'filled'); xlabel('omega/omega_0'); title('Phase of ck');N = 10，A = 1，T=1tt = 0.5Magnitude of ck0.80.60.40.20-10-8-6-4-20w/w0Phase of ck24681043210-10-8-6-4-20w/w0246810tt = 0.3Magnitude of ck0.80.60.40.20-10-8-6-4-20w/w0Phase of ck24681043210-

22、10-8-6-4-20w/w0246810tt = 0.7Magnitude of ck0.80.60.40.20-10-8-6-4-20w/w0Phase of ck24681043210-10-8-6-4-20w/w0246810思考題：1.什么是吉伯斯現(xiàn)象？產(chǎn)生吉伯斯現(xiàn)象的是什么？答：將具有不連續(xù)點的周期函數(shù)（如矩形脈沖）進(jìn)行傅立葉級數(shù)展開后，選取有限項進(jìn)行。當(dāng)選取的項數(shù)越多，在所的波形中出現(xiàn)的峰起越靠近原信號的不連續(xù)點。當(dāng)選取的項數(shù)很大時，該峰起值趨于一個常數(shù)，大約等于總跳變值的 9%。這種現(xiàn)象稱為吉伯斯現(xiàn)象。是在不連續(xù)點附近所有的正弦信號均具有相同的變化趨勢，該趨勢在有限項內(nèi)無法被消

23、除。2. 以周期矩形脈沖信號為例，說明周期信號的頻譜有什么特點？ 答：周期信號的頻譜是具有周期性的一系列的脈沖信號。3. 周期矩形脈沖信號的有效頻帶寬度與信號的時域?qū)挾戎g有什么關(guān)系？答：時域?qū)挾仍酱?，有效頻帶寬度越小。4. 隨著矩形脈沖信號參數(shù)的變化，其頻譜結(jié)構(gòu)如何變化？答：頻譜包絡(luò)形狀不變，過零點不變，普賢間隔隨著 T 變大而縮小。2.2.2syms t A tt; tt=input('tt=');A=input('A=');x=A*( heaviside(t+ (tt/2) )-heaviside(t- (tt/2) ); X=fourier(x)ezpl

24、ot(abs(X),-20,20);A=1;tt=1;A=1;tt=2;A=1;tt=5;思考題：1).比較矩形脈沖信號和周期矩形脈沖信號的頻譜，兩者之間有何異同？答：矩形脈沖信號頻譜為連續(xù)函數(shù)，周期矩形脈沖信號的頻譜為一系列的脈沖。2).根據(jù)矩形脈沖寬度變化時頻譜的變化規(guī)律，說明信號的有效頻帶寬度與其時域?qū)挾戎g有什么關(guān)系？答：信號有效頻帶寬度越大，時域?qū)挾仍叫　?.3.3Double N; double N1; N=input('N=');N1=input('N1=');n=-N+1:N-1;x(n>=-N1 & n<= N1)=1; x

25、(n>N1 & n<(N-N1)=0; x(n>-N+N1 & n<-N1)=0; x(n>=N-N1 & n<=N+N1)=1;subplot(211); stem(n,x,'filled');xlabel('n');title('x(n)'); subplot(212); X=fft(x); stem(n,X,'filled');xlabel('k');title('X(k)');N1=5,N=20 的運行結(jié)果x(n)10.50-20-

26、15-10-50n X(k)510152020100-10-20-15-10-50k5101520x(n)10.50-20-15-10-50n X(k)5101520151050-5-20-15-10-50k5101520思考題：1).以周期方波序列為例，說明周期序列與連續(xù)周期信號的頻譜有何異同。答：周期序列的頻譜向外越來越大，連續(xù)周期信號頻譜則是中間向兩邊越來越小。2).隨著周期方波序列占空比的變化，其頻譜如何隨之變化？ 答：隨著占空比越來越大，頻譜密度也越來越大。4.4.1double N1; N1=input('N1='); w=-pi:0.01*pi:pi; n=-N1

27、:N1;x(n>=(-N1) & n<=N1)=1;X=x*exp(-j*(n'*w) ); subplot(211); stem(n,x,'filled');xlabel('n');title('x(n)'); subplot(212); plot(w/pi,abs(X);N=10x(n)10.50-10-8-6-4-20n2468103020100-1-0.8-0.6-0.4-0.200 20.40.60.81N=5x(n)10.50-5-4-3-2-10n12345151050-1-0.8-0.6-0.4-0.2

28、00 20.40.60.81思考題：1).隨著矩形脈沖序列寬度的變化，其頻譜如何隨之變化？其寬度與頻譜的有效寬度有何關(guān)系？答：寬度越大，頻譜變化就越密集。矩形脈沖信號序列寬度越大，頻譜有效寬度越大。試驗中所遇到的問題及解決方法1. 此次實驗難度與前兩次相比較有明顯的提升，花費的時間也明顯增加。很多程序都是在與同學(xué)討論之后才寫出來的。本次實驗中我遇到的問題就是容易誤把頻譜也做了累加，得到了錯誤的波形和結(jié)果， 在經(jīng)過和上的范例對比之后，我找出了錯誤，得到了正確的波形和結(jié)果。心得體會四五1.主要收獲就是從實驗的角度了解了吉伯斯現(xiàn)象，同時也掌握了對連續(xù)時間信號以及離散時間信號的頻域分析方法，很使用多時

29、域上難以分析的信號在頻域上是非常直觀和間接的。2. 信號的頻域分析方法對于信號與系統(tǒng)這門課程的學(xué)習(xí)乃至是今后要從事的有關(guān)信號處理的工作而言都是關(guān)鍵、有用的。利用將十分抽象的表達(dá)式轉(zhuǎn)化為了直觀的圖像，并且從側(cè)面印證了很多課上學(xué)習(xí)到的東西，加深了對課程的理解?；蛟S這個實驗的意義就是印證課上知識加深理解，并掌握分析方法，為以后工作與學(xué)習(xí)打下堅實的基礎(chǔ)。信號與系統(tǒng)實驗報告實驗名稱：LTI 系統(tǒng)的頻域分析課程名稱：信號與系統(tǒng)實驗實驗時間：10 月 28 日任課教師：王群實驗地點：4-427實驗教師：何冰松, 范哲意實驗類型： 原理驗證 綜合設(shè)計 創(chuàng)新學(xué)生尹迪學(xué)號/班級：1120131225/051113

30、61組號：學(xué)院：信息與電子學(xué)院同組搭檔：專業(yè)：信息工程成績：一、1.2.二、1.2.三、1.實驗?zāi)康募由顚?LTI 系統(tǒng)頻率響應(yīng)概念的掌握和理解學(xué)習(xí)和掌握 LTI 系統(tǒng)頻率特性的分析方法實驗原理與方法利用利用計算連續(xù)時間系統(tǒng)的頻率響應(yīng)計算離散時間系統(tǒng)的頻率響應(yīng)相應(yīng)程序的代碼及實驗結(jié)果1.2b=0.04 0 0;a=0.04 0.4 2;freqs(b,a); H,w=freqs(b,a); subplot(211); plot(w,abs(H);axis(0,30,0,1);set(gca,'xtick',0:2:30);set(gca,'ytick',0 0.

31、4 0.707 1); xlabel('omega(rad/s)'); ylabel('Magnitude'); title('|H(jomega)|');grid on; subplot(212); plot(w,angle(H);xlabel('omega(rad/s)'); ylabel('Phase');title('phi(omega)'); grid on;|H(jw)|10.7070.400 2468 10 1214 16 18 20w(rad/s) phi(w)22 24 26 28

32、 30432100102030405060w(rad/s)7080901002.2.1R=input('R=');MagnitudePhaseC=input('C=');b=1;a=R*C 1;H,w=freqs(b,a); subplot(211); plot(w,abs(H); xlabel('omega(rad/s)'); ylabel('Magnitude'); title('|H(jomega)|'); grid onsubplot(212); plot(w,angle(H); xlabel('o

33、mega(rad/s)'); ylabel('Phase');title('phi(omega)'); grid on；RC=1 時|H(jw)|10.500123456w(rad/s) phi(w)789100-0.5-1-1.50123456w(rad/s)78910RC=0.1 時|H(jw)|10.500102030405060w(rad/s) phi(w)7080901000-0.5-1-1.50102030405060w(rad/s)708090100PhasePhaseMagnitudeMagnitudeRC=0.01 時|H(jw)|10

34、.500100200300400500600w(rad/s) phi(w)70080090010000-0.5-1-1.50100200300400500600w(rad/s)7008009001000RC=0.001 時|H(jw)|10.500100020003000 4000 5000 6000w(rad/s) phi(w)700080009000100000-0.5-1-1.50100020003000 4000 5000 6000w(rad/s)700080009000100002.2運用 RC=0.005 時的濾波系統(tǒng)t=0:0.001:0.2; b=1;a=0.001 1;sys

35、=tf(b,a); x=cos(100*t)+cos(3000*t); lsim(sys,x,t);PhaseMagnitudePhaseMagnitudeLinear Simulation Results21.510.50-0.5-1-1.5-200.020.040.060.080.10.120.140.160.180.2Time (seconds)3.3.1xn+ xn-1+ xn-2+ xn-3+ xn-4+ xn-5+ xn-6+xn-7+ xn-8=yn系統(tǒng)函數(shù)：H(z)=1+z-1+z-2+z-3+z-4+z-5+z-6+z-7+z-83.2b=1 1 1 1 1 1 1 1 1;

36、 a=1;impz(b,a);3.3.1b=1 1 1 1 1 1 1 1 1; a=1;H,w=freqz(b,a); plot(w/pi,abs(H)幅頻特性曲線Amplitudeb=1 1 1 1 1 1 1 1 1;a=1;H,w=freqz(b,a); plot(w/pi,angle(H);相頻特性b=1 1 1 1 1 1 1 1 1;a=1;H,w=freqz(b,a); zplane(b,a);零級點圖4.n=-10:10;x=cos(0.3*pi*n)+0.5*cos(0.8*pi*n); stem(n,x,'filled');xlabel('n

37、9;);title('x(n)');x(n)圖像n=-10:10;y=2*cos(0.3*pi*n); stem(n,y,'filled');xlabel('n');title('y(n)');stem(n,y,'filled');y(n)圖像四、1.2.五、1.2.實驗中所遇到的問題及解決方法在實驗 3 中，錯誤的以為輸入只有 8 項，導(dǎo)致圖像相似但是確。在實驗 4 中，對題目理解有一些問題，和心得體會交流后解決問題。第一次發(fā)現(xiàn)原來還可以畫零極點圖，對其有了更深的認(rèn)識。在使用的時候，由于其所保有的函數(shù)庫非常大，基

38、本涵蓋了我們所有的應(yīng)用，所以要勤于查閱，多使用 HELP 等幫助性功能，有助于我們學(xué)習(xí)。信號與系統(tǒng)實驗報告實驗名稱：連續(xù)時間系統(tǒng)的復(fù)頻域分析課程名稱：信號與系統(tǒng)實驗實驗時間：10 月 30 日任課教師：王群實驗地點：4-427實驗教師：何冰松, 范哲意實驗類型： 原理驗證 綜合設(shè)計 創(chuàng)新學(xué)生尹迪學(xué)號/班級：1120131225/05111361組號：學(xué)院：信息與電子學(xué)院同組搭檔：專業(yè)：信息工程成績：一、實驗?zāi)康?.掌握拉斯變換及其反變換的定義，并掌握實現(xiàn)方法2. 學(xué)習(xí)和掌握連續(xù)時間系統(tǒng)系統(tǒng)函數(shù)的定義及復(fù)頻域分析方法3. 掌握系統(tǒng)零極點的定義，加深理解系統(tǒng)零極點分布于系統(tǒng)特性的關(guān)系二、實驗原理及

39、方法1. 利用2. 利用3. 利用實現(xiàn)拉斯變換實現(xiàn)連續(xù)時間系統(tǒng)的系統(tǒng)函數(shù)表示實現(xiàn)連續(xù)時間系統(tǒng)的零極點分析三、相應(yīng)程序代碼及實驗結(jié)果1.h=sym('heaviside(t)-heaviside(t-2)'); h1=laplace(h); x=sym('heaviside(t)');h2=laplace(x); Y=h2.*h1; y=ilaplace(Y)ezplot(y,0,5);2.1.b=1;a=1 2 2 1;sys=tf(b,a); pzmap(sys); grid on;b=1;a=1 2 2 1;H,w=freqs(b,a); subplot(2

40、11); plot(w,abs(H);xlabel('omega(pi)');ylabel('幅頻特性'); grid on; subplot(212); plot(w,angle(H);xlabel('omega(pi)');ylabel('相頻特性'); grid on;沖激響應(yīng)：b=1;a=1 2 2 1;r,p,k=residue(b,a) r =1.0000 + 0.0000i-0.5000 - 0.2887i-0.5000 + 0.2887ip =-1.0000 + 0.0000i-0.5000 + 0.8660i-0

41、.5000 - 0.8660ik =X(t)=(exp(-1)+(-0.5-0.2887i)*(exp(-0.5+0.866i)+(-0.5+0.2887i)*(exp(-0.5 - 0.8660i);2.b=1 0 1;a=1 2 -3 3 3 2;sys=tf(b,a); pzmap(sys); grid on;b=1 0 1;a=1 2 -3 3 3 2;H,w=freqs(b,a); subplot(211); plot(w,abs(H);xlabel('omega(pi)');ylabel('幅頻特性'); grid on ; subplot(212)

42、; plot(w,angle(H);xlabel('omega(pi)');ylabel('相頻特性'); grid on ;沖激響應(yīng)b=1 0 1;a=1 2 -3 3 3 2;r,p,k=residue(b,a)r =0.0769 + 0.0000i-0.0300 - 0.0881i-0.0300 + 0.0881i-0.0085 - 0.1436i-0.0085 + 0.1436ip =-3.1704 + 0.0000i0.9669 + 0.9540i0.9669 - 0.9540i-0.3817 + 0.4430i-0.3817 - 0.4430ik =

43、X(t)=(0.0769)*exp(-3.1704)+(-0.03-0.0881i)*(exp(0.9669+0.9540i)+(-0.03+0.0881i)*(exp(0.9669-0.9540i)+(-0.0085 - 0.1436i)*(exp(-0.3817+0.4430i)+(-0.0085+0.1436i)*(exp(-0.3817-0.4430i);3.1.t=0:0.01:10;b=1;a=1 0;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=heaviside(t);plot(t,x);2.t=0:0.01:10;

44、b=1;a=1 2;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=exp(-2*t).*heaviside(t);plot(t,x);3.t=0:0.01:10;b=1;a=1 -2;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=exp(2*t).*heaviside(t);plot(t,x);4.t=0:0.01:10;b=1;a=1 0 4;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=(sin(2*t)/

45、2).*heaviside(t);plot(t,x);5.t=0:0.01:10;b=1;a=1 2 17;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=0.25*exp(-t).*sin(4*t).*heaviside(t);plot(t,x);6.t=0:0.01:10;b=1;a=1 -2 17;sys=tf(b,a); subplot(211); pzmap(sys); subplot(212);x=0.25*exp(t).*sin(4*t).*heaviside(t);plot(t,x);4.1.t=0:0.01:10;b

46、=1;a=1 2 17;sys=tf(b,a); subplot(211); pzmap(sys);subplot(212);x=0.25*exp(-t).*sin(4*t).*heaviside(t);plot(t,x);2.t=0:0.01:10;b=1 8;a=1 2 17;sys=tf(b,a); subplot(211); pzmap(sys); grid on ; subplot(212);h2=exp(-t).*(cos(4*t)+(7*sin(4*t)/4).*heaviside(t); plot(t,h2);grid on;3.t=0:0.01:10;b=1 -8;a=1 2 17;sys=tf(b,a); subplot(211);pzmap(sys);subplot(212);h3=(exp(-t).*(cos(4*t