長(zhǎng)沙理工大二第二學(xué)期信號(hào)與系統(tǒng)實(shí)驗(yàn)總結(jié)

實(shí)驗(yàn)一(1)程序以下:f=sym('((2-3*exp(-t))*Heaviside(t))');ezplot(f,[-3,10]);gridon;xlabel('時(shí)間');ylabel('幅度');title('f(t)=(2-3e^-^t)ε(t)')結(jié)果以下：(3)程序以下:f=sym('(sin(pi*t)*Heaviside(t))');ezplot(f,[-3,10]);gridon;xlabel('時(shí)間');ylabel('幅度');title('f(t)=sin(πt)ε(t)')結(jié)果以下：(7)程序以下:k=-2:5;x=2.^k;y=stepfun(k,0);f=x.*y;stem(k,f)gridon;xlabel('時(shí)間');ylabel('幅度')title('f(t)=2^kε(t)')結(jié)果以下：程序以下:k=-2:17;x=sin(k*pi/4);y=stepfun(k,0);f=x.*y;stem(k,f)gridon;xlabel('時(shí)間');ylabel('幅度')title('f(t)=2^kε(t)')結(jié)果以下：(3)程序以下

:t=-4::4;

？u1=stepfun(t,-2);u2=stepfun(t,2);g1=2*(u1-u2);u3=stepfun(t,0);g2=t.*u3-t.*u2;g3=g1+g2;t2=2-t;plot(t2,g3)

？axis([-4,8,,6])gridon;

