| 5 | 1/1 | 返回列表 |
| 查看: 1208 | 回復: 6 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[交流]
1stopt1.5解四元多項式方程組無穩(wěn)定解,求高版本幫忙運行!
|
|||
|
大家好,我在算一個四個參數的非線性多項式方程組如下,得到的目標函數大概0.13,求高版本的運行結果,看是否有更準確的解。 Parameters x2[0.0001,0.00025],K1[-3.3,3.3],K2[-3.3,3.3],K3[0.1,9] ; Function -0.0002483857122054368* K1^2 - 0.00016214005342224397* K3^2 +K1* (0.11832038801446097 + 0.0004143201485611828* K3 -2.597716388604742* x2) -6791.987324581259*(-0.09254221017610412 + x2) *(0.0014465166404352989 + x2) +(-0.07341654746178977 + 2.166562300403883* x2) = 0.8738; 0.9126034755805909 - 7.408480476916429*10^(-7)* K1^2 -1.8521201192291073*10^(-7)* K2^2 - 0.0036164973644227326* K3 +0.0018405191448191696* K3^2 + 4.83607293527163*10^(-7)* K3^3 -1.3692495660106428*10^(-7)* K3^4 + K1 *(0.006461909187108292 + 7.408480476916429*10^(-7)* K2 + 1.2357726636344952*10^(-6)* K3 - 6.178863318172476*10^(-7)* K3^2 ) + K2 *(-0.003230954593554146 - 6.178863318172476*10^(-7)* K3 + 3.34494856322119*10^(-7)* K3^2 ) + 0.004893093285395663* K3^2* x2 = 0.9034; 0.9119390763671137 - 0.000014217569378577079* K1^2 - 3.5543923446442698*10^(-6)* K2^2 - 0.0158429576267302* K3 + 0.008055685207646872* K3^2 + 9.280877865753863*10^(-6)* K3^3 - 2.62771843187007*10^(-6)* K3^4 + K1* (0.02830798508143708+ 0.000014217569378577079* K2 + 0.0000237156372836732* K3 - 0.0000118578186418366* K3^2 ) + K2* (-0.01415399254071854 - 0.0000118578186418366* K3 + 6.419270242197259*10^(-6)* K3^2 ) +0.02143540054716326* K3^2 *x2= 0.8912000000000001; 0.9108915129935684 - 0.00006986738872089048*K1^2 - 0.00001746684718022262*K2^2 - 0.03512051199327993* K3 + 0.017832729459977437* K3^2 + 0.00004560770440092129* K3^3 - 0.000012913024740019884* K3^4 + K1 *(0.06275286174349187 + 0.00006986738872089048* K2 + 0.0001165423993892175* K3 - 0.00005827119969460875* K3^2 ) + K2* (-0.031376430871745936- 0.00005827119969460875* K3 + 0.0000315453111128709* K3^2 ) +0.04751778422529179* K3^2 *x2 = 0.8808; 另外,我原方程是等式兩百乘以500,得到目標函數是3.05,兩個情況的誤差差不多,但是解卻差不少,這是問什么呢? 我得到一個解,另它作為牛頓迭代法的初值去迭代,卻無解,這說明了什么呢? 急求!謝謝了 |
新蟲 (小有名氣)
|
1、樓上dingd的回答可能有欠缺,4個方程4個未知數,一般情況下(不加任何服務約束)也就只有一組解,如果是線性滿秩方程組的話,這是沒有錯的。如果是非線性方程組的話,存在多組解(均滿足設定的誤差規(guī)則)的情況就顯得很正常了。下面給出一組在所設定的范圍內找到的一組解,具體誤差是否滿足你的要求,可以自己測試看 x2=0.00018602746110 k1=-0.65443662426355 k2= -0.35539381460637 k3=0.10000000000929 2、用那個做初值進行牛頓迭代反而無解,只能說明你那個初值進去,采用那種牛頓迭代算法無法找到滿足條件的解而已。 |

鐵桿木蟲 (著名寫手)
鐵桿木蟲 (職業(yè)作家)
鐵桿木蟲 (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調劑 +4 | Koxui 2026-03-24 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 278求調劑 +4 | 我可以上岸的對?/a> 2026-03-19 | 4/200 |
|
|
[考研] 340求調劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 資源與環(huán)境 調劑申請(333分) +7 | holy J 2026-03-21 | 7/350 |
|
|
[考研] 求調劑 +5 | 林之夕 2026-03-24 | 5/250 |
|
|
[考研] 一志愿211 初試270分 求調劑 +5 | 谷雨上岸 2026-03-23 | 6/300 |
|
|
[考研] 321求調劑 +4 | Ymlll 2026-03-24 | 4/200 |
|
|
[考研] 【雙一流院校新能源、環(huán)境材料,材料加工與模擬招收大量調劑】 +4 | Higraduate 2026-03-22 | 7/350 |
|
|
[考研] 一志愿北京化工大學 070300 學碩 336分 求調劑 +7 | vv迷 2026-03-22 | 7/350 |
|
|
[考研] 工科0856求調劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
[考研] 350求調劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 317求調劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
[考研] 308求調劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 尋找調劑 +4 | 倔強芒? 2026-03-21 | 4/200 |
|
|
[考研] 285求調劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 085601調劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 330求調劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 329求調劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
[考研] 一志愿西安交通大學 學碩 354求調劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|