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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

alice919

金蟲 (小有名氣)

應助: 2 (幼兒園)
金幣: 1083.6
紅花: 2
帖子: 52
在線: 59.7小時
蟲號: 439553
注冊: 2007-10-15
性別: MM
專業(yè): 化學反應工程

[求助] 幫我看看問題出在哪？

options =

               Display: 'iter'
            MaxFunEvals: []
               MaxIter: []
                  TolFun: []
                  TolX: []
            FunValCheck: []
               OutputFcn: []
               PlotFcns: []
         ActiveConstrTol: []
               Algorithm: []
AlwaysHonorConstraints: []
         BranchStrategy: []
         DerivativeCheck: []
            Diagnostics: []
         DiffMaxChange: []
         DiffMinChange: []
            FinDiffType: 'central'
      GoalsExactAchieve: []
            GradConstr: []
               GradObj: []
               HessFcn: []
               Hessian: []
               HessMult: []
            HessPattern: []
            HessUpdate: []
         InitialHessType: []
      InitialHessMatrix: []
      InitBarrierParam: []
   InitTrustRegionRadius: []
               Jacobian: []
               JacobMult: []
            JacobPattern: []
            LargeScale: 'off'
      LevenbergMarquardt: []
         LineSearchType: []
               MaxNodes: []
            MaxPCGIter: []
         MaxProjCGIter: []
            MaxRLPIter: []
            MaxSQPIter: []
               MaxTime: []
         MeritFunction: []
               MinAbsMax: []
   NodeDisplayInterval: []
      NodeSearchStrategy: []
      NonlEqnAlgorithm: []
      NoStopIfFlatInfeas: []
         ObjectiveLimit: []
   PhaseOneTotalScaling: []
         Preconditioner: []
      PrecondBandWidth: []
         RelLineSrchBnd: []
RelLineSrchBndDuration: []
            ScaleProblem: []
               Simplex: []
   SubproblemAlgorithm: []
                  TolCon: []
               TolConSQP: []
            TolGradCon: []
                  TolPCG: []
               TolProjCG: []
            TolProjCGAbs: []
               TolRLPFun: []
            TolXInteger: []
               TypicalX: []
            UseParallel: []

                                       Norm of    First-order Trust-region
Iteration  Func-count    f(x)       step       optimality radius
   0       3 2.06697e+014                   2.49e+014             1
   1       6 1.03581e+012       0.92278    2.86e+013             1
   2       9 4.38753e+007    0.0759773    5.33e+011          2.31
   3       12    0.0824035 0.000500989    3.85e+007          2.31
   4       15 1.18085e-013 2.17124e-008    3.77e+003          2.31
   5       16 1.18085e-013 3.37394e-017    3.77e+003          2.31

Equation solved, fsolve stalled.

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

<stopping criteria details>

x =

1.0000
0.9993

fval =

  1.0e-006 *

-0.3436 0.0006

extiflag =

   2

output =

   iterations: 5
      funcCount: 16
      algorithm: 'trust-region-dogleg'
firstorderopt: 3.7653e+003
      message: [1x774 char]

JAC =

  1.0e+010 *

1.0957 -0.0000
-0.0013 0.0013
幫我看看我的問題出在哪，我該怎么解決？這是一個用來計算化學平衡組成

[ 來自小組化工工程工藝 ]

回復此樓

» 猜你喜歡

304求調(diào)劑已經(jīng)有6人回復
材料工程專碩調(diào)劑已經(jīng)有6人回復
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復
求材料調(diào)劑已經(jīng)有8人回復
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復
一志愿華中科技大學，080502，354分求調(diào)劑已經(jīng)有4人回復
一志愿吉林大學材料學碩321求調(diào)劑已經(jīng)有6人回復
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復
330求調(diào)劑已經(jīng)有3人回復

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

請問哪位可以幫幫我，看看問題出在哪里？已經(jīng)有6人回復
關(guān)于留學期限的問題，大家?guī)臀铱纯?/a> 已經(jīng)有5人回復

