- 注册时间
- 2010-7-23
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 851
- 在线时间
- 小时
|
发表于 2011-2-24 21:37:07
|
显示全部楼层
楼主的题目,我用Forcal算了个似乎是更优的值(实际也差不多):- !using["fcopt","math"];
- y(p1,p2,p3,p4,x)=(p1/(ln(p2/x)))^(x/p3)*p4;
- f(p1,p2,p3,p4) = [y(p1,p2,p3,p4,0.25)-4.28]^2 + [y(p1,p2,p3,p4,0.5)-6.63]^2 + [y(p1,p2,p3,p4,0.6)-7.61]^2 + [y(p1,p2,p3,p4,0.8)-14.55]^2 + [y(p1,p2,p3,p4,0.9)-22.77]^2;
- Opt[HFor("f"), optmax,1000, optmode,100, optdeep,50, optwaysimdeep, optwayconfra];
复制代码 结果(可能要求解几次):
-2.665064291582835 0.1469425133157456 -3.353692523962099e-002 710981.3329454004 2.387625927609502e-002 |
|