24小时热门版块排行榜    

查看: 1706  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 让我中一个面上吧! +8 大萍1987 2026-08-20 10/500 2026-08-20 21:48 by 大萍1987
[基金申请] 言之凿凿,明天官网会通知放榜? +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见