24小时热门版块排行榜    

查看: 1440  |  回复: 11

阿斯顿求爱者

新虫 (初入文坛)

[求助] 求帮忙用遗传算法优化下这个目标函数,帮帮忙啦 已有2人参与

目标函数是下面这个,使Z最小,bcdef的范围分别为(0,0.1)(0,0.1)(0,0.1)(1550,1750),(0,8),求帮帮忙了 没有金币了,给程序的可以给你冲十块话费
function Z=fanyan5(b,c,d,e,f)
f1=-0.61454+30./(1500*sin(acos(1500*cos(b)./e)))+cos(pi/180)*(30*tan(pi/180)./tan(acos(1500*cos(b)./e))+f)./(e.*sin(b-pi/180))+cos(pi/180)*(1000*tan(pi/180)+f)./(e.*sin(b-pi/180));
f2=-1000+30./(tan(acos(1500*cos(b)./e)))+cos(pi/180)*cos(b).*(30*tan(pi/180)./tan(acos(1500*cos(b)./e))+f)./(sin(b-pi/180))+cos(pi/180).*cos(b-pi/90).*(1000*tan(pi/180)+f)./(sin(b-pi/180));
f3=-0.62666+30./(1500*sin(acos(1500*cos(c)./e)))+cos(pi/180)*(30*tan(pi/180)./tan(acos(1500*cos(c)./e))+f)./(e.*sin(c-pi/180))+cos(pi/180)*(1020*tan(pi/180)+f)./(e.*sin(c-pi/180));
f4=-1020+30./(tan(acos(1500*cos(c)./e)))+cos(pi/180)*cos(c).*(30*tan(pi/180)./tan(acos(1500*cos(c)./e))+f)./(sin(c-pi/180))+cos(pi/180).*cos(c-pi/90).*(1020*tan(pi/180)+f)./(sin(c-pi/180));
f5=-0.63878+30./(1500*sin(acos(1500*cos(d)./e)))+cos(pi/180)*(30*tan(pi/180)./tan(acos(1500*cos(d)./e))+f)./(e.*sin(d-pi/180))+cos(pi/180)*(1040*tan(pi/180)+f)./(e.*sin(d-pi/180));
f6=-1040+30./(tan(acos(1500*cos(d)./e)))+cos(pi/180)*cos(d).*(30*tan(pi/180)./tan(acos(1500*cos(d)./e))+f)./(sin(d-pi/180))+cos(pi/180).*cos(d-pi/90).*(1040*tan(pi/180)+f)./(sin(d-pi/180));
Z=0.5*(f1.^2+f2.^2+f3.^2+f4.^2+f5.^2+f6.^2);
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
阿斯顿求爱者: 金币+5, ★★★★★最佳答案 2015-07-19 18:28:59
目标函数值(最小): 1.79624201608983E-16
b: 0.0471509580783501
c: 0.0468948188165174
d: 0.0466491725546775
e: 1650.04770675958
f: 4.59200659099393
3楼2015-06-17 21:56:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

醇熟的先森

新虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
试试看  我运行一下    你这没有其他约束条件了吗
2楼2015-06-17 08:18:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿斯顿求爱者

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by dingd at 2015-06-17 21:56:07
目标函数值(最小): 1.79624201608983E-16
b: 0.0471509580783501
c: 0.0468948188165174
d: 0.0466491725546775
e: 1650.04770675958
f: 4.59200659099393

你好 想问下你这是用matlab 的程序处理的吗 可以私聊下
4楼2015-07-16 21:19:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
4楼: Originally posted by 阿斯顿求爱者 at 2015-07-16 21:19:17
你好 想问下你这是用matlab 的程序处理的吗 可以私聊下...

是用1stOpt算的,自己可以试试!
5楼2015-07-17 11:01:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿斯顿求爱者

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by dingd at 2015-07-17 11:01:20
是用1stOpt算的,自己可以试试!...

请问用的1stopt哪个版本呢,我用的1.5, 老是提示输入的不正确 麻烦你了
6楼2015-07-17 11:53:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿斯顿求爱者

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by dingd at 2015-07-17 11:01:20
是用1stOpt算的,自己可以试试!...

