求高手给下面非线性规划的lingo的代码
来源:学生作业帮 编辑:搜狗做题网作业帮 分类:综合作业 时间:2024/07/03 15:03:53
求高手给下面非线性规划的lingo的代码
![](http://img.wesiedu.com/upload/c/43/c43cee3587d639059ce47a682ee2c69f.jpg)
![](http://img.wesiedu.com/upload/c/43/c43cee3587d639059ce47a682ee2c69f.jpg)
![求高手给下面非线性规划的lingo的代码](/uploads/image/z/8497701-45-1.jpg?t=%E6%B1%82%E9%AB%98%E6%89%8B%E7%BB%99%E4%B8%8B%E9%9D%A2%E9%9D%9E%E7%BA%BF%E6%80%A7%E8%A7%84%E5%88%92%E7%9A%84lingo%E7%9A%84%E4%BB%A3%E7%A0%81)
sets:
row/1..7/:x,p,a,b;
col/1..3/:r;
endsets
max=@sum(row(i):p(i)*(1-@exp(-(x(i)/(6.30589*10^5))^2)));
@for(col(j):x(j)=(1-@sum(row(i)|i#ge#4:p(i)*x(i)))*r(j)/p(j));
@sum(col:r)=1;
@bnd(0.5,r(1),0.8);
@bnd(6*10^5,x(1),5*10^6);
@for(row(i)|i#lt#@size(row):a(i)=是一样的,都是大于等于的意思,所以用这种方式来表示rj大于0,在合理的情况下也可以用一个充分小的正数来限制,例如r>=10^-8;
row/1..7/:x,p,a,b;
col/1..3/:r;
endsets
max=@sum(row(i):p(i)*(1-@exp(-(x(i)/(6.30589*10^5))^2)));
@for(col(j):x(j)=(1-@sum(row(i)|i#ge#4:p(i)*x(i)))*r(j)/p(j));
@sum(col:r)=1;
@bnd(0.5,r(1),0.8);
@bnd(6*10^5,x(1),5*10^6);
@for(row(i)|i#lt#@size(row):a(i)=是一样的,都是大于等于的意思,所以用这种方式来表示rj大于0,在合理的情况下也可以用一个充分小的正数来限制,例如r>=10^-8;