| 查看: 1483 | 回复: 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); |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有4人回复
-
2025年遐想
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
-
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
-
求个博导看看
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
rita128
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 96.9
- 帖子: 21
- 在线: 10.3小时
- 虫号: 1129918
- 注册: 2010-10-23
- 性别: MM
- 专业: 无机合成和制备化学
3楼2013-04-21 11:44:41