| 查看: 1882 | 回复: 21 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
yiyao0828铁虫 (初入文坛)
|
[求助]
求助回归方程一个
|
||||
|
求助一个回归模型： X1=68.6， 66.8， 66， 64.3， 62.8， 62， 61.3， 60.5， 59.8， 59.1， 58.4， 57.1， 55.8， 55.1， 53.9 X2=25.6， 25， 29， 28.2， 27.1， 31.4， 25.8， 30.1， 34.6， 28.7， 33.2， 31.6， 29.8， 34.7， 32.7 X3=5.75， 8.14， 5.01， 7.48， 10.1， 6.61， 13， 9.34， 5.56， 12.2， 8.35， 11.3， 14.4， 10.2， 13.4 Y =0.34， 0.30， 0.39， 0.37， 0.31， 0.37， 0.45， 0.37， 0.38， 0.33， 0.43， 0.35， 0.37， 0.36， 0.32 其中X1+X2+X3=100 要求回归模型具有极大值，每列值是相对应的，比如：当X1=68.6 X2=25.6 X3=5.75 Y=0.34。能给出极值就更好了，回归方程中能尽量不出现X3！ |
回复此楼» 猜你喜欢
filecode=后面第一个是大写字母
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有81人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
感觉是下周放榜了
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
|
本帖小结： 1) 7楼 MaxFunction y/Abs(p1*x1+p2*x2+p3*x4+p4-y); 应更正为：MaxFunction y/Abs(p1*x1+p2*x2+p3*x3+p4-y); 2) 7楼 或 12楼之后，得到拟合公式 y=x1*p1+x2*p2+x3*p3+p4 =x1*0.00149252373465775+x2*0.00113838998918746+x3*0.0275918372945969+0.126862232328834 ************* 这才成功了一半，关键要看第三步 *************** 3) Mont carlo寻优 代码如下 Constant p1=0.00149252373465775, p2=0.00113838998918746, p3=0.0275918372945969, p4=0.126862232328834; Parameter x1[50, 60], x2[20, 35], x3[10, 20], y; MaxFunction y=p1*x1 + x2*p2 + p3*x3 +p4; -0.1<x1+x2+x3-100<0.1; ---------------------------------------------------------------------------- 运行结果： x1=59.4199134740019 ; x2=25.2562755846884 ; x3=15.2252918648651 y= 0.664393130687968 the max ! |
13楼2013-11-24 16:37:19
peterflyer 
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 146666
- 红花: 1374
- 帖子: 93091
- 在线: 7695小时
- 虫号: 1482829
- 注册: 2011-11-08
- 性别: GG
- 专业: 功能陶瓷
2楼2013-11-24 10:35:30
yiyao0828
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 83.5
- 散金: 12
- 帖子: 19
- 在线: 7.9小时
- 虫号: 2754556
- 注册: 2013-10-26
- 性别: MM
- 专业: 有机高分子功能材料
3楼2013-11-24 12:26:39
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
yiyao0828: 金币+8 2013-11-24 21:43:32
yiyao0828: 金币+7 2013-11-24 22:07:43
感谢参与，应助指数 +1
yiyao0828: 金币+8 2013-11-24 21:43:32
yiyao0828: 金币+7 2013-11-24 22:07:43
|
1stop代码 Parameters p1, p2, p3, p4; Variable x1,x2,x3,y ; MinFunction p1*x1+p2*x2+p3*x4+p4-y; x1+x2+x3=100; Data; 68.6 25.6 5.75 0.34 66.8 25 8.14 0.3 66 29 5.01 0.39 64.3 28.2 7.48 0.37 62.8 27.1 10.1 0.31 62 31.4 6.61 0.37 61.3 25.8 13 0.45 60.5 30.1 9.34 0.37 59.8 34.6 5.56 0.38 59.1 28.7 12.2 0.33 58.4 33.2 8.35 0.43 57.1 31.6 11.3 0.35 55.8 29.8 14.4 0.37 55.1 34.7 10.2 0.36 53.9 32.7 13.4 0.32 |
4楼2013-11-24 12:41:45
