| 查看: 1620 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
alice919金虫 (小有名气)
|
[求助]
帮我看看问题出在哪？
|
||
|
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 帮我看看我的问题出在哪，我该怎么解决？这是一个用来计算化学平衡组成 [ 来自小组 化工工程工艺 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问哪位可以帮帮我，看看问题出在哪里？
已经有6人回复
- 关于留学期限的问题，大家帮我看看
已经有5人回复
- 关于克隆技术——pcr产物测序问题
已经有16人回复
- 请大家帮我看看我的southern blot 图片，不知道问题出在哪里了，谢谢
已经有16人回复
- 各位大侠，帮我看看老外审稿人提出的一个问题。谢谢你们了！
已经有7人回复
- 求助高斯中出现错误 希望帮忙解决
已经有9人回复
- 帮我看看这个极化曲线为什么这样子？问题出在哪里？
已经有5人回复
- 大家爱帮我看看这段MATLAB程序为啥出错了
已经有5人回复
- PCR好多次了，一直没结果，麻烦各位帮我看看是哪出问题了？
已经有15人回复
- 哪位大神帮我看看这个从牛肚子里出来的东西是什么
已经有14人回复
- 请帮我看看提取质粒的图是出了什么问题
已经有19人回复
- 求高手帮看看我的笔记本问题出在哪？
已经有6人回复
- 大家帮我看一个fortran的程序，我总是计算不出正确的结果
已经有20人回复
- 麻烦塔菲尔曲线高手帮我看看这几个图，说明什么问题，或者能从中得出什么信息？
已经有12人回复
- 【求助】锂电池三电极结果，帮我看看哪出问题了
已经有4人回复
- 【求助】水热法制备磷酸铁锂，大家帮我看看 问题出在哪
已经有35人回复
- 【求助/交流】急!有关PCR中出现的一些问题,大家帮我看看啊!
已经有5人回复
alice919
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1083.6
- 红花: 2
- 帖子: 52
- 在线: 59.7小时
- 虫号: 439553
- 注册: 2007-10-15
- 性别: MM
- 专业: 化学反应工程
3楼2013-07-26 15:04:32
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2013-07-26 14:51:35
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2013-07-26 15:06:35
alice919
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1083.6
- 红花: 2
- 帖子: 52
- 在线: 59.7小时
- 虫号: 439553
- 注册: 2007-10-15
- 性别: MM
- 专业: 化学反应工程
|
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
投票:
10