怎么对文本中的数字进行分组然后运算

Herowinter

翻了楼主的另一个帖子，找到计算E的这些参数的意义了。

1. 如果我想要继续对得到的数据进行处理，比如说对以上得到的数组 [ Yi ]  以L个为单元进行分组， 即第一组为 [  Y1 Y2 Y3 ... YL].    对这Ｌ个数字进行最小二乘法拟合，得 Y=a1 x + a0, (这里的x为数组对应的下标数，即1～ L).
   
2. 计算偏差量  E = sqrt( sum (Yi - (a1*xi +a0) )^2)这就很麻烦了是不是?能帮忙解决一下吗？

复制代码

LikeLx

回复 10# Herowinter
你瘦的不够彻底啊，你看我，完美身材

   

Herowinter

回复 1# wpd0508
不知道结果对不对，有错请指出，把大神提供的最小2乘法封装了下，
用来求数组a中a[1],..a[n]对应的系数a0,a1

1. #!/bin/bash
   
2. 
   
3. awk -vl=5 '
   
4.     function least_square(a,n,a0,a1){
   
5.         for(x=1;x<=n;x++)
   
6.             s += a[x];
   
7.         avg = s/n;
   
8.         for(x=1; x<=n; x++)
   
9.         {
   
10.                 y = a[x] - avg
    
11.                 sx += x
    
12.                 sy += y
    
13.                 sxy += x * y
    
14.                 sx2 += x * x
    
15.         }
    
16.         a0 = (sx2 * sy - sx * sxy) / (n * sx2 - sx * sx)
    
17.         a1 = (n * sxy - sx * sy) / (n * sx2 - sx * sx)
    
18.     }
    
19.       
    
20.     {
    
21.       sum = 0
    
22.       count = 0
    
23.       k = 0
    
24.       for(i=1;i<=NF;i++){
    
25.          a[++count]=$i
    
26.          if(count==l){
    
27.              least_square(a,l,a0,a1);
    
28.              sum=0;
    
29.              count=0;
    
30.              for(j=1;j<=l;j++){
    
31.                  sum += a[j]-(a1*j+a0)^2
    
32.              }
    
33.              E[++k]=sqrt(sum);
    
34.             
    
35.          }
    
36.       }
    
37.       if(count>0){
    
38.           least_square(a,count,a0,a1);
    
39.           sum=0;
    
40.           for(j=1;j<=count;j++){
    
41.                  sum += a[j]-(a1*j+a0)^2
    
42.              }
    
43.           E[++k]=sqrt(sum);
    
44.       }
    
45.       sum=0;
    
46.       for(i in E)
    
47.           sum += E[i]
    
48.       print sum/k
    
49.     }
    
50.    
    
51.    ' i

复制代码

1. cat i
   
2. 1 2 4 5 6 13 141 3 11 13 44 11
   
3. 1 2 4 5 6 7 8 9 10

复制代码

结果

1. ./test.sh
   
2. 8.37082
   
3. 8.7449

复制代码

wpd0508

大大，非常感谢你教我了在LINUX里如何封装函数，sum += a[j]-(a1*j+a0) 这里少了一个平方，还有就是为什么得到的是a1,a0，而不是E的平均数吗？回复 13# Herowinter


   

Herowinter

回复 14# wpd0508
13楼的代码我改了，加了平方，但是计算结果没影响。。。
不知到哪里有问题。按照你的需求，要先算出a0,a1才能代入
公式计算E的吧，我这里最后打印的就是每行的所有E的平均值。

有个数学知识请教下，我输入Y=[1,2,3,4,5]，最小2乘法算出来的
不是a1=1,a0=0,Y=x吗？
   

wpd0508

如果a1=1,a0=1 ,Y=1+x大大还有就是awk 'BEGIN{getline t < "b"}{for(i=0;i++<NFprintf （ "%d\n",i,X=X[i-1]+$i-t）>c}' a  这个代码能产生的数是以一行1个数的形式放入文件c中的,如何改为以空格为分隔符放入文件c中呢
回复 15# Herowinter


   

Herowinter

回复 16# wpd0508
上面我写错是a1=1,a0=0，你的原始文本不是1行
10几个数吗？怎么又变成1行1个了。。。，你重新
描述需求吧，输入数据格式，期望结果。
   

wpd0508

对的，我是想让你帮我看看另外一个问题，
cat test2
1 3 3 5 4 2
test3
3
awk 'BEGIN{getline t < "test3"}{for(i=0;i++<NFprintf ("X%d=%d\n",i,X=X[i-1]+$i-t)> "test4"}' test2
X1=-2
X2=-2
X3=-2
X4=0
X5=1
X6=0
如何将-2 -2 -2 0 1 0 这几个数以这种形式放入test4 ,前面提到的问题就是从这一步开始的。回复 17# Herowinter


   

Herowinter

回复 18# wpd0508
你从test2 => test4 仅仅是为了把
1 2 3 4 5变为

1
2
3
4
5  ???
我是有点被你搞糊涂了，我13楼的代码不需要这种转换，
你直接把里面的文件名改为test2应该可以的。

还有，想要得到正确的结果，必须要保证你提供的最小2乘法
算法能计算出正确的a0,a1，为什么我输入5个点(1,1)  (2,2)
(3,3)  (4,4)  (5,5)，计算的结果不是a0=0,a1=1,  y=x呢？


   

ly5066113

回复 15# Herowinter


把 y = a[x] - avg 改为
y = a[x]
就是你想要的结果了。

楼主原来的需求里是把样本数据减去均值后的数据作为样本的，所以我的代码里减去了均值。