774099895

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 87.9
帖子: 152
在線: 112.4小時
蟲號: 1482396
注冊: 2011-11-08
性別: GG
專業(yè): 能源化工

[求助] 關(guān)于非線性方程組的問題已有2人參與

大神，麻煩你們幫我看一下這個非線性方程組怎么求解，我自己編的程，但是計算出來后始終等于初值，請指出問題所在。。。。。。
function NOX
x0=[0 0]'; x=fsolve(@concentration,x0)
function f=concentration(x)
a1=5.28E-05;
a2=-6.22E-06;
a3=-1.05E-05;
b1=2.96E-05;
b2=-3.79E-06;
b3=-3.79E-06;
b4=-1.05E-05;
b5=-8.81E-06;
f(1)=a1+a2*x(1)+a3*x(1)*x(2);
f(2)=b1+b2*x(2)+b3*x(2)^(3/2)+b4*x(1)*x(2)+b5*x(2)^2;

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有4人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有4人回復(fù)
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？已經(jīng)有9人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有7人回復(fù)
286求調(diào)劑已經(jīng)有6人回復(fù)
280求調(diào)劑已經(jīng)有10人回復(fù)
材料學(xué)碩301分求調(diào)劑已經(jīng)有4人回復(fù)
278求調(diào)劑已經(jīng)有8人回復(fù)
求助已經(jīng)有5人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

請教一個matlab求解非線性方程組的問題已經(jīng)有9人回復(fù)
用matlab解非線性方程組編程已經(jīng)有2人回復(fù)
origin 非線性自定義擬合問題~~方程組編寫及運行~~ 已經(jīng)有8人回復(fù)
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
【求助】非線性方程組的求解問題已經(jīng)有6人回復(fù)
擬牛頓法求解復(fù)雜（帶根號的）非線性方程組已經(jīng)有17人回復(fù)
非線性方程組問題已經(jīng)有4人回復(fù)
matlab求解非線性方程組，錯誤提示怎么解決已經(jīng)有5人回復(fù)
急切求助！一個mathematica求解一個簡單一元非線性方程的問題已經(jīng)有4人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
關(guān)于matlab線性規(guī)劃背包問題的解法求救呀~~方程知道但不知道如何解答···· 已經(jīng)有4人回復(fù)
matlab求解非線性方程組，求助！已經(jīng)有6人回復(fù)
非線性方程組matlab擬合已經(jīng)有9人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
fslove解非線性方程組解不出來已經(jīng)有3人回復(fù)
【求助】用mathematica 5.0求解一個非線性方程組失敗，特發(fā)帖求助！已經(jīng)有5人回復(fù)
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！已經(jīng)有18人回復(fù)
【求助】求一個用牛頓迭代法解非線性方程組的程序已經(jīng)有10人回復(fù)

世界上有木有另一個我，過著我想過的生活

1樓 2013-12-18 20:16:17

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

hppdyx

木蟲 (知名作家)

應(yīng)助: 60 (初中生)
金幣: 3242.5
散金: 1145
紅花: 48
帖子: 6010
在線: 297.6小時
蟲號: 2428506
注冊: 2013-04-21
性別: GG
專業(yè): 巖土與基礎(chǔ)工程

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
774099895: 金幣+2, ★★★很有幫助, 謝謝 2013-12-19 09:26:12

我試過，應(yīng)該是你的torlerance太小了，默認(rèn)是1e-4，但是你的參數(shù)都太小。所以可以把torlerance減小。我試過,結(jié)果是x= 4.5759和0.4939，不知道對不對。

options = optimset('TolFun', 1e-10);
x=fsolve(@concentration,x0, options)

贊一下(1人)

回復(fù)此樓

不以風(fēng)騷驚天下，但求淫蕩動世人

2樓2013-12-19 09:21:58

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

774099895

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 87.9
帖子: 152
在線: 112.4小時
蟲號: 1482396
注冊: 2011-11-08
性別: GG
專業(yè): 能源化工

引用回帖:

2樓: Originally posted by hppdyx at 2013-12-19 09:21:58
我試過，應(yīng)該是你的torlerance太小了，默認(rèn)是1e-4，但是你的參數(shù)都太小。所以可以把torlerance減小。我試過,結(jié)果是x= 4.5759和0.4939，不知道對不對。

options = optimset('TolFun', 1e-10);
x=fsolve(@concen ...

大神你好，我這個方程組初值是不知道的，我是隨便寫了一個初值，請問大神如果在不知道初值的前提下能不能解得這樣的方程，如果可以麻煩您能賜教一下，，，，，，，

贊一下

回復(fù)此樓

世界上有木有另一個我，過著我想過的生活

3樓2013-12-19 09:25:45

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

dingd

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

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

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

★
感謝參與，應(yīng)助指數(shù) +1
774099895: 金幣+1, ★有幫助 2013-12-19 09:58:44

試試1stOpt，不需初值，簡單易用：

CODE:

Constant a1=5.28E-05, a2=-6.22E-06, a3=-1.05E-05, b1=2.96E-05, b2=-3.79E-06,

         b3=-3.79E-06, b4=-1.05E-05, b5=-8.81E-06;

Function a1+a2*x1+a3*x1*x2;

         b1+b2*x2+b3*x2^(3/2)+b4*x1*x2+b5*x2^2;

x1: 4.60379123777867
x2: 0.499886522153301

贊一下(1人)

回復(fù)此樓

4樓2013-12-19 09:45:09

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

774099895

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 87.9
帖子: 152
在線: 112.4小時
蟲號: 1482396
注冊: 2011-11-08
性別: GG
專業(yè): 能源化工

引用回帖:

4樓: Originally posted by dingd at 2013-12-19 09:45:09
試試1stOpt，不需初值，簡單易用：

Constant a1=5.28E-05, a2=-6.22E-06, a3=-1.05E-05, b1=2.96E-05, b2=-3.79E-06,
b3=-3.79E-06, b4=-1.05E-05, b5=-8.81E-06;
Function a1+a2*x1+a3*x1*x2;
...

大哥，我?guī)Щ胤匠讨序炞C了一下，結(jié)果不等于零啊，不知道是我的問題沒表達(dá)清楚，我重新給您說明一下吧，我是要解
a1+a2*x1+a3*x1*x2=0;
b1+b2*x2+b3*x2^(3/2)+b4*x1*x2+b5*x2^2=0;
這兩個方程組，我是新手，對matlab不是很了解，希望大哥能教小弟一下，怎么解這兩個方程，謝謝您了，，，，，，

