Bash循环中变量传递给grep文件名失败，求解决方案

我有一批带统一前缀的SDF文件，需要提取含坐标的特定行，后续用awk计算均方根偏差（RMSD）。单文件下正则表达式验证有效，但在for循环中使用bash变量时，grep无法正确识别变量，尝试了多种引号、转义、大括号组合都未解决。

单文件正常执行示例 #

单文件下执行以下命令可正确提取坐标行：

grep " [A-Z] " ../glide-dock_SP_8CHM/ligands/3579831839.sdf | grep " 'F\|S\|C\|O\|N "

输出结果：

2.7118    4.0281   21.0125 S   0  0  0  0  0  0
    4.0921    3.8708   21.4967 O   0  0  0  0  0  0
    1.8208    2.8602   20.9648 O   0  0  0  0  0  0
    1.9954    5.2598   21.8979 N   0  0  1  0  0  0
    2.8079    4.6789   19.3978 N   0  0  0  0  0  0
    2.3264    6.6613   21.7269 C   0  0  0  0  0  0
    0.5679    5.5328   21.8880 C   0  0  0  0  0  0
    3.8805    4.1396   18.5518 C   0  0  0  0  0  0
    1.5591    4.8798   18.6409 C   0  0  0  0  0  0
    0.9197    6.9883   22.3042 C   0  0  0  0  0  0
    0.8840    7.1094   23.8539 C   0  0  0  0  0  0
    0.1007    8.1332   21.6760 C   0  0  0  0  0  0
   -0.1257    6.7465   24.4613 O   0  0  0  0  0  0
    1.9559    7.6325   24.4759 N   0  0  0  0  0  0
    0.9264    9.3760   21.3740 C   0  0  0  0  0  0
    2.1242    7.8085   25.9208 C   0  0  1  0  0  0
    1.5512    9.5358   20.1194 C   0  0  0  0  0  0
    1.0887   10.3742   22.3570 C   0  0  0  0  0  0
    3.4722    7.2469   26.3919 C   0  0  0  0  0  0
    1.7964    9.2425   26.3963 C   0  0  0  0  0  0
    2.3648   10.6567   19.8673 C   0  0  0  0  0  0
    1.8942   11.5002   22.1018 C   0  0  0  0  0  0
    3.4888    6.2142   27.0561 O   0  0  0  0  0  0
    4.6188    7.8580   26.0568 N   0  0  0  0  0  0
    2.6369   10.3760   25.7768 C   0  0  0  0  0  0
    2.5397   11.6366   20.8599 C   0  0  0  0  0  0
    4.8138    9.1019   25.3276 C   0  0  0  0  0  0
    4.1501   10.2988   26.0247 C   0  0  0  0  0  0

循环执行时的错误情况 #

初始尝试命令 #

尝试遍历ligand-list.txt中的前缀，命令如下：

for i in $(cat ligand-list.txt)
do
grep " [A-Z] " ../glide-dock_SP_8CHM/ligands/$i.sdf | grep " 'F\|S\|C\|O\|N " | awk '{print $1,$2,$3}' > $i_8CHM.xyz
done

错误输出：

grep: ../glide-dock_SP_8CHM/ligands/ligand-list.txt.sdf: No such file or directory

转义双引号尝试 #

给变量加转义双引号后，变量被当作字面量：

for i in $(cat ligand-list.txt)
do
grep " [A-Z] " ../glide-dock_SP_8CHM/ligands/"\$i".sdf | grep " 'F\|S\|C\|O\|N " | awk '{print $1,$2,$3}' > $i_8CHM.xyz
done

错误输出：

grep: ../glide-dock_SP_8CHM/ligands/$i.sdf: No such file or directory

大括号转义尝试 #

使用转义大括号后，符号被当作文件名的一部分：

for i in $(cat ligand-list.txt)
do
grep " [A-Z] " ../glide-dock_SP_8CHM/ligands/"\{$i}".sdf | grep " 'F\|S\|C\|O\|N " | awk '{print $1,$2,$3}' > $i_8CHM.xyz
done

错误输出：

grep: ../glide-dock_SP_8CHM/ligands/\{3579831839}.sdf: No such file or directory

问题澄清与解决方案 #

核心问题分析 #

1. for i in $(cat ...)的隐患：这种写法会把文件内容按空格拆分，若前缀包含空格或特殊字符，会导致变量解析错误；极端情况下甚至会把文件名本身当作变量值，出现找ligand-list.txt.sdf的错误。
2. 变量边界不明确：$i_8CHM.xyz会被bash解析为变量i_8CHM，而非i加上_8CHM，导致输出文件名错误。
3. 引号使用错误：双引号内的变量会被bash解析，你之前转义$i反而让它变成了字面量；文件名路径应直接用双引号包裹，避免空格拆分参数。

正确循环写法 #

使用while read安全读取文件，并明确变量边界：

while IFS= read -r i; do
  grep " [A-Z] " "../glide-dock_SP_8CHM/ligands/${i}.sdf" | grep " 'F\|S\|C\|O\|N " | awk '{print $1,$2,$3}' > "${i}_8CHM.xyz"
done < ligand-list.txt

优化版本（合并grep） #

将两次grep合并为一次，减少进程调用：

while IFS= read -r i; do
  grep -E " [A-Z] .* '[FSCON] " "../glide-dock_SP_8CHM/ligands/${i}.sdf" | awk '{print $1,$2,$3}' > "${i}_8CHM.xyz"
done < ligand-list.txt

关键说明 #

• while IFS= read -r i：安全读取每行内容，IFS=保留行首行尾空格，-r禁止反斜杠转义，避免解析错误。
• ${i}.sdf：用大括号明确变量i的边界，防止bash把后续字符当作变量名的一部分。
• 双引号包裹路径/文件名：确保带空格的路径或文件名被当作单个参数传递，避免被拆分。

内容的提问来源于stack exchange，提问作者blakem