信号与系统课后matlab作业.

(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]) 图象为：

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])

图象为：

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)

图象为：

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])

图象为：

（1）function yt=x2_3(t)

yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); （2）function yt=x2_3(t)

yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=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])

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]')

（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.6649

-0.0239

-0.0630

-0.0314

-0.0127

从（1）（2）对比我们当抽样间隔越小时数值精度越高。

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); stem(k,y,'file');

M3-8

k=0:30;

a=[1 0.7 -0.45 -0.6];

b=[0.8 -0.44 0.36 0.02];

h=impz(b,a,k);

stem(k,h,'file')