非线性规划模型Python代码.pdfVIP

fromscipyimportoptimizeasopt

importnumpyasnp

fromscipy.optimizeimportminimize

#目标函数

defobjective(x):

returnx[0]**2+x[1]**2+x[2]**2+8

#约束条件

defconstraint1(x):

returnx[0]**2-x[1]+x[2]**2#不等约束

defconstraint2(x):

return-(x[0]+x[1]**2+x[2]**2-20)#不等约束

defconstraint3(x):

return-x[0]-x[1]**2+2

defconstraint4(x):

returnx[1]+2*x[2]**2-3#不等约束

#边界约束

b=(0.0,None)

bnds=(b,b,b)

con1={type:ineq,fun:constraint1}

con2={type:ineq,fun:cons