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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

crazyiyy

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1696.5
紅花: 1
帖子: 78
在線: 85.6小時
蟲號: 2483453
注冊: 2013-05-27
性別: GG
專業(yè): 通信理論與系統(tǒng)

[求助] MATLAB解非線性方程組已有1人參與

本人有個非線性方程組，已用1stOpt解出，但需要將結(jié)果導(dǎo)入到Matlab中，十分麻煩。因此想用Matlab將此方程組解出來，肯定路過大神幫忙！題目如下
1/((4.5+(8.5)*(x+y-x*y)+(33/2)*(x+y-x*y)^2+(33/2)*(x+y-x*y)^3+(33/2)*(x+y-x*y)^4)*(1-(x+y-x*y))/(1-(x+y-x*y)^5)+(1-(12/(48*2^3-4)))*(1-x)+6*(12/(48*2^3-4))+(1-(12/(48*2^3-4)))*(1-x)*(12*(1-y)+1))-z=0;
(13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)=x;
(1-(1-z)^8)/(2-(1-z)^9)=y;
注：其中x,y,z都是[0,1]數(shù)，一定需要用matlab解出來，感激不盡！

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有144人回復(fù)
哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學工程學碩調(diào)劑信息已經(jīng)有26人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有34人回復(fù)

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

matlab求解非線性方程組已經(jīng)有3人回復(fù)
matlab 解非線性偏微分方程組已經(jīng)有6人回復(fù)
MATLAB求解五元非線性方程組，在線等已經(jīng)有7人回復(fù)
非線性二階微分方程組求解，matlab 已經(jīng)有12人回復(fù)
求助解復(fù)雜非線性方程組的好的方法已經(jīng)有24人回復(fù)
matlab非線性方程組該如何編程。求助已經(jīng)有10人回復(fù)
用matlab解非線性方程組已經(jīng)有4人回復(fù)
matlab 解超定非線性方程組已經(jīng)有8人回復(fù)
用matlab求解一個非線性方程組的解已經(jīng)有4人回復(fù)
求助用matlab或者（1st0pt)編程解方程組已經(jīng)有3人回復(fù)
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復(fù)
matlab如何求解一個非線性微分方程組已經(jīng)有8人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
求Matlab解非線性方程代碼已經(jīng)有12人回復(fù)
matlab解非線性方程組解不開已經(jīng)有3人回復(fù)
求助matlab 非線性方程組全部解！已經(jīng)有6人回復(fù)
Matlab：Newton-Raphson 解線非線性性方程組的問題已經(jīng)有12人回復(fù)
matlab求解非線性方程組，錯誤提示怎么解決已經(jīng)有5人回復(fù)
求助 MATLAB解方程組-fslove 已經(jīng)有7人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
matlab求解非線性方程組，求助！已經(jīng)有6人回復(fù)
如何使用matlab求解非線性方程組的所有整數(shù)解？已經(jīng)有9人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
請教Matlab高手解多元非線性方程組中的參數(shù) 已經(jīng)有6人回復(fù)
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！已經(jīng)有18人回復(fù)

用心做好每一件事

1樓 2014-07-14 15:43:54

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

月只藍

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

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
crazyiyy: 金幣+100, ★★★★★最佳答案, 衷心地感謝您對我的幫助！ 2014-07-14 16:20:17

CODE:

function jie_fun

u0=[0.5 0.5 0.5];

u=fsolve(@myfun,u0)

f_val=myfun(u)

function f=myfun(u)

x=u(1);

y=u(2);

z=u(3);

f(1)=1/((4.5+(8.5)*(x+y-x*y)+(33/2)*(x+y-x*y)^2+(33/2)*(x+y-x*y)^3+(33/2)*...

    (x+y-x*y)^4)*(1-(x+y-x*y))/(1-(x+y-x*y)^5)+(1-(12/(48*2^3-4)))*(1-x)+6*...

    (12/(48*2^3-4))+(1-(12/(48*2^3-4)))*(1-x)*(12*(1-y)+1))-z;

f(2)= (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)-x;

f(3)=(1-(1-z)^8)/(2-(1-z)^9)-y;

計算結(jié)果：
u =

