matlab simulink汽车优化设计遍历法

75 阅读 0 评论 50 点赞

我是靠谱客的博主灵巧战斗机，最近开发中收集的这篇文章主要介绍matlab simulink汽车优化设计遍历法，觉得挺不错的，现在分享给大家，希望可以做个参考。

概述

1、内容简介

略
455-可以交流、咨询、答疑

2、内容说明

略

3、仿真分析

clc
close all
clear
m1 = 55;
m2 = 400;
k1 = 300000;
k2 = 40000;
b2 = 4000;
g = 9.8;
t = 0:0.0001:2;
v = 20;%车速
h = 0.05;
L = 3;
R = (L^2/4+h^2)/2/h;%路面不平顺等效的圆半径
delta = asin(L/2/R);
degre = -delta:0.0001:delta;
x0 = R*sin(degre)+L/2;
y0 = R*cos(degre)-R+h;
figure
plot(x0,y0)
xlabel 水平位移
ylabel 路面不平顺
x = t*v;
y = zeros(1,length(x));
for i = 1:length(x)
if x(i)<L
y(i) = interp1(x0,y0,x(i));
end
end
figure
plot(x,y)
xlabel 水平距离
ylabel 路面不平顺
title 完整路面的路面
simin = [t' y']; %路面数据输入到simulink

sim model1.mdl
figure
plot(tout,yout(:,4))
xlabel 时间/s
ylabel 车体位移/m
fop_init = trapz(tout,abs(yout(:,1))); %vehicle comfort
k_1 = 200000:50000:400000;
k_2 = 30000:5000:50000;
b_2 = 500:500:5000;
temp = [];
for i = 1:length(k_1)
for j = 1:length(k_2)
for m = 1:length(b_2)
k1 = k_1(i);
k2 = k_2(j);
b2 = b_2(m);
sim model1.mdl
fop_1 = trapz(tout,abs(yout(:,1)));
% max(-yout(:,3))*300000<=0.8*(m1+m2)*g
if max(yout(:,2))<=0.08 && max(-yout(:,3))*300000<=0.8*(m1+m2)*g %判断结果是否满足限制条件
aa = 1
if fop_1<fop_init %如果本次结果比上次的值小就更新目标函数值
bb = 1
fop_init = fop_1;
temp = [temp; i j m];
end
end
end
end
end
k11 = 300000;
k21 = 40000;
b21 = 4000;
if isempty(temp)
disp('未能得到最优的结果')
disp(['k1 = ',num2str(k11)])
disp(['k2 = ',num2str(k21)])
disp(['b2 = ',num2str(b21)])
else
disp('优化得到的结果为')
disp(['k1 = ',num2str(k_1(temp(end,1)))])
disp(['k2 = ',num2str(k_2(temp(end,2)))])
disp(['b2 = ',num2str(b_2(temp(end,3)))])
end