？xlabel('時(shí)間');ylabel('

幅度')title('f(2-t)')結(jié)果以下：(6)程序以下:t=-4::4;

？u1=stepfun(t,-2);u2=stepfun(t,2);g1=2*(u1-u2);u3=stepfun(t,0);g2=t.*u3-t.*u2;g3=g1+g2;t2=*t-2;plot(t2,g3)

？axis([-8,4,,6])gridon;

？xlabel('

時(shí)間');ylabel('

幅度')title('f')結(jié)果以下：實(shí)驗(yàn)二第一題程序以下：subplot(3,2,1)t=-20::20;y=square*t,25);%矩形脈沖y(y<0)=0;plot(t-1,y)gridon;xlabel('時(shí)間t');ylabel('幅度');title('矩形脈沖信號(hào)τ/T=1/4');subplot(3,2,2)t=-20::20;y=square*t,25);yy=fft(y);N=size(y);N=N(2);z=abs(fft(y,64*32))/N*2;stem(z);gridon;xlabel('角頻次w');ylabel('幅度');title('矩形脈沖信號(hào)τ/T=1/4的頻譜圖');subplot(3,2,3);t=-50::50;y=square*t,;y(y<0)=0;plot(t-1,y)gridon;axis([-50500]);xlabel('時(shí)間t');ylabel('幅度');title('矩形脈沖信號(hào)τ/T=1/16');subplot(3,2,4)t=-50::50;y=square*t,;yy=fft(y);N=size(y);N=N(2);z=abs(fft(y,128*64))/N*2;stem(z);gridon;axis([-50001]);xlabel('角頻次w');ylabel('幅度');title('矩形脈沖信號(hào)τ/T=1/16的頻譜圖');subplot(3,2,5);t=-200::200;y=square*t,;y(y<0)=0;plot(t-1,y)gridon;axis([-2002000]);xlabel('時(shí)間t');ylabel('幅度');title('矩形脈沖信號(hào)τ/T=1/64');subplot(3,2,6)t=-200::200;y=square*t,;yy=fft(y);N=size(y);N=N(2);z=abs(fft(y,128*64))/N*2;stem(z);gridon;axis([-200-19801]);xlabel('角頻次w');ylabel('幅度');title('矩形脈沖信號(hào)τ/T=1/64的頻譜圖');結(jié)果以下：實(shí)驗(yàn)2第二題程序以下：subplot(2,1,1)t=-2::2;g=rectpuls(t,2);%周期方波y=*(1+cos(pi*t)).*g;plot(t,y);axis([])gridon;xlabel('t');ylabel('幅度');title('升余弦脈沖圖形');subplot(2,1,2)t=-4::4;q=rectpuls(t,2*pi);f=((1+cos(t*pi))/2).*q;k=-500:500;w=pi*k/(500*;Fw=*f*exp(-j*t'*w);stem(w,abs(Fw))axis([-10104])gridon;xlabel('w');ylabel('幅度');title('升余弦脈沖的頻譜圖形');結(jié)果以下：實(shí)驗(yàn)三.習(xí)題P80，(1)程序以下：functionshiyan1a=[1,4,3];b=[1];sys=tf(b,a);subplot(221);t=-1::4;t0=0;ut=stepfun(t,t0);plot(t,ut);axis([-1,4,,]);gridon;xlabel('時(shí)間');ylabel('f(t)');title('激勵(lì)的時(shí)域波形')subplot(222);t=0::8;y1=impulse(sys,t);plot(t,y1);gridon;xlabel('時(shí)間');ylabel('h(t)');title('沖激響應(yīng)的時(shí)域波形')subplot(223);t=0::8;y2=step(sys,t);plot(t,y2);gridon;xlabel('時(shí)間');ylabel('g(t)');title('階躍響應(yīng)的時(shí)域波形')圖形以下：.習(xí)題P111，（4）程序以下：functionshiyansan2a=[1,3,2];b=[1];subplot(221);t=-4::8;t0=0;ut=stepfun(t,t0);plot(t,ut);axis([-4,8,,]);gridon;xlabel('時(shí)間');ylabel('f(k');title('激勵(lì)的時(shí)域波形')subplot(222);k=-2:18;y1=impz(b,a,k);stem(k,y1);gridon;xlabel('時(shí)間');ylabel('h(k)');title('單位序列響應(yīng)的時(shí)域波形')subplot(223);B=18;n=[0:B-1];y2=stepz(b,a,B);stem(n,y2);gridon;xlabel('時(shí)間');ylabel('g(k)');title('階躍響應(yīng)的時(shí)域波形')圖形以下：實(shí)驗(yàn)4程序以下：functionshiyan4w=-28::28;a=[1,3,2];b=[0,1,0];H=freqs(b,a,w);subplot(211);plot(w,abs(H));title('幅頻響應(yīng)曲線');xlabel('w');ylabel('幅頻');subplot(212);plot(w,angle(H));title('相頻響應(yīng)曲線');xlabel('w');ylabel('相頻');結(jié)果以下：實(shí)驗(yàn)五實(shí)驗(yàn)五（1）（2）程序以下：symsabts%syms函數(shù)用于創(chuàng)立符號(hào)對(duì)象,定義變量的laplace(sin(2*t-pi/4))？laplace(t^2*exp(-2*t))？laplace(t*exp(-a*t)*cos(b*t))？ilaplace((s^2+4*s+5)/(s^2+3*s+2))ilaplace(((s+1)*(s+4))/(s*(s+2)*(s+3)))？ilaplace((s^2+4*s)/((s+1)*(s^2-4)))運(yùn)轉(zhuǎn)結(jié)果以下：ans=-1/2*2^(1/2)*(s-2)/(s^2+4)ans=2/(s+2)^3ans=(s^2+2*s*a+a^2-b^2)/(s^2+2*s*a+a^2+b^2)^2ans=dirac(t)-exp(-2*t)+2*exp(-t)ans=2/3+exp(-2*t)-2/3*exp(-3*t)ans=2*sinh(2*t)+exp(-t)實(shí)驗(yàn)五（3）程序以下：functionshiyan5a=[1,4];b=[1,3,2,0];sys=tf(b,a);figure(1);pzmap(sys);title('零極點(diǎn)散布圖');c=[1,3,2,0];d=[1,4];sys=tf(d,c);subplot(2,2,1);pzmap(sys);F=sym('(s+4)/(s*(s^2+3*s+2))');f=ilaplace(F)b=[1,4];a=[1,3,2];w=linspace(-5,5,