lingo求救!!!!!!高手来啊!!!!

来源:百度知道 编辑:UC知道 时间:2024/07/02 04:51:42
想用lingo求最优值,非线性的,例如y=x*exp(-a*x-d*x2-b)/(exp(-a*x-d*x2-b)+c);,
a,b,c都是常数,这很容易实现,就是我有一连串对应的x等变量,求出对应的一连串y;用数组或什么存起来,由于对应的变量太多``一个一个的求最优值难啊,```高手救啊啊```或者有没有其它软件能办到的``重谢!!!!

liu.txt
4
5
8
7
1
2
47
7
5
3
1
2
1
4
45

model:
sets:
liu/1..15/:x,y;
endsets
data:
x=@text('f:\liu.txt');
a=1;
b=-452;
c=103000000000;
d=5;
enddata
submodel ll:
@for(liu:y-x*@exp(-a*x-d*x^2-b)/(@exp(-a*x-d*x^2-b)+c)=0);
endsubmodel
calc:
@divert('f:\liu1.txt');
@solve(ll);
@writefor(liu(j):y(j),@newline(1));
@divert();
endcalc
end

liu1.txt

4
5
8
7
1
2
0
7
5
3
1
2
1
4
0