關(guān)于克隆技術(shù)——pcr產(chǎn)物測序問題已經(jīng)有16人回復
請大家?guī)臀铱纯次业膕outhern blot 圖片，不知道問題出在哪里了，謝謝已經(jīng)有16人回復
各位大俠，幫我看看老外審稿人提出的一個問題。謝謝你們了！已經(jīng)有7人回復
求助高斯中出現(xiàn)錯誤希望幫忙解決已經(jīng)有9人回復
幫我看看這個極化曲線為什么這樣子？問題出在哪里？已經(jīng)有5人回復
大家愛幫我看看這段MATLAB程序為啥出錯了已經(jīng)有5人回復
PCR好多次了，一直沒結(jié)果，麻煩各位幫我看看是哪出問題了？已經(jīng)有15人回復
哪位大神幫我看看這個從牛肚子里出來的東西是什么已經(jīng)有14人回復
請幫我看看提取質(zhì)粒的圖是出了什么問題已經(jīng)有19人回復
求高手幫看看我的筆記本問題出在哪？已經(jīng)有6人回復
大家?guī)臀铱匆粋€fortran的程序，我總是計算不出正確的結(jié)果已經(jīng)有20人回復
麻煩塔菲爾曲線高手幫我看看這幾個圖，說明什么問題，或者能從中得出什么信息？已經(jīng)有12人回復
【求助】鋰電池三電極結(jié)果，幫我看看哪出問題了已經(jīng)有4人回復
【求助】水熱法制備磷酸鐵鋰，大家?guī)臀铱纯?問題出在哪已經(jīng)有35人回復
【求助/交流】急!有關(guān)PCR中出現(xiàn)的一些問題,大家?guī)臀铱纯窗? 已經(jīng)有5人回復

1樓 2013-07-26 14:31:37

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

月只藍

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

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

【答案】應助回帖

★ ★
感謝參與，應助指數(shù) +1
xiegangmai: 金幣+2, 謝謝參與！ 2013-08-01 22:33:51

“方程被求解，fsolve函數(shù)停止計算”
迭代最后一步與上一步的差值已經(jīng)小于options中的預設(shè)值。

試著把options中的第4和第5項的TolFun和TolX調(diào)得更小一些，試試。

贊一下

回復此樓

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

2樓2013-07-26 14:51:35

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

alice919

金蟲 (小有名氣)

應助: 2 (幼兒園)
金幣: 1083.6
紅花: 2
帖子: 52
在線: 59.7小時
蟲號: 439553
注冊: 2007-10-15
性別: MM
專業(yè): 化學反應工程

引用回帖:

2樓: Originally posted by 月只藍 at 2013-07-26 14:51:35
“方程被求解，fsolve函數(shù)停止計算”
迭代最后一步與上一步的差值已經(jīng)小于options中的預設(shè)值。

試著把options中的第4和第5項的TolFun和TolX調(diào)得更小一些，試試。

改了，還是不行

回復此樓

3樓2013-07-26 15:04:32

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

月只藍

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

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

引用回帖:

3樓: Originally posted by alice919 at 2013-07-26 15:04:32
改了，還是不行...

是解非線性方程組嗎？再改一下初值看看，同時看看目前的解代入原方程有多接近0�？刹豢梢园汛蠓匠探M發(fā)上來，我們一起看看。

贊一下

回復此樓

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

4樓2013-07-26 15:06:35

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

alice919

金蟲 (小有名氣)

應助: 2 (幼兒園)
金幣: 1083.6
紅花: 2
帖子: 52
在線: 59.7小時
蟲號: 439553
注冊: 2007-10-15
性別: MM
專業(yè): 化學反應工程

引用回帖:

function ChemEquil
clear all
clc
global P
P=600; %壓力單位：kpa
%m=50; %氫氣：草酸二甲酯摩爾比
%KP1=[23.4 21.8 20.3 18.9 17.6 16.4 15.3 14.2 13.3 12.4 11.5 10.8];
%KP2=[0.0393 0.0375 0.0358 0.0342 0.0328 0.0314 0.0301 0.0289 0.0277 0.0267...
% 0.0257 0.0247];
%KP3=[4.4346E05 3.6476E05 3.0202E05 2.5168E05 2.1101E05 1.7796E05 1.5094E05...
%1.2872E05 1.1036E05 0.9509E05 0.8233E05 0.7163E05];
%x0=0.001;
%x=fsolve(@equil1,x0,[])
x0=[1 0 ]';
options=optimset('Display','iter', 'LargeScale','off','TolFun', 1.0e-22,'TolX',1.0e-22)
[x,fval,extiflag,output,JAC]=fsolve(@equil1,x0,options)

%...................................................
%平衡組成關(guān)系式
function f=equil1(x)
global P
%f=x^2*(51-x)-21.8*((50-2*x)^2*(1-x)*P)
tmp1=50-2*x(1)-2*x(2);
tmp2=51-x(1)-x(2);
f(1)=(x(1)-x(2))*(x(1)+x(2))*tmp2*100-8630*((1-x(1))*tmp1^2*P);
f(2)=(x(1)+x(2))*x(2)*tmp2*100-10.4*((x(1)-x(2))*tmp1^2*P);
這是源程序，可以幫我看看錯在哪？是不是要改算法

