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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 318求調(diào)劑,計(jì)算材料方向 +7 | 吸喵有害笙命 2026-04-01 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +9 | 太容易1018 2026-04-01 | 9/450 |
|
|
[考研] 324分 085600材料與化工 +17 | 呆鵝oor 2026-03-27 | 17/850 |
|
|
[考研] 285求調(diào)劑 +7 | AZMK 2026-03-30 | 13/650 |
|
|
[考研] 350求調(diào)劑 +7 | 阿佳~ 2026-03-31 | 7/350 |
|
|
[考研] 0703求調(diào)劑 +4 | zizimo 2026-03-31 | 4/200 |
|
|
[考研] 材料0856 英一數(shù)二 323 求調(diào)劑 +9 | 袁sy 2026-04-01 | 9/450 |
|
|
[考研] 0710生物學(xué)求調(diào)劑! +6 | 敘述文 2026-03-31 | 6/300 |
|
|
[考研] 311(085601)求調(diào)劑 +12 | liziyeyeye 2026-03-28 | 13/650 |
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 4/200 |
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
[考研] 求調(diào)劑 +9 | akdhjs 2026-03-31 | 10/500 |
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
[考研] 求化學(xué)調(diào)劑 +12 | wulanna 2026-03-28 | 12/600 |
|
|
[考博] 材料專業(yè)申博 +5 | 杜雨婷dyt 2026-03-29 | 5/250 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|