将混合动力电动重型车辆具有速度约束的生态驾驶问题的全局最优解与恒速轮廓驾驶进行了比较（Matlab代码实现）

 👨‍🎓个人主页：研学社的博客 

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🌈3 Matlab代码实现

🎉4 参考文献

---

💥1 概述

在本文中，将混合动力电动重型车辆具有速度约束的生态驾驶问题的全局最优解与恒速轮廓驾驶进行了比较。

• paramterets.m 包含问题中使用的参数集。在此文件中，可以修改速度和加速度约束。
• problem_config.m包括与问题相关的一般参数，例如采样时间。

📚2 运行结果

 部分代码：

function [ par ] = parameters( par )

%Electric Motor - parameters
par.em.beta2 = 2*2.526044703610551e-04;
par.em.beta1 = 2*1.005;
par.em.beta0 = 2*0.292803996942460;

%Vehicle - parameters
par.veh.Af = 7.54; %frontal area
par.veh.Cdrag = 0.7; %drag coefficient
par.veh.Croll = 0.007; %rolling coefficient
par.veh.rhoa = 1.184; 
par.veh.m = 15.95*1000;
par.veh.Cd = 0.5*par.veh.Cdrag*par.veh.rhoa*par.veh.Af; %drag constant
par.veh.Cr =par.veh.Croll*9.81*par.veh.m; %rolling resistance force
par.veh.g = 9.81;
par.veh.gamma1 = (1.5*2*pi)/(21000);
par.veh.gamma0 = -pi/4;
par.veh.beta = -((450)/2)*par.veh.gamma1;
par.veh.s0 = 0;
par.veh.sN = 21000;
par.veh.v0 = 70/3.6; % %[m/s]
par.veh.vN = 70/3.6; % %[m/s]
par.veh.vMax = 80/3.6; %[m/s]
par.veh.vMin = 60/3.6; %60/3.6; %[m/s]
par.veh.aMax=0.2*9.8;  %[m/s^2]
par.veh.aMin=-0.2*9.8; %[m/s^2]

end

、

fig = figure;
    left_color = [39, 75, 130]/255;
    right_color =  [36, 96, 63]/255;
    set(fig,'defaultAxesColorOrder',[left_color; right_color]);
    
    subplot(2,1,1)%---------------------
    
    yyaxis left
    axis([0, max(s/1000),min([v_bl;v]*3.6) ,max([v_bl;v]*3.6)])
    
    plot(s_bl/1000,v_bl*3.6,'LineWidth',4)
    hold on
    plot(s/1000,v*3.6,'-.','LineWidth',4,'Color',[67, 137, 198]/255)

    xlabel('Displacement \sl s \rm [km]')
    ylabel('Velocity \sl v \rm [km/h]')
    hold on
    
    yyaxis right
    ar=area(s_bl/1000,h_bl);
    ar.FaceAlpha =0.3;
    ar.FaceColor = [36, 96, 63]/255;
    ar.LineStyle  =':';
    grid on
    legend('Constant speed driving', 'Eco-driving');
    
    
    
   subplot(2,1,2)%------------------------------------
    
    yyaxis left
    axis([0, max(s/1000) ,min(u/1000) ,max(u/1000)])
  
    plot(s/1000,u/1000,'-.','LineWidth',4,'Color',[67, 137, 198]/255)
    xlabel('Displacement \sl s \rm [km]')
    ylabel('Traction Force \sl u \rm [kN]')
    
    yyaxis right
    ar=area(s_bl/1000,h_bl);
    ar.FaceAlpha =0.3;
    ar.FaceColor = [36, 96, 63]/255;
    ar.LineStyle  =':';
    ylabel('Elevation \sl h \rm [m]')
    grid on

   saveas(gcf, 'results/ExampleB.fig')
   print('results/plot', '-dpdf')
    
    

🌈3 Matlab代码实现

🎉4 参考文献

部分理论来源于网络，如有侵权请联系删除。