信号与系统课后matlab作业.doc

M2-1 (1) t=-2:0.001:4; T=2; xt=rectpuls(t-1,T); plot(t,xt) axis([-2,4,-0.5,1.5]) 图象为： (2) t=sym(t); y=Heaviside(t); ezplot(y,[-1,1]); grid on axis([-1 1 -0.1 1.1]) 图象为： (3) A=10;a=-1;B=5;b=-2; t=0:0.001:10; xt=A*exp(a*t)-B*exp(b*t); plot(t,xt) 图象为： (4) t=sym(t); y=t*Heaviside(t); ezplot(y,[-1,3]); grid on axis([-1 3 -0.1 3.1]) 图象为： (5) A=2;w0=10*pi;phi=pi/6; t=0:0.001:0.5; xt=abs(A*sin(w0*t+phi)); plot(t,xt) 图象为： (6) A=1;w0=1;B=1;w1=2*pi; t=0:0.001:20; xt=A*cos(w0*t)+B*sin(w1*t); plot(t,xt) 图象为： (7) A=4;a=-0.5;w0=2*pi; t=0:0.001:10; xt=A*exp(a*t).*cos(w0*t); plot(t,xt) 图象为： (8) w0=30; t=-15:0.001:15; xt=cos(w0*t).*sinc(t/pi); plot(t,xt) axis([-15,15,-1.1,1.1]) 图象为： M2-3 function yt=x2_3(t) yt=(t).*(t=0t=2)+2*(t=2t=3)-1*(t=3t=5); （2） function yt=x2_3(t) yt=(t).*(t=0t=2)+2*(t=2t=3)-1*(t=3t=5); t=0:0.001:6; subplot(3,1,1) plot(t,x2_3(t)) title(x(t)) axis([0,6,-2,3]) subplot(3,1,2) plot(t,x2_3(0.5*t)) title(x(0.5t)) axis([0,11,-2,3]) subplot(3,1,3) plot(t,x2_3(2-0.5*t)) title(x(2-0.5t)) axis([-6,5,-2,3]) 图像为： M2-9 （1） k=-4:7; xk=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; stem(k,xk,file) （2） k=-12:21; x=[-3,0,0,-2,0,0,3,0,0,1,0,0,-2,0,0,-3,0,0,-4,0,0,2,0,0,-1,0,0,4,0,0,1,0,0,-1]; subplot(2,1,1) stem(k,x,file) title(3倍内插) t=-1:2; y=[-2,-2,2,1]; subplot(2,1,2) stem(t,y,file) title(3倍抽取) axis([-3,4,-4,4]) （3） k=-4:7; x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; subplot(2,1,1) stem(k+2,x,file) title(x[k+2]) subplot(2,1,2) stem(k-4,x,file) title(x[k-4]) （4） k=-4:7; x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; stem(-fliplr(k),fliplr(x),file) title(x[-k]) M3-1 （1） ts=0;te=5;dt=0.01; sys=tf([2 1],[1 3 2]); t=ts:dt:te; x=exp(-3*t); y=lsim(sys,x,t); plot(t,y) xlabel(Time(sec)) ylabel(y(t)) （2） ts=0;te=5;dt=1; sys=tf([2 1],[1 3 2]); t=ts:dt:te; x=exp(-3*t); y=lsim(sys,x,t) y = 0 0.6649 -0.0239 -0.0630 -0.0314 -0.0127 从（1）（2）对比我们当抽样间隔越小时数值精度越高。 M3-4 x=[0.85,0.53,0.21,0.67,0.84,0.12]; kx=-2:3; h=[0.68,0.37,0.83,0.52,0.71]; kh=-1:3; y=conv(x,h); k=kx(1)+kh(1):kx(end)+kh(end); st