while循环执行awk批量提取生成空输出文件的问题

你的脚本生成空文件的核心原因是：awk代码被单引号包裹时，shell变量$specie_name不会被解析，awk实际处理的是字面量字符串"$specie_name"，自然匹配不到任何包含目标物种名的行。

下面提供几种可行的修复方案：

方案1：用awk的-v选项传递变量（推荐） #

这是最安全、不易出错的方式，通过-v参数将shell变量传递给awk内部使用：

while IFS= read -r specie_name
do 
    awk -F'\t' -v target="$specie_name" '$3 ~ target {print}' "$krakenfile" > "${specie_name}_lines.txt"
done < "$fungalnames"

补充说明：

• 给$krakenfile和$fungalnames加上双引号，避免文件名包含空格、特殊字符时出现解析错误
• -v target="$specie_name"将shell的specie_name变量赋值给awk的target变量，awk内部直接引用target即可完成匹配

方案2：转义双引号中的awk变量 #

如果坚持用双引号包裹awk代码，需要转义awk里的$符号（防止被shell解析为shell变量），同时保留物种名的变量替换：

while IFS= read -r specie_name
do 
    awk -F'\t' "\$3 ~ \"$specie_name\" {print}" "$krakenfile" > "${specie_name}_lines.txt"
done < "$fungalnames"

注意：这种方式如果物种名包含*、&等shell特殊字符，可能会出现意外问题，因此仅适合物种名无特殊字符的场景。

优化方案：单次遍历大文件（适合大型$krakenfile） #

如果你的$krakenfile体积很大，多次遍历会浪费IO资源，可以用awk一次性完成所有物种的匹配与输出，只需要读一次源文件：

awk -F'\t' '
    BEGIN {
        # 先读取所有目标物种名存入数组
        while ((getline specie < ARGV[2]) > 0) {
            species[specie] = 1
        }
        close(ARGV[2])
        ARGV[2] = ""  # 标记第二个参数已处理，避免awk重复读取
    }
    {
        # 检查当前行是否匹配任一目标物种
        for (s in species) {
            if ($3 ~ s) {
                print > (s "_lines.txt")
                break  # 匹配到一个物种后停止检查，避免重复写入
            }
        }
    }
' "$krakenfile" "$fungalnames"

内容的提问来源于stack exchange，提问作者Mohamed Aziz Hadj Nasr