小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
帮我看看问题出在哪？
& && && && && && & 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.0& && && && && && && &2.49e+014& && && && && &1
& &&&1& && && & 6& & 1.0& && &&&0.92278& && &2.86e+013& && && && && &1
& &&&2& && && & 9& & 4.3& && &0.0759773& && &5.33e+011& && && && &2.31
& &&&3& && && &12& && & 0.0824035& & 0.& && &3.85e+007& && && && &2.31
& &&&4& && && &15& & 1.1& &2.1& && &3.77e+003& && && && &2.31
& &&&5& && && &16& & 1.1& &3.3& && &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&
& & 1.0000
& & 0.9993
&&1.0e-006 *
& &-0.3436& & 0.0006
extiflag =
& && & iterations: 5
& && &&&funcCount: 16
& && &&&algorithm: 'trust-region-dogleg'
& & firstorderopt: 3.
& && && & message:
&&1.0e+010 *
& & 1.0957& &-0.0000
& &-0.0013& & 0.0013
帮我看看我的问题出在哪，我该怎么解决？这是一个用来计算化学平衡组成
改了，还是不行
是解非线性方程组吗？再改一下初值看看，同时看看目前的解代入原方程有多接近0。可不可以把待求方程组发上来，我们一起看看。
function ChemEquil& &
P=600;& && && && && & %压力单位：kpa
%m=50;& && && && &&&%氢气：草酸二甲酯摩尔比
%x0=0.001;
%x=fsolve(@equil1,x0,)
options=optimset('Display','iter', 'LargeScale','off','TolFun', 1.0e-22,'TolX',1.0e-22)
=fsolve(@equil1,x0,options)
%...................................................
%平衡组成关系式
function f=equil1(x)
%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);
这是源程序，可以帮我看看错在哪？是不是要改算法
目标函数值fval =
&&1.0e-006 *
& &-0.3436& & 0.0006
& & 1.0000
& &&&0.9993
这一组解，已经使fval达到1e-6的高精度了，已经算是成功地解方程了，你觉得有哪里不对吗？
你的结果出现“Equation solved, fsolve stalled.”这个吗，以前我算过别的平衡组成，都没有出现这种情况
这一组x已经使得原方程相当接近于0了，如果这组解符合实际的话，就可以采用。
MATLAB的Help中Equation Solved, fsolve Stalled的意思是：
fsolve found a point where the sum of squares of function values is less than the square root of the TolFun tolerance. However, the last step was less than the TolX tolerance, indicating the function may be changing rapidly, or that the function is not smooth near the final point. This is the meaning of stalled.
问题是，不符合
在确认方程没问题的情况下，可以根据x的预估区间调节初值。你希望两个x各落在什么范围内？
算出来的结果选择性不对
x1=0.9~1;X2=0.05~0.2
那你检查一下方程是否写对了。比如方程中各参数的单位是否一致的，比如温度是K还是℃，压力是Pa、atm还是kPa等等。
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研→ 帮我看看我得问题？？？？？
帮我看看我得问题？？？？？
健康咨询描述：
正常人倒立是一定会还是可能会出现头部血管跳动？
曾经的治疗情况和效果：
想得到怎样的帮助：谢谢谢谢谢谢谢谢
wolaifanni
其他类似问题
医生回复区
擅长: 牙体牙髓病,颌面外科
帮助网友：27453称赞：2021
&&&&&&病情分析：&&&&&&你好，从你叙述的情况看是会的。&&&&&&指导意见：&&&&&&因为这时候血液需要回流到心脏，需要一定的压力会出现你说的情况。
您可能关注的问题
用药指导/吃什么药好
本品用于湿热瘀滞所致的带下病。...
参考价格：￥28
本品用于跌打损伤，跖骨、趾骨骨折，瘀血肿痛，吐血...
参考价格：￥17.5
下载APP，免费快速问医生推荐到广播
290371 人聚集在这个小组
(格桑梅朵)
(格桑梅朵)
(走近你 离开你)
第三方登录：