实验2截图及结果
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
要求：输入示例程序1：
示例程序1计算数据
Eigenvalue Damping Freq. (rad/s)
0.00e+000 + 1.00e+001i -6.12e-017 1.00e+001 0.00e+000 - 1.00e+001i -6.12e-017 1.00e+001
Eigenvalue Damping Freq. (rad/s)
-2.50e+000 + 9.68e+000i 2.50e-001 1.00e+001 -2.50e+000 - 9.68e+000i 2.50e-001 1.00e+001
Eigenvalue Damping Freq. (rad/s)
-5.00e+000 + 8.66e+000i 5.00e-001 1.00e+001 -5.00e+000 - 8.66e+000i 5.00e-001 1.00e+001
Eigenvalue Damping Freq. (rad/s)
-7.50e+000 + 6.61e+000i 7.50e-001 1.00e+001 -7.50e+000 - 6.61e+000i 7.50e-001 1.00e+001
Eigenvalue Damping Freq. (rad/s)
-1.00e+001 1.00e+000 1.00e+001 -1.00e+001 1.00e+000 1.00e+001
Eigenvalue Damping Freq. (rad/s)
-5.00e+000 1.00e+000 5.00e+000 -2.00e+001 1.00e+000 2.00e+001
要求：输入示例程序2：
要求：<2>求系统G(s)=10/(s2+2s+10)的特征根、闭环极点和单位阶跃响应。分别改变系统阻尼比和自振角频率，观察系统单位阶跃响应的变化规律。
程序：
num=10; i=0;
for sigma=0: 0.25: 1.25
den=[1 2*sigma*sqrt(10) 10];
damp(den)
sys=tf(num, den);
i=i+1;
step(sys,2)
hold on
end
grid
hold off
title('阻尼比不同时阶跃响应曲线')
lab1='ζ=0'; text(0.3,0.9,lab1),
lab2='ζ=0.25'; text(0.3,1.5,lab2),
lab3='ζ=0.5'; text(0.3,1.2,lab3),
lab4='ζ=0.75'; text(0.3,1.05,lab4),
lab5='ζ=1'; text(0.35,0.9,lab5),
lab6='ζ=1.25'; text(0.35,0.8,lab6)
程序图像：
程序计算数据：
Eigenvalue Damping Freq. (rad/s)
0.00e+000 + 3.16e+000i -6.12e-017 3.16e+000 0.00e+000 - 3.16e+000i -6.12e-017 3.16e+000
Eigenvalue Damping Freq. (rad/s)
-7.91e-001 + 3.06e+000i 2.50e-001 3.16e+000
-7.91e-001 - 3.06e+000i 2.50e-001 3.16e+000 Eigenvalue Damping Freq. (rad/s)
-1.58e+000 + 2.74e+000i 5.00e-001 3.16e+000 -1.58e+000 - 2.74e+000i 5.00e-001 3.16e+000
Eigenvalue Damping Freq. (rad/s)
-2.37e+000 + 2.09e+000i 7.50e-001 3.16e+000 -2.37e+000 - 2.09e+000i 7.50e-001 3.16e+000
Eigenvalue Damping Freq. (rad/s)
-3.16e+000 1.00e+000 3.16e+000
-3.16e+000 1.00e+000 3.16e+000
Eigenvalue Damping Freq. (rad/s)
-1.58e+000 1.00e+000 1.58e+000
-6.32e+000 1.00e+000 6.32e+000
要求：<3>已知系统结构如图，使用梅森公式和SIMULINK结构图模型求系统的闭环传递函数。
将系统结构图存为smg.mdl文件，可以下面程序转化为传递函数模型：
[A，B，C，D]=linmod(‘smg’); %将结果图转化为状态空间模型
[num,den]=ss2tf(A, B, C, D); %将状态空间模型转化为多项式模型
printsys(num, den, ‘s’) %输出传递函数
程序：
程序计算结果：
num/den =
-1 s^4 - 10 s^3 - 23 s^2 + 10 s + 24 ---------------------------------------------
s^5 + 15 s^4 + 85 s^3 + 269 s^2 + 598 s + 600
>>