谁可以用matlab或lingo解决下列非线性规划问题

min=5.6*x+6.0*y+1.2*x^2+1.2*y^2-1.2*41-36*m
st x+y>=15.2971
6x+6y+x^2+y^2-41-30m>=sqrt(m^2+5.0981^2+2*5.0981m+19.4904^2
-10<=m<=5
x>0
y>0
x^2+y^2-41-30m>0
要求有完整的程序和过程,写的好追加分
第二个类似的
min=5.6*sqrt((y-5)^2+(x+10)^2)+6*sqrt((x-5)^2+(y-8)^2)+7.2*y+135
st
-10<=x<=5
0<=y<=8
2xy+5y-55-13x>=0
帮我看看下面程序有什么错误
%目标函数fun1.m
function f=fun1(x);
f=5.6*x(1)+6.0*x(2)+1.2*x(1)^2+1.2*x(2)^2-1.2*41-36*x(3);
%建立非线性约束条件mycon1.m
function[g,ceq]=mycon(x);
g=[sqrt(x(3)^2+5.0981^2+2*5.0981*x(3)+19.4904^2)+30*x(3)+41-6*x(1)-6*x(2)-x(1)^2-x(2)^2,30*x(3)+41-x(1)^2-x(2)^2]
ceq=[];
x0=[1;1;0];A=[-1,-1];
b=[-15.2971];Aeq=[];beq=[];VLB=[0;0;-10]; VUB=[Inf;Inf;5];
options=optimset;[x,fval]=fmincon(@fun1,x0,A,b,Aeq,beq,VLB,VUB,@mycon1,options)

第一问很简单,楼下都会做,第二问程序如下:
model:
min=5.6*@sqrt((y-5)^2+(x+10)^2)+6*@sqrt((x-5)^2+(y-8)^2)+7.2*y+135;
@bnd(-10,x,5);
@bnd(0,y,8);
2*x*y+5*y-55-13*x>=0 ;
end
运行结果如下:
Local optimal solution found.
Objective value: 250.9731
Extended solver steps: 5
Total solver iterations: 52

Variable Value Reduced Cost
Y 0.4335663 0.000000
X -4.354467 0.000000

Row Slack or Surplus Dual Price
1 250.9731 -1.000000
2 0.000000 -0.2563279E-01
以后有什么问题找我啊,我很乐意帮助新手,想当年我学lingo的时候都没人可以请教的!
温馨提示:答案为网友推荐,仅供参考
第1个回答  2010-09-13
第一问解出来了…
程序如下:
model:
min=5.6*x+6.0*y+1.2*x^2+1.2*y^2-1.2*41-36*m;
x+y>=15.2971;
6*x+6*y+x^2+y^2-41-30*m>=(m^2+5.0981^2+2*5.0981*m+19.4904^2)^1/2;
x^2+y^2-41-30*m>0;
@free(m);
@bnd(-10,m,5);
end
第二问恐怕是目标函数有错,我用lingo一解马上就程序关机了…(⊙o⊙)
第一问答案如下:
Local optimal solution found at iteration: 96
Objective value: 216.3679

Variable Value Reduced Cost
X 8.340160 0.000000
Y 6.956940 0.000000
M -0.9881297 -0.3322783E-07

Row Slack or Surplus Dual Price
1 216.3679 -1.000000
2 0.000000 -1.679351
3 0.000000 -1.055410
4 106.6012 0.000000
相似回答