找回密码
 注册
关于网站域名变更的通知
查看: 505|回复: 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);/ y! I, R7 f- u: B/ W$ T

# p; U9 c" s* y$ k- 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], 劳烦大神指点看怎么编程,能求出自变量范围内的结果 ,谢谢
$ I! `5 r6 S! q& N7 [3 O5 i
( T7 Z" _; S2 V9 s! o0 K8 q3 Z% y( w6 G5 C# o" G1 ?

该用户从未签到

2#
 楼主| 发表于 2020-7-8 13:57 | 只看该作者
fitnessFunction=@minf_5;nvars=5;( f- m0 \  _0 @% D9 A$ \- M& ]9 M
options=gaoptimset;& x3 w. t" }+ Y- c
options=gaoptimset(options,'PopInitRange',[0;1]);
- K9 ]- n1 O1 d+ eoptions = gaoptimset(options,'PopulationSize',100);1 w5 }8 L& [6 m( G1 d! m
options=gaoptimset(options,'CrossoverFraction',0.9);
5 {5 _* e  y. S. p  ~( xoptions=gaoptimset(options,'Generations',1000);
, T( Q5 \( `' coptions=gaoptimset(options,'CrossoverFcn',{@crossoverheuristic,0.9});
) o( I. D. x$ @7 `( q/ M, ]/ Zoptions=gaoptimset(options,'MutationFcn',{@mutationuniform,0.02});
1 I$ H& s' E8 D( Y  T[x fval reason output population scores]=ga(@minf_5,nvars,options)
0 [4 {8 ]" o/ r- bOptimization terminated: average change in the fitness value less than options.TolFun.
( T' g* d+ h. \9 Y, N; u
6 c4 Z1 J- u8 r. ]x =
% D( U; O& K2 u4 {- V. a
9 n! Y- J3 m* B- T  n    0.6026    0.0979    0.2421    0.9712    0.9449
# u( h' A) u( z, G, w9 l9 u1 {
" C* T/ o6 K4 G4 ?; T/ P
+ z% r1 I/ v* \fval =
; m3 @5 L% S& b0 n
( O; F% ?$ q7 ^; R  -1.6983e+03 - 1.4575e+01i& A9 j7 F/ f* g9 v, K$ C+ o
3 u2 S5 U7 p, D# [1 l; q9 a
附上我自己的遗传算法编程,首先是第三行的自变量范围,只能定义一个区间,不知道怎么分开定义五个自变量的范围,再就是这样计算结果竟然是个复数,还都是负号,跟我想要的正实数实在是相去甚远,

该用户从未签到

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

该用户从未签到

4#
发表于 2020-7-9 10:50 | 只看该作者
& d: v3 L% L# N) k+ X5 u* O
理论上最小值就是-2.83599462789046,如果非要为正,那就是0了。

点评

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

该用户从未签到

5#
 楼主| 发表于 2020-7-9 14:14 | 只看该作者
zaiyiaaaa 发表于 2020-7-9 10:50
4 b7 b8 w( {$ e( [3 [9 Q8 `理论上最小值就是-2.83599462789046,如果非要为正,那就是0了。
; P, x7 w) v. p( ^5 J! y6 i
好吧,我再看看+ y% h: W0 t6 h; k" Z$ i3 J
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-24 17:25 , Processed in 0.171875 second(s), 24 queries , Gzip On.

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

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

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