有谁用了线性拟合的，进来交流一下呢

p141592653

y=a+bx，x我是用行来代入，输出的y就是改行所对应的列，但问题在于，利用计算后得出的a和b代入后，算出的y不是很大就是很小，擦，求大神

vic_hong

p141592653 发表于 2013-6-19 13:44
你的是说讲Exy乘上n，E(x^2)乘上n，其他都不用乘n是不是这样子。然后你说x,y对调，是不是将x作为列输出， ...

对的。反过来就是为了避免垂直时的无限大嘛。公式可以从你的公式全部乘以n^2转换而来。

p141592653

vic_hong 发表于 2013-6-19 13:15
我觉得b放大了10倍，后面的a求解就不对了。不要用这个平均的公式，因为每个和先平均会有精度问题。用(n*E ...

x,y对调，是不是将x作为列输入，然后y是行输出

p141592653

vic_hong 发表于 2013-6-19 13:15
我觉得b放大了10倍，后面的a求解就不对了。不要用这个平均的公式，因为每个和先平均会有精度问题。用(n*E ...

你的是说讲Exy乘上n，E(x^2)乘上n，其他都不用乘n是不是这样子。然后你说x,y对调，是不是将x作为列输出，然后y是行输出

vic_hong

p141592653 发表于 2013-6-18 23:35
额，是因为输出的值不在数组范围才加的。。。

我觉得b放大了10倍，后面的a求解就不对了。不要用这个平均的公式，因为每个和先平均会有精度问题。用(n*Exy-Ex*Ey)/(n*E(x^2)-(Ex)^2)。还有用int存斜率是不是有精度问题。 还有你要注意的是垂直线的问题，一般来说你要求斜率的倒数，这时候你要把x，y对调。

p141592653

CQQ 发表于 2013-6-18 22:32
用高斯曲线求解就行了

是不是这个方程y=1/(0.4sqrt(2pi))e^(-0.5((x-1)/0.4)^2)

p141592653

vic_hong 发表于 2013-6-18 22:23
b=10*(x_y_pj-y_pj*x_pj)/(x_pf_pj-x_pj*x_pj);
  a=y_pj-b*x_pj;
乘以10？

额，是因为输出的值不在数组范围才加的。。。

登高使者

请问楼主是用在哪里的呢？是优化小s的吧？我试过最小二乘法。 我的情况跟你一样！

CQQ

用高斯曲线求解就行了

vic_hong

b=10*(x_y_pj-y_pj*x_pj)/(x_pf_pj-x_pj*x_pj);
  a=y_pj-b*x_pj;
乘以10？