求帮忙看为什么这个if循环工作不正常

alphafighter

<code>
while(defined(my $m_z_value = <$ff>)) {
        $m_z_value =~ s/^\s+//;
    chomp $m_z_value;
        if ($tmp_to <= $to and $m_z_value <= $to){
                if ($m_z_value < $from){
                goto MIRACLE;
                }
                else{
                        if($m_z_value = $from){
                                $mass = $mass + $m_z_value;
                                $whole_mass = $whole_mass +$m_z_value;
                                ++$bin_pepnumber;
                                ++$whole_pepnumber;
                       
                                print "$bin_number\t\t $tmp_from\t $tmp_to\t $bin_pepnumber\t\t $mass/$bin_pepnumber\n";
                        }
                        else{
                                if ($m_z_value >= $tmp_from and $m_z_value <= $tmp_to) {       
                                $mass = $mass + $m_z_value;
                                $whole_mass = $whole_mass +$m_z_value;
                                ++$bin_pepnumber;
                                ++$whole_pepnumber;       
                                }
                                else{
                                        if ($m_z_value>$tmp_to){
                                        print "$bin_number\t\t $tmp_from\t $tmp_to\t $bin_pepnumber\t\t $mass/$bin_pepnumber\n";
                                        ++$bin_number;
                                        $mass = $m_z_value;
                                        $whole_mass = $whole_mass + $m_z_value;
                                        $bin_pepnumber = 1;
                                        ++$whole_pepnumber;
                                        $tmp_from = $tmp_from + $intervals;
                                        $tmp_to = $tmp_to + $intervals;
                                        }
                                        }
                                }
                        }
                }               
MIRACLE:
        }
<\code>

在这个循环里，我想表达的其实是，每个if基本都是并列的，阅读一个文件的每一行，判断这一行分别在大于某区间，在区间内和大于区间的时候该做什么。
但是这个程序运行到第二个if，假如这一行的值正好等于$from的时候，在debug中会，$m_z_value这个变量会直接变成$from的值，不论怎么改它都会在运行到这里的时候变成这个值。

求救！这个问题已经至少放倒了俩人了。

815138698

