版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

M.bySong

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 419.1
散金: 26
帖子: 145
在线: 10.3小时
虫号: 4403056
注册: 2016-02-15
专业: 零件成形制造

[求助] 1stopt高版本代跑，很急，先谢谢大神们了！已有1人参与

大神上午好，最近碰到一个多目标优化的问题，想用理想点法与平方和加权法求多目标最优。三个目标函数，分别跑出最优值以后，加权求解。希望大神能够帮忙代跑一下，不胜感激。就这点金币了，完成后加v有大红包，现在很急，先谢谢大神们了！！！

求第一个目标目标函数最优值op1：
parameter x1=[6422,10000],x2=[0,6422],x3=[0,6422],h1=[0,2280],h2=[1600,10000],h3=[1241,10000] ;
minfunction 6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3));
x1-10000<=0
x3-x1<0;
x2-x3<0;
h2-h1<0;
h3-h2<0;
h3-x2<0;
sqrt((x1-x3)^2+(h1-h3)^2)-3000<=0;
-sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h<0;
sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-8150<=0;
(sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-sqrt((x1-x3)^2+(h1-h3)^2))/(2*sqrt((x1-x3)^2+(h1-h3)^2))-0.8<=0;
x1-(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6000<=0;
-x1+(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)<=0;
sqrt((x1-x2)^2+(-h3+h2)^2)-7000<0;
-6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3))<0;
——————————————————————————————————————————————————————————————————————————————————第2个目标目标函数最优值op2：
parameter x1=[6422,10000],x2=[0,6422],x3=[0,6422],h1=[0,2280],h2=[1600,10000],h3=[1241,10000] ;
minfunction sqrt((x1-x2)^2+(-h3+h2)^2);
x1-10000<=0
x3-x1<0;
x2-x3<0;
h2-h1<0;
h3-h2<0;
h3-x2<0;
sqrt((x1-x3)^2+(h1-h3)^2)-3000<=0;       //收拢长度
-sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h<0;          //H
sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-8150<=0;             //展开长度
(sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-sqrt((x1-x3)^2+(h1-h3)^2))/(2*sqrt((x1-x3)^2+(h1-h3)^2))-0.8<=0;    //导向比
x1-(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6000<=0;       //导向行程
-x1+(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)<=0;       //导向行程
sqrt((x1-x2)^2+(-h3+h2)^2)-7000<0;                         //拉杆长度
-6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3))<0;

——————————————————————————————————————————————————————————————————————————————————第3个目标目标函数最优值op3：
parameter x1=[6422,10000],x2=[0,6422],x3=[0,6422],h1=[0,2280],h2=[1600,10000],h3=[1241,10000] ;
minfunction sqrt((h3-269)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6422)^2);
x1-10000<=0
x3-x1<0;
x2-x3<0;
h2-h1<0;
h3-h2<0;
h3-x2<0;
sqrt((x1-x3)^2+(h1-h3)^2)-3000<=0;       //收拢长度
-sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h<0;          //H
sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-8150<=0;             //展开长度
(sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-sqrt((x1-x3)^2+(h1-h3)^2))/(2*sqrt((x1-x3)^2+(h1-h3)^2))-0.8<=0;    //导向比
x1-(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6000<=0;       //导向行程
-x1+(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)<=0;       //导向行程
sqrt((x1-x2)^2+(-h3+h2)^2)-7000<0;                         //拉杆长度
-6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3))<0;

————————————————————————————————————————————————————————————————————————
将下面代码中的op1，op2，op3替换为上面计算出的结果，求多目标加权结果
parameter x1=[6422,10000],x2=[0,6422],x3=[0,6422],h1=[0,2280],h2=[1600,10000],h3=[1241,10000] ;
ConstStr f1=6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3));
conststr f2=minfunction sqrt((x1-x2)^2+(-h3+h2)^2);
conststr f3=sqrt((h3-269)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6422)^2);
minfunction ((f1-op1)/op1)^2+((f2-op2)/op2)^2+((f3-op3)/op3)^2;
//minfunction 6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3));
//minfunction sqrt((x1-x2)^2+(-h3+h2)^2);
//minfunction sqrt((h3-269)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6422)^2);
x1-10000<=0
x3-x1<0;
x2-x3<0;
h2-h1<0;
h3-h2<0;
h3-x2<0;
sqrt((x1-x3)^2+(h1-h3)^2)-3000<=0;       //收拢长度
-sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h<0;          //H
sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-8150<=0;             //展开长度
(sqrt((x3-h3)^2+((sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)+h1)^2)-sqrt((x1-x3)^2+(h1-h3)^2))/(2*sqrt((x1-x3)^2+(h1-h3)^2))-0.8<=0;    //导向比
x1-(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)-6000<=0;       //导向行程
-x1+(sqrt((x1-x2)^2+(h2-h3)^2-(x2-h3)^2)-h2)<=0;       //导向行程
sqrt((x1-x2)^2+(-h3+h2)^2)-7000<0;                         //拉杆长度
-6422*600000*sqrt((h1-h3)^2+(x1-x3)^2)*(x1-x2)/((-x3*h3+x1*h1)*(x1-x2)-(-x2*h3+x1*h2)*(x1-x3))<0;

回复此楼

» 猜你喜欢

重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。已经有0人回复
UJN物理学专业调剂已经有12人回复
物理学I论文润色/翻译怎么收费? 已经有121人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名已经有0人回复
山西大同大学物理学硕士研究生招收调剂生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
一志愿物理所凝聚态物理352求调剂已经有0人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有1人回复

1楼 2019-07-11 11:09:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31016.3
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 856.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
小红豆: 金币+20, 应助指数+1, 回帖应助 2019-07-11 14:39:04
M.bySong: 金币+38 2019-07-11 14:53:51

目标函数值(最小): 1.87232999181561
x1: 6422
x2: 2282.96873670372
x3: 4485.76040752828
h1: 2279.99154310903
h2: 1600.03227522495
h3: 1241.37911564834
h: -1006.24586353898

赞一下(2人)

回复此楼

数值计算

2楼2019-07-11 14:11:41