| 查看: 2182 | 回复: 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呀，怎么会是一个复数范围内的数值呢？为什么会出现呢？是因为在实数范围内找不到解吗？求计算的大神求助！ |
回复此楼» 猜你喜欢
时间戳又变了
已经有13人回复
-
科研孤儿太难了
已经有4人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，8月21号(明天)放榜
已经有16人回复
-
时间戳变了，能看出什么问题？
已经有21人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
只有每年这种时候来逛逛小木虫
已经有21人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 二元函数极值求解
已经有5人回复
- C++中实现复杂二元函数曲线
已经有8人回复
- 请教函数求解问题
已经有11人回复
- 金币求助--求一个二元函数的最大值和其他极值
已经有5人回复
- 二元一次函数求值问题
已经有8人回复
- 【求助】二元函数拟合
已经有12人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2014-08-20 21:37:49
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 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
xiashaoyan
新虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 2574.9
- 散金: 30
- 红花: 10
- 帖子: 498
- 在线: 184.1小时
- 虫号: 2507802
- 注册: 2013-06-14
- 专业: 计算机软件
|
好的， 参数： 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; b3=2.760513*10^-5; t3=1.064; z=31.5; 其实方程原来是5个，现在我编程2个了 原来的5个方程为： 1/x(1)+1/x(2)-1/t3=0 x(5)/t3-x(4)/x(2)-x(3)/x(1)-1/z=0 x(3)^2-(a1+b1*f+(a2+b2*f)/(x(1)^2-a3^2)+(a4+b3*f)/(x(1)^2-a5^2)-a6*x(1)^2)=0 x(4)^2-(a1+b1*f+(a2+b2*f)/(x(2)^2-a3^2)+(a4+b3*f)/(x(2)^2-a5^2)-a6*x(2)^2)=0 x(5)^2-(a1+b1*f+(a2+b2*f)/(t3^2-a3^2)+(a4+b3*f)/(t3^2-a5^2)-a6*t3^2)=0 x(1)-x(5)的取值大概在：1.4-1.8，1.4-1.8，1.5-2.5,1.5-2.5（这些范围也可能不准，以求到的解的范围为准） |
5楼2014-08-20 21:45:40
40