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

p141592653

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

p141592653

附上代码：其中p是我指向计算出来的中线的数组，count是有效行，最近向远处计算
void ni_he(int *p,int count)
{
  //x--行    y--列
  int x_pj=0;
  int y_pj=0;
  int x_y_pj=0;
  int x_pf_pj=0;
  int i;
  int a,b;
  for(i=count;i>0;i--)  
  {
    y_pj+=p[h-i];
     x_pj+=i;
     x_pf_pj+=i*i;
     x_y_pj+=i*p[h-i];
  }
  y_pj/=count;
  x_pj/=count;
  x_pf_pj/=count;
  x_y_pj/=count;
  
  b=10*(x_y_pj-y_pj*x_pj)/(x_pf_pj-x_pj*x_pj);
  a=y_pj-b*x_pj;
  
  for(i=h-1;i>=0;i--)     
  {
    if(a+b*i>l-1 || a+b*i<0)
      break;
    mid_path[i]=a+b*i;
  }
}

p141592653

x_pj是行的平均   y_pj是列的平均  x_y_pj是行*列的和的平均  x_pf_pj是行的平方和的平均

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？

CQQ

用高斯曲线求解就行了

登高使者

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

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？

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

p141592653

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

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

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

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

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