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

Herowinter

翻了樓主的另一個(gè)帖子，找到計(jì)算E的這些參數(shù)的意義了。

1. 如果我想要繼續(xù)對(duì)得到的數(shù)據(jù)進(jìn)行處理，比如說(shuō)對(duì)以上得到的數(shù)組 [ Yi ]  以L個(gè)為單元進(jìn)行分組， 即第一組為 [  Y1 Y2 Y3 ... YL].    對(duì)這Ｌ個(gè)數(shù)字進(jìn)行最小二乘法擬合，得 Y=a1 x + a0, (這里的x為數(shù)組對(duì)應(yīng)的下標(biāo)數(shù)，即1～ L).
   
2. 計(jì)算偏差量  E = sqrt( sum (Yi - (a1*xi +a0) )^2)這就很麻煩了是不是?能幫忙解決一下嗎？

復(fù)制代碼

LikeLx

回復(fù) 10# Herowinter
你瘦的不夠徹底啊，你看我，完美身材

   

Herowinter

回復(fù) 1# wpd0508
不知道結(jié)果對(duì)不對(duì)，有錯(cuò)請(qǐng)指出，把大神提供的最小2乘法封裝了下，
用來(lái)求數(shù)組a中a[1],..a[n]對(duì)應(yīng)的系數(shù)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

復(fù)制代碼

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

復(fù)制代碼

結(jié)果

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

復(fù)制代碼

wpd0508

大大，非常感謝你教我了在LINUX里如何封裝函數(shù)，sum += a[j]-(a1*j+a0) 這里少了一個(gè)平方，還有就是為什么得到的是a1,a0，而不是E的平均數(shù)嗎？回復(fù) 13# Herowinter


   

Herowinter

回復(fù) 14# wpd0508
13樓的代碼我改了，加了平方，但是計(jì)算結(jié)果沒(méi)影響。。。
不知到哪里有問(wèn)題。按照你的需求，要先算出a0,a1才能代入
公式計(jì)算E的吧，我這里最后打印的就是每行的所有E的平均值。

有個(gè)數(shù)學(xué)知識(shí)請(qǐng)教下，我輸入Y=[1,2,3,4,5]，最小2乘法算出來(lái)的
不是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  這個(gè)代碼能產(chǎn)生的數(shù)是以一行1個(gè)數(shù)的形式放入文件c中的,如何改為以空格為分隔符放入文件c中呢
回復(fù) 15# Herowinter


   

Herowinter

回復(fù) 16# wpd0508
上面我寫錯(cuò)是a1=1,a0=0，你的原始文本不是1行
10幾個(gè)數(shù)嗎？怎么又變成1行1個(gè)了。。。，你重新
描述需求吧，輸入數(shù)據(jù)格式，期望結(jié)果。
   

wpd0508

對(duì)的，我是想讓你幫我看看另外一個(gè)問(wèn)題，
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 這幾個(gè)數(shù)以這種形式放入test4 ,前面提到的問(wèn)題就是從這一步開始的。回復(fù) 17# Herowinter


   

Herowinter

回復(fù) 18# wpd0508
你從test2 => test4 僅僅是為了把
1 2 3 4 5變?yōu)?br />
1
2
3
4
5  ???
我是有點(diǎn)被你搞糊涂了，我13樓的代碼不需要這種轉(zhuǎn)換，
你直接把里面的文件名改為test2應(yīng)該可以的。

還有，想要得到正確的結(jié)果，必須要保證你提供的最小2乘法
算法能計(jì)算出正確的a0,a1，為什么我輸入5個(gè)點(diǎn)(1,1)  (2,2)
(3,3)  (4,4)  (5,5)，計(jì)算的結(jié)果不是a0=0,a1=1,  y=x呢？


   

ly5066113

回復(fù) 15# Herowinter


把 y = a[x] - avg 改為
y = a[x]
就是你想要的結(jié)果了。

樓主原來(lái)的需求里是把樣本數(shù)據(jù)減去均值后的數(shù)據(jù)作為樣本的，所以我的代碼里減去了均值。