| 查看: 1926 | 回復(fù): 17 | |||
xiashaoyan新蟲 (正式寫手)
|
[求助]
用遺傳算法求一個(gè)二元函數(shù)的最大值是一個(gè)復(fù)數(shù)范圍的值,為什么會(huì)這樣呢? 已有2人參與
|
|
求兩個(gè)二元非線性方程組的解; 方程組為: 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)在我將兩個(gè)方程轉(zhuǎn)換成求一個(gè)函數(shù)的極大值 1-f^2-g^2,f為第一個(gè)方程組的左邊,g為第二個(gè)方程組的左邊,現(xiàn)在我用的是多種群遺傳算法求解上述函數(shù)的最大值,但是運(yùn)行的結(jié)果是 Warning: Imaginary parts of complex X and/or Y arguments ignored > In MPGA at 51 最優(yōu)值為:449451.5954+2663.266383i 對(duì)應(yīng)的自變量取值:49.9593 0.316498 從結(jié)果來(lái)看,函數(shù)的最大值是449451.5954+2663.266383i,我很納悶,最大值按理來(lái)說(shuō)是1呀,怎么會(huì)是一個(gè)復(fù)數(shù)范圍內(nèi)的數(shù)值呢?為什么會(huì)出現(xiàn)呢?是因?yàn)樵趯?shí)數(shù)范圍內(nèi)找不到解嗎?求計(jì)算的大神求助! |
榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
新蟲 (正式寫手)
榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +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; 其實(shí)方程原來(lái)是5個(gè),現(xiàn)在我編程2個(gè)了 原來(lái)的5個(gè)方程為: 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(這些范圍也可能不準(zhǔn),以求到的解的范圍為準(zhǔn)) |
木蟲之王 (文學(xué)泰斗)
peterflyer
| 問(wèn)題應(yīng)該出在第二個(gè)函數(shù)式子里,因?yàn)槔锩嬗衅椒礁母?hào)。這可能是在某些x1和x2進(jìn)行取值時(shí),根號(hào)里的結(jié)果為負(fù)導(dǎo)致了復(fù)數(shù)的產(chǎn)生。樓主在遺傳算法優(yōu)化過(guò)程中應(yīng)該嚴(yán)格注意各個(gè)函數(shù)的實(shí)數(shù)定義域,保證只在各個(gè)函數(shù)的實(shí)數(shù)定義域中進(jìn)行搜索優(yōu)化就不會(huì)出現(xiàn)這個(gè)問(wèn)題了。 |
榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
新蟲 (正式寫手)
榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
新蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
[考研] 考研調(diào)劑 +4 | 芬達(dá)46 2026-03-12 | 4/200 |
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123! 2026-03-11 | 4/200 |
|
|
[考研] 化工學(xué)碩306求調(diào)劑 +9 | 42838695 2026-03-12 | 9/450 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|