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

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

回復此樓

» 猜你喜歡

求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復
335求調(diào)劑已經(jīng)有5人回復
求調(diào)劑已經(jīng)有7人回復
一志愿吉大化學322求調(diào)劑已經(jīng)有4人回復
環(huán)境學碩288求調(diào)劑已經(jīng)有8人回復
341求調(diào)劑(一志愿湖南大學070300) 已經(jīng)有6人回復
一志愿北京化工大學 070300 學碩 336分求調(diào)劑已經(jīng)有7人回復
299求調(diào)劑已經(jīng)有5人回復
306求0703調(diào)劑一志愿華中師范已經(jīng)有8人回復
一志愿河北工業(yè)大學0817化工278分求調(diào)劑已經(jīng)有11人回復

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

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

關于克隆技術(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人回復
【求助/交流】急!有關PCR中出現(xiàn)的一些問題,大家?guī)臀铱纯窗? 已經(jīng)有5人回復

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

已閱回復此樓關注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中的預設值。

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

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)

%...................................................
%平衡組成關系式
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

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

查看全部 13 個回答

月只藍

主管區(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中的預設值。

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

贊一下

回復此樓

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

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

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

alice919

金蟲 (小有名氣)

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

引用回帖:

改了，還是不行

回復此樓

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

已閱回復此樓關注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

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

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

[考研] 284求調(diào)劑 +9	Zhao anqi 2026-03-22	9/450	2026-03-23 23:11 by peike
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-23	4/200	2026-03-23 21:25 by pswait
[考研] 一志愿重慶大學085700資源與環(huán)境，總分308求調(diào)劑 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta
[考研] 上海電力大學材料防護與新材料重點實驗室招收調(diào)劑研究生（材料、化學、電化學，環(huán)境） +3	我愛學電池 2026-03-23	3/150	2026-03-23 17:16 by AZMK
[考研] 08工學調(diào)劑 +7	用戶573181 2026-03-20	11/550	2026-03-23 15:47 by 我愛學習學習使?/a>
[考研] 350求調(diào)劑 +6	weudhdk 2026-03-19	6/300	2026-03-23 15:47 by tangyuan0840221
[考研] 298求調(diào)劑一志愿211 +3	上岸6666@ 2026-03-20	3/150	2026-03-22 15:50 by ColorlessPI
[考研] 一志愿華中科技大學071000，求調(diào)劑 +4	沿岸有貝殼6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學員8dgXkO
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3	李軼男003 2026-03-20	3/150	2026-03-21 13:28 by 搏擊518
[考研] 求調(diào)劑 +3	白QF 2026-03-21	3/150	2026-03-21 13:12 by zhukairuo
[考研] 化學求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 華東師范大學-071000生物學-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] A區(qū)線材料學調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 材料學碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 一志愿中國海洋大學，生物學，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo