24小时热门版块排行榜

>论坛更新日志 (3900)
>文献求助 (405)
>虫友互识 (340)
>导师招生 (316)
>博后之家 (160)
>硕博家园 (154)
>考博 (142)
>休闲灌水 (141)
>论文投稿 (114)
>基金申请 (90)
>招聘信息布告栏 (54)
>考研 (54)
>教师之家 (51)
>公派出国 (43)
>绿色求助(高悬赏) (37)
>找工作 (37)

返回列表

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;

回复此楼

» 猜你喜欢

求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有206人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

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

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

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.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

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

相关版块跳转我要订阅楼主 M.bySong 的主题更新

返回列表