0.8138 0.3049 0.0733

目標函數(shù)值f_val =

1.0e-008 *

0.0008 0.1291 -0.0216

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2014-07-14 16:03:41

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

月只藍

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

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

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

以下程序第對2樓中程序的改進，提高解的小數(shù)點后的位數(shù)。

CODE:

function jie_fun

clear a;clc

format long

u0=[0.5 0.5 0.5];

u=fsolve(@myfun,u0)

f_val=myfun(u)

function f=myfun(u)

x=u(1);

y=u(2);

z=u(3);

f(1)=1/((4.5+(8.5)*(x+y-x*y)+(33/2)*(x+y-x*y)^2+(33/2)*(x+y-x*y)^3+(33/2)*...

    (x+y-x*y)^4)*(1-(x+y-x*y))/(1-(x+y-x*y)^5)+(1-(12/(48*2^3-4)))*(1-x)+6*...

    (12/(48*2^3-4))+(1-(12/(48*2^3-4)))*(1-x)*(12*(1-y)+1))-z;

f(2)= (13*(1-(1-z)^8)+9*z*(1-z)^8)*(1-x)*(1-y)-x;

f(3)=(1-(1-z)^8)/(2-(1-z)^9)-y;

計算結(jié)果：
u =

0.813803728151198 0.304867313866725 0.073289552279229

f_val =

1.0e-008 *

0.000793234922192 0.129125732328816 -0.021619173118381

已驗證，上述結(jié)果與1stopt結(jié)果一致。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

crazyiyy

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

3樓2014-07-14 16:14:10

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

crazyiyy

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1696.5
紅花: 1
帖子: 78
在線: 85.6小時
蟲號: 2483453
注冊: 2013-05-27
性別: GG
專業(yè): 通信理論與系統(tǒng)

送紅花一朵

引用回帖:

3樓: Originally posted by 月只藍 at 2014-07-14 16:14:10
以下程序第對2樓中程序的改進，提高解的小數(shù)點后的位數(shù)。

function jie_fun
clear a;clc
format long
u0=;
u=fsolve(@myfun,u0)
f_val=myfun(u)

function f=myfun(u)
x=u(1);
y=u(2);
z=u(3);
...

您好！您的答案的確跟1stOpt一致，還有個問題希望您能為我解惑。如果我需要對方程組中的某些系數(shù)進行變化，進而得到一組解，我應(yīng)該如何操作。
例如我需要對方程1中的48*2^3 變?yōu)?8*2^x 其中x=1,2,3,4 這樣我應(yīng)該得到4組解，

贊一下

回復(fù)此樓

用心做好每一件事

4樓2014-07-14 16:39:24

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

月只藍

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

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

引用回帖:

4樓: Originally posted by crazyiyy at 2014-07-14 16:39:24
您好！您的答案的確跟1stOpt一致，還有個問題希望您能為我解惑。如果我需要對方程組中的某些系數(shù)進行變化，進而得到一組解，我應(yīng)該如何操作。
例如我需要對方程1中的48*2^3 變?yōu)?8*2^x 其中x=1,2,3,4 這樣我應(yīng) ...

我的處理辦法是把x定義為全局變量。
以前已經(jīng)有過相關(guān)代碼了，參見：
http://www.gaoyang168.com/bbs/viewthread.php?tid=6267100&authorid=1122189

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

crazyiyy

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

5樓2014-07-14 16:45:00

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

月只藍

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

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

引用回帖:

以及：
http://www.gaoyang168.com/bbs/viewthread.php?tid=6087640

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

6樓2014-07-14 16:47:27

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

crazyiyy

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1696.5
紅花: 1
帖子: 78
在線: 85.6小時
蟲號: 2483453
注冊: 2013-05-27
性別: GG
專業(yè): 通信理論與系統(tǒng)

送紅花一朵

引用回帖:

5樓: Originally posted by 月只藍 at 2014-07-14 16:45:00
我的處理辦法是把x定義為全局變量。
以前已經(jīng)有過相關(guān)代碼了，參見：
http://www.gaoyang168.com/bbs/viewthread.php?tid=6267100&authorid=1122189...

