 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 1200 | 回復: 6 | |||
[交流]
1stopt1.5解四元多項式方程組無穩(wěn)定解，求高版本幫忙運行！！
|
|
大家好，我在算一個四個參數(shù)的非線性多項式方程組如下，得到的目標函數(shù)大概0.13，求高版本的運行結(jié)果，看是否有更準確的解。 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，得到目標函數(shù)是3.05，兩個情況的誤差差不多，但是解卻差不少，這是問什么呢？ 我得到一個解，另它作為牛頓迭代法的初值去迭代，卻無解，這說明了什么呢? 急求！謝謝了 |
回復此樓» 猜你喜歡
哈爾濱理工大學物理系招收物理學考研調(diào)劑
已經(jīng)有6人回復
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有86人回復
-
求調(diào)劑
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復
-
光學工程學碩調(diào)劑信息
已經(jīng)有31人回復
-
歡迎加入課題組
已經(jīng)有0人回復
-
散金幣，求好運，祝面上順利！
已經(jīng)有65人回復
-
華東師范大學芯片設計徐瓏真實水平咋樣
已經(jīng)有1人回復
-
調(diào)劑
已經(jīng)有0人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 1stopt解非線性方程組結(jié)果不穩(wěn)定
已經(jīng)有2人回復
- 求大神幫忙用1stopt 5.0解一個方程組，原先的版本太低，每次結(jié)果不一樣。
已經(jīng)有11人回復
- 用1stOpt解出非線性方程組的解作為初值用Matlab求解，解不出來
已經(jīng)有4人回復
- 誰能幫忙用1Stopt運行下程序呀？謝謝，我的1.5版本的運行不了啊
已經(jīng)有6人回復
- 求助1stopt擬合動力學參數(shù)
已經(jīng)有4人回復
- 求助：以下非線性擬合方程可以用1stopt求解么？用什么方法
已經(jīng)有16人回復
- 1stopt擬合問題求助
已經(jīng)有7人回復
- 求大神用1stopt擬合數(shù)據(jù) 緊急
已經(jīng)有6人回復
- 1stopt 解非線性方程組
已經(jīng)有3人回復
- 求助，1stopt 解方程組為什么總是出現(xiàn)常數(shù)設置錯誤呢
已經(jīng)有8人回復
- 求助，1stopt解非線性方程組。
已經(jīng)有2人回復
- sos!!1stopt解非線性方程組，一直出現(xiàn)常數(shù)定義錯誤的問題！求大神幫助！
已經(jīng)有12人回復
WanderingHeart
鐵桿木蟲 (著名寫手)
- 應助: 90 (初中生)
- 金幣: 13890.8
- 散金: 5
- 紅花: 18
- 帖子: 1885
- 在線: 467.1小時
- 蟲號: 257464
- 注冊: 2006-06-04
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2015-03-30 12:59:27
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
3樓2015-03-30 13:53:41
shikang999
新蟲 (小有名氣)
- 計算強帖: 1
- 應助: 27 (小學生)
- 金幣: 1454.7
- 紅花: 7
- 帖子: 182
- 在線: 206.4小時
- 蟲號: 2201981
- 注冊: 2012-12-23
- 性別: GG
- 專業(yè): 固體力學
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
1、樓上dingd的回答可能有欠缺，4個方程4個未知數(shù)，一般情況下（不加任何服務約束）也就只有一組解，如果是線性滿秩方程組的話，這是沒有錯的。如果是非線性方程組的話，存在多組解（均滿足設定的誤差規(guī)則）的情況就顯得很正常了。下面給出一組在所設定的范圍內(nèi)找到的一組解，具體誤差是否滿足你的要求，可以自己測試看 x2=0.00018602746110 k1=-0.65443662426355 k2= -0.35539381460637 k3=0.10000000000929 2、用那個做初值進行牛頓迭代反而無解，只能說明你那個初值進去，采用那種牛頓迭代算法無法找到滿足條件的解而已。 |
4樓2015-03-30 16:57:05
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
5樓2015-03-30 19:21:52
6樓2015-03-31 13:59:21
|
特別謝謝你，你的說法也是我在想的，這組解的誤差完全滿足我的需求，唉，只是我的一些參數(shù)值又變了一下，所以方程組又變了，依然用1.5解，得到的誤差更大了。能否再幫我運行一下，謝謝了。 這里x2這個參數(shù)值最好在5.8*10^-5左右，目標函數(shù)誤差在0.1這個order的都可以接受。若有多個相似誤差的解，方便的話一塊貼上來。thanks again. Parameters x2[0.000056,0.000065],K1[-3.3,3.3],K2[-3.3,3.3],K3[1,9] ; Function -0.1176618358 *K1^2 - 0.08556995644 *K3^2 + K1* (55.94823395 + 0.2071600743* K3 - 1230.554191 *x2) - 3.217406066*10^6 *(-0.09245993440 + x2)* (0.001528166642 + x2) + K3 *(-36.61948488 + 1083.281150 x2)= 436.9; 456.3017378 - 0.0003509442655 *K1^2 - 0.00008773606637* K2^2 - 1.803399598 *K3 + 0.9178216086* K3^2 + 0.0002552253695* K3^3 - 0.00007226260461* K3^4 + K1 *(3.055537706 + 0.0003509442655* K2 + 0.0006178863318* K3 - 0.0003089431659* K3^2) + K2 *(-1.527768853 - 0.0003089431659* K3 + 0.0001672474282 *K3^2) + 2.446546643 *K3^2* x2= 451.7; 455.9695382 - 0.006734949843* K1^2 - 0.001683737461* K2^2 - 7.900236206* K3 + 4.016963725* K3^2 + 0.004898014389* K3^3 - 0.001386787206* K3^4 + K1* (13.38553565 + 0.006734949843* K2 + 0.01185781864 *K3 - 0.005928909321* K3^2) + K2 *(-6.692767825 - 0.005928909321* K3 + 0.003209635121* K3^2) + 10.71770027 *K3^2* x2 = 445.6; 455.4457565 - 0.03309661069* K1^2 - 0.008274152672* K2^2 - 17.51316559 *K3 + 8.891553762* K3^2 + 0.02406961880* K3^3 - 0.006814892071 *K3^4 + K1*(29.67292323 + 0.03309661069* K2 + 0.05827119969 *K3 - 0.02913559985* K3^2) + K2* (-14.83646161 - 0.02913559985* K3 + 0.01577265556* K3^2) + 23.75889211* K3^2 *x2 = 440.4 ; |
7樓2015-03-31 14:06:27
相關版塊跳轉(zhuǎn)
我要訂閱樓主 dollapym 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +4 | 晨昏線與星海 2026-03-19 | 5/250 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +6 | z1z2z3879 2026-03-18 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿福大288有機化學，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 東南大學364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|