matlab在信號與系統(tǒng)系統(tǒng)中

Matlab在信號與系統(tǒng)中的應(yīng)用1、連續(xù)時間信號的波形內(nèi)容：1、繪制以下信號的波形：a.3e4tplot(t,y)606050403020109x109x105876543210b.5cos(2t)3 4t=-pi:0.01:pi;plot(t,y)543210-1-2-3-4-5-4 -3 -2 -1 0 1 2 3 4Sa(t)y=sin(t)./t;plot(t,y)10.80.60.40.20-0.2-0.4-10 -8 -6 -4 -2 0 2 4 6 8 10t=-5:0.1:5;x=(t>=0);plot(t,x)10.90.80.70.60.50.40.30.20.10-5 -4 -3 -2 -1 0 1 2 3 4 5y=(t==0);plot(t,y)10.90.80.70.60.50.40.30.20.10-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 12、連續(xù)時間信號的時域分析內(nèi)容：1、繪制以下信號的波形：ysincost 3 3t=-5:0.01:5;y=sin(5*pi*t)+cos(pi/3*t);plot(t,y)21.510.50-0.5-1-1.5-2-5 -4 -3 -2 -1 0 1 2 3 4 54ye2tsin2t4 t=-5:0.01:5;x104x10421.510.50-0.5-1-5 -4 -3 -2 -1 0 1 2 3 4 52)

u(t)x(t)y(t)的圖形，t-35，0.01。t=-3:0.01:5;z=(t>0);x=exp(-0.5*t).*z;y=x.*(1.5*t+3);plot(t,x,t,y)32.521.510.50-3 -2 -1 0 1 2 3 4 53

f(t)eatu(t)(a=1f2t)sintut)f(t)f1(tf2(t)1的時間范圍取為0-10，步長值取為0.1。繪制三個信號的波形。1function[w,tw]=conv1(u,tu,v,tv)T=tu(2)-tu(1);w=T*conv(u,v);tw=tu(1)+tv(1)+T*[0:length(u)+length(v)-2]';t=[0.1:0.1:10]';a=1;u=heaviside(t);f1=exp((-a).*t).*u;f2=sin(t).*u;[f,t1]=conv1(f1,t,f2,t);plot(t,f1,'r',t,f2,'g',t1,f)10.80.60.40.20-0.2-0.4-0.6-0.8-10 2 4 6 8 10 12 14 16 18 203、連續(xù)時間系統(tǒng)的時域分析內(nèi)容：1、 已知系統(tǒng)的微分方程為r''(t)4r'(t)3r(t)e(t)當(dāng)系統(tǒng)的輸入信號為e(t)4cos(2t)，求系統(tǒng)的零狀態(tài)響應(yīng)，并畫出圖形。a=[1,4,3];b=[1];sys=tf(b,a);lsim(sys,e,t)LinearSimulationResultsInput:In(1)LinearSimulationResultsInput:In(1)Time(sec):0.755Amplitude:0.111321p0pmA-1-2-3-40 1 2 3 4 5Time(sec)

6 7 8 9 10r''(t4r'(t3r(te'(t2e(t，求系統(tǒng)的沖激響應(yīng)，同時求出理論解，繪圖比較。a=[1,4,3];b=[1,2];sys=tf(b,a);t=[0:0.1:5]';[h1,t1]=impulse(sys,t);plot(t1,h1)10.90.80.70.60.50.40.30.20.100 0.5 1 1.5 2 2.5 3 3.5 4 4.5 53.求下列微分方程描述的系統(tǒng)的階躍響應(yīng)，同時求出理論解，繪圖比較。r''(t)'(t)2r(t)e'(t)3e(t)a=[1,3,2];b=[1,3];sys=tf(b,a);t=[0:0.1:5]';step(sys)StepResponse1.5StepResponse1ppmA0.500 1 2 3Time

4 5 64、連續(xù)時間信號的頻譜分析內(nèi)容：1、求以下信號的傅里葉變換及繪制頻譜圖a.f(t)cost 2 2symstF=cos((pi*t)/2);f=fourier(F);w=[0:0.1:5];f=pi*(dirac(-pi/2-w)+dirac(pi/2-w));x=abs(f);y=angle(f);plot(w,x,'*',w,y,'p')10.80.60.40.20-0.2-0.4-0.6-0.8-10 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5b.f(t)2e0.4tu(t)symstu=heaviside(t);F=2*exp(-2/5*t).*u;f=fourier(F);w=[0:0.1:3];f=2./(i*w+2/5);x=abs(f);y=angle(f);plot(w,x,w,y)543210-1-200.511.522.53c.f(t)u(t2)u(t2)symstf=fourier(heaviside(t+2)-heaviside(t-2));f=(i./w-pi.*dirac(-w))./exp(2*i*w)-exp(2*i.*w).*(i./w-pi*dirac(-w));abs(f);w=0:0.1:5;angle(f);plot(w,abs(f),w,angle(f))43210-1-2-3-40 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5d.f(t)etsymstF=exp(-abs(t));f=fourier(F);f=2./(w.^2+1);abs(f);w=0:0.1:5;angle(f);plot(w,abs(f),'*',w,angle(f),'r')21.81.61.41.210.80.60.40.200 0.5 1 1.5 2 2.5 3 3.54 4.555、連續(xù)時間系統(tǒng)的s域分析內(nèi)容：1．a(chǎn). 3s(s4)(s2)symssS=(3*s)/((s+4)*(s+2));f=ilaplace(S)f=6/exp(4*t)-3/exp(2*t)b.symss

4s5s25s6S=(4*s+5)/(s^2+5*s+6);f=ilaplace(S)f=7/exp(3*t)-3/exp(2*t)2.LTIa.r'(t3r(t)3e(t)a=[1,3];b=[3];freqs(b,a)uM

e(hP

101010100-20-40-60-8010

100Frequency(rad/s)100Frequency(rad/s)

101101b.3r''(t)4r'(t)r(t)e''(t)5e(t)a=[3,4,1];b=[1,0,5];freqs(b,a)uM

102100101010-2 10-1 100 101Frequency(rad/s)e(eaP

500-50-100-150

10

Frequency

100

101分別畫出它們的系統(tǒng)頻率響應(yīng)的幅值和相位特性曲線。3、有一模擬濾波器，其傳遞函數(shù)為：H(s)

0.2s20.3s1s20.4s1用freqs函數(shù)畫出它的幅頻特性和相位特性曲線。b=[0.2,0.3,1];a=[1,0.4,1];freqs(b,a)uM10-1 Frequency(rad/s)e(hP

0-20-40-60-80-100-12010-1

Frequency(rad/s)4.已知系統(tǒng)函數(shù)H（s）為：s1H(s)

s27s10求該系統(tǒng)的零極點，并畫出零極點分布圖。b=[1,1];a=[1,7,10];pzmap(sys)MapMap0.80.6iAmI

0.40.20-0.2-0.4-0.6-0.8-1-5 -4.5 -4 -3.5 -3 -2.5RealAxis

-2 -1.5 -1 -0.5 06、離散時間信號的產(chǎn)生與運算內(nèi)容：n=[-10:10]范圍內(nèi)產(chǎn)生離散信號：2nf(n)0

3n3其余n=[0:30]范圍內(nèi)畫出下列信號：naf(n)cos2( )2 43、畫出如下信號的波形 a.ysinn cos 2 3 4b.y(2)nsinn4 4、兩個有限長序列分別如下圖所示，求兩個信號的卷積和。x(x(n)13221121 2nx(n)243211101 23n