第三次实验课课堂作业(4学时)
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
MATLAB 第三次实验课课堂作业（4学时）
学号 姓名
注：1）此课堂作业作为本课程结业成绩的重要依据，请同学们认真、独立完成，不得抄袭。
2）请在授课教师规定的时间内完成；
3）完成作业后，请以word 格式保存，文件名为：学号+姓名
已知系统的差分方程为：
（1）给出该系统的零点向量q 和极点向量p 值，并在z 平面上绘制零极点图
（给出命令语句）；
a=[0.2 0.1 0];b=[1 -0.4 -0.5];
[q,p,k]=tf2zp(a,b)
q =
0 -0.5000 p =
0.9348 -0.5348 k =
0.2000
>> zplane(q,p)
[]0.4[1]0.5[2]0.2[]0.1[1]
y n y n y n x n x n ----=+-
-1
-0.5
00.5
1
-1
-0.8-0.6-0.4-0.200.20.4
0.60.8
1Real Part
I m a g i n a r y P a r t
（2）给出求该系统的单位冲激响应，并作出冲激响应的时域图（给出命令语句）； a=[0.2 0.1 0];b=[1 -0.4 -0.5]; [h,t]=impz(a,b,7,1000) h =
0.2000 0.1800 0.1720 0.1588 0.1495 0.1392 0.1304 t =
0 0.0010 0.0020 0.0030 0.0040 0.0050 0.0060
>> subplot(1,2,1);stem(t,h) >> subplot(1,2,2);impz(a,b)
2
4
6
x 10
-3
n (samples)
A m p l i t u d e
（3）绘制该系统的频响特性图（给出命令语句）; >> num=[0.2 0.1 0] num =
0.2000 0.1000 0
>> den=[1 -0.4 -0.5] den =
1.0000 -0.4000 -0.5000
>> num=[0.2 0.1 0] num =
0.2000 0.1000 0
>> freqz(num,den)
00.10.2
0.30.40.50.60.70.80.91
-80
-60-40-20
Normalized Frequency (⨯π rad/sample)
P h a s e (d e g r e e s )
00.10.2
0.30.40.50.60.70.80.91
-20
-10
10
Normalized Frequency (⨯π rad/sample)
M a g n i t u d e (d B )
（4）试在Simulink 下绘制系统的结构图，并给出该系统的单位阶跃响应波形图。
Gain3
plot(tout,yout)
0102030405060708090100
0.5
1
1.5
2
2.5
3