找回密码
 注册
关于网站域名变更的通知
查看: 498|回复: 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);' }; r  u! P( M1 D! c3 `. a3 d+ W& e

4 s9 s8 k9 a* b6 Z2 P' |. d函数中有五个自变量,自变量之间没有任何的约束关系,但是自变量有范围,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], 劳烦大神指点看怎么编程,能求出自变量范围内的结果 ,谢谢0 L3 i! W+ K. x( u

, W* r! S! D# Q2 C; r
8 o1 b; ~8 ~2 Y  B

该用户从未签到

2#
 楼主| 发表于 2020-7-8 13:57 | 只看该作者
fitnessFunction=@minf_5;nvars=5;  P, ^% ]# r7 Q2 ^) Y) r, D0 B# h
options=gaoptimset;- d. l$ C' }+ D
options=gaoptimset(options,'PopInitRange',[0;1]);, H, o* J4 F  P/ v2 N+ A
options = gaoptimset(options,'PopulationSize',100);; U  e9 C0 D3 P
options=gaoptimset(options,'CrossoverFraction',0.9);
9 U) R; c" y) F4 y: j. I8 H! L% Zoptions=gaoptimset(options,'Generations',1000);9 y$ ]8 s' `8 v) m. a
options=gaoptimset(options,'CrossoverFcn',{@crossoverheuristic,0.9});2 ~5 C7 K; w- w1 p2 E: a
options=gaoptimset(options,'MutationFcn',{@mutationuniform,0.02});
& o" ]6 S/ q* W9 N[x fval reason output population scores]=ga(@minf_5,nvars,options)
: z, q' v# o, QOptimization terminated: average change in the fitness value less than options.TolFun.
5 `. g  T4 r" I/ H+ d1 J0 T0 p* R2 B( T" q: H  P
x =
8 z9 F0 M  R/ q  d- E& f& S( A, ~2 w7 B
    0.6026    0.0979    0.2421    0.9712    0.9449. s- K. Q, e7 I2 \% r7 {6 ]

: o/ i$ c1 L& q# G* Y0 ]7 n6 B4 J- P/ v
fval =1 P2 |1 @) m# X" E. A

8 W( l5 T3 k" \# u" M  W# }  -1.6983e+03 - 1.4575e+01i
8 Q+ j+ O8 ^6 d. Z( h' ~4 X8 W( p& o, J' |2 \
附上我自己的遗传算法编程,首先是第三行的自变量范围,只能定义一个区间,不知道怎么分开定义五个自变量的范围,再就是这样计算结果竟然是个复数,还都是负号,跟我想要的正实数实在是相去甚远,

该用户从未签到

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

该用户从未签到

4#
发表于 2020-7-9 10:50 | 只看该作者
' d3 v* j, S& M- L+ w+ m6 H2 s
理论上最小值就是-2.83599462789046,如果非要为正,那就是0了。

点评

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

该用户从未签到

5#
 楼主| 发表于 2020-7-9 14:14 | 只看该作者
zaiyiaaaa 发表于 2020-7-9 10:50! Y% f3 y* }  k# D- ^9 k7 T
理论上最小值就是-2.83599462789046,如果非要为正,那就是0了。
4 J& q/ p. V9 C6 e4 e7 i8 W
好吧,我再看看# e1 l! B& H. F0 s# r/ n
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-5 19:42 , Processed in 0.140625 second(s), 24 queries , Gzip On.

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

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

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