24小时热门版块排行榜    

查看: 1681  |  回复: 3

rita128

新虫 (初入文坛)

[求助] 1stopt做最小二乘法回归参数,请求大家帮忙指导下

第一次用这个做最小二乘法回归参数。一共六个参数。有八组数据,目标方程一个。应该是最简单的那种吧。代码如下,运行不出来,不知道怎么定义常量或者输入数据。求助懂的人帮忙指点一下。谢谢
Title "Type your title here"i;
Parameters x1,x2,x3,x4,x5,x6 ;
Constant a(1:8)= [0.593492611, 0.615383905,0.605448133,0.692841215,0.655361944,0.757020028,0.671381587,0.735878492];
Constant b(1:8)=[ 5.1612, 4.6650 ,3.7857,6.7947,2.3167 ,1.4693,7.0426 ,10.0069];


Conststr I=3*a;

Conststr f1=-0.76*(I^0.5/(1+1.2*I^0.5)+2/1.2*ln(1+1.2*I^0.5));//f'(i)
Conststr B1=-0.00804+1.4354/(2*I)*(1-(1+2*I^0.5)*EXP((-2)*I^0.5));//B
Conststr B2=1.4354/(2*I^2)*(-1+(1+2*I^0.5+2*I)*exp((-2)*I^0.5));//B'
Conststr F1=1/(2*I)*(1-(1+2*I^0.5)*exp((-2)*I^0.5));//F(i)
Conststr F2=1/(2*I^2)*(-1+(1+2*I^0.5+2*I)*EXP((-2)*I^0.5));//F'()i
MinFunction F=sum(i=1:8)(((f1+8/3*a*(B1+2*a*0.001838)+2*a^2*(2*B2+4/3*0.001838)+b*x3/3+4/3*b*a*(x1+F1*x2)+8/3*a*b*(x1+F1*x2)+4*a^2*b*F2*x2)+ln(a)+1.5*ln(b)+1.5*(a*x3+4*a^2*(x1+F1*x2)+2*b*x4+3*b^2*x5)-lnx6)^2);
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
rita128: 金币+20, ★★★★★最佳答案, 哦,非常感谢你,让我学了好多东西。要不真的没有接触这个软件。遇到好多问题,真的太谢谢了。 2013-04-21 00:02:42
ben_ladeng: 金币+2, 谢谢应助 2013-04-21 20:51:23
ben_ladeng: 专家考核 2013-04-21 20:51:32
几点错误:
1:“lnx6”应该是“ln(x6)”吧
2:1stOpt不区分大小写
3:“Sum”的写法有问题

正确代码如下:
CODE:
Title "Type your title here"i;
Parameters x1,x2,x3,x4,x5,x6 ;
Constant a(1:8)= [0.593492611, 0.615383905,0.605448133,0.692841215,0.655361944,0.757020028,0.671381587,0.735878492];
Constant b(1:8)=[ 5.1612, 4.6650 ,3.7857,6.7947,2.3167 ,1.4693,7.0426 ,10.0069];


Conststr I=3*a;

Conststr f1=-0.76*(I^0.5/(1+1.2*I^0.5)+2/1.2*ln(1+1.2*I^0.5));//f'(i)
Conststr BB1=-0.00804+1.4354/(2*I)*(1-(1+2*I^0.5)*EXP((-2)*I^0.5));//B
Conststr BB2=1.4354/(2*I^2)*(-1+(1+2*I^0.5+2*I)*exp((-2)*I^0.5));//B'
Conststr FF1=1/(2*I)*(1-(1+2*I^0.5)*exp((-2)*I^0.5));//F(i)
Conststr F2=1/(2*I^2)*(-1+(1+2*I^0.5+2*I)*EXP((-2)*I^0.5));//F'()i
MinFunction sum(j=1:8,a,b)(((f1+8/3*a*(BB1+2*a*0.001838)+2*a^2*(2*BB2+4/3*0.001838)+b*x3/3+4/3*b*a*(x1+FF1*x2)+8/3*a*b*(x1+FF1*x2)+4*a^2*b*F2*x2)+ln(a)+1.5*ln(b)+1.5*(a*x3+4*a^2*(x1+FF1*x2)+2*b*x4+3*b^2*x5)-ln(x6))^2);

结果:
目标函数值(最小): 0.000360743527320632
x1: -0.140392701334565
x2: 0.216320938632417
x3: 2.29529578782098
x4: -0.370075186018573
x5: 0.00657741865607194
x6: 1.24296439011043
2楼2013-04-20 20:24:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rita128

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2013-04-20 20:24:12
几点错误:
1:“lnx6”应该是“ln(x6)”吧
2:1stOpt不区分大小写
3:“Sum”的写法有问题

正确代码如下:

Title "Type your title here"i;
Parameters x1,x2,x3,x4,x5,x6 ;
Constant a(1:8) ...

请问我这个代码是任何版本都能运行吗,我把这个代码重新运行点了计算之后没有任何反应。。。。是什么原因呢?
3楼2013-04-21 11:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

不能保证你用破解版的能运行成功。
4楼2013-04-21 14:08:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rita128 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 u1xkze8ybe 2026-07-30 7/350 2026-07-31 12:10 by syhhffkrym
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 3/150 2026-07-31 12:08 by syhhffkrym
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 6/300 2026-07-31 12:05 by syhhffkrym
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +6 锦衣卫寒战 2026-07-28 6/300 2026-07-31 11:27 by syl200707
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 小木虫看见有人已经知道结果了 +15 1234567wang 2026-07-25 19/950 2026-07-31 08:30 by ZJTJZ
[基金申请] 2026年国自然面上资助率 +8 布布和一二 2026-07-30 8/400 2026-07-31 01:51 by alongwaytogo
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 17/850 2026-07-30 20:10 by kissu88
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 没消息就是被刷了呗 +13 lambert2014 2026-07-24 14/700 2026-07-29 10:38 by msh145800
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见