($m_z_value = $from）？？

alphafighter

回复 2# 815138698

没看明白你的意思。
帖子里说了，这一句想表达的是如果等于的话执行下面区块的文件。
我在后面的修改中将“=”改成了“==”就好使了，但是这个查不到原因。而且表述感觉也并不合法。

所以想知道问题的根源所在。

   

jason680

回复 3# alphafighter

不是语法(syntax)错误
而是语义错误,原因是误用==跟=语法

解决方法打开warnings跟strict
use strict;
use warnings;
   

alphafighter

1. #! /usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. # get the mz value and make it easier to use{
   
6. my $dummyfile = "dummy.pepmasses"; #filename defined here
   
7. my @mzco;
   
8. 
   
9. 
   
10. open my $fh, "<" , $dummyfile or die "Unable to open $dummyfile:$!";
    
11. readline $fh;
    
12. while (my $line = <$fh>){
    
13.         #read each line in file
    
14.         chomp $line;
    
15.         my $mz_value = (split/\s+/,$line)[2]; #pick column 2th at everyline
    
16.        
    
17.         push (@mzco,$mz_value);  #add them all in one array @mzco
    
18. }
    
19. open (FILE, ">tmp.data");
    
20. print FILE join("\n",@mzco),"\n";
    
21. close $fh;
    
22. close FILE;
    
23. 
    
24. my $tmpdata = "tmp.data";
    
25. open my $fa, "<", $tmpdata or die "Unable to open $tmpdata:$!";
    
26. my @mvalue = <$fa>;
    
27. my @somvalue = sort { $a <=> $b } @mvalue;
    
28. 
    
29. open (FI,">tmpf.data");
    
30. print FI"@somvalue";
    
31. close FI;
    
32. 
    
33. # }
    
34. 
    
35. my $tmpf = "tmpf.data";
    
36. open my $ff, "<" , $tmpf or die "Unable to open $tmpf:$!";
    
37. 
    
38. my $text = ();
    
39. 
    
40. # get rage and intervals
    
41. print "what range of information do you want to know?\n";
    
42. # sample text : 1000-4000-250-result
    
43. 
    
44. chomp($text = <STDIN>);
    
45. 
    
46. my ($from,$to,$intervals) = split (/\s+/,$text);
    
47. 
    
48. my $whole_pepnumber = 0;
    
49. my $whole_mass = 0;
    
50. my $mass = 0;
    
51. my $bin_pepnumber = 0;
    
52. my $bin_number = 0;
    
53. my $tmp_from = $from;
    
54. my $tmp_to = $tmp_from + $intervals;
    
55. 
    
56. # head line of final result
    
57. print "Bin number\t From\t To\t peptide number\t average m\/z value\n";
    
58. # get the number of bins
    
59. readline $ff;
    
60. while(defined(my $m_z_value = <$ff>)) {
    
61.         $m_z_value =~ s/^\s+//;
    
62.         GOD:
    
63.     chomp $m_z_value;
    
64.         if ($tmp_to <= $to and $m_z_value <= $to){
    
65.                        
    
66.                 if($m_z_value >> $tmp_to){
    
67.                         ++$bin_number;
    
68.                         print "$bin_number\t\t $tmp_from\t $tmp_to\t 0\t\t -\n";
    
69.                         $tmp_from = $tmp_from + $intervals;
    
70.                         $tmp_to = $tmp_to + $intervals;
    
71.                         goto GOD;
    
72.                 }else{
    
73.                 if ($m_z_value < $from){
    
74.                 goto MIRACLE;
    
75.                 }
    
76.                 else{
    
77.                         if($m_z_value == $from){
    
78.                                 $mass = $mass + $m_z_value;
    
79.                                 $whole_mass = $whole_mass +$m_z_value;
    
80.                                 ++$bin_pepnumber;
    
81.                                 ++$whole_pepnumber;
    
82.                        
    
83.                                 print "$bin_number\t\t $tmp_from\t $tmp_to\t $bin_pepnumber\t\t $mass/$bin_pepnumber\n";
    
84.                         }
    
85.                         else{
    
86.                                 if ($m_z_value >= $tmp_from and $m_z_value <= $tmp_to) {       
    
87.                                 $mass = $mass + $m_z_value;
    
88.                                 $whole_mass = $whole_mass +$m_z_value;
    
89.                                 ++$bin_pepnumber;
    
90.                                 ++$whole_pepnumber;       
    
91.                                 }
    
92.                                 else{
    
93.                                         if ($m_z_value > $tmp_to){
    
94.                                         print "$bin_number\t\t $tmp_from\t $tmp_to\t $bin_pepnumber\t\t $mass/$bin_pepnumber\n";
    
95.                                         ++$bin_number;
    
96.                                         $mass = $m_z_value;
    
97.                                         $whole_mass = $whole_mass + $m_z_value;
    
98.                                         $bin_pepnumber = 0;
    
99.                                         ++$whole_pepnumber;
    
100.                                         $tmp_from = $tmp_from + $intervals;
     
101.                                         $tmp_to = $tmp_to + $intervals;
     
102.                                         goto GOD;
     
103.                                         }
     
104.                                         }
     
105.                                 }
     
106.                         }
     
107.                 }               
     
108. MIRACLE:
     
109.        
     
110.         }
     
111. }
     
112. print "The peptide number in the range is $whole_pepnumber\n";
     
113. my $average_mass = $whole_mass / $whole_pepnumber;
     
114. print "The average m\/z value is $whole_mass/$whole_pepnumber";

复制代码

上面的是完整的代码。是有use strict和use warning的。但是源代码在运行中并没有反应任何有关后面的if跟else的部分。
因为不能直接插入图片，所以我把异常的情况更详细地说一下。
100到300范围，区间大小为20的时候。
其他数据都正常，只有240到260是0，而原处理数据中这个部分很明显不是0。经检查，发现，在其他行运行的时候，
if($m_z_value >> $tmp_to){ （这一行是后加的判定）
这一行的>>都是作大于号的，只有运行到240到260的时候，这一行会突然失效，我用了debug，当运行到这个区间的时候，$mz_v_value最小的值为245，最大的为259，$tmp_to为260，明明这个if在这里是false，但是程序就是会执行下面的代码。然后到260到280就会又恢复正常。

然后每次的范围总会漏掉最后一个区间。

我试验了更大的区间，1000到5000，区间大小为50，结果发现，之前那个莫名其妙失效的效果的出现完全没有规律，3,4,12,6,5,3,4,3，4. 前面的是它会出现的大概间隔。
@jason680

程序猿zy

你以为 >> 是干嘛的？？  ！！！！！！！！！！！！！！！！！！！！

zhlong8

代码排版很漂亮，但是不是你写的唉，感觉你什么也不懂

>> 是位移操作，相当于除以 2 的 N 次方。if () 的括号里面是赋值表达式，相当于判断赋的那个值是不是 false (包括 0, "", "0" 这个语法是抄C语言的，C语言中也是常见错误，平时 use strict; use warnings; 会提醒在 if 中用了赋值表达式，你需要的可能是 ==