ode45_solver.asv

自己编的四节龙格琨塔程序

%ode45执行函数
%解微分方程
clear all
%x=zeros(8,1);
ui=[5;-5;5;-5];
x0=[1;0;3e+7;6e+5;1;0;3e+7;6e+5];
[t,dxdt]=ode45(@valve_motor_pinion,[0 10],x0,[],ui);
%x=rk4(@valve_motor_pinion,0,10,x0,100);
subplot(2,2,1);plot(dxdt(:,1),'-');title('小齿轮1的角位移')；
subplot(2,2,2);plot(dxdt(:,2),'-.');title('小齿轮1的角速度')；
subplot（2,2,3）;plot(dxdt(:,3),'--');title('马达1的油液输入压力PA1');
subplot(2,2,4);plot(dxdt(:,4),':');title('马达1的油液输出压力PA2');