直线制动优化设计MATLAB程序

直线制动优化设计MATLAB程序
%质量300KG，轴距1.6M，质心到前轴距离0.848M，质心高度0.353M,首先建立目标函数
function m=zhidong_m(x);
l=1.6;a=0.848;b=1-a;h=0.353;G=300*9.8;%赛车原始数据
z=0.2:0.05:0.8;ii=length(z);
for i=1:ii
faif=x.*z.*l./(b+z.*h)
fair=(1-x).*z.*l./(a-z.h)
dm=(faif-z).^2+(fair-z).2;
end
m=max(dm);
%建立约束函数
function[c,ceq]=zhidong_y(x);
for z=0.2:0.05:0.8;
l=1.6;a=0.848;b=1-a;h=0.353;G=3009.8;%赛车原始数据
faif=x.*z.*l./(b+z.*h);
fair=(1-x).*z.l./(a-z.h);
ceq=[];
if z>0.3&z<0.45
c(1)=fair-z-0.05;
end
c(2)=fair-(z+0.07)/0.85;
c(3)=0.85(faif-0.2)+0.1-z;
c(4)=0.85(fair-0.2)+0.1-z;
c(5)=fair-faif;
end

%主程序
clear all
clc
lb=0.4;
ub=0.6;
x0=0.4; %制动力分配比初值
l=1.6;a=0.848;b=1-a;h=0.353;G=3009.8;%赛车原始数据
[beta,y]=fmincon(‘zhidong_m’,x0,[],[],[],[],lb,ub,‘zhidong_y’)
z=0.2:0.1:0.8;
x=beta;%前后制动力分配比
fai0=(lbeta-b)/h;%同步附着系数
faif=x.z.l./(b+z.h);
fair=(1-x).z.l./(a-z.h);
Er=a./(l.(1-beta)+fair.h);%后轴制动效率
Ef=b./(l.beta-faif.h);%前轴制动效率
f1=0.1:10:3000;
f2=0.5(G/hsqrt(b.^2+4hlf1/G)-(Gb/h+2f1));%I曲线
ff2=(1-beta)./beta.f1;%实际制动力分配曲线
plot(f1,f2,f1,ff2)
grid on
fai=0.7;%路面附着系数
v=80;%制动初速度
am=bfai9.8/(b+(fai0-fai)h);%最大制动减速度计算
s=0.4v./3.6+v^2/25.92/abs(am);%制动距离
fprintf(‘汽车最佳制动力分配比beta=%3.3f \n’,beta)
fprintf(‘汽车最大减速度am=%3.2fm/s2 \n’,am)
fprintf(‘汽车同步附着系数fai0=%3.3f \n’,fai0)
fprintf(‘汽车车速为80公里每小时的制动距离distance=%3.3f m \n’,s)

这个程序有点小问题，按照自己的参数修改一下

按照这个程序执行的结果：
汽车最佳制动力分配比beta=0.468
汽车最大减速度am=2.08m/s2
汽车同步附着系数fai0=1.690
汽车车速为80公里每小时的制动距离distance=127.664 m