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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

monixia

木蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 2535.6
紅花: 3
帖子: 80
在線: 24.2小時
蟲號: 1817862
注冊: 2012-05-15
性別: GG
專業(yè): 化學反應工程

[求助] 求解決matlab解非線性方程組出現(xiàn)的問題已有1人參與

各位matlab大神，我計算的方程組如下，
>> m=0.00040;
f=@(x)([3.14159*m^3+1.9635e-7*x(1)-0.002*3.14159*0.0005^2;
6*(m+x(2))*1.9635e-7-3.14159*m^3+1.9635e-7*x(3)+1.9635e-7*x(6)*(m+x(2))/x(4)-0.002*3.14159*0.0005^2;
3.14159*m^3*x(4)/(m+x(2))+1.9635e-7*x(1)*x(5)/(x(1)+x(3))-1/6*1e-6;
(6*1.9635e-7*(m+x(2))-3.14159*m^3)*x(4)/(m+x(2))+1.9635e-7*x(3)*x(5)/(x(1)+x(3))+1.9635e-7*x(6)-1/6*1e-6;
fix(0.05/(x(1)+x(3)))*(32*0.0010087*x(5)/0.0005^2*x(3)+7.16*0.07275*(3*(1.0087*x(5)/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x(5)/72.75)^(2/3))))-6455.68*x(6)
0.05/(m+x(2))*(32*0.0010087*x(4)/0.0005^2*x(2)+7.16*0.07275*(3*(1.0087*x(4)/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x(4)/72.75)^(2/3))))-6455.68*x(6)]);
options = optimset('MaxIter',1e30,'MaxFunEvals',1e50,'TolFun',1e-50);
x=fsolve(f,[0.005 0.0005 0.001 0.15 0.9 0.3],options)

No solution found.

fsolve stopped because the relative size of the current step is less than the
default value of the step size tolerance squared, but the vector of function values
is not near zero as measured by the selected value of the function tolerance.

<stopping criteria details>

x =

0.0070 -0.0001 0.0030 0.0344 1.0563 0.4352

fsolve stopped because the relative norm of the current step, 2.848645e-15, is less than
max(options.TolX^2,eps) = 1.000000e-12. However, the sum of squared function values,
r = 3.833839e-15, exceeds sqrt(options.TolFun) = 1.000000e-25.

Optimization Metric                                              Options
relative norm(step) = 2.85e-15             max(TolX^2,eps) = 1e-12 (default)
r = 3.83e-15                                  sqrt(TolFun) =  1.0e-25 (selected)

如上，初值已經取的很接近實驗值了，還是沒有滿意的結果。將matlab算的帶入六個方程，
ans =
4.7168e-12
ans =
  -8.4160e-11
ans =
1.5716e-09
ans =
  -1.5222e-09
ans =
15.7955
ans =
  232.7567
后兩項收斂不是很好。我已經嘗試把optimset調到很小了，可是還是沒有滿意的結果。
求助各位大牛，這個該怎么解？還是方程需要修改？

如果圓滿解決，可以加金幣~

回復此樓

» 猜你喜歡

【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經有0人回復
中德博士后2026項目5月30日截止- 國內博士后入站、德國兩年全職已經有0人回復
物理學I論文潤色/翻譯怎么收費? 已經有235人回復
面向全球申請人的中德博士后2026項目5月30日截止- 國內博后入站德國兩年全職已經有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經有2人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。已經有0人回復
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！已經有1人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人已經有0人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人已經有0人回復
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8 已經有0人回復

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