贊一下

回復此樓

5樓2013-07-26 15:07:09

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

月只藍

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

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

【答案】應助回帖

★ ★
xiegangmai: 金幣+2, 謝謝參與！ 2013-08-01 22:34:01

引用回帖:

5樓: Originally posted by alice919 at 2013-07-26 15:07:09
function ChemEquil
clear all
clc
global P
P=600; %壓力單位：kpa
%m=50; %氫氣：草酸二甲酯摩爾比
%KP1=;
%KP2=;
%KP3=;
%x0=0.001;
%x=fsolve(@equil1,x0,[])
x0= ...

目標函數(shù)值fval =

1.0e-006 *

-0.3436 0.0006

x =

1.0000
0.9993
這一組解，已經(jīng)使fval達到1e-6的高精度了，已經(jīng)算是成功地解方程了，你覺得有哪里不對嗎？

贊一下

回復此樓

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

6樓2013-07-26 15:23:01

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

alice919

金蟲 (小有名氣)

應助: 2 (幼兒園)
金幣: 1083.6
紅花: 2
帖子: 52
在線: 59.7小時
蟲號: 439553
注冊: 2007-10-15
性別: MM
專業(yè): 化學反應工程

引用回帖:

6樓: Originally posted by 月只藍 at 2013-07-26 15:23:01
目標函數(shù)值fval =

1.0e-006 *

-0.3436 0.0006

x =

1.0000
0.9993
這一組解，已經(jīng)使fval達到1e-6的高精度了，已經(jīng)算是成功地解方程了，你覺得有哪里不對嗎？...

你的結(jié)果出現(xiàn)“Equation solved, fsolve stalled.”這個嗎，以前我算過別的平衡組成，都沒有出現(xiàn)這種情況

贊一下

回復此樓

7樓2013-07-26 15:31:25

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

月只藍

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

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

【答案】應助回帖

★ ★
xiegangmai: 金幣+2, 鼓勵討論交流！ 2013-08-01 22:34:08

引用回帖:

7樓: Originally posted by alice919 at 2013-07-26 15:31:25
你的結(jié)果出現(xiàn)“Equation solved, fsolve stalled.”這個嗎，以前我算過別的平衡組成，都沒有出現(xiàn)這種情況...

出現(xiàn)了。
這一組x已經(jīng)使得原方程相當接近于0了，如果這組解符合實際的話，就可以采用。

MATLAB的Help中Equation Solved, fsolve Stalled的意思是：
fsolve found a point where the sum of squares of function values is less than the square root of the TolFun tolerance. However, the last step was less than the TolX tolerance, indicating the function may be changing rapidly, or that the function is not smooth near the final point. This is the meaning of stalled.

贊一下

回復此樓

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

8樓2013-07-26 15:40:42

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

alice919

金蟲 (小有名氣)

應助: 2 (幼兒園)
金幣: 1083.6
紅花: 2
帖子: 52
在線: 59.7小時
蟲號: 439553
注冊: 2007-10-15
性別: MM
專業(yè): 化學反應工程

引用回帖:

8樓: Originally posted by 月只藍 at 2013-07-26 15:40:42
出現(xiàn)了。
這一組x已經(jīng)使得原方程相當接近于0了，如果這組解符合實際的話，就可以采用。

MATLAB的Help中Equation Solved, fsolve Stalled的意思是：
fsolve found a point where the sum of squares of func ...

問題是，不符合

回復此樓

9樓2013-07-26 15:44:18

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

月只藍

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

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

引用回帖:

9樓: Originally posted by alice919 at 2013-07-26 15:44:18
問題是，不符合...

在確認方程沒問題的情況下，可以根據(jù)x的預估區(qū)間調(diào)節(jié)初值。你希望兩個x各落在什么范圍內(nèi)？

贊一下

回復此樓

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

10樓2013-07-26 15:46:34

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 alice919 的主題更新

返回列表

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

[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +4	Ymlll 2026-03-18	6/300	2026-03-18 22:15 by li123456789.
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 085600材料與化工 +5	安全上岸！ 2026-03-16	5/250	2026-03-18 15:33 by cmz0325
[考研] 化學工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 311求調(diào)劑 +6	26研0 2026-03-15	6/300	2026-03-18 14:43 by haxia
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習 +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 274求調(diào)劑 +5	時間點 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky