资源

解决大规模NLP问题(nlps4)

/*****************************************************************//*                                                               *//*S A S S A M P L E L I B R A R Y*//*                                                               *//*名称:nlps4*//*标题:解决大规模NLP问题(nlps4)*//*产品:或*//*系统:所有*//*关键字:或*//*过程:OPTMODEL*//*数据:*//*                                                               *//*支持:更新:*//*参考:*//*MISC:非线性规划求解器示例4*//*数学程序设计一章*//*                                                               *//*****************************************************************/过程操作模型;num nx=100;数量ny=100;变量x{1..nx,0..ny+1}>=0;变量y{0..nx+1,1..ny}>=0;最小f=(和{i在1..nx-1中,j在1..ny-1}中(x[i,j]-1)^2+和{i在1..nx-1中,j在1..ny-1}(y[i,j]-1)^2+和{i在1..nx-1}(x[i,ny]-1)^2+1..ny-1}中的和{j(y[nx,j]-1)^2) / 2;con con1{i在2..nx-1中,j在2..ny-1}中:(x[i,j]-x[i-1,j])+(y[i,j]-y[i、j-1])=1;con con2{i in 2..nx-1}:x[i,0]+(x[i、1]-x[i-1,1])+y[i,1]=1;con con3{i in 2..nx-1}:x[i,ny+1]+(x[i、ny]-x[i-1,ny])-y[i,ny-1]=1;concon4{jin2..ny-1}:y[0,j]+(y[1,j]-y[1,j-1])+x[1,j]=1;con5{j in 2..ny-1}:y[nx+1,j]+(y[nx,j]-y[nx,j-1])-x[nx-1,j]=1;对于{i在1..nx-1}x[i,ny].lb=1;对于1..ny-1}y[nx,j].lb=1中的{j;用nlp求解;退出;