贊一下

回復(fù)此樓

世界上有木有另一個我，過著我想過的生活

5樓2013-12-19 09:58:31

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

dingd

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

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

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

★ ★
774099895: 金幣+2, ★★★★★最佳答案 2013-12-19 10:13:54

驗證對嗎？我這邊都接近于0�。�

x1: 4.60379123780314
x2: 0.499886522148953

5.28e-5+(-6.22e-6)*x1+(-1.05e-5)*x1*x2=-7.049685813E-17
2.96e-5+(-3.79e-6)*x2+(-3.79e-6)*x2^(3/2)+(-1.05e-5)*x1*x2+(-8.81e-6)*x2^2=1.539867782E-16

贊一下(1人)

回復(fù)此樓

6樓2013-12-19 10:10:09

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

774099895

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 87.9
帖子: 152
在線: 112.4小時
蟲號: 1482396
注冊: 2011-11-08
性別: GG
專業(yè): 能源化工

引用回帖:

6樓: Originally posted by dingd at 2013-12-19 10:10:09
驗證對嗎？我這邊都接近于0�。�

x1: 4.60379123780314
x2: 0.499886522148953

5.28e-5+(-6.22e-6)*x1+(-1.05e-5)*x1*x2=-7.049685813E-17
2.96e-5+(-3.79e-6)*x2+(-3.79e-6)*x2^(3/2)+(-1.05e-5)*x1*x2+(-8 ...

小弟太笨了，不太了解，真是麻煩您了，謝謝您啊，能不能順便幫我看一下下面這個編程，為啥我運行后就什么也沒有，這樣編程是不是不對，希望您有時間給予指導(dǎo)，如果方便真心希望
function non
a1=5.28E-05;
a2=-6.22E-06;
a3=-1.05E-05;
b1=2.96E-05;
b2=-3.79E-06;
b3=-3.79E-06;
b4=-1.05E-05;
b5=-8.81E-06;
[x y]=solve('a1+a2*x+a3*x*y=0','b1+b2*y+b3*y^(3/2)+b4*x*y+b5*y^2=0')

贊一下

回復(fù)此樓

世界上有木有另一個我，過著我想過的生活

7樓2013-12-19 10:16:57

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

dingd

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

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

Matlab數(shù)值解的話用fsolve。Solve是求解析解，不一定有。

贊一下

回復(fù)此樓

8樓2013-12-19 10:39:00

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

hppdyx

木蟲 (知名作家)

應(yīng)助: 60 (初中生)
金幣: 3242.5
散金: 1145
紅花: 48
帖子: 6010
在線: 297.6小時
蟲號: 2428506
注冊: 2013-04-21
性別: GG
專業(yè): 巖土與基礎(chǔ)工程

引用回帖:

3樓: Originally posted by 774099895 at 2013-12-19 09:25:45
大神你好，我這個方程組初值是不知道的，我是隨便寫了一個初值，請問大神如果在不知道初值的前提下能不能解得這樣的方程，如果可以麻煩您能賜教一下，，，，，，，...

可以檢驗一下是否接近與零，還有就是可以自己寫程序求解，不用內(nèi)置函數(shù)，來檢驗一下。。。

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

贊一下(2人)

回復(fù)此樓

不以風(fēng)騷驚天下，但求淫蕩動世人

9樓2013-12-19 11:10:59

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

萌兮萌兮___

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
帖子: 1
在線: 2.4小時
蟲號: 4353111
注冊: 2016-01-13

引用回帖:

你好，我按照你的方法改了但是還是運行結(jié)果和初值始終一樣呀

贊一下

回復(fù)此樓

10樓2016-01-15 01:08:28

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

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

[考研] 求調(diào)劑 +4	十三加油 2026-03-21	4/200	2026-03-21 16:20 by barlinike
[考研] 材料學(xué)碩301分求調(diào)劑 +4	Liyouyumairs 2026-03-21	4/200	2026-03-21 16:01 by 18337163304
[考研] 299求調(diào)劑 +4	shxchem 2026-03-20	6/300	2026-03-21 11:54 by 無際的草原
[考研] 0856材料專碩353求調(diào)劑 +3	NIFFFfff 2026-03-20	3/150	2026-03-21 10:23 by luoyongfeng
[考研] 材料 271求調(diào)劑 +3	展信悅_ 2026-03-21	3/150	2026-03-21 10:17 by 王笑宇一定是研?/a>
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-17	8/400	2026-03-21 02:12 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 北科281學(xué)碩材料求調(diào)劑 +5	tcxiaoxx 2026-03-20	5/250	2026-03-20 21:35 by laoshidan
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 生物學(xué)調(diào)劑招人�。�！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm