西电电院自动控制技术上机报告40页.doc

自动控制技术上机实验报告班级：021215 学号：021214 姓名：时域分析程序源代码：close all;clear all;ft = 30;M=1;B=5;K=20; %系统参数t0=0;tfinal = 5;tspan = t0 tfinal; %设置仿真开始和结束时间x0 = 0,0; %系统初始值，零初始条件options = odeset(AbsTol,1e-6;1e-6); %设置仿真计算精度t,x = ode113(xt4odefile,tspan,x0,options);%微分方程求解，计算位移x(:,1)和速度x(:,2) a = 1/M*(ft-B*x(:,2)-K*x(:,1); %计算加速度 i = 1; while (abs(a(i)0.0001|(abs(x