找回密码
 注册
关于网站域名变更的通知
查看: 508|回复: 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);
0 K* D+ u1 y# [5 ]% L2 X
" D3 I; t" Q* f& 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], 劳烦大神指点看怎么编程,能求出自变量范围内的结果 ,谢谢
$ m: d! Y$ u2 ]  a4 ^5 t1 a/ y* S" e0 }: Z% X% E  E$ A/ ?- l8 f# _

8 E8 u' x% O3 {2 a

该用户从未签到

2#
 楼主| 发表于 2020-7-8 13:57 | 只看该作者
fitnessFunction=@minf_5;nvars=5;
3 |) B) c4 }2 A  ^2 Yoptions=gaoptimset;6 c+ O3 D: t8 k1 R; `
options=gaoptimset(options,'PopInitRange',[0;1]);
3 n4 M% |% Y; U3 N8 foptions = gaoptimset(options,'PopulationSize',100);
$ ~: {4 X# f" I# Xoptions=gaoptimset(options,'CrossoverFraction',0.9);, s! f) l0 N/ l6 F7 ]$ v
options=gaoptimset(options,'Generations',1000);
! A, F! V9 v. _* r+ [options=gaoptimset(options,'CrossoverFcn',{@crossoverheuristic,0.9});
3 r+ h" p9 B9 Yoptions=gaoptimset(options,'MutationFcn',{@mutationuniform,0.02});0 ?# j- f% \! t% f+ c
[x fval reason output population scores]=ga(@minf_5,nvars,options)
) q1 z3 E2 C: e; R1 |) YOptimization terminated: average change in the fitness value less than options.TolFun.
$ G5 l* h! x* k
( m! J# T& \" @% Lx =
, I9 ?: ]/ }, T; t, T, G% u- W3 v8 A) L$ [$ s
    0.6026    0.0979    0.2421    0.9712    0.9449
( ~! J0 D! {- B3 p8 t  r* t% u4 |0 d( d1 {' U6 t) I

% Q* f6 F% N' jfval =
. n$ g. B" I; V" @, Q/ X: Y
0 h) ~8 ~/ J3 Z! G3 v4 E+ M4 I# j7 O8 N  -1.6983e+03 - 1.4575e+01i
( _& Y  u) K, c' Y( D6 y8 C) t5 j5 W. Z' ^9 T% g' m( q# {
附上我自己的遗传算法编程,首先是第三行的自变量范围,只能定义一个区间,不知道怎么分开定义五个自变量的范围,再就是这样计算结果竟然是个复数,还都是负号,跟我想要的正实数实在是相去甚远,

该用户从未签到

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

该用户从未签到

4#
发表于 2020-7-9 10:50 | 只看该作者

, I% O2 f" M$ J8 Q. D- Q9 D理论上最小值就是-2.83599462789046,如果非要为正,那就是0了。

点评

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

该用户从未签到

5#
 楼主| 发表于 2020-7-9 14:14 | 只看该作者
zaiyiaaaa 发表于 2020-7-9 10:50( z, ^0 \- K. L9 c; W. Y9 J
理论上最小值就是-2.83599462789046,如果非要为正,那就是0了。
, E- i4 g% y& N4 w& Q* v. r
好吧,我再看看
8 Z8 P- {1 D$ y) v+ s! C8 H
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-25 00:18 , Processed in 0.187500 second(s), 25 queries , Gzip On.

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

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

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