| 查看: 2359 | 回复: 17 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
xiashaoyan新虫 (正式写手)
|
[求助]
用遗传算法求一个二元函数的最大值是一个复数范围的值，为什么会这样呢？ 已有2人参与
|
||||
|
求两个二元非线性方程组的解； 方程组为： 1/x(1)+1/x(2)-1/t3=0， sqrt(a1+b1*f+(a2+b2*f)/(t3^2-a3^2)+(a4+b3*f)/(t3^2-a5^2)-a6*t3^2)/t3-sqrt(a1+b1*f+(a2+b2*f)/(X(i,2)^2-a3^2)+(a4+b3*f)/(X(i,2)^2-a5^2)-a6*X(i,2)^2)/X(i,2)-sqrt(a1+b1*f+(a2+b2*f)/(X(i,1)^2-a3^2)+(a4+b3*f)/(X(i,1)^2-a5^2)-a6*X(i,1)^2)/X(i,1)-1/z=0， 方程中除了x(1),x(2)（是变量）,其他参数都是已经的常数，现在我将两个方程转换成求一个函数的极大值 1-f^2-g^2,f为第一个方程组的左边，g为第二个方程组的左边，现在我用的是多种群遗传算法求解上述函数的最大值，但是运行的结果是 Warning: Imaginary parts of complex X and/or Y arguments ignored > In MPGA at 51 最优值为：449451.5954+2663.266383i 对应的自变量取值：49.9593 0.316498 从结果来看，函数的最大值是449451.5954+2663.266383i，我很纳闷，最大值按理来说是1呀，怎么会是一个复数范围内的数值呢？为什么会出现呢？是因为在实数范围内找不到解吗？求计算的大神求助！ |
回复此楼» 猜你喜欢
27博士申请
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
厉娜
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有12人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有15人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有16人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有15人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3916
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24222
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiashaoyan: 金币+5, ★★★很有帮助 2014-08-21 16:53:14
感谢参与，应助指数 +1
xiashaoyan: 金币+5, ★★★很有帮助 2014-08-21 16:53:14
|
Constant T=50, f=(T-24.5)/(T+570.82), a1=5.319725, a2=0.09147285, a3=0.3165008, a4=100.2028, a5=11.37639, a6=0.01497046, //b1=4.753469*10^-7, //b2=3.310965*10^-8, b1=0, b2=0, b3=2.760513*10^-5, t3=1.064, z=31.5; Parameters x1, x2, x3, x4, x5; Function 1/x1+1/x2-1/t3=0; x5/t3-x4/x2-x3/x1-1/z=0; x3^2-(a1+b1*f+(a2+b2*f)/(x1^2-a3^2)+(a4+b3*f)/(x1^2-a5^2)-a6*x1^2)=0; x4^2-(a1+b1*f+(a2+b2*f)/(x2^2-a3^2)+(a4+b3*f)/(x2^2-a5^2)-a6*x2^2)=0; x5^2-(a1+b1*f+(a2+b2*f)/(t3^2-a3^2)+(a4+b3*f)/(t3^2-a5^2)-a6*t3^2)=0; 优化算法: 粒子群算法CFPSO 1 函数表达式 1: 1/x1+1/x2-1/1.064-(0) 2: x5/1.064-x4/x2-x3/x1-1/31.5-(0) 3: x3^2-(5.319725+0*0.0410747076447279+(0.09147285+0*0.0410747076447279)/(x1^2-0.3165008^2)+(100.2028+2.760513e -5*0.0410747076447279)/(x1^2-11.37639^2)-0.01497046*x1^2)-(0) 4: x4^2-(5.319725+0*0.0410747076447279+(0.09147285+0*0.0410747076447279)/(x2^2-0.3165008^2)+(100.2028+2.760513e -5*0.0410747076447279)/(x2^2-11.37639^2)-0.01497046*x2^2)-(0) 5: x5^2-(5.319725+0*0.0410747076447279+(0.09147285+0*0.0410747076447279)/(1.064^2-0.3165008^2)+(100.2028 +2.760513e-5*0.0410747076447279)/(1.064^2-11.37639^2)-0.01497046*1.064^2)-(0) 目标函数值: 0.000393372471519715 x1: 2.14627793576053 x2: 2.11178358874413 x3: 2.11382505280947 x4: 2.11472993815495 x5: 2.1471740096698 ====== 计算结束 ====== |
11楼2014-08-21 15:36:11
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3916
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24222
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2014-08-20 20:58:24
xiashaoyan
新虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 2574.9
- 散金: 30
- 红花: 10
- 帖子: 498
- 在线: 184.1小时
- 虫号: 2507802
- 注册: 2013-06-14
- 专业: 计算机软件
3楼2014-08-20 21:25:00
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3916
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24222
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2014-08-20 21:37:49
