alice919

金蟲 (小有名氣)

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

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

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
幫我看看我的問題出在哪，我該怎么解決？這是一個用來計算化學(xué)平衡組成

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

回復(fù)此樓

» 猜你喜歡

【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑已經(jīng)有4人回復(fù)
317求調(diào)劑已經(jīng)有5人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有10人回復(fù)
一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。已經(jīng)有3人回復(fù)
一志愿福大288有機化學(xué)，求調(diào)劑已經(jīng)有3人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有6人回復(fù)
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有6人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有10人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有11人回復(fù)
一志愿南昌大學(xué)，327分，材料與化工085600 已經(jīng)有3人回復(fù)

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

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

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

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

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

alice919

金蟲 (小有名氣)

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

引用回帖:

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

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

改了，還是不行

回復(fù)此樓

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

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

查看全部 13 個回答

月只藍(lán)

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

專家經(jīng)驗: +1059
應(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
xiegangmai: 金幣+2, 謝謝參與！ 2013-08-01 22:33:51

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

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

贊一下

回復(fù)此樓

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

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

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

月只藍(lán)

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

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

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

alice919

金蟲 (小有名氣)

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

引用回帖:

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);
這是源程序，可以幫我看看錯在哪？是不是要改算法

贊一下

回復(fù)此樓

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

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

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

[考研] 294求調(diào)劑材料與化工專碩 +11	陌の森林 2026-03-18	11/550	2026-03-19 13:22 by houyaoxu
[考研] 招收調(diào)劑碩士 +3	lidianxing 2026-03-19	7/350	2026-03-19 12:31 by 山東小li
[考研] 材料080500調(diào)劑求收留 +4	一顆meteor 2026-03-13	4/200	2026-03-19 10:32 by 30660438
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 331求調(diào)劑（0703有機化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 274求調(diào)劑 +5	時間點 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 307求調(diào)劑 +5	超級伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty