查看单个帖子
旧 2009-01-10, 19:41   #3
jAc
初级会员
 
注册日期: 2007-12-08
帖子: 6
声望力: 0
jAc 正向着好的方向发展
默认 回复: [求助]12.4037+1065/T-lnT+x1+0.6159/(x2*(x3+T))=0

利用下面的代码,取500组数据{x1,x2,x3},然后去平均值
t = 1; p = {}; While[t < 500, s = {}; n = 1;
While[n <= 3, T = RandomReal[{318, 388}];
s = Insert[s,
12.4037 + 1065/T - Log[T] + x1 + 0.6159/(x2*(x3 + T)) == 0, -1];
n++]; p =
Insert[p, {x1, x2, x3} /. Flatten[Solve[s, {x1, x2, x3}]], -1];
t++]
Mean[p]
这样就得到了平均值{-4.96599, -0.000332984, 50.2058}
做回归要知道拟合类型,你没有给,我就求不出了
jAc 当前离线   回复时引用此帖