应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 求助回归方程一个
222/3返回列表上一页123下一页
查看: 1881  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
10楼: Originally posted by yiyao0828 at 2013-11-24 14:55:21
能详细的阐述下优选后如何进行拟合吗？谢谢...

取y数值大的。
回复此楼
11楼2013-11-24 15:25:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
Parameters p1[0, ], p2[0, ], p3[0, 0.01], p4;
Variable x1, x2, x3, y ;
Function  z=p1*x1+p2*x2+p3*x3+p4;

data;
66        29        5.01        0.39
64.3        28.2        7.48        0.37
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
58.4        33.2        8.35        0.43
55.8        29.8        14.4        0.37
55.1        34.7        10.2        0.36
求助回归方程一个
拟合效果.JPG
赞一下 回复此楼
12楼2013-11-24 15:30:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
本帖小结：
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiyao0828

铁虫 (初入文坛)
引用回帖:
13楼: Originally posted by feixiaolin at 2013-11-24 16:37:19
本帖小结：
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.00149252373 ...

X1 X2 X3 的取值是自己定的吗？如果我想选取50<X1<70,25<X2<35,
5<X3<15可以吗？这样的取值范围运行结果是多少！
赞一下 回复此楼
14楼2013-11-24 21:42:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
14楼: Originally posted by yiyao0828 at 2013-11-24 21:42:21
X1 X2 X3 的取值是自己定的吗？如果我想选取50<X1<70,25<X2<35,
5<X3<15可以吗？这样的取值范围运行结果是多少！...

看12楼data中数据。
想选取50<X1<70,25<X2<35, 5<X3<15也可以的。
用13楼中代码即可。
赞一下 回复此楼
15楼2013-11-24 21:49:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiyao0828

铁虫 (初入文坛)
引用回帖:
15楼: Originally posted by feixiaolin at 2013-11-24 21:49:00
看12楼data中数据。
想选取50<X1<70,25<X2<35, 5<X3<15也可以的。
用13楼中代码即可。...

能帮忙运算一下吗？我实在是不懂，谢谢！非常感谢！
赞一下 回复此楼
16楼2013-11-24 21:57:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiyao0828

铁虫 (初入文坛)
引用回帖:
15楼: Originally posted by feixiaolin at 2013-11-24 21:49:00
看12楼data中数据。
想选取50<X1<70,25<X2<35, 5<X3<15也可以的。
用13楼中代码即可。...

大哥，帮帮忙吧，就最后一步了！非常的感谢你！
赞一下 回复此楼
17楼2013-11-24 22:19:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
17楼: Originally posted by yiyao0828 at 2013-11-24 22:19:08
大哥，帮帮忙吧，就最后一步了！非常的感谢你！...

54.0424666
33.16403282
14.46068189

0.644272281
赞一下 回复此楼
18楼2013-11-24 22:33:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiyao0828

铁虫 (初入文坛)
引用回帖:
18楼: Originally posted by feixiaolin at 2013-11-24 22:33:21
54.0424666
33.16403282
14.46068189

0.644272281...

X1+X2+X3的总和超出总量100好多呀！
赞一下 回复此楼
19楼2013-11-24 22:46:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiyao0828

铁虫 (初入文坛)
引用回帖:
18楼: Originally posted by feixiaolin at 2013-11-24 22:33:21
54.0424666
33.16403282
14.46068189

0.644272281...

能告诉我你用的什么软件吗？好像挺好用的啊
赞一下 回复此楼
20楼2013-11-24 22:48:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yiyao0828 的主题更新
222/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] filecode=后面第一个是大写字母 +7 wangze12014 2026-08-14 9/450 2026-08-17 14:59 by wangze12014
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +33 医学老男孩 2026-08-13 80/4000 2026-08-17 14:40 by yq199112088
[基金申请] 感觉是下周放榜了 +3 angus9576 2026-08-17 6/300 2026-08-17 14:32 by 启萌科技
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 EWi2p09MOOv8 2026-08-16 6/300 2026-08-17 14:02 by EYK67XJMj64E
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6GojgJvkDudM 2026-08-16 6/300 2026-08-17 13:58 by EYK67XJMj64E
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 11:46 by 4GBAYCdVQoK3
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
关闭
请填处理意见
关闭