数维杯节能列车运行优化问题一思路及代码

一、问题背景

 

二、问题描述

假设一辆列车在水平轨道上运行，从站台A运行至站台B，其间距为5144.7m，运行的速度上限为100km/h，列车质量为176.3t，列车旋转部件惯性的旋转质量 因数p = 1.08 ，列车电机的最大牵引力为310KN，机械制动部件的最大制动力为 760KN。列车受到的阻力满足Davis阻力方程fv = 2.0895 + 0.0098v + 0.006v 2 ，该 公式中的速度单位为m/s，阻力单位为KN。 你如何通过建模方法编写程序以获得列车运行过程的速度-距离曲线、牵引 制动力-距离曲线、时间-距离曲线与能量消耗-距离曲线？程序的运行时间是多 长？需要获取列车以最短时间到达站台B、在最短运行时间上分别增加10s、20s、50s、150s、300s到达站台B总共六组曲线。

三、问题分析

当时在选择题目时，初步观察来看，感觉B题更适合自己，而且感觉上是要比A题更简单一些，所以选择了B题，但是真正做起来才发现。。。。。

B题的第一问需要求解的东西：最短运行时间、最短运行时间下的四条曲线、改变时间后的 四条曲线。

1.首先是对于最短时间求解。

求最短时间，首先要明确，应该将整个路程分为三个阶段（问题1不考虑惰行），牵引阶段、匀速阶段、制动阶段。只要保证牵引阶段以最大牵引力加速，制动阶段以最大制动力减速，就能保证总的时间为最短。

 可以看出这是前两个是微分方程，于是想到求解微分方程数值解最常用的龙格库塔方法。

主程序

clear;clc;close all
%主程序：
t0=0;
tf=28;   %时间区间
V0=0;    %初始速度为0
[t,Vt]=ode45(@jiasu,[t0,tf],V0);

%绘图牵引阶段
figure(1)
plot(t,Vt)
xlabel('时间')
ylabel('速度')
title('牵引阶段的V-T图')


t0=0;
tf=7;   %重新给定时间区间
V0=27.78;    %制动开始速度为27.78

[t1,Vt1]=ode45(@zhidong,[t0,tf],V0);

%绘图制动阶段
figure(2)
plot(t1,Vt1)
xlabel('时间')
ylabel('速度')
title('制动阶段的V-T图')

%绘图牵引路程
t2=(0:0.1:12);
[t2,y]=ode45('lucheng1',t,[0;0]);%t2是时间，y是路程ds;
%路程1函数返回值是[ds,d2s];给定时间区间，和初始值，用龙格库塔求数值解解
figure(3)
plot(t2,y(:,1))
xlabel('时间')
ylabel('路程')
title('牵引阶段路程时间曲线')


%绘图制动路程
t3=(0:1:4)';
[t3,y1]=ode45('lucheng2',t,[100;27.78]);%t2是时间，y是路程ds;
%路程1函数返回值是[ds,d2s];给定时间区间，和初始值，用龙求数值解
figure(4)
plot(t3,y1(:,1))
xlabel('时间')
ylabel('路程')
title('制动阶段路程时间曲线')

加速方程函数代码

function xdot=jiasu(t,v)
%牵引阶段微分方程
Fq=310000;fv=(2089.5-9.8*v-6*v*v);mp=190404;
xdot=(Fq-fv)/mp;
end

制动方程函数代码

function xdot=zhidong(t,v)
%制动阶段微分方程
B=760000;fv=(2089.5-9.8*v-6*v*v);mp=190404;
xdot=(B+fv)/-mp;
end

牵引阶段运动方程

function dx = lucheng1(t,x)%返回值 函数名 输入值
dx=zeros(2,1);
dx(1)=x(2);
dx(2)=(310000-2089.5-9.8*x(1)-6*x(1)*x(1))/190404;   %dx(2)是加速度

end

制动阶段运动方程

function dx = lucheng2(t,x)%返回值 函数名 输入值
dx=zeros(2,1);
dx(1)=x(2);
dx(2)=(-760000-2089.5-9.8*x(1)-6*x(1)*x(1))/190404;   %dx(2)是加速度

end

得到两个路程的距离以及时间，之后在计算得出匀速阶段的路程，时间。三个时间相加得到最短距离时间。大概为197.3s。