怎么對文本中的數(shù)字進行分組然后運算

Herowinter

回復 20# ly5066113
多謝大神指點，現(xiàn)在(i i) i=1,2,...5,得到的結果是y=x了。

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

復制代碼

1. ./test.sh
   
2. ###y = -0.3 + 1.3x
   
3. ###y = 75.2 + -13x
   
4. ###y = 77 + -33x
   
5. 8.37082
   
6. ###y = 0 + 1x
   
7. 3.872988
   

復制代碼

1. cat test2
   
2. 1 2 4 5 6 13 141 3 11 13 44 11
   
3. 1 2 3 4 5

復制代碼

wpd0508

是為了將1 3 3 5 4 2  先求均值，然后求出與均值的偏差X1，X2，X3，X4，X5，X6，最后累加 Y1=X1 ,Y2=Y1+X2,...但是我用awk 'BEGIN{getline t < "test3"}{for(i=0;i++<NFprintf ("X%d=%d\n",i,Y=Y[i-1]+$i-t) > "test4"}' test2得到的是
Y1=-2
Y2=-2



這樣的形式，而我想的是將-2  -2  1  。以這種形式放入test4中，然后進行最小二乘法回復 19# Herowinter [/

   

wpd0508

解決了 ，我去原來/N在shell中也表示換行，，
我真是愚，謝謝各位大大們這幾天的幫忙。。 回復 22# wpd0508


   

wpd0508

解決了 ，我去原來/N在shell中也表示換行，，
我真是愚，謝謝各位大大們這幾天的幫忙。。 回復 22# wpd0508


   

wpd0508

解決了 ，我去原來/N在shell中也表示換行，，
我真是愚，謝謝各位大大們這幾天的幫忙。。 回復 22# wpd0508


   

Herowinter

回復 22# wpd0508

1. [masonzhu@rat122 zzz]$ awk 'BEGIN{getline t < "test3"}{for(i=1;i<=NF;i++){Y[i]=Y[i-1]+$i-t;s=s?s" "Y[i]:Y[i]};print s > "test4"}' test2
   
2. [masonzhu@rat122 zzz]$ cat test4
   
3. -2 -2 -2 0 1 0
   
4. 
   
5. [masonzhu@rat122 zzz]$ cat test2
   
6. 1 3 3 5 4 2
   
7. [masonzhu@rat122 zzz]$ cat test3
   
8. 3

復制代碼