| 5 | 1/1 | 返回列表 |
| 查看: 1207 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
1stopt1.5解四元多項(xiàng)式方程組無穩(wěn)定解,求高版本幫忙運(yùn)行!
|
|||
|
大家好,我在算一個四個參數(shù)的非線性多項(xiàng)式方程組如下,得到的目標(biāo)函數(shù)大概0.13,求高版本的運(yùn)行結(jié)果,看是否有更準(zhǔn)確的解。 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,得到目標(biāo)函數(shù)是3.05,兩個情況的誤差差不多,但是解卻差不少,這是問什么呢? 我得到一個解,另它作為牛頓迭代法的初值去迭代,卻無解,這說明了什么呢? 急求!謝謝了 |
鐵桿木蟲 (職業(yè)作家)
鐵桿木蟲 (著名寫手)
新蟲 (小有名氣)
|
1、樓上dingd的回答可能有欠缺,4個方程4個未知數(shù),一般情況下(不加任何服務(wù)約束)也就只有一組解,如果是線性滿秩方程組的話,這是沒有錯的。如果是非線性方程組的話,存在多組解(均滿足設(shè)定的誤差規(guī)則)的情況就顯得很正常了。下面給出一組在所設(shè)定的范圍內(nèi)找到的一組解,具體誤差是否滿足你的要求,可以自己測試看 x2=0.00018602746110 k1=-0.65443662426355 k2= -0.35539381460637 k3=0.10000000000929 2、用那個做初值進(jìn)行牛頓迭代反而無解,只能說明你那個初值進(jìn)去,采用那種牛頓迭代算法無法找到滿足條件的解而已。 |

鐵桿木蟲 (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +5 | 加一一九 2026-03-24 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +5 | kotoko_ik 2026-03-23 | 6/300 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +4 | 722865 2026-03-23 | 4/200 |
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +7 | 生物工程調(diào)劑 2026-03-17 | 12/600 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +3 | WW.' 2026-03-23 | 5/250 |
|
|
[考研] 招08考數(shù)學(xué) +6 | laoshidan 2026-03-20 | 14/700 |
|
|
[考研] 333求調(diào)劑 +6 | 87639 2026-03-21 | 10/500 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 268求調(diào)劑 +9 | 簡單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|