ode45_solver.m

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

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