信号实验2matlab第二次

1、MATLAB第二次实验项目程序图形实验四1（1）syms t wGt1=sym(Heaviside(t+1)-Heaviside(t-1);Gt2=sym(Heaviside(2*t+1)-Heaviside(2*t-1);Fw1=fourier(Gt1,t,w);Fw2=fourier(Gt2,t,w);FFw1=maple(convert,Fw1,piecewise);FFw2=maple(convert,Fw2,piecewise);FFP1=abs(FFw1);FFP2=abs(FFw2);subplot(2,1,1);ezplot(FFP1,-10*pi 10*pi);grid;ax

2、is(-10*pi 10*pi 0 2.2)subplot(2,1,2);ezplot(FFP2,-10*pi 10*pi);grid;axis(-10*pi 10*pi 0 2.2)实验四1（3）syms t mGt=sym(Heaviside(t);f=exp(-t)*Gt;Fw=fourier(f,t,w);FFw=maple(convert,Fw,piecewise);FFP=abs(FFw);ezplot(FFP,-10*pi 10*pi);grid;axis(-10*pi 10*pi 0 1.2)实验四2(2)syms t wFw=sym(i*w)2+5*i*w-8)/(i*w)2

3、+6*i*w+5);ft=ifourier(Fw,w,t)运算结果：ft =dirac(t)+(2*exp(-5*t)-3*exp(-t)*heaviside(t)实验五1（1）num=1 0 1;den=1 2 5;z,p=tf2zp(num,den), zplane(z,p)z = 0 + 1.0000i 0 - 1.0000ip = -1.0000 + 2.0000i -1.0000 - 2.0000i实验五2（1）num=0 2 0;den=1 20.5 1;w=logspace(-2,2);freqs(num,den,w),实验六1（1）实现单位冲激响应：a=1 4 4;b=1 3;

4、impulse(b,a,0:0.05:4);实现实现零状态响应的程序：a=1 4 4;b=1 3;p1=0.01;t1=0:p1:4;x1=exp(-t1);lsim(b,a,x1,t1);hold on;p2=0.5;t2=0:p2:4;x2=exp(-t2);lsim(b,a,x2,t2);hold off实验六（2）实现单位函数响应及阶跃响应的程序a=1 -0.7 0.1;b=7 -2 0;y1=impz(b,a,0:20);y1=y1y2=dstep(b,a,0:20);y2=y2subplot(2,1,1);impz(b,a,0:20);grid;subplot(2,1,2);dst

5、ep(b,a,0:20);grid;实现零状态响应及阶跃响应的程序a=1,-0.7,0.1;b=7,-2,0;k=0:20;x=1.k;y=filter(b,a,x)；y=ysubplot(2,1,1)stem(k,x)subplot(2,1,2)stem(k,y) y1 =y1 = Columns 1 through 10 7.0000 2.9000 1.3300 0.6410 0.3157 0.1569 0.0783 0.0391 0.0195 0.0098 Columns 11 through 20 0.0049 0.0024 0.0012 0.0006 0.0003 0.0002 0.

6、0001 0.0000 0.0000 0.0000 Column 21 0.0000y2 = Columns 1 through 10 7.0000 9.9000 11.2300 11.8710 12.1867 12.3436 12.4218 12.4609 12.4805 12.4902 Columns 11 through 20 12.4951 12.4976 12.4988 12.4994 12.4997 12.4998 12.4999 12.5000 12.5000 12.5000 Column 21 12.5000y = Columns 1 through 9y = Columns 1 through 9 7.0000 9.9000 11.2300 11.8710 12.1867 12.3436 12.4218 12.4609 12.4805 Columns 10 through 18 12.4902 12.4951