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ù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有109人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(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的回帖

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的回帖

查看全部 10 個回答

月只藍

主管區(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è): 宇宙學(xué)
管轄: 計算模擬區(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

目標(biāo)函數(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è): 宇宙學(xué)
管轄: 計算模擬區(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的回帖

月只藍

主管區(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è): 宇宙學(xué)
管轄: 計算模擬區(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的回帖

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

[考研] 本人考085602 化學(xué)工程專碩 +21	不知道叫什么！ 2026-03-15	23/1150	2026-03-22 00:12 by BruceLiu320
[考研] 269專碩求調(diào)劑 +5	金恩貝 2026-03-21	5/250	2026-03-21 22:37 by zhyzzh
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[基金申請] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 0856材料專碩353求調(diào)劑 +3	NIFFFfff 2026-03-20	3/150	2026-03-21 10:23 by luoyongfeng
[考研] 316求調(diào)劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 086500 325 求調(diào)劑 +3	領(lǐng)帶小熊 2026-03-19	3/150	2026-03-20 18:38 by 盡舜堯1
[考研] 08工學(xué)調(diào)劑 +5	用戶573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機用戶 2026-03-17	9/450	2026-03-20 14:24 by 無懈可擊111
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82