|
% ********************************************************************
F=61005*((1+x(1))*(1.02+7.25*x(1))*(120-x(3))*x(2)*log(100/(x(3)-40))/(((1+x(1))*(2.93+1.92*x(1))/x(2)^2)^16*(140-x(3))))^66+...
84177+6785*x(1)*(120-x(3))+336*(1+x(1))*(1.02+7.25*x(1))*(120-x(3))*log(100/(x(3)-40))*(40+1/2*x(3))/(((1+x(1))*(2.93+1.92*x(1))/x(2)^2)^16*(140-x(3)))+47.11*x(1);
================================================
求解结果:
Optimization terminated: first-order optimality measure less
than options.TolFun and maximum constraint violation is less
than options.TolCon.
Active inequalities (to within options.TolCon = 1e-006):
lower upper ineqlin ineqnonlin
1 3
x =
0 0.45 120
fval =
84177
exitflag =
1
output =
iterations: 2
funcCount: 12
stepsize: 1
algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
firstorderopt: 0
cgiterations: []
message: [1x144 char]
===============================================
1stopt代码:
Title "Type your title here";
Parameters x1[0,1],x2[.4,5],x3[40,120];
//Variable ;
minFunction 61005*((1+x1)*(1.02+7.25*x1)*(120-x3)*x2*log(100/(x3-40))/(((1+x1)*(2.93+1.92*x1)/x2^2)^16*(140-x3)))^66+84177+6785*x1*(120-x3)+336*(1+x1)*(1.02+7.25*x1)*(120-x3)*log(100/(x3-40))*(40+1/2*x3)/(((1+x1)*(2.93+1.92*x1)/x2^2)^16*(140-x3))+47.11*x1;
=================================================
求解结果:
====== 结果 ======
迭代数: 17
计算用时(时:分:秒:毫秒): 00:00:00:156
计算中止原因: 达到收敛判定标准
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
函数表达式: 61005*((1+x1)*(1.02+7.25*x1)*(120-x3)*x2*log(100/(x3-40))/(((1+x1)*(2.93+1.92*x1)/x2^2)^16*(140-x3)))^66
+84177+6785*x1*(120-x3)+336*(1+x1)*(1.02+7.25*x1)*(120-x3)*log(100/(x3-40))*(40+1/2*x3)/(((1+x1)*(2.93
+1.92*x1)/x2^2)^16*(140-x3))+47.11*x1
目标函数值(最小): 84177
x1: 0
x2: 0.401332868776814
x3: 45.1969769224524
====== 计算结束 ======
================================================
至于什么固定常量画3D图什么的,自己看看surf命令怎么用,整个问题是很简单的,看看书就可以了。
另外这个问题用我所取的常量(固定x或者固定t2)感觉绘图意义不大,多处数值inf,固定t2只有当其等于120才能画出一个接近40度的倾侧平面。另外亦有多点处可取在84177值附近,极值意义同样不大。 |
|