當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

crazyiyy

木蟲 (小有名氣)

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

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

本人有個(gè)非線性方程組，已用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)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有116人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有79人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)

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

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求解一個(gè)非線性方程組的解已經(jīng)有4人回復(fù)
求助用matlab或者（1st0pt)編程解方程組已經(jīng)有3人回復(fù)
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復(fù)
matlab如何求解一個(gè)非線性微分方程組已經(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求解非線性方程組，錯(cuò)誤提示怎么解決已經(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ù)
求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
請教Matlab高手解多元非線性方程組中的參數(shù) 已經(jīng)有6人回復(fù)
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實(shí)數(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小時(shí)
蟲號(hào): 2483453
注冊: 2013-05-27
性別: GG
專業(yè): 通信理論與系統(tǒng)

送紅花一朵

引用回帖:

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

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

贊一下

回復(fù)此樓

用心做好每一件事

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

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

查看全部 10 個(gè)回答

月只藍(lán)

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

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

計(jì)算結(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ā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

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

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

月只藍(lán)

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

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

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

以下程序第對2樓中程序的改進(jìn)，提高解的小數(shù)點(diǎn)后的位數(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;

計(jì)算結(jié)果：
u =

0.813803728151198 0.304867313866725 0.073289552279229

f_val =

1.0e-008 *

0.000793234922192 0.129125732328816 -0.021619173118381

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

贊一下

回復(fù)此樓

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

crazyiyy

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

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小時(shí)
蟲號(hào): 2483453
注冊: 2013-05-27
性別: GG
專業(yè): 通信理論與系統(tǒng)

送紅花一朵

引用回帖:

3樓: Originally posted by 月只藍(lán) at 2014-07-14 16:14:10
以下程序第對2樓中程序的改進(jìn)，提高解的小數(shù)點(diǎn)后的位數(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一致，還有個(gè)問題希望您能為我解惑。如果我需要對方程組中的某些系數(shù)進(jìn)行變化，進(jìn)而得到一組解，我應(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的回帖

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

[考研] 299求調(diào)劑 +6	某某某某位 2026-03-21	6/300	2026-03-24 12:22 by 544594351
[考研] 279分求調(diào)劑一志愿211 +18	chaojifeixia 2026-03-19	20/1000	2026-03-24 10:34 by dolphin_ycj
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +5	niko- 2026-03-22	5/250	2026-03-23 22:01 by fuyu_
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-23	4/200	2026-03-23 21:25 by pswait
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +3	WW.' 2026-03-23	5/250	2026-03-23 17:18 by 楓翼ljj
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 306求調(diào)劑 +9	chuanzhu川燭 2026-03-18	9/450	2026-03-23 13:17 by luoyongfeng
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 311求調(diào)劑 +6	冬十三 2026-03-18	6/300	2026-03-22 20:18 by edmund7
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-20	3/150	2026-03-22 16:00 by ColorlessPI
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 求調(diào)劑 +7	Auroracx 2026-03-22	7/350	2026-03-22 12:38 by 素顏傾城1988
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne