xk3=cos(13*pi/5*k3); subplot(2,3,6) stem(k3,xk3) title('x3[k]')

运行结果如右上图：

2-8分别用square函数和 Swtooth画出题M2-8图所 示的离散周期序列。 程序如下：

%program2-81 k=0:1:40;

A=1;T=20;w0=2*pi/T;

xk=A*sawtooth(w0*k,0.5); stem(k,xk)

运行结果如右图所示：

%program2-82 k=0:1:40;

A=-1;T=10;w0=2*pi/T; xk=A*square(w0*k,50); stem(k,xk)

运行结果如右图：

实验二

3-1一个连续时间LTI 系统满足的微分方程为： y″(t)+3y′(t)+2y(t) =2x′(t)+x(t) %program3_1

ts=0;te=10;dt=0.001; sys=tf([2 1],[1 3 2]); t=ts:dt:te;

x=1*exp(-3*t)*1; y=lsim(sys,x,t); plot(t,y);

xlabel('Time(sec)') ylabel('y(t)')

3-2程序如下： %program3-2

ts=0;te=5;dt=0.01; sys=tf([6],[1 3]); t=ts:dt:te;

y=impulse(sys,t); plot(t,y);

xlabel('Time(sec)') ylabel('h(t)')

运行结果如有图所示：

%program3-2

ts=0;te=5;dt=0.01; sys=tf([6],[1 3]); t=ts:dt:te; y=step(sys,t); plot(t,y);

xlabel('Time(sec)') ylabel('h(t)')

运行结果如有图所示：

3-4程序如下： %program3_4

x=[0.85,0.53,0.21,0.67,0.84,0.12]; y=[0.68,0.37,0.83,0.52,0.71]; z=conv(x,y); n=length(z); stem(-3:n-4,z);

运行结果如有图所示：

3-7某离散时间LTI系统 满足差分方程

Y[k]+0.42y[k-1]-0.19y [k-2]=0.31x[k]+0.68x [k-1]，若x(t)=0.6ku(k) 试求系统的零状态响应 Yzs(k),并画出前20点图。 程序如下 %program3_7 k=0:100;

a=[1,0.42,-0.19]; b=[0.31,0.68]; A=1;a=0.6; xk=A*a.^k; stem(k,xk)

y=filter(b,a,xk); axis([0,20,0,1]) xlabel('k'); ylabel('Yzs[k]')

运行结果如右图所示：

3-8某离散时间LTI系统 满足差分方程

Y[k]+0.7y[k-1]-0.45y [k-2]-0.6y[k-3]=0.8x[k] -0.44x[k-1]+0.36x[k-2] +0.02x[k]，

试利用impz函数求其单位

脉冲响应，,并画出前30点图。 程序如下： %program3_8 k=0:100;

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)

axis([0,30,-1,0.8]) 运行结果如右图所示：

实验三

M4-1 试求题M4-1图所示周期矩形信号和周期三角波信号的频谱，并画出频谱图，取A=1，T=2。

得如下程序 4_1_1程序如下 k=-4:0.01:4;

A=1;T=2;w0=2*pi/T;

xk=A*square(w0*k,50); k1=-4:4;N=length(k);k1=(0:N)*0.5*T; x=A*square(w0*k1,50); X=fft(x);

subplot(3,1,1) stem(k1,real(X)) xlabel('k');

title('X[k]的实部'); axis([0,40,-20,60]) subplot(3,1,2) stem(k1,imag(X)); xlabel('k');

title('X[k]的虚部'); axis([0,40,-20,60]) xr=ifft(X); subplot(3,1,3)