有2个公式的文本数据处理问题？

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行)
   

复制代码

根据上面一个文件1.dat，对第4列和第5列经过公式转换成下面一个文件，前三列值不变，第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来实现这个文本处理功能，菜鸟我不知道该怎么处理，谢谢各位！

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://bbs.chinaunix.net/thread-2309494-1-1.html
awk初学之常见问题
1. awk ‘{code}1’ 中的“1”是干什么的？
...
   

tengfei0311

回复 3# jason680

您好，我想问一下，0.5*atan2(a,b) 用这个atan求出来的值是不是一个弧度值？  这里我想要一个角度值。用arctan求，就是求反函数 。这里这样写是不是不对？


   

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本质区别。

我看见别人用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.

...