版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

rita128

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 96.9
帖子: 21
在线: 10.3小时
虫号: 1129918
注册: 2010-10-23
性别: MM
专业: 无机合成和制备化学

[求助] 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);

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有49人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
三无产品还有机会吗已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

麻烦高人指点一下1stopt和spss的非线性回归中weibull的两参数的公式应该怎么表达啊？已经有9人回复
关于用lstopt进行非线性最小二乘法回归问题已经有5人回复

1楼 2013-04-20 17:08:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 96.9
帖子: 21
在线: 10.3小时
虫号: 1129918
注册: 2010-10-23
性别: MM
专业: 无机合成和制备化学

引用回帖:

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

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

不能保证你用破解版的能运行成功。

赞一下

回复此楼

4楼2013-04-21 14:08:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 rita128 的主题更新

返回列表