自动控制上机作业..docxVIP

自动控制上机大作业 班级： 学号： 姓名： 1.1 设质量——阻尼——弹簧系统的微分运动方程为 式中，x(t)为位移输出信号，f(t)为输入的力信号。质量为M=1kg，粘性摩擦系数为 B = 5N / m? s?1，弹簧的弹性系数为K=20N/m。当t=0 时，施加外力f(t)=30N，试问系统 何时达到稳定？并画出该机械系统位移、速度随时间变化的曲线以及速度与位移的关系曲 线。 提示：龙格-库塔法求解微分方程数值解的函数：odel13(), 调用方式：[T,Y] = ODE113(ODEFUN,TSPAN,Y0,OPTIONS)。 其中ODEFUN 为用户自定义的系统微分方程的描述，本题中可使用xt4odefile.m 文件定义的函数； TSPAN 表示计算开始和结束的时间；Y0 表示微分方程的初始条件；OPTION 为计算精度的可选参数，由 odese()函数设置。 odel13()函数只接受一阶微分方程的形式，使用时需要先将高阶方程化为若干个一阶微分方程； 绘图函数：plot()，subplot(); 程序： ft = 30; M=1; B=5; K=20; %系统参数 tspan = [0 20]; %设置仿真开始和结束时间 x0 = [0,0]; %系统初始值，零初始条件 options = odeset