| 查看: 1501 | 回复: 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 帮我看看我的问题出在哪，我该怎么解决？这是一个用来计算化学平衡组成 [ 来自小组 化工工程工艺 ] |
回复此楼» 猜你喜欢
312求调剂
已经有5人回复
-
网络空间安全0839招调剂
已经有4人回复
-
329求调剂
已经有3人回复
-
【2026考研调剂】制药工程 284分 求相关专业调剂名额
已经有8人回复
-
085602 289分求调剂
已经有7人回复
-
材料学硕333求调剂
已经有4人回复
-
各位老师您好：本人初试372分
已经有6人回复
-
334分 一志愿武理-080500 材料求调剂
已经有3人回复
-
0854电子信息求调剂
已经有9人回复
-
347求调剂
已经有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人回复
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
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2013-07-26 14:51:35
alice919
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1083.6
- 红花: 2
- 帖子: 52
- 在线: 59.7小时
- 虫号: 439553
- 注册: 2007-10-15
- 性别: MM
- 专业: 化学反应工程
3楼2013-07-26 15:04:32
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2013-07-26 15:06:35
6