版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

RK1224

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 0.5
帖子: 7
在線: 2.2小時(shí)
蟲號(hào): 1898291
注冊(cè): 2012-07-18

[求助] 非線性方程組的系數(shù)的數(shù)量級(jí)差很大，結(jié)果發(fā)散已有3人參與

小生最近做點(diǎn)東西，模型簡(jiǎn)化為非線性方程組，輸入一組參數(shù)后，發(fā)現(xiàn)方程的系數(shù)數(shù)量級(jí)相差很大，用牛頓迭代結(jié)果發(fā)散。請(qǐng)問(wèn)各位大神有遇到過(guò)這種情況嗎？有什么好的方法嗎？
簡(jiǎn)化后的方程代入一組參數(shù)如下：
(4.1018e-07)*x1^4+38716e*x1-38708*x2-(9.2333e-08)*x3^4-297.7296)=0;

38708*x1-38943*x2+1174.2e=0;

-23679*x2+47828e*x3-24149e*x4=0;

(9.2333e-08)*x1^4+23914e*x2-47834*x3-(9.2333e-08)*x3^4+23914*x4 -5.3514*x5+4613.8=0;

5.3514*x3-83.6514*x5+1566=0;

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

matlab編程繪圖

» 猜你喜歡

329求調(diào)劑已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有11人回復(fù)
一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一）已經(jīng)有3人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有14人回復(fù)
320求調(diào)劑0856 已經(jīng)有3人回復(fù)
307求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有14人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有18人回復(fù)

1樓 2015-08-10 19:50:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hsdwlxhx

金蟲 (正式寫手)

應(yīng)助: 20 (小學(xué)生)
金幣: 2986.8
散金: 204
紅花: 6
帖子: 430
在線: 318.5小時(shí)
蟲號(hào): 1537165
注冊(cè): 2011-12-14
性別: GG
專業(yè): 配位化學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

5維迭代，你也太bug了

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

BringLightonChip

2樓2015-08-10 20:09:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

RK1224

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 0.5
帖子: 7
在線: 2.2小時(shí)
蟲號(hào): 1898291
注冊(cè): 2012-07-18

引用回帖:

2樓: Originally posted by hsdwlxhx at 2015-08-10 20:09:31
5維迭代，你也太bug了

我是菜鳥。。。希望你指點(diǎn)一下

贊一下

回復(fù)此樓

3樓2015-08-10 20:43:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hsdwlxhx

金蟲 (正式寫手)

應(yīng)助: 20 (小學(xué)生)
金幣: 2986.8
散金: 204
紅花: 6
帖子: 430
在線: 318.5小時(shí)
蟲號(hào): 1537165
注冊(cè): 2011-12-14
性別: GG
專業(yè): 配位化學(xué)

【答案】應(yīng)助回帖

引用回帖:

3樓: Originally posted by RK1224 at 2015-08-10 20:43:33

我是菜鳥。。。希望你指點(diǎn)一下...

指點(diǎn)談不上，你的方程可以簡(jiǎn)化，先用數(shù)值方法跑一跑，看看有沒(méi)有解，大概在什么區(qū)間，然后再用迭代計(jì)算，matlab有fsolve函數(shù)可以直接用

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

BringLightonChip

4樓2015-08-10 22:02:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

1：第一個(gè)公式左右括號(hào)不匹配；
2：第2、第3公式中的“e”是什么？應(yīng)該是“e+0”吧！

贊一下

回復(fù)此樓

5樓2015-08-10 23:28:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

RK1224

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 0.5
帖子: 7
在線: 2.2小時(shí)
蟲號(hào): 1898291
注冊(cè): 2012-07-18

引用回帖:

5樓: Originally posted by dingd at 2015-08-10 23:28:57
1：第一個(gè)公式左右括號(hào)不匹配；
2：第2、第3公式中的“e”是什么？應(yīng)該是“e+0”吧！

謝謝你的回復(fù)，確實(shí)是我沒(méi)寫清楚，我改了一下如下：
(4.1018e-07)*x1^4+38716*x1-38708*x2-(9.2333e-08)*x3^4-297.7296=0;

38708*x1-38943*x2+1174.2=0;

-23679*x2+47828*x3-24149*x4=0;

(9.2333e-08)*x1^4+23914*x2-47834*x3-(9.2333e-08)*x3^4+23914*x4 -5.3514*x5+4613.8=0;

5.3514*x3-83.6514*x5+1566=0;
請(qǐng)問(wèn)你遇到過(guò)這種問(wèn)題么？主要的問(wèn)題就是矩陣奇異

贊一下

回復(fù)此樓

6樓2015-08-11 13:32:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

引用回帖:

6樓: Originally posted by RK1224 at 2015-08-11 13:32:07
謝謝你的回復(fù)，確實(shí)是我沒(méi)寫清楚，我改了一下如下：
(4.1018e-07)*x1^4+38716*x1-38708*x2-(9.2333e-08)*x3^4-297.7296=0;

38708*x1-38943*x2+1174.2=0;

-23679*x2+47828*x3-24149*x4=0;

(9.2333e-08)*x ...

解：
x1: 6.06356999532981
x2: 6.05713138122965
x3: 15.5431308164225
x4: 24.8445089532369
x5: 19.7148823600203

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

7樓2015-08-11 14:29:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

引用回帖:

MATLAB代碼：

CODE:

function solve_eq

clear all;clc

format long

x0=[0 0 0 0 0];

options=optimset('TolFun',1e-9);

[x fval]=fsolve(@fun,x0,options);

[x;fval]'

function f=fun(u)

x1=u(1);

x2=u(2);

x3=u(3);

x4=u(4);

x5=u(5);

f(1)=(4.1018e-07)*x1^4+38716*x1-38708*x2-(9.2333e-08)*x3^4-297.7296;

f(2)=38708*x1-38943*x2+1174.2;

f(3)=-23679*x2+47828*x3-24149*x4;

f(4)=(9.2333e-08)*x1^4+23914*x2-47834*x3-(9.2333e-08)*x3^4+23914*x4 -5.3514*x5+4613.8;

f(5)=5.3514*x3-83.6514*x5+1566;

計(jì)算結(jié)果：
         x1~x5                      方程殘差
6.063569995329909    -0.000000000031662
6.057131381229750    -0.000000000011596
  15.543130816422794    -0.000000000116415
  24.844508953237408    0.000000000042746
  19.714882360020336  -0.000000000000227

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

8樓2015-08-11 14:34:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 RK1224 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	5/250	2026-03-19 18:03 by zcl123
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +8	薛云鵬 2026-03-13	8/400	2026-03-19 15:36 by haoshis
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 0817調(diào)劑 +3	沒(méi)有答案_ 2026-03-14	3/150	2026-03-19 09:51 by Xu de nuo
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無(wú)際的草原
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7	手機(jī)用戶 2026-03-17	7/350	2026-03-18 14:34 by vgtyfty
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu

24小時(shí)熱門版塊排行榜

[求助] 非線性方程組的系數(shù)的數(shù)量級(jí)差很大，結(jié)果發(fā)散 已有3人參與

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] 非線性方程組的系數(shù)的數(shù)量級(jí)差很大，結(jié)果發(fā)散已有3人參與