有2個公式的文本數(shù)據(jù)處理問題？

tengfei0311

如何用bash shell腳本處理一下文件：

1. 
   
2.        第1列        第2列      第3列(C0)  第4列(A1)   第5列(A2)
   
3. 文件:1.dat
   
4.         90.0000   33.0000    3.0841      0.0000       -0.0052   
   
5.         90.5000   33.0000    3.0787      0.0003       -0.0071   
   
6.         91.0000   33.0000    3.0694      0.0005       -0.0094   
   
7.         91.5000   33.0000    3.0564      0.0006       -0.0116   
   
8.         (共2000行)
   

復制代碼

根據(jù)上面一個文件1.dat，對第4列和第5列經(jīng)過公式轉換成下面一個文件，前三列值不變，第4,5列的A和θ計算公式為：
要得到下面文件：

1. 
   
2. 第1列    第2列    第3列(C0)              第4列                 第5列
   
3. 90.0000   33.0000    3.0841           A                        θ   
   
4. 90.5000   33.0000    3.0787           A                        θ
   
5. 91.0000   33.0000    3.0694           A                        θ  
   
6. 91.5000   33.0000    3.0564           A                        θ
   
7. 92.0000   33.0000    3.0430           A                        θ
   

復制代碼

想用bash shell來實現(xiàn)這個文本處理功能，菜鳥我不知道該怎么處理，謝謝各位！

bikkuri

1. [root@hp8:~]# cat file
   
2.         90.0000   33.0000    3.0841      0.0000       -0.0052
   
3.         90.5000   33.0000    3.0787      0.0003       -0.0071
   
4.         91.0000   33.0000    3.0694      0.0005       -0.0094
   
5.         91.5000   33.0000    3.0564      0.0006       -0.0116
   
6. [root@hp8:~]# cat file|awk '{a1=$1;a2=$2;m=a1*a1+a2*a2;n=a1/a2;cmd1="echo \"scale=4; sqrt("m")\"|bc -l";cmd1|getline a;cmd2="echo \"scale=4; 0.5*a("n")\"|bc -l";cmd2|getline t;print$1,$2,$3,a,t}'
   
7. 90.0000 33.0000 3.0841 95.8592 .6096
   
8. 90.5000 33.0000 3.0787 96.3288 .6105
   
9. 91.0000 33.0000 3.0694 96.7987 .6114
   
10. 91.5000 33.0000 3.0564 97.2689 .6123
    
11. [root@hp8:~]#
    

復制代碼

回復 1# tengfei0311


   

jason680

回復 1# tengfei0311

$ awk '{OFS="\t";a=$4;b=$5;$4=sqrt(a*a+b*b);$5=0.5*atan2(a,b)}1' 1.dat
90.0000        33.0000        3.0841        0.0052        1.5708
90.5000        33.0000        3.0787        0.00710634        1.54968
91.0000        33.0000        3.0694        0.00941329        1.54423
91.5000        33.0000        3.0564        0.0116155        1.54496

   

tengfei0311

回復 2# bikkuri

謝謝您的幫助 好像你求出的四五列值不對啊   

tengfei0311

回復 3# jason680


謝謝您的幫助，中間有一點不懂。

1. $5=0.5*atan2(a,b)}1

復制代碼

這里的1是什么意思呢？
   

jason680

回復 5# tengfei0311

http://72891.cn/thread-2309494-1-1.html
awk初學之常見問題
1. awk ‘{code}1’ 中的“1”是干什么的？
...
   

tengfei0311

回復 3# jason680

您好，我想問一下，0.5*atan2(a,b) 用這個atan求出來的值是不是一個弧度值？  這里我想要一個角度值。用arctan求，就是求反函數(shù) 。這里這樣寫是不是不對？


   

jason680

回復 7# tengfei0311

弧度與角度

2 * pi = 2 * 3.1415926 = 圓 = 360度
      pi =  3.1415926 = 半圓 = 180度

弧度 = 角度 / 180 * pi
角度 = 弧度 / pi * 180

tengfei0311

回復 8# jason680


你這樣就是對的嗎？我只是看見別人論文里面都是arctan或者tan-1  我不知道atan跟arctan本質區(qū)別。

我看見別人用matlab寫的code是這樣的形式

1.     MagAzim(i,2) = 0.5*atan2(cdata(i,6),cdata(i,5));

復制代碼

jason680

回復 9# tengfei0311


https://www.gnu.org/software/gawk/manual/gawk.html#Numeric-Functions

9.1.2 Numeric Functions

The following list describes all of the built-in functions that work with numbers. Optional parameters are enclosed in square brackets ([ ]):

atan2(y, x)

    Return the arctangent of y / x in radians. You can use ‘pi = atan2(0, -1)’ to retrieve the value of pi.

cos(x)

    Return the cosine of x, with x in radians.

...