找回密码
 注册
关于网站域名变更的通知
查看: 472|回复: 4
打印 上一主题 下一主题

多维函数

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-7-8 13:36 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
函数如下:y=(0.0206*x(2)^2/(2*9.81*x(1))+(0.5*(1-(x(2)/(4*x(3)*9.81*(x(4)-1040)/(3*0.43*1040))^0.5*(1-(x(3)/x(1))^2)*2.7183^(-(2.65*x(5)-3.32*x(5)^2.2))))+((0.25*(x(2)/(4*x(3)*9.81*(x(4)-1040)/(3*0.43*1040))^0.5*(1-(x(3)/x(1))^2)*2.7183^(-(2.65*x(5)-3.32*x(5)^2.2)))-1)^2+x(5)*x(2)/(4*x(3)*9.81*(x(4)-1040)/(3*0.43*1040))^0.5*(1-(x(3)/x(1))^2)*2.7183^(-(2.65*x(5)-3.32*x(5)^2.2)))^0.5)*(x(4)/1040-1)+1.8413*(x(2)/(9.81*x(1)))^2.7736*(0.5*(1-(x(2)/(4*x(3)*9.81*(x(4)-1040)/(3*0.43*1040))^0.5*(1-(x(3)/x(1))^2)*2.7183^(-(2.65*x(5)-3.32*x(5)^2.2))))+((0.25*(x(2)/(4*x(3)*9.81*(x(4)-1040)/(3*0.43*1040))^0.5*(1-(x(3)/x(1))^2)*2.7183^(-(2.65*x(5)-3.32*x(5)^2.2)))-1)^2+x(5)*x(2)/(4*x(3)*9.81*(x(4)-1040)/(3*0.43*1040))^0.5*(1-(x(3)/x(1))^2)*2.7183^(-(2.65*x(5)-3.32*x(5)^2.2)))^0.5)*0.43*x(3)/x(1)*(x(4)/1040-1))/(x(4)*x(5)/1040);5 q2 u. f* F& D1 r6 z

+ U  v  X5 y3 s/ S函数中有五个自变量,自变量之间没有任何的约束关系,但是自变量有范围,x(1)=[0.2,0.4],x(2)=[1.234,3],x(3)=[0.01,0.04],x(4)=[1190,1398],x(5)=[0.15,0.4], 劳烦大神指点看怎么编程,能求出自变量范围内的结果 ,谢谢" G$ \1 n* |; F6 J3 C3 a
* h' H( _5 V# u& A* }* h) i* Z
7 x% r+ l; R' @  f$ L* m+ y

该用户从未签到

2#
 楼主| 发表于 2020-7-8 13:57 | 只看该作者
fitnessFunction=@minf_5;nvars=5;
2 K5 {5 Y, q  ]  T0 N: R' [5 Yoptions=gaoptimset;" ?4 l% K8 h1 x8 i6 b: K" R8 D
options=gaoptimset(options,'PopInitRange',[0;1]);
4 h  l9 F% H8 L  coptions = gaoptimset(options,'PopulationSize',100);( p9 Q; G5 k) f9 r( [1 Z+ \( E. t
options=gaoptimset(options,'CrossoverFraction',0.9);
4 K+ e' g% G, E& _7 @8 Voptions=gaoptimset(options,'Generations',1000);
  Z% A2 D: O2 N, ^  k3 m: Yoptions=gaoptimset(options,'CrossoverFcn',{@crossoverheuristic,0.9});
9 k' P7 f& ~% J8 _4 y0 Y0 u3 l' moptions=gaoptimset(options,'MutationFcn',{@mutationuniform,0.02});7 Z& s' g9 o+ V8 f% h5 y
[x fval reason output population scores]=ga(@minf_5,nvars,options)
2 k$ v+ P% A+ l) H# jOptimization terminated: average change in the fitness value less than options.TolFun.
% w+ H0 m. U$ a" N  A; f+ x/ m6 i- A: d- s& W
x =5 M, j& F% E! r& x% n7 [% q! V
" m7 ]( _2 X; a" }
    0.6026    0.0979    0.2421    0.9712    0.9449: Z  N2 ]9 s4 Y/ V
/ W+ R  O0 F3 s, c; Z  U

8 H; B% j9 }( @! S% }+ f' r% xfval =/ p/ g: F& j2 s/ m) l* M

1 |2 w0 R6 G1 Y+ J  -1.6983e+03 - 1.4575e+01i
3 @* i7 X, L/ y2 k8 k+ K: x$ Q/ B, `1 p1 m
附上我自己的遗传算法编程,首先是第三行的自变量范围,只能定义一个区间,不知道怎么分开定义五个自变量的范围,再就是这样计算结果竟然是个复数,还都是负号,跟我想要的正实数实在是相去甚远,

该用户从未签到

3#
发表于 2020-7-8 14:02 | 只看该作者
帮你顶一下

该用户从未签到

4#
发表于 2020-7-9 10:50 | 只看该作者
* g  Q6 _) y8 d6 r* ^
理论上最小值就是-2.83599462789046,如果非要为正,那就是0了。

点评

好吧,我再看看  详情 回复 发表于 2020-7-9 14:14

该用户从未签到

5#
 楼主| 发表于 2020-7-9 14:14 | 只看该作者
zaiyiaaaa 发表于 2020-7-9 10:50! t# n9 L4 C6 k5 C+ z
理论上最小值就是-2.83599462789046,如果非要为正,那就是0了。
' D2 b3 i0 t4 S- P( I
好吧,我再看看
! b! w# L8 v: b8 C, c+ V6 B
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-7-26 09:36 , Processed in 0.109375 second(s), 27 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表