求一組數(shù)的偏差和累加

wpd0508

已知文本1中有1 3 3 5 4 2 這組數(shù) 文本2存放著他們的平均數(shù)3 然后通過SHELL腳本如何實現(xiàn)求出各數(shù)的偏差和累加。X1=$1-3 ,X2=X1+$2-3......

yestreenstars

醬紫？

1. awk 'BEGIN{getline t < "b"}{for(i=0;i++<NF;)printf "X%d=%d\n",i,X[i]=X[i-1]+$i-t}' a
   
2. X1=-2
   
3. X2=-2
   
4. X3=-2
   
5. X4=0
   
6. X5=1
   
7. X6=0

復(fù)制代碼

wpd0508

大大厲害，謝了，順便問一下其中的a,b 是指的文件名嗎？
回復(fù) 2# yestreenstars


   

wpd0508

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


   

yestreenstars

回復(fù) 3# wpd0508
是的~{:2_172:}

   

yestreenstars

回復(fù) 4# wpd0508
沒看懂~{:2_169:}

   

q1208c

這么復(fù)雜的計算, shell可能處理不了.

就算是能處理, 你以后維護的成本也會非常高.

建議使用 perl 或 python 吧.

yestreenstars

回復(fù) 7# q1208c
你是來抬杠的~{:2_172:}

   

wpd0508

但是我已經(jīng)得到了YI這一組數(shù)，a1和a0是隨便選的2個數(shù)，只需要進行最后一步的運算E = sqrt( sum (Yi - (a1*xi +a0) )^2)比如說E = sqrt( sum (Yi - (1*xi +2) )^2) 這也不能實現(xiàn)嗎，真的很需要解決這個問題，求大大們點撥一下。

wpd0508

a0 = （∑Yi) / n - a1（∑Xi) / n （式1-
a1 = [n∑Xi Yi - （∑Xi ∑Yi)] / [n∑Xi2 - （∑Xi)2 )]   我錯了，a0 和a1居然是2個值，我去他妹的。