您好！我還有最后一個問題，就是matlab種一定要給出初值，請問初值的調(diào)試有沒有什么方法或者經(jīng)驗啊，或者說如何判斷當前初值得到的解是正解。

贊一下

回復(fù)此樓

用心做好每一件事

7樓2014-07-14 17:06:19

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

月只藍

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

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

引用回帖:

7樓: Originally posted by crazyiyy at 2014-07-14 17:06:19
您好！我還有最后一個問題，就是matlab種一定要給出初值，請問初值的調(diào)試有沒有什么方法或者經(jīng)驗啊，或者說如何判斷當前初值得到的解是正解。...

把算出來的解代回原方程，看能有多接近于0，3樓正是給出了這些數(shù)值。
此外，對于多解的情況，最終解的如果有物理意義要求，舉例說，比如未知數(shù)是開爾文溫度T，是不可能小于0的，在解里面就應(yīng)當取滿足實際物理意義的。
在直接調(diào)用fsolve的情況下，初值基本靠自己嘗試，大概判斷一下。
究其本質(zhì)而言，解方程組問題是最優(yōu)化問題，調(diào)用具有全局優(yōu)化能力的算法，比如說遺傳算法，可以降低對初值的要求。
1stopt軟件解決最優(yōu)化問題的強大之處其一在于在操作者層面上不需要賦予初值，即可獲得相當優(yōu)秀的結(jié)果。

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

8樓2014-07-14 18:37:07

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

mjjackey

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46
帖子: 9
在線: 4小時
蟲號: 3696153
注冊: 2015-02-25
性別: GG
專業(yè): 模式識別

引用回帖:

8樓: Originally posted by 月只藍 at 2014-07-14 18:37:07
把算出來的解代回原方程，看能有多接近于0，3樓正是給出了這些數(shù)值。
此外，對于多解的情況，最終解的如果有物理意義要求，舉例說，比如未知數(shù)是開爾文溫度T，是不可能小于0的，在解里面就應(yīng)當取滿足實際物理意義 ...

用fsolve，貌似也只能設(shè)置三種算法，沒辦法設(shè)置更多的算法“調(diào)用具有全局優(yōu)化能力的算法，比如說遺傳算法，”

贊一下

回復(fù)此樓

9樓2015-03-02 16:40:35

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

月只藍

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

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

引用回帖:

9樓: Originally posted by mjjackey at 2015-03-02 16:40:35
用fsolve，貌似也只能設(shè)置三種算法，沒辦法設(shè)置更多的算法“調(diào)用具有全局優(yōu)化能力的算法，比如說遺傳算法，”...

“究其本質(zhì)而言，解方程組問題是最優(yōu)化問題，調(diào)用具有全局優(yōu)化能力的算法，比如說遺傳算法，可以降低對初值的要求。
”
先將原問題轉(zhuǎn)化為最優(yōu)化問題，調(diào)用遺傳算法算出初值，該初值代入fsolve函數(shù)。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

10樓2015-03-02 17:20:20

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

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

返回列表

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

[考研] 302求調(diào)劑 +6	呼呼呼。。。。 2026-03-17	6/300	2026-03-17 23:56 by 星空星月
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3	Eember. 2026-03-17	7/350	2026-03-17 20:20 by 花125533
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風里
[考研] 工科材料085601 279求調(diào)劑 +3	困于星晨 2026-03-17	3/150	2026-03-17 14:08 by ms629
[考研] 283求調(diào)劑 +3	聽風就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 274求調(diào)劑 +5	時間點 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] [導(dǎo)師推薦]西南科技大學國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 一志愿985，本科211，0817化學工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +5	d如愿上岸 2026-03-12	8/400	2026-03-16 15:19 by peike
[考研] 277材料科學與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考博] 東華理工大學化材專業(yè)26屆碩士博士申請 +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6	步川酷紫123 2026-03-13	6/300	2026-03-13 21:59 by 星空星月
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao

24小時熱門版塊排行榜

[求助] MATLAB解非線性方程組 已有1人參與

» 猜你喜歡

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

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

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

» 本帖已獲得的紅花（最新10朵）

» 本帖已獲得的紅花（最新10朵）

[求助] MATLAB解非線性方程組已有1人參與