大神 1stopt里的代码怎么写的呢 这样不行吗 求教了
ConstStr  -0.6145354547+30/(1500*sin(acos(1500*cos(b)/e)))+cos(pi/180)*(30*tan(pi/180)/tan(acos(1500*cos(b)/e))+f)/(e*sin(b-pi/180))+cos(pi/180)*(1000*tan(pi/180)+f)/(e*sin(b-pi/180))=g;
ConstStr  -1000+30/(tan(acos(1500*cos(b)/e)))+cos(pi/180)*cos(b)*(30*tan(pi/180)/tan(acos(1500*cos(b)/e))+f)/(sin(b-pi/180))+cos(pi/180)*cos(b-pi/90)*(1000*tan(pi/180)+f)/(sin(b-pi/180))=h;
ConstStr  -0.6266559814+30/(1500*sin(acos(1500*cos(c)/e)))+cos(pi/180)*(30*tan(pi/180)/tan(acos(1500*cos(c)/e))+f)/(e*sin(c-pi/180))+cos(pi/180)*(1020*tan(pi/180)+f)/(e*sin(c-pi/180))=i;
ConstStr  -1020+30/(tan(acos(1500*cos(c)/e)))+cos(pi/180)*cos(c)*(30*tan(pi/180)/tan(acos(1500*cos(c)/e))+f)/(sin(c-pi/180))+cos(pi/180)*cos(c-pi/90)*(1020*tan(pi/180)+f)/(sin(c-pi/180))=j;
ConstStr  -0.6387765073+30/(1500*sin(acos(1500*cos(d)/e)))+cos(pi/180)*(30*tan(pi/180)/tan(acos(1500*cos(d)/e))+f)/(e*sin(d-pi/180))+cos(pi/180)*(1040*tan(pi/180)+f)/(e*sin(d-pi/180))=k;
ConstStr  -1040+30/(tan(acos(1500*cos(d)/e)))+cos(pi/180)*cos(d)*(30*tan(pi/180)/tan(acos(1500*cos(d)/e))+f)/(sin(d-pi/180))+cos(pi/180)*cos(d-pi/90)*(1040*tan(pi/180)+f)/(sin(d-pi/180))=l;
Parameter b[0,1],c[0,1],d[0,1],e[1550,1750],f[0,8];
MinFunction 0.5*(g^2+h^2+i^2+j^2+k^2+l^2);
7楼2015-07-17 20:17:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

CODE:
ConstStr  g=-0.6145354547+30/(1500*sin(arccos(1500*cos(b)/e)))+cos(pi/180)*(30*tan(pi/180)/tan(arccos(1500*cos(b)/e))+f)/(e*sin(b-pi/180))+cos(pi/180)*(1000*tan(pi/180)+f)/(e*sin(b-pi/180));
ConstStr  h=-1000+30/(tan(arccos(1500*cos(b)/e)))+cos(pi/180)*cos(b)*(30*tan(pi/180)/tan(arccos(1500*cos(b)/e))+f)/(sin(b-pi/180))+cos(pi/180)*cos(b-pi/90)*(1000*tan(pi/180)+f)/(sin(b-pi/180));
ConstStr  i=-0.6266559814+30/(1500*sin(arccos(1500*cos(c)/e)))+cos(pi/180)*(30*tan(pi/180)/tan(arccos(1500*cos(c)/e))+f)/(e*sin(c-pi/180))+cos(pi/180)*(1020*tan(pi/180)+f)/(e*sin(c-pi/180));
ConstStr  j=-1020+30/(tan(arccos(1500*cos(c)/e)))+cos(pi/180)*cos(c)*(30*tan(pi/180)/tan(arccos(1500*cos(c)/e))+f)/(sin(c-pi/180))+cos(pi/180)*cos(c-pi/90)*(1020*tan(pi/180)+f)/(sin(c-pi/180));
ConstStr  k=-0.6387765073+30/(1500*sin(arccos(1500*cos(d)/e)))+cos(pi/180)*(30*tan(pi/180)/tan(arccos(1500*cos(d)/e))+f)/(e*sin(d-pi/180))+cos(pi/180)*(1040*tan(pi/180)+f)/(e*sin(d-pi/180));
ConstStr  l=-1040+30/(tan(arccos(1500*cos(d)/e)))+cos(pi/180)*cos(d)*(30*tan(pi/180)/tan(arccos(1500*cos(d)/e))+f)/(sin(d-pi/180))+cos(pi/180)*cos(d-pi/90)*(1040*tan(pi/180)+f)/(sin(d-pi/180));
Parameter b=[0,1],c=[0,1],d=[0,1],e=[1550,1750],f=[0,8];
MinFunction 0.5*(g^2+h^2+i^2+j^2+k^2+l^2);

b: 0.0458544215030493
c: 0.045625376869945
d: 0.0454057168442035
e: 1650.00336089837
f: 3.98464708386331
8楼2015-07-17 23:27:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿斯顿求爱者

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by dingd at 2015-07-17 23:27:53
ConstStr  g=-0.6145354547+30/(1500*sin(arccos(1500*cos(b)/e)))+cos(pi/180)*(30*tan(pi/180)/tan(arccos(1500*cos(b)/e))+f)/(e*sin(b-pi/180))+cos(pi/180)*(1000*tan(pi/180)+f)/(e*sin(b-pi/180));
Const ...

你好,我想知道我写的那个程序可以运行吗 我在1。5版本上老提示有误,需要修改哪儿吗
9楼2015-07-19 11:50:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
9楼: Originally posted by 阿斯顿求爱者 at 2015-07-19 11:50:51
你好,我想知道我写的那个程序可以运行吗 我在1。5版本上老提示有误,需要修改哪儿吗...

改过后的代码上面已给出,老版本能否运行自己试试。
10楼2015-07-19 13:08:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿斯顿求爱者 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 X5bHxOULxFqb 2026-10-30 23/1150 2026-11-01 23:49 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 23:45 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 21:52 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 21:45 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 21:26 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 14:24 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见