Perl脚本计算蛋白配体原子距离后无法移动文件问题排查

看起来你已经把距离计算的逻辑搞定了，但文件移动的问题出在几个关键的小细节上，我帮你拆解并修复：

核心错误点 #

1. last语句直接阻断了移动命令
   你在判断距离≤5时，先写了last;，这会直接跳出当前的foreach $atom循环，后面的移动代码根本没机会执行！必须调整顺序，或者用标记来控制循环跳出。

2. 移动命令用错了参数
   你写的是system ("mv $fh mp");，但$fh是文件句柄（代表打开的文件流），不是文件名！正确的应该用遍历的文件名变量$txt，也就是system("mv '$txt' mp/");。给文件名加单引号还能避免文件名含空格、特殊字符时出错。

3. 缺少错误检查逻辑

   • mkdir "mp"没检查目录是否已存在，重复创建会抛出警告，改成mkdir "mp" unless -d "mp";更稳妥
   • 打开文件时没加错误处理，万一文件打不开你根本不知道，应该加上or die "Can't open $txt: $!";
   • system命令的执行结果没检查，就算移动失败也不会有提示，建议判断返回值

4. 循环跳出逻辑不完善
   找到符合条件的距离后，你需要停止处理当前文件的所有剩余行，否则脚本会继续读取后续内容，甚至可能重复触发移动操作。

修复后的完整脚本 #

#!/usr/local/bin/perl
use strict;
use warnings;
use File::Glob;

# 仅当mp目录不存在时创建
mkdir "mp" unless -d "mp";

for my $txt ( glob '*.txt' ) {
    # 打开文件并添加错误捕获
    open my $fh, '<', $txt or die "无法打开文件 $txt: $!";
    my (@refer, @points);
    my $part = 0;
    my $found_close_pair = 0; # 标记是否找到符合条件的原子对

    while (my $line = <$fh>) {
        chomp($line);
        # 遇到HETATM时切换到配体原子部分
        $part++ if ($line =~ /^HETATM/);
        
        # 提取坐标并转换为数值（避免字符串运算隐患）
        my @array = map { $_ + 0 } (substr($line, 30, 8), substr($line,38,8), substr($line,46,8));

        if ($part == 0) {
            push @refer, [ @array ]; # 存储蛋白原子坐标
        } else {
            push @points, [ @array ]; # 存储配体原子坐标
        }

        # 没有配体原子时跳过计算
        next unless @points;

        foreach my $ref_atom (@refer) {
            my ($x1, $y1, $z1) = @{$ref_atom};
            foreach my $lig_atom (@points) {
                my ($x, $y, $z) = @{$lig_atom};
                my $dist = sqrt( ($x-$x1)**2 + ($y-$y1)**2 + ($z-$z1)**2 );
                
                if ($dist <= 5) {
                    print "文件 $txt 中检测到近距离原子对，距离: $dist\n";
                    $found_close_pair = 1;
                    # 跳出两层foreach循环
                    last foreach 1..2;
                }
            }
            last if $found_close_pair;
        }
        last if $found_close_pair; # 跳出文件行读取循环
    }

    # 找到符合条件的原子对后移动文件
    if ($found_close_pair) {
        my $move_cmd = qq{mv '$txt' mp/};
        if (system($move_cmd) == 0) {
            print "成功将 $txt 移动到 mp/ 目录\n";
        } else {
            warn "移动 $txt 失败，退出码: $?";
        }
    }
    close $fh;
}

额外优化说明 #

• 把坐标字符串显式转换为数值，避免Perl自动转换时可能出现的意外问题
• 用$found_close_pair标记来统一控制循环跳出，逻辑更清晰
• 移动命令用qq{}包裹，方便处理文件名中的特殊字符
• 增加了移动操作的成功/失败反馈，方便调试

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