Bash动态变量赋值问题：无法创建atom_1等带索引变量

我明白你遇到的问题了——在bash里直接动态生成带索引的变量名（比如atom_1、numberofatoms_1）的写法不对，这才导致了bash: atom_1=x: command not found的报错。咱们来一步步解决这个问题，顺便优化一下你的脚本。

错误原因解析 #

你尝试用numberofatoms_$i=$(( $numberofatoms + $noa ))这种方式赋值时，bash会先解析$i得到具体数字（比如1），然后把整个表达式变成numberofatoms_1=xxx。但bash的变量赋值语法要求左边必须是一个固定的合法变量名，它会把numberofatoms_1当成要执行的命令，而不是变量，所以才会报“command not found”。

最佳解决方案：使用数组 #

bash的数组天生就是为这种带索引的变量场景设计的，我们可以定义两个数组：

• atom_names：存储每个原子的类型（对应你的atom_1、atom_2）
• atom_counts：存储对应原子的数量（对应你的numberofatoms_1、numberofatoms_2）

这样既符合bash的语法规范，又能轻松实现你想要的索引式变量效果。

修改后的完整脚本 #

cd ~/A/
# 定义两个数组，用来存储原子类型和对应数量
declare -a atom_names
declare -a atom_counts

while read line
do
    cd "folder_$line" || continue  # 加上||continue，避免目录不存在时脚本崩溃
    
    # 直接读取input的第6、7行，不用临时文件
    temp1=$(sed -n 6p input)
    temp2=$(sed -n 7p input)
    
    # 统计原子类型的数量，不用中间文件
    atoms=$(wc -w <<< "$temp1")
    echo "原子类型总数: $atoms"
    
    # 重置数组索引（每次处理一个folder时清空）
    unset atom_names atom_counts
    declare -a atom_names
    declare -a atom_counts
    
    for ((i=1; i<=atoms; i++))
    do
        # 获取第i个原子类型
        atom=$(cut -d' ' -f$i <<< "$temp1")
        atom_names[$i]="$atom"
        
        # 获取第i个原子的数量并累加
        noa=$(cut -d' ' -f$i <<< "$temp2")
        atom_counts[$i]=$noa  # 若需累加可改为：atom_counts[$i]=$(( ${atom_counts[$i]:-0} + $noa ))
    done
    
    # 输出你想要的格式
    echo "=== 文件夹folder_$line的原子统计结果 ==="
    for ((i=1; i<=atoms; i++))
    do
        echo "atom_$i=${atom_names[$i]}、numberofatoms_$i=${atom_counts[$i]}"
    done
    
    cd ..
done < ~/A/B.txt  # 用重定向代替cat|while，避免子shell导致的变量丢失问题

关键改动说明 #

1. 替换临时文件：用$(sed -n 6p input)直接捕获命令输出，不用生成temp1.txt、junk等临时文件，脚本更简洁高效。
2. 使用数组存储：用atom_names[$i]和atom_counts[$i]代替动态命名的变量，完全符合bash语法规范。
3. 改用for循环：把原有的while循环改成for ((i=1; i<=atoms; i++))，写法更简洁，避免expr和bc的不必要使用。
4. 避免子shell问题：把cat ~/A/B.txt | while read line改成while read line done < ~/A/B.txt，这样循环内的变量修改不会因为子shell环境丢失。
5. 增加错误处理：cd "folder_$line" || continue确保如果目录不存在时，脚本不会崩溃，而是跳过继续处理下一行。

额外优化建议 #

• 如果input文件的第6、7行格式固定，还可以用awk替代sed+cut，比如直接用awk -v i=$i 'NR==6{print $i}' input获取第6行第i列，执行效率更高。
• 如果需要跨文件夹累加同一原子类型的总数量，可以把数组定义放在循环外，每次处理时直接累加对应索引的数值。

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