内容发布更新时间 : 2024/11/16 9:21:19星期一 下面是文章的全部内容请认真阅读。

5*x^4+8*(1+tan(4*x^2)^2)*x 第6题：

f1=sym('x^3+3*x^2+5*x+7'); f2=sym('8*x^3-6*x^2+4*x-2'); f=f1*f2 collect(f) (f-f1)/f2 collect(ans) 结果： f =

(x^3+3*x^2+5*x+7)*(8*x^3-6*x^2+4*x-2) f =

8*x^6+18*x^5+26*x^4+36*x^3-28*x^2+18*x-14 h =

(-21+8*x^6+18*x^5+26*x^4+35*x^3-31*x^2+13*x)/(8*x^3-6*x^2+4*x-2) h =

(-21+8*x^6+18*x^5+26*x^4+35*x^3-31*x^2+13*x)/(8*x^3-6*x^2+4*x-2)

实验八：线性时不变系统的时域响应

一、实验目的

1．掌握线性时不变系统的三种描述形式——传递函数描述法、零极点增益描述法、状态空间描述法。

2．掌握三种描述形式之间的转换。

3．掌握连续和离散系统频率响应的求解。

二、实验内容（边做实验，边将生成结果和图形拷贝到Word文档中）

1．生成20个点的单位脉冲信号、单位阶跃信号，并记录下函数命令和波形。 2．生成占空比为30％的矩形波。 3．将连续系统H(s)?0.5达式。

3?5z?1?2z?24．将离散系统H(z)?转化为零极点增益的描述?1?2?3?41?1.6z?1.3z?0.9z?0.5z(s?1)(s?3)转化为传递函数的形式，并显示其表

(s?1)(s?2)(s?4)形式，并显示其表达式。

5．分别求实验内容3和4的频率响应（对离散系统取256样点，采样频率取8000Hz）。

6．分别求实验内容3和4的单位冲激响应（对离散系统，作60样点图）。

三、设计提示

1．显示传递函数模型用tf(b,a)；显示零极点增益模型用zpk(z,p,k)。注意：z、p为列向量。

2．连续系统频率响应用freqs函数；离散系统用freqz函数。 3．连续系统冲激响应用impulse函数；离散系统用impz函数。

四、实验报告要求

1．编写实现第二节实验内容中所使用的函数文件，并记录相应的生成结果。 2．书写实验报告时要结构合理，层次分明，在分析描述的时候，需要注意语言的流畅。

五、参考答案 第一题：

k=-4:15

x =[zeros(1,4),1,zeros(1,15)]

subplot (1,2,1) stem(k,x)

y=[zeros(1,4),ones(1,16)] subplot (1,2,2) stem(k,y)

结果：

第二题

f=2;

t=0:0.001:2 w=2*pi*f*t

y=square(w,30) axis([0 2 0 2]) hold on plot(t,y)

结果：

第3题： k=0.5;

z=[1,-3]';

p=[-1,-2,-4]'; sys_zpk=zpk(z,p,k) [b,a]=zp2tf(z,p,k); sys_tf=tf(b,a) 结果：

Zero/pole/gain: 0.5 (s-1) (s+3) ----------------- (s+1) (s+2) (s+4)

Transfer function: 0.5 s^2 + s - 1.5 ---------------------- s^3 + 7 s^2 + 14 s + 8 第4题：

b=[3 5 2];

a=[1 -1.6 1.3 -.9 .5];

systf=tf(b,a,'variable','z^-1')%显示传递函数 [z,p,k]=tf2zp(b,a);

syszpk=zpk(z,p,k,'variable','z^-1')%显示零点增益

结果：

Transfer function:

3 + 5 z^-1 + 2 z^-2

--------------------------------------------- 1 - 1.6 z^-1 + 1.3 z^-2 - 0.9 z^-3 + 0.5 z^-4 Sampling time: unspecified Zero/pole/gain:

3 z^-2 (1+z^-1) (1+0.6667z^-1)

------------------------------------------------------------- (1 - 1.685z^-1 + 0.8654z^-2) (1 + 0.08497z^-1 + 0.5778z^-2) Sampling time: unspecified 第5题： k=0.5;

z=[1,-3]';

p=[-1,-2,-4]';

[b,a]=zp2tf(z,p,k); freqs(b,a) 结果：

b=[3 5 2];

a=[1 -1.6 1.3 -.9 .5];

freqz(b,a,256,8000) % Fs=8000 结果： 第6题： k=0.5;

z=[1,-3]';

p=[-1,-2,-4]'; sys3=zpk(z,p,k);

impulse(sys3)%频率响应 结果：

b=[3 5 2];

a=[1 -1.6 1.3 -.9 .5]; impz(b,a,60)