实验一 计算机仿真研究
B06050614 张叶锋
1) 实验内容
考虑如图所示的质量-弹簧-阻尼系统,当质量系数 m=1,弹簧刚性系数 k=4 时,
为了使系统的单位阶跃响应不发生振荡,阻尼系数 f(0<=f<=10)应在什么范围
内取值。
2)实验代码:
>> m=1;k=4;
>> f=input('请输入阻尼系数 f:');
请输入阻尼系数 f:1
>> t=0;T=0.01;
>> A=[0 1;-k/m -f/m];
>> B=[0 1/m]';
>> tmax=10;
>> x=[0 0]';
>> Y=0;
>> H=t;
>> while(t<tmax)
xs=x+(A*x+B)*T;
y=xs(1);
t=t+T;
Y=[Y;y];H=[H;t];
x=xs;
end
>> plot(H,Y,'k');
>> grid
3)实验结果
得到最终的图如下,可以看出系统不稳定,发生振荡。
评论0