| 查看: 1924 | 回復: 17 | |||
xiashaoyan新蟲 (正式寫手)
|
[求助]
用遺傳算法求一個二元函數(shù)的最大值是一個復數(shù)范圍的值,為什么會這樣呢? 已有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)(是變量),其他參數(shù)都是已經(jīng)的常數(shù),現(xiàn)在我將兩個方程轉(zhuǎn)換成求一個函數(shù)的極大值 1-f^2-g^2,f為第一個方程組的左邊,g為第二個方程組的左邊,現(xiàn)在我用的是多種群遺傳算法求解上述函數(shù)的最大值,但是運行的結(jié)果是 Warning: Imaginary parts of complex X and/or Y arguments ignored > In MPGA at 51 最優(yōu)值為:449451.5954+2663.266383i 對應(yīng)的自變量取值:49.9593 0.316498 從結(jié)果來看,函數(shù)的最大值是449451.5954+2663.266383i,我很納悶,最大值按理來說是1呀,怎么會是一個復數(shù)范圍內(nèi)的數(shù)值呢?為什么會出現(xiàn)呢?是因為在實數(shù)范圍內(nèi)找不到解嗎?求計算的大神求助! |
榮譽版主 (文壇精英)
![]() |
專家經(jīng)驗: +518 |
新蟲 (正式寫手)
榮譽版主 (文壇精英)
![]() |
專家經(jīng)驗: +518 |
新蟲 (正式寫手)
|
好的, 參數(shù): 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個,現(xiàn)在我編程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(這些范圍也可能不準,以求到的解的范圍為準) |
木蟲之王 (文學泰斗)
peterflyer
| 問題應(yīng)該出在第二個函數(shù)式子里,因為里面有平方根的根號。這可能是在某些x1和x2進行取值時,根號里的結(jié)果為負導致了復數(shù)的產(chǎn)生。樓主在遺傳算法優(yōu)化過程中應(yīng)該嚴格注意各個函數(shù)的實數(shù)定義域,保證只在各個函數(shù)的實數(shù)定義域中進行搜索優(yōu)化就不會出現(xiàn)這個問題了。 |
榮譽版主 (文壇精英)
![]() |
專家經(jīng)驗: +518 |
新蟲 (正式寫手)
榮譽版主 (文壇精英)
![]() |
專家經(jīng)驗: +518 |
新蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 考研化學學碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 085600材料與化工 +4 | 安全上岸! 2026-03-16 | 4/200 |
|
|
[考研] 267一志愿南京工業(yè)大學0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 藥學383 求調(diào)劑 +3 | 藥學chy 2026-03-15 | 4/200 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 學碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學上- 2026-03-12 | 5/250 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|