matlab解二元方程組，錯誤提示，求解決方法已經有11人回復
MATLAB解非線性方程組已經有9人回復
matlab 解非線性偏微分方程組已經有6人回復
MATLAB求解五元非線性方程組，在線等已經有7人回復
非線性二階微分方程組求解，matlab 已經有12人回復
Matlab中ode45出現(xiàn)大誤差，關于用ODE45解偏微分方程的疑問，有待解決，大家共同討論！已經有4人回復
Matlab用牛頓法求解非線性方程組問題已經有8人回復
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經有11人回復
matlab如何求解一個非線性微分方程組已經有8人回復
sos!!1stopt解非線性方程組，一直出現(xiàn)常數(shù)定義錯誤的問題！求大神幫助！已經有12人回復
matlab中利用LM算法解非線性方程組的問題已經有5人回復
MATLAB求解非線性方程組已經有5人回復
關于matlab求解非線性指數(shù)方程組出現(xiàn)問題已經有3人回復
【求助】一些A、B類錯誤求解已經有15人回復
求Matlab高手解決線性方程組的迭代求解問題已經有17人回復
matlab求解非線性方程組，錯誤提示怎么解決已經有5人回復
求助matlab---fsolve解非線性方程組已經有6人回復
求助！matlab用fsolve函數(shù)求解非線性方程組的問題！已經有19人回復
matlab報錯，求解已經有5人回復
matlab求解非線性方程組已經有16人回復
matlab問題求助已經有16人回復
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經有11人回復
請教Matlab高手解多元非線性方程組中的參數(shù) 已經有6人回復
matlab的fsove 命令求解非線性方程組已經有6人回復
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！已經有18人回復
【求助】求教matlab解非線性方程組已經有9人回復

1樓 2014-11-05 11:36:41

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

monixia

木蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 2535.6
紅花: 3
帖子: 80
在線: 24.2小時
蟲號: 1817862
注冊: 2012-05-15
性別: GG
專業(yè): 化學反應工程

引用回帖:

2樓: Originally posted by 月只藍 at 2014-11-05 13:40:03
1stOpt代碼（注意不是MATLAB軟件代碼）:
Constant m=0.00040;
Parameters x(1:6);
Function 3.14159*m^3+1.9635e-7*x1-0.002*3.14159*0.0005^2;
6*(m+x2)*1.9635e-7-3.14159*m^3+1.9635e-7*x3+1.9635e-7*x6*(m+ ...

非常感謝你的解答，1stOpt用起來真方便。不過，我有個問題就是，我希望找到初值附近的零點，但是幾乎每次算出來的值都差很多（至少應該沒有負數(shù)出現(xiàn)）。是不是不太容易實現(xiàn)在初值附近尋找呢？

贊一下

回復此樓

3樓2014-11-05 19:09:09

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

查看全部 7 個回答

月只藍

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

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

【答案】應助回帖

感謝參與，應助指數(shù) +1

1stOpt代碼（注意不是MATLAB軟件代碼）:

CODE:

Constant m=0.00040;

Parameters x(1:6);

Function 3.14159*m^3+1.9635e-7*x1-0.002*3.14159*0.0005^2;

6*(m+x2)*1.9635e-7-3.14159*m^3+1.9635e-7*x3+1.9635e-7*x6*(m+x2)/x4-0.002*3.14159*0.0005^2;

3.14159*m^3*x4/(m+x2)+1.9635e-7*x1*x5/(x1+x3)-1/6*1e-6;

(6*1.9635e-7*(m+x2)-3.14159*m^3)*x4/(m+x2)+1.9635e-7*x3*x5/(x1+x3)+1.9635e-7*x6-1/6*1e-6;

(0.05/(x1+x3))*(32*0.0010087*x5/0.0005^2*x3+7.16*0.07275*(3*(1.0087*x5/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x5/72.75)^(2/3))))-6455.68*x6;

0.05/(m+x2)*(32*0.0010087*x4/0.0005^2*x2+7.16*0.07275*(3*(1.0087*x4/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x4/72.75)^(2/3))))-6455.68*x6;

計算結果：

CODE:

