24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

alice919

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 1083.6
红花: 2
帖子: 52
在线: 59.7小时
虫号: 439553
注册: 2007-10-15
性别: MM
专业: 化学反应工程

[求助] 帮我看看问题出在哪？

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
帮我看看我的问题出在哪，我该怎么解决？这是一个用来计算化学平衡组成

[ 来自小组化工工程工艺 ]

回复此楼

» 猜你喜欢

假如你的研究生提出不合理要求已经有4人回复
论文终于录用啦！满足毕业条件了已经有27人回复
所感已经有3人回复
要不要辞职读博？已经有7人回复
不自信的我已经有11人回复
北核录用已经有3人回复
实验室接单子已经有3人回复
磺酰氟产物，毕不了业了！已经有8人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有10人回复
26申博（荧光探针方向，有机合成）已经有4人回复

» 本主题相关价值贴推荐，对您同样有帮助:

请问哪位可以帮帮我，看看问题出在哪里？已经有6人回复
关于留学期限的问题，大家帮我看看已经有5人回复
关于克隆技术——pcr产物测序问题已经有16人回复
请大家帮我看看我的southern blot 图片，不知道问题出在哪里了，谢谢已经有16人回复
各位大侠，帮我看看老外审稿人提出的一个问题。谢谢你们了！已经有7人回复
求助高斯中出现错误希望帮忙解决已经有9人回复
帮我看看这个极化曲线为什么这样子？问题出在哪里？已经有5人回复
大家爱帮我看看这段MATLAB程序为啥出错了已经有5人回复
PCR好多次了，一直没结果，麻烦各位帮我看看是哪出问题了？已经有15人回复
哪位大神帮我看看这个从牛肚子里出来的东西是什么已经有14人回复
请帮我看看提取质粒的图是出了什么问题已经有19人回复
求高手帮看看我的笔记本问题出在哪？已经有6人回复
大家帮我看一个fortran的程序，我总是计算不出正确的结果已经有20人回复
麻烦塔菲尔曲线高手帮我看看这几个图，说明什么问题，或者能从中得出什么信息？已经有12人回复
【求助】锂电池三电极结果，帮我看看哪出问题了已经有4人回复
【求助】水热法制备磷酸铁锂，大家帮我看看问题出在哪已经有35人回复
【求助/交流】急!有关PCR中出现的一些问题,大家帮我看看啊! 已经有5人回复

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

alice919

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 1083.6
红花: 2
帖子: 52
在线: 59.7小时
虫号: 439553
注册: 2007-10-15
性别: MM
专业: 化学反应工程

引用回帖:

2楼: Originally posted by 月只蓝 at 2013-07-26 14:51:35
“方程被求解，fsolve函数停止计算”
迭代最后一步与上一步的差值已经小于options中的预设值。

试着把options中的第4和第5项的TolFun和TolX调得更小一些，试试。

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