西电电院自动控制技术上机报告.docVIP

自动控制技术上机实验报告 班级：021215 学号4 姓名 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(i,2))0.0001)) i = i+1; end disp(稳态时系统的位移、速度和加速度及对应的时间分别为： ); %显示计算结果 result = sprintf(位移 d=%6.4f\n,x(i,1)); disp(result); result = sprintf(速度 v=%8.6f\n,x(i,2)); disp(resul