函數(shù)表達式 1: 3.14159*0.0004^3+1.9635e-7*x1-0.002*3.14159*0.0005^2

         2: 6*(0.0004+x2)*1.9635e-7-3.14159*0.0004^3+1.9635e-7*x3+1.9635e-7*x6*(0.0004+x2)/x4-0.002*3.14159*0.0005^

            2

         3: 3.14159*0.0004^3*x4/(0.0004+x2)+1.9635e-7*x1*x5/(x1+x3)-1/6*1e-6

         4: (6*1.9635e-7*(0.0004+x2)-3.14159*0.0004^3)*x4/(0.0004+x2)+1.9635e-7*x3*x5/(x1+x3)+1.9635e-7*x6-1/6*1e

            -6

         5: (0.05/(x1+x3))*(32*0.0010087*x5/0.0005^2*x3+7.16*0.07275*(3*(1.0087*x5/72.75)^(2/3)/0.0005/(1+3.34*(1.0087

            *x5/72.75)^(2/3))))-6455.68*x6

         6: 0.05/(0.0004+x2)*(32*0.0010087*x4/0.0005^2*x2+7.16*0.07275*(3*(1.0087*x4/72.75)^(2/3)/0.0005/(1+3.34

            *(1.0087*x4/72.75)^(2/3))))-6455.68*x6

目標函數(shù)值: 8.44256812628831E-7

x1: 1.38116851215046

x2: 0.213584946434081

x3: 0.106474657467166

x4: 0.0176943028978422

x5: 0.2480594437683

x6: 0.0180988189842738

6個方程目標函數(shù)值：

ans =

  2.6982e-007

ans =

  3.1421e-007

ans =

 -1.2143e-007

ans =

 -1.3880e-007

ans =

  1.4211e-013

ans =

 -8.5265e-014

1.5版本計算結果不唯一，以下再例舉6組解：

CODE:

目標函數(shù)值: 8.81354127610376E-9

x1: 0.0270233957846693

x2: -4.08878886642551E-5

x3: 0.0126231922603729

x4: 0.0323838969507018

x5: 1.10773562695781

x6: 0.383966561657027

目標函數(shù)值: 5.88138553710961E-7

x1: 2.08840573452721

x2: -0.0811262956407749

x3: 0.416992754705403

x4: 0.19145905588163

x5: 1.1205806122223

x6: 0.187004870338678

目標函數(shù)值: 6.80817986394196E-7

x1: 2.10149925984767

x2: -0.218796314718131

x3: 0.718745532551135

x4: 0.109052962294469

x5: 0.42221573528721

x6: 0.10785410685195

目標函數(shù)值: 1.10547660307994E-6

x1: 2.53429273796369

x2: -0.218841716843157

x3: 1.13544861896696

x4: 0.431086688290742

x5: 1.38395310620881

x6: 0.42858810590596

目標函數(shù)值: 1.32464424554872E-6

x1: 4.35394653830248

x2: -0.286471171270745

x3: 2.04705166768062

x4: 0.401695771492555

x5: 1.24966197787108

x6: 0.39985126552927

目標函數(shù)值: 7.67093610988707E-7

x1: 1.6997802237498

x2: -0.00709436222448859

x3: -0.505554887761384

x4: 0.00623112985186669

x5: 0.00101373791179563

x6: -0.000417366283335469

贊一下

回復此樓

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

2樓2014-11-05 13:40:03

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

月只藍

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

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

【答案】應助回帖

引用回帖:

3樓: Originally posted by monixia at 2014-11-05 19:09:09
非常感謝你的解答，1stOpt用起來真方便。不過，我有個問題就是，我希望找到初值附近的零點，但是幾乎每次算出來的值都差很多（至少應該沒有負數(shù)出現(xiàn)）。是不是不太容易實現(xiàn)在初值附近尋找呢？...

可以設定x的取值范圍

贊一下

回復此樓

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

4樓2014-11-05 20:52:13

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

月只藍

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

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

【答案】應助回帖

引用回帖:

比如：
讓所有x都大于0：

CODE:

Constant m=0.00040;

Parameters x(1:6)[0,];

Function 3.14159*m^3+1.9635e-7*x1-0.002*3.14159*0.0005^2;

