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

crazyiyy

木蟲 (小有名氣)

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

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

本人有個(gè)非線性方程組，已用1stOpt解出，但需要將結(jié)果導(dǎo)入到Matlab中，十分麻煩。因此想用Matlab將此方程組解出來(lái)，肯定路過(guò)大神幫忙！題目如下
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解出來(lái)，感激不盡！

回復(fù)此樓

» 猜你喜歡

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

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

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求解非線性方程組說(shuō)無(wú)解，一定是方程組本身無(wú)解，還是有可能程序有問(wèn)題呢？已經(jīng)有11人回復(fù)
matlab如何求解一個(gè)非線性微分方程組已經(jīng)有8人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
求Matlab解非線性方程代碼已經(jīng)有12人回復(fù)
matlab解非線性方程組解不開(kāi) 已經(jīng)有3人回復(fù)
求助matlab 非線性方程組全部解！已經(jīng)有6人回復(fù)
Matlab：Newton-Raphson 解線非線性性方程組的問(wèn)題已經(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ù)
請(qǐng)教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的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(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
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

以下程序第對(duì)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小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

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

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

月只藍(lán)

主管區(qū)長(zhǎng) (職業(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
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
crazyiyy: 金幣+100, ★★★★★最佳答案, 衷心地感謝您對(duì)我的幫助！ 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小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

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

crazyiyy

木蟲 (小有名氣)

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

送紅花一朵

引用回帖:

3樓: Originally posted by 月只藍(lán) at 2014-07-14 16:14:10
以下程序第對(duì)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è)問(wèn)題希望您能為我解惑。如果我需要對(duì)方程組中的某些系數(shù)進(jìn)行變化，進(jìn)而得到一組解，我應(yīng)該如何操作。
例如我需要對(duì)方程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的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(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
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

引用回帖:

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

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

贊一下(1人)

回復(fù)此樓

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

crazyiyy

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

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

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

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

[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +7	楓橋ZL 2026-03-18	9/450	2026-03-21 09:56 by JineShine
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +6	墨墨漠 2026-03-20	6/300	2026-03-21 09:23 by 天青等煙雨
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 354求調(diào)劑 +5	Tyoumou 2026-03-18	8/400	2026-03-21 00:35 by JourneyLucky
[考研] 299求調(diào)劑 +4	shxchem 2026-03-20	5/250	2026-03-21 00:24 by 星空星月
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 319求調(diào)劑 +3	小力氣珂珂 2026-03-20	3/150	2026-03-20 19:47 by JourneyLucky
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong