设为首页收藏本站

LinuxSir.cn，穿越时空的Linuxsir!

 找回密码
 注册
搜索
热搜: shell linux mysql
LinuxSir.cn，穿越时空的Linuxsir!»论坛 › Linux 综合讨论区 —— LinuxSir.cn › shell进阶应用、shell编程 › 请问如下一段 sed或 awk代码怎么写？对我来说，比较难哦 ...
123下一页
返回列表 发新帖
查看: 2315|回复: 24

请问如下一段 sed或 awk代码怎么写？对我来说，比较难哦

[复制链接]
发表于 2007-1-19 09:58:07 | 显示全部楼层 |阅读模式
如题，我的目录结构如下：


  2. w_3.9039nm/
  3.           ----->p_2.50e+01kpa/result.txt
  4.                 p_4.00e+01kpa/result.txt
  5.                 p_5.50e+01kpa/result.txt
  6.                 p_7.00e+01kpa/result.txt

  8. w_4.2588nm/
  9.           ----->p_2.50e+01kpa/result.txt
  10.                 p_4.00e+01kpa/result.txt
  11.                 p_5.50e+01kpa/result.txt
  12.                 p_7.00e+01kpa/result.txt
复制代码


其中的result.txt为不同的w和p下的计算数据，格式相同，但数字不同，举例如下：


  3.       rhosa=  0.2879                  drhos=  0.0014
  4.        nksa=  453.44                   dnks=  2.1915
  5.       ucmsa=   14.02                  ducms=  0.0476
  6.       ucfsa=    1.79                  ducfs=  0.0124
  7.       ucpsa=   12.23                  ducps=  0.0519
  8.     dhuffsa=  6.3425 kj/mol         ddhuffs=  0.2167 kj/mol
  9.     dhufpsa=  6.3603 kj/mol         ddhufps=  0.7861 kj/mol
  10.        dhsa= 12.7028 kj/mol            ddhs=  0.9143 kj/mol
  11.       adswt= 11.0945 wt%
  12.      adsmgg=  110.94 mg/g
  13.    adsmmolg=  6.9168 mmol/g
  14.     adsgcm3=  0.1387 g/cm3
  15. adsmmolcm3=  8.6443 mmol/cm3
  16.    fracfill=  0.4065

复制代码


我想提取每个p条件下的数据，以及对应的每个result.txt文件里的第1列变量里adswt的数据，然后输出到一个名为adswt.dat的文件中，其文件格式如下：


  3.                  w,3.9039nm          w,4.2588nm
  4. p,kpa            adswt, wt%          adswt,wt%         
  5. 2.50e+01         xxxxxx              xxxxxx
  6. 4.00e+01         xxxxxx              xxxxxx
  7. 5.50e+01         xxxxxx              xxxxxx
  8. 7.00e+01         xxxxxx              xxxxxx

复制代码


其中的，xxxxxx为各个p和w条件下的result.txt文件里的adswt数据，其值不相同。

请问shell中该如何用awk或者sed命令写该段代码？？
回复

使用道具 举报

lastart
发表于 2007-1-19 13:06:35 | 显示全部楼层
试试, 我没试过的.

  2. #!/bin/bash
  3. wdir=(w_*nm)
  4. echo -n >adswt.dat
  5. for ((i=0;i<${#wdir[*]};i++))
  6. do
  7. echo -ne "\t"${wdir[i]/_/,}>>adswt.dat
  8. done
  9. echo >>adswt.dat
  10. echo -ne 'p,kpa'>>adswt.dat
  11. for ((i=0;i<${#wdir[*]};i++))
  12. do
  13. echo -ne '\tadswt,wt%'>>adswt.dat
  14. done
  15. echo >>adswt.dat
  16. cd ${wdir[1]}
  17. pdir=(p_*kpa)
  18. cd ..
  19. for ((i=0;i<${#pdir[*]};i++))
  20. do
  21.         a=$(echo ${pdir[i]}|sed 's/p_\|kpa//g')
  22.         echo -ne "$a">>adswt.txt
  23.         for ((j=0;j<${#wdir[*]};j++))
  24.         do
  25.                 echo -ne "\t">>adswt.txt
  26.                 awk '/adswt/ {print $2}' "${wdir[j]}/${pdir[i]}/result.txt" >>adswt.txt
  27.         done
  28.         echo >> adswt.txt
  29. done
复制代码
回复 支持 反对

使用道具 举报

发表于 2007-1-19 20:14:45 | 显示全部楼层
有点晕了。呵呵。
lz试试看吧：
  2. #!/bin/bash
  4. TITLE="title"
  5. CONTENT="content"
  6. ADSWT="adswt"
  8. wlist=""
  10. echo | awk '{printf("%10-s%10-s"," "," ")}'>$TITLE
  11. for i in `find . -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f2 | uniq`;
  12. do
  13.         w=`echo $i | sed -n -e 's/_/,/;p'`
  14.         wlist=${wlist}" "$w
  15.         echo $w | awk '{printf("%10-s%10-s",$1," ")}'>>$TITLE;
  16. done
  17. echo | awk '{printf("\n")}'>>$TITLE
  19. echo | awk '{printf("%10-s%10-s","p,kpa"," ")}'>>$TITLE
  20. for i in ${wlist};
  21. do
  22.         echo | awk '{printf("%10-s%10-s","adswt, wt%"," ")}'>>$TITLE
  23. done
  24. echo | awk '{printf("\n")}'>>$TITLE
  26. echo "" > $CONTENT
  27. for i in `find . -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f2 | uniq`;
  28. do
  29.         w=`echo $i | sed -n -e 's/_/,/;p'`
  30.         for j in `find ./$i -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f3 | uniq`;
  31.         do
  32.                 p=`echo $j | cut -d'_' -f2 | sed -n -e 's/kpa//;p'`
  33.                 a=`cat ./$i/$j/result.txt | grep "adswt" | awk '{print $2}'`
  35.                 line=`grep $p $CONTENT|cut -d' ' -f1`
  36.                 if [ "$line" = "$p" ]
  37.                 then
  38.                         awk '{if($1 != pinfo) print $0}' pinfo=$p ${CONTENT} >> ${CONTENT}.t
  39.                         awk '{if($1 == pinfo) printf("%s%10-s%10-s",$0," ",ainfo)}' pinfo=$p ainfo=$a ${CONTENT} >> ${CONTENT}.t
  40.                         mv -f ${CONTENT}.t ${CONTENT}
  41.                 else
  42.                         echo $p | awk '{printf("\n%10-s%10-s",$1," ")}' >> $CONTENT
  43.                         echo $a | awk '{printf("%10-s",$1)}' >> $CONTENT
  44.                 fi
  45.         done
  46. done
  48. echo | awk '{printf("\n")}' >> $CONTENT
  50. cat $TITLE $CONTENT > $ADSWT
  51. rm $CONTENT $TITLE
  53. cat $ADSWT
复制代码

俺机器上的结果（w_4.2588nm目录下没有p＝2.50e+01的子目录）
adswt的数值都是我自己瞎写的。
~/programming/ksh/example1$ mysh
                    w,3.9039nm          w,4.2588nm         
p,kpa               adswt, wt%          adswt, wt%         


2.50e+01            08.0945   
4.00e+01            09.0945             11.0945   
5.50e+01            10.0945             12.0945
回复 支持 反对

使用道具 举报

 楼主| 发表于 2007-1-19 21:16:53 | 显示全部楼层
能否把如下输出的结果 做成code输出

w,3.9039nm w,4.2588nm
p,kpa adswt, wt% adswt, wt%


2.50e+01 08.0945
4.00e+01 09.0945 11.0945
5.50e+01 10.0945 12.0945


那样我才能看出真正是否对齐了
回复 支持 反对

使用道具 举报

 楼主| 发表于 2007-1-19 21:27:44 | 显示全部楼层
Post by peter_huang
有点晕了。呵呵。
lz试试看吧：


  2. #!/bin/bash

  4. TITLE="title"
  5. CONTENT="content"
  6. ADSWT="adswt"

  8. wlist=""

  10. echo | awk '{printf("%10-s%10-s"," "," ")}'>$TITLE
  11. for i in `find . -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f2 | uniq`;
  12. do
  13.         w=`echo $i | sed -n -e 's/_/,/;p'`
  14.         wlist=${wlist}" "$w
  15.         echo $w | awk '{printf("%10-s%10-s",$1," ")}'>>$TITLE;
  16. done
  17. echo | awk '{printf("\n")}'>>$TITLE

  19. echo | awk '{printf("%10-s%10-s","p,kpa"," ")}'>>$TITLE
  20. for i in ${wlist};
  21. do
  22.         echo | awk '{printf("%10-s%10-s","adswt, wt%"," ")}'>>$TITLE
  23. done
  24. echo | awk '{printf("\n")}'>>$TITLE

  26. echo "" > $CONTENT
  27. for i in `find . -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f2 | uniq`;
  28. do
  29.         w=`echo $i | sed -n -e 's/_/,/;p'`
  30.         for j in `find ./$i -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f3 | uniq`;
  31.         do
  32.                 p=`echo $j | cut -d'_' -f2 | sed -n -e 's/kpa//;p'`
  33.                 a=`cat ./$i/$j/result.txt | grep "adswt" | awk '{print $2}'`

  35.                 line=`grep $p $CONTENT|cut -d' ' -f1`
  36.                 if [ "$line" = "$p" ]
  37.                 then
  38.                         awk '{if($1 != pinfo) print $0}' pinfo=$p ${CONTENT} >> ${CONTENT}.t
  39.                         awk '{if($1 == pinfo) printf("%s%10-s%10-s",$0," ",ainfo)}' pinfo=$p ainfo=$a ${CONTENT} >> ${CONTENT}.t
  40.                         mv -f ${CONTENT}.t ${CONTENT}
  41.                 else
  42.                         echo $p | awk '{printf("\n%10-s%10-s",$1," ")}' >> $CONTENT
  43.                         echo $a | awk '{printf("%10-s",$1)}' >> $CONTENT
  44.                 fi
  45.         done
  46. done

  48. echo | awk '{printf("\n")}' >> $CONTENT

  50. cat $TITLE $CONTENT > $ADSWT
  51. rm $CONTENT $TITLE

  53. cat $ADSWT
复制代码


俺机器上的结果（w_4.2588nm目录下没有p＝2.50e+01的子目录）
adswt的数值都是我自己瞎写的。
~/programming/ksh/example1$ mysh
                    w,3.9039nm          w,4.2588nm         
p,kpa               adswt, wt%          adswt, wt%         


2.50e+01            08.0945   
4.00e+01            09.0945             11.0945   
5.50e+01            10.0945             12.0945



  1. echo | awk '{printf("%10-s%10-s"," "," ")}'>$TITLE
复制代码

这段代码不太好吧，因为实际上我是有14个w目录和40个p目录，你这样做就只能输出2个w，最好将w设置成变量，然后做个循环打印TITLE。
回复 支持 反对

使用道具 举报

发表于 2007-1-19 22:15:18 | 显示全部楼层
呵呵，请仔细分析一下，echo | awk '{printf("%10-s%10-s"," "," ")}'>$TITLE
这个只是打一个最初的开头。
打印w是在这里：
  2. for i in `find . -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f2 | uniq`;
  3. do
  4.         w=`echo $i | sed -n -e 's/_/,/;p'`
  5.         wlist=${wlist}" "$w
  6.         echo $w | awk '{printf("%10-s%10-s",$1," ")}'>>$TITLE;
  7. done
  8. echo | awk '{printf("\n")}'>>$TITLE
复制代码


我机器上的输出：
  2.                     w,3.9039nm          w,4.2588nm         
  3. p,kpa               adswt, wt%          adswt, wt%         
  6. 2.50e+01            08.0945   
  7. 4.00e+01            09.0945             11.0945   
  8. 5.50e+01            10.0945             12.0945
复制代码
回复 支持 反对

使用道具 举报

 楼主| 发表于 2007-1-19 23:26:45 | 显示全部楼层
echo | awk '{printf("%10-s%10-s"," "," ")}'>$TITLE

这句是什么意思？执行出来是什么结果？

还有，%10-s还是 %-10s啊 ？？
回复 支持 反对

使用道具 举报

lastart
发表于 2007-1-19 23:35:27 | 显示全部楼层
果然有点小问题, 文件名前后用的不一样, 格式也不是很合要求
改了一下.

  2. #!/bin/bash
  3. exec >adswt.dat
  4. wdir=(w_*nm)
  5. printf '%20s' " "
  6. for ((i=0;i<${#wdir[*]};i++))
  7. do
  8. printf '%-20s' ${wdir[i]/_/,}
  9. done
  10. echo
  11. printf '%-20s' 'p,kpa'
  12. for ((i=0;i<${#wdir[*]};i++))
  13. do
  14. printf '%-20s' 'adswt,wt%'
  15. done
  16. echo
  17. pdir=($(ls w_*nm|sed '/^w_.*:$/d'|sort -u))
  18. for ((i=0;i<${#pdir[*]};i++))
  19. do
  20.     a=$(echo ${pdir[i]}|sed 's/p_\|kpa//g')
  21.     printf '%-20s' "$a"
  22.     for ((j=0;j<${#wdir[*]};j++))
  23.     do
  24.         if [ -f "${wdir[j]}/${pdir[i]}/result.txt" ]
  25.         then
  26.             awk '/adswt/ {printf("%-20s",$2)}' "${wdir[j]}/${pdir[i]}/result.txt"
  27.         else
  28.             printf "%-20s" " "
  29.         fi
  30.     done
  31.     echo
  32. done
复制代码

应该没问题, 考虑并不是所有的w目录里的p目录都一样多时测试效果:

  2.                     w,3.2nm             w,3.3nm
  3. p,kpa               adswt,wt%           adswt,wt%
  4. 2.3e+01             11.0945             13.0945
  5. 2.4e+01             12.0945             14.0945
  6. 2.5e+01                                 18.0945
复制代码
回复 支持 反对

使用道具 举报

发表于 2007-1-20 11:25:00 | 显示全部楼层
发现俺的也有问题，这题真不好做。
如果w目录下的p目录数量不一致，很容易出错。

改了一下：
  2. #!/bin/bash
  4. TITLE="title"
  5. CONTENT="content"
  6. ADSWT="adswt"
  7. NA="-------"
  9. wlist=""
  11. echo | awk '{printf("%10-s%10-s"," "," ")}'>$TITLE
  12. for i in `find . -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f2 | uniq`;
  13. do
  14.         w=`echo $i | sed -n -e 's/_/,/;p'`
  15.         wlist=${wlist}" "$w
  16.         echo $w | awk '{printf("%10-s%10-s",$1," ")}'>>$TITLE;
  17. done
  18. echo | awk '{printf("\n")}'>>$TITLE
  20. echo | awk '{printf("%10-s%10-s","p,kpa"," ")}'>>$TITLE
  21. for i in ${wlist};
  22. do
  23.         echo | awk '{printf("%10-s%10-s","adswt, wt%"," ")}'>>$TITLE
  24. done
  25. echo | awk '{printf("\n")}'>>$TITLE
  27. echo "" > $CONTENT
  28. counter=0
  29. for i in `find . -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f2 | uniq`;
  30. do
  31.         counter=`expr $counter + 1`
  32.         for j in `find ./$i -name "result.txt" -exec dirname "{}" \; | cut -d'/' -f3 | uniq`;
  33.         do
  34.                 p=`echo $j | cut -d'_' -f2 | sed -n -e 's/kpa//;p'`
  35.                 a=`cat ./$i/$j/result.txt | grep "adswt" | awk '{print $2}'`
  37.                 line=`grep $p $CONTENT|cut -d' ' -f1`
  38.                 if [ "$line" = "$p" ]
  39.                 then
  40.                         awk '
  41.                         {
  42.                                 if($1 == pinfo) {
  43.                                         y=0
  44.                                         printf("%s",$0)
  45.                                         if(NF == winfo)
  46.                                                 printf("%10-s"," ")
  47.                                         else
  48.                                                 for(x=NF;x<=winfo+1;x++) {
  49.                                                         if(y%2==0)
  50.                                                                 printf("%10-s"," ")
  51.                                                         else
  52.                                                                 printf("%10-s",na)
  53.                                                         y++
  54.                                                 }
  55.                                         printf("%10-s\n",ainfo);
  56.                                 }
  57.                                 else
  58.                                         printf("%s\n",$0)
  59.                         }' pinfo=$p ainfo=$a winfo=$counter na=$NA ${CONTENT} >> ${CONTENT}.t
  60.                         mv -f ${CONTENT}.t ${CONTENT}
  61.                 else
  62.                         echo | awk '{printf("\n%10-s",pinfo)}' pinfo=$p >> ${CONTENT}
  63.                         grep $p ${CONTENT} | awk '
  64.                         {
  65.                                 y=0
  66.                                 if(NF == winfo)
  67.                                         printf("%10-s"," ")
  68.                                 else
  69.                                         for(x=NF;x<=winfo+1;x++) {
  70.                                                 if(y%2==0)
  71.                                                         printf("%10-s"," ")
  72.                                                 else
  73.                                                         printf("%10-s",na)
  74.                                                 y++
  75.                                         }
  76.                                 printf("%10-s",ainfo)
  77.                         }' ainfo=$a winfo=$counter na=$NA >> ${CONTENT}
  78.                 fi
  79.         done
  80. done
  82. echo | awk '{printf("\n")}' >> $CONTENT
  84. cat $TITLE $CONTENT > $ADSWT
  85. rm $CONTENT $TITLE
  87. cat $ADSWT
复制代码

俺的目录结构如下：
  2. |-----w_3.9039nm
  3. |     |-----p_2.50e+01kpa
  4. |     |     |-----result.txt
  5. |     |-----p_5.50e+01kpa
  6. |     |     |-----result.txt
  7. |-----readme
  8. |-----mysh
  9. |-----w_4.2588nm
  10. |     |-----p_4.00e+01kpa
  11. |     |     |-----result.txt
  12. |     |-----p_5.50e+01kpa
  13. |     |     |-----result.txt
  14. |-----w_5.8079nm
  15. |     |-----p_2.50e+01kpa
  16. |     |     |-----result.txt
  17. |     |-----p_5.50e+01kpa
  18. |     |     |-----result.txt
  19. |-----w_6.7903nm
  20. |     |-----p_4.00e+01kpa
  21. |     |     |-----result.txt
  22. |-----adswt
复制代码

最后的输出如下：
  2.                     w,3.9039nm          w,4.2588nm          w,5.8079nm          w,6.7903nm         
  3. p,kpa               adswt, wt%          adswt, wt%          adswt, wt%          adswt, wt%         
  6. 2.50e+01            08.0945             -------             13.0945   
  7. 5.50e+01            10.0945             12.0945             14.0945   
  8. 4.00e+01            -------             11.0945             -------             15.0945
复制代码

如果lz不需要结果中那些个-------，那么一个简单的sed就可以了。
回复 支持 反对

使用道具 举报

 楼主| 发表于 2007-1-22 10:41:55 | 显示全部楼层
我的代码如下：


  3.     10      ROOT_DIRECTORY=/export/home/pengx/testsh
  4.     11      WORK_DIRECTORY=$ROOT_DIRECTORY/work
  5.     12      TASKS_DIRECTORY=$ROOT_DIRECTORY/tasks
  6.     13      RESULTS_DIRECTORY=$ROOT_DIRECTORY/results
  7.     14  
  8.     15      PORE_WIDTH_DIRECTORY=(0.3549 0.7098)
  9.     16  
  10.     17      PRESSURE_DIRECTORY=(1.00e-08 2.50e-08 5.00e-08)
  11.     18  
  12.     19      RESULT_FILES=(rho nk ucm ucf ucp dh dhuff dhufp adswt adsmgg adsmmolg adsgcm3 adsmmolcm3 fracfill)


  15.     99      echo -ne "Processing final data ........................\t"
  16.    100      cd $RESULTS_DIRECTORY
  17.    101      x=0
  18.    102      for i in ${RESULT_FILES[@]}
  19.    103      do
  20.    104          if [ -f ${i}.txt ]; then
  21.    105              mv ${i}.txt _${i}.old
  22.    106          fi
  23.    107          touch ${i}.txt
  24.    108          printf "%-20s" " " >> ${i}.txt
  25.    109          for j in ${PORE_WIDTH_DIRECTORY[@]}
  26.    110          do
  27.    111              printf "%-20s" "w,${j}nm" >> ${i}.txt
  28.    112          done
  29.    113          printf "\n" >> ${i}.txt
  30.    114          printf "%-20s" "p,kpa" >> ${i}.txt
  31.    115          x=`expr $x+1`
  32.    116          for j in ${PORE_WIDTH_DIRECTORY[@]}
  33.    117          do
  34.    118              for k in ${PRESSURE_DIRECTORY[@]}
  35.    119              do
  36.    120                  while [ $x -lt 6 ]
  37.    121                  do
  38.    122                      printf "%-20s" "${i}" >> ${i}.txt
  39.    123                  done
  40.    124                  while [ $x -ge 6 -a $x -lt 14 ]
  41.    125                  do
  42.    126                      a=`cat $WORK_DIRECTORY/w_${j}nm/p_${k}kpa/result.dat | grep "${i}" | awk '{print $3}'`
  43.    127                      printf "%-20s" "${i}," "$a" >> ${i}.txt
  44.    128                  done
  45.    129                  while [ $x -eq 14 ]
  46.    130                  do
  47.    131                      printf  "%-20s" "${i}" >> ${i}.txt
  48.    132                  done
  49.    133              done
  50.    134          done
  51.    135          printf "\n" >> ${i}.txt
  52.    136          for j in ${PORE_WIDTH_DIRECTORY[@]}
  53.    137          do
  54.    138              for k in ${PRESSURE_DIRECTORY[@]}
  55.    139              do
  56.    140                  printf "%-20s\n" "${k}" >> ${i}.txt
  57.    141              done
  58.    142              for k in ${PRESSURE_DIRECTORY[@]}
  59.    143              do
  60.    144                  a=`cat $WORK_DIRECTORY/w_${j}nm/p_${k}kpa/result.dat | grep "${i}" | awk '{print $2}'`
  61.    145                  printf "%-20s" "$a" >> ${i}.txt
  62.    146              done
  63.    147          done
  64.    148          printf "\n" >> ${i}.txt
  65.    149      done
  66.    150      echo "Done."
复制代码


结果执行时，出现如下的错误：
./test.sh: line 120: [: 0+1: integer expression expected
./test.sh: line 124: [: 0+1: integer expression expected
./test.sh: line 129: [: 0+1: integer expression expected

请问是什么意思？
回复 支持 反对

使用道具 举报

下一页 »
123下一页
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表