6*(m+x2)*1.9635e-7-3.14159*m^3+1.9635e-7*x3+1.9635e-7*x6*(m+x2)/x4-0.002*3.14159*0.0005^2;

3.14159*m^3*x4/(m+x2)+1.9635e-7*x1*x5/(x1+x3)-1/6*1e-6;

(6*1.9635e-7*(m+x2)-3.14159*m^3)*x4/(m+x2)+1.9635e-7*x3*x5/(x1+x3)+1.9635e-7*x6-1/6*1e-6;

(0.05/(x1+x3))*(32*0.0010087*x5/0.0005^2*x3+7.16*0.07275*(3*(1.0087*x5/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x5/72.75)^(2/3))))-6455.68*x6;

0.05/(m+x2)*(32*0.0010087*x4/0.0005^2*x2+7.16*0.07275*(3*(1.0087*x4/72.75)^(2/3)/0.0005/(1+3.34*(1.0087*x4/72.75)^(2/3))))-6455.68*x6;

計算結果：
目標函數(shù)值: 2.37290707692599E-8
x1: 0.0105848644277172
x2: 0.0161346934843435
x3: 0.000420183576622565
x4: 0.115861612534128
x5: 0.876247555934376
x6: 0.132259599167984

贊一下

回復此樓

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

5樓2014-11-05 20:57:57

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

查看全部 7 個回答

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

[考研] 求收留 +4	1943443204 2026-03-28	4/200	2026-03-29 10:35 by 唐沐兒
[考研] 348求調劑 +5	小懶蟲不懶了 2026-03-28	5/250	2026-03-29 10:34 by 唐沐兒
[考研] 279求調劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 一志愿北京化工大學材料與化工（085600）296求調劑 +13	稻妻小編 2026-03-26	13/650	2026-03-29 08:20 by fmesaito
[考研] 291求調劑 +4	Y-cap 2026-03-29	5/250	2026-03-29 08:08 by 無際的草原
[考研] 求調劑 +7	爭取九點睡 2026-03-28	8/400	2026-03-28 21:07 by 爭取九點睡
[考研] 394求調劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 一志愿廈門大學化學學碩307求調劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 【求調劑】085601材料工程專碩 \| 總分272 \| +6	腳滑的守法公民 2026-03-27	6/300	2026-03-28 11:02 by gjlllb
[考研] 085404求調劑，總分309，本科經歷較為豐富 +4	來財aa 2026-03-25	4/200	2026-03-28 07:41 by 棒棒球手
[考研] 274求調劑 +17	顧九笙要謙虛 2026-03-24	23/1150	2026-03-27 15:16 by caszguilin
[考研] 322求調劑 +4	我真的很想學習 2026-03-23	4/200	2026-03-27 13:51 by 楊楊楊紫
[考研] 333求調劑 +3	question挽風 2026-03-23	3/150	2026-03-27 11:29 by 不吃魚的貓
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調劑】 +4	Higraduate 2026-03-22	8/400	2026-03-26 20:34 by Higraduate
[考研] 材料與化工304求B區(qū)調劑 +3	邱gl 2026-03-25	3/150	2026-03-25 19:03 by Ainin_
[考研] 生物技術與工程 +3	1294608413 2026-03-25	4/200	2026-03-25 18:02 by 1294608413
[考研] 【2026考研調劑】制藥工程 284分求相關專業(yè)調劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 0854電子信息求調劑 +7	α____ 2026-03-22	9/450	2026-03-25 13:37 by α____
[考研] 一志愿河北工業(yè)大學0817化工278分求調劑 +7	jhybd 2026-03-23	12/600	2026-03-24 09:03 by jhybd
[考研] 一志愿國科過程所081700，274求調劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

monixia

[求助] 求解決matlab解非線性方程組出現(xiàn)的問題 已有1人參與

» 猜你喜歡

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

monixia

月只藍

【答案】應助回帖

月只藍

【答案】應助回帖

月只藍

【答案】應助回帖

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 求解決matlab解非線性方程組出現(xiàn)的問題已有1人參與

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