版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

>論壇更新日志 (6620)
>導師招生 (1380)
>考研 (1280)
>蟲友互識 (588)
>文獻求助 (479)
>休閑灌水 (243)
>基金申請 (189)
>考博 (185)
>招聘信息布告欄 (145)
>碩博家園 (137)
>論文投稿 (117)
>博后之家 (96)
>綠色求助(高懸賞) (46)
>教師之家 (41)
>SciFinder/Reaxys (39)
>找工作 (35)

返回列表

crazyiyy

木蟲 (小有名氣)

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

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

本人有個非線性方程組，已用1stOpt解出，但需要將結果導入到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]數，一定需要用matlab解出來，感激不盡！

回復此樓

» 猜你喜歡

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

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

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

用心做好每一件事

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

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

月只藍

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

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

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應助指數 +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;

計算結果：
u =

0.8138 0.3049 0.0733

目標函數值f_val =

1.0e-008 *

0.0008 0.1291 -0.0216

贊一下

回復此樓

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

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

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

月只藍

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

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

【答案】應助回帖

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

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;

計算結果：
u =

0.813803728151198 0.304867313866725 0.073289552279229

f_val =

1.0e-008 *

0.000793234922192 0.129125732328816 -0.021619173118381

已驗證，上述結果與1stopt結果一致。

贊一下

回復此樓

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

crazyiyy

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

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

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

crazyiyy

木蟲 (小有名氣)

應助: 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樓中程序的改進，提高解的小數點后的位數。

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一致，還有個問題希望您能為我解惑。如果我需要對方程組中的某些系數進行變化，進而得到一組解，我應該如何操作。
例如我需要對方程1中的48*2^3 變?yōu)?8*2^x 其中x=1,2,3,4 這樣我應該得到4組解，

贊一下

回復此樓

用心做好每一件事

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

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

月只藍

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

專家經驗: +1059
計算強帖: 8
應助: 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一致，還有個問題希望您能為我解惑。如果我需要對方程組中的某些系數進行變化，進而得到一組解，我應該如何操作。
例如我需要對方程1中的48*2^3 變?yōu)?8*2^x 其中x=1,2,3,4 這樣我應 ...

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

贊一下(1人)

回復此樓

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

crazyiyy

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

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

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

月只藍

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

專家經驗: +1059
計算強帖: 8
應助: 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

贊一下

回復此樓

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

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

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

crazyiyy

木蟲 (小有名氣)

應助: 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定義為全局變量。
以前已經有過相關代碼了，參見：
http://www.gaoyang168.com/bbs/viewthread.php?tid=6267100&authorid=1122189...

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

贊一下

回復此樓

用心做好每一件事

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

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

月只藍

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

專家經驗: +1059
計算強帖: 8
應助: 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種一定要給出初值，請問初值的調試有沒有什么方法或者經驗啊，或者說如何判斷當前初值得到的解是正解。...

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

贊一下(1人)

回復此樓

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

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

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

mjjackey

新蟲 (初入文壇)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

月只藍

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

專家經驗: +1059
計算強帖: 8
應助: 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，貌似也只能設置三種算法，沒辦法設置更多的算法“調用具有全局優(yōu)化能力的算法，比如說遺傳算法，”...

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

贊一下

回復此樓

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

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

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

相關版塊跳轉我要訂閱樓主 crazyiyy 的主題更新

返回列表

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

[考研] 302求調劑 +5	呼呼呼。。。。 2026-03-17	5/250	2026-03-17 22:54 by lbsjt
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調劑 +6	yangfz 2026-03-17	6/300	2026-03-17 21:30 by hmn_wj
[考研] 341求調劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 070300化學319求調劑 +4	錦鯉0909 2026-03-17	4/200	2026-03-17 18:21 by 重科小霸王
[考研] 本人考085602 化學工程專碩 +16	不知道叫什么！ 2026-03-15	18/900	2026-03-17 17:05 by ruiyingmiao
[考研] 290求調劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 工科材料085601 279求調劑 +3	困于星晨 2026-03-17	3/150	2026-03-17 14:08 by ms629
[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[考研] 275求調劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 311求調劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 070300化學學碩求調劑 +6	太想進步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 0703化學調劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 22408總分284求調劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 中科大材料專碩319求調劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 招收0805（材料）調劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 材料371求調劑 +9	鱷魚? 2026-03-11	11/550	2026-03-13 22:53 by JourneyLucky
[考研] 求材料調劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 一志愿211化學學碩310分求調劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 295求調劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 283求調劑，材料、化工皆可 +8	蘇打水7777 2026-03-11	10/500	2026-03-13 09:06 by Linda Hu

24小時熱門版塊排行榜

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

» 猜你喜歡

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

【答案】應助回帖

【答案】應助回帖

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

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

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