电子科技大学Perl语言编程各章课后作业提示

附录各章课后作业提示第一章2. say输出时默认换行（say Just like "print", but implicitly appends a newline.）。

如果看不出来，请分别选中两个脚本中的say与print一词，然后Ctrl+Shift+H，并查看文档。

3. Linux上可用：perl -e 'print "Hello, Strawberry Perl!\n";'Windows上可用：perl -e "print \"Hello, Strawberry Perl!\n\";"4. 假设日志文件名为error.log，且目前已在该文件所在目录，那么可以使用以下命令查看最近的100行日志信息：tail -100 error.log5. 假设在Linux系统Shell中已转到Chapter01文件夹里了。

输入：diff -u muscle skin | sort > venn，打开生成的结果文件venn，首先是在肌肉与皮肤都有的多肽（443条），中间+开头的多肽是只见于皮肤的（670条），最后-开头的多肽只见于肌肉（704条）。

6. 假设在Linux系统Shell中已转到Chapter01文件夹里了。

（1）先解压laf.tgz和update.tgz。

tar –xf laf.tgztar –xf update.tgz（2）生成laf和update文件夹里的文献列表ls ./laf > lafiles ls ./update updatefiles（3）比较lafiles与updatefiles的差异diff -u updatefiles lafiles | sort > diffiles打开diffiles文件，首先罗列的是laf.tgz中已被处理并已更新的文献（68篇），接着以+开头的文件就是laf.tgz中尚未处理的文献（214个）。

第二章2. 假设在Linux系统Shell中已转到Chapter02文件夹里了。

（1）解压缩all_mimodb_peptide.gz文件gzip -d all_mimodb_peptide.gz （17427条）（2）找出所有出现且仅出现过一次的多肽sort < all_mimodb_peptides | uniq -u > appear_only_once （15687条）（3）找出所有出现过多次的多肽sort < all_mimodb_peptides | uniq -d > appear_more_times （672条）（4）出现多次的仅收录一条，建立非冗余的多肽数据集sort < all_mimodb_peptides | uniq > nonredundant （16359条）3. 办法之一#!D:\Perl\perl\bin\perl.exemy $comment= '>gi|224589802:c2182439-2181009 Homo sapiens chromosome 11, GRCh37.p10 Primary Assembly';my $insulin= " AGCCCTCCAGGACAGGCTGCATCAGAAGAGGCCATCAAGCAGGTCTGTTCCAAGGGCCTTTGCGTCAGGT GGGCTCAGGATTCCAGGGTGGCTGGACCCCAGGCCCCAGCTCTGCAGCAGGGAGGACGTGGCTGGGCTCG TGAAGCATGTGGGGGTGAGCCCAGGGGCCCCAAGGCAGGGCACCTGGCCTTCAGCCTGCCTCAGCCCTGC CTGTCTCCCAGATCACTGTCCTTCTGCCATGGCCCTGTGGATGCGCCTCCTGCCCCTGCTGGCGCTGCTG GCCCTCTGGGGACCTGACCCAGCCGCAGCCTTTGTGAACCAACACCTGTGCGGCTCACACCTGGTGGAAG CTCTCTACCTAGTGTGCGGGGAACGAGGCTTCTTCTACACACCCAAGACCCGCCGGGAGGCAGAGGACCT GCAGGGTGAGCCAACTGCCCATTGCTGCCCCTGGCCGCCCCCAGCCACCCCCTGCTCCTGGCGCTCCCAC CCAGCATGGGCAGAAGGGGGCAGGAGGCTGCCACCCAGCAGGGGGTCAGGTGCACTTTTTTAAAAAGAAG TTCTCTTGGTCACGTCCTAAAAGTGACCAGCTCCCTGTGGCCCAGTCAGAATCTCAGCCTGAGGACGGTG TTGGCTTCGGCAGCCCCGAGATACATCAGAGGGTGGGCACGCTCCTCCCTCCACTCGCCCCTCAAACAAA TGCCCCGCAGCCCATTTCTCCACCCTCATTTGATGACCGCAGATTCAAGTGTTTTGTTAAGTAAAGTCCT GGGTGACCTGGGGTCACAGGGTGCCCCACGCTGCCTGCCTCTGGGCGAACACCCCATCACGCCCGGAGGA GGGCGTGGCTGCCTGCCTGAGTGGGCCAGACCCCTGTCGCCAGGCCTCACGGCAGCTCCATAGTCAGGAG ATGGGGAAGATGCTGGGGACAGGCCCTGGGGAGAAGTACTGGGATCACCTGTTCAGGCTCCCACTGTGAC GCTGCCCCGGGGCGGGGGAAGGAGGTGGGACATGTGGGCGTTGGGGCCTGTAGGTCCACACCCAGTGTGG GTGACCCTCCCTCTAACCTGGGTCCAGCCCGGCTGGAGATGGGTGGGAGTGCGACCTAGGGCTGGCGGGC AGGCGGGCACTGTGTCTCCCTGACTGTGTCCTCCTGTGTCCCTCTGCCTCGCCGCTGTTCCGGAACCTGC TCTGCGCGGCACGTCCTGGCAGTGGGGCAGGTGGAGCTGGGCGGGGGCCCTGGTGCAGGCAGCCTGCAGC CCTTGGCCCTGGAGGGGTCCCTGCAGAAGCGTGGCATTGTGGAACAATGCTGTACCAGCATCTGCTCCCT CTACCAGCTGGAGAACTACTGCAACTAGACGCAGCCCGCAGGCAGCCCCACACCCGCCGCCTCCTGCACC GAGAGAGATGGAATAAAGCCCTTGAACCAGC";my @reverseq= reverse (split //, $insulin);foreach my $base(@reverseq){if ($base eq 'T'){$base= 'A'; next; };if ($base eq 'A'){$base= 'T'; next; };if ($base eq 'C'){$base= 'G'; next; }if ($base eq 'G'){$base= 'C'; next; }}print "$comment\n";print @reverseq;4. 办法之一#!D:\Perl\perl\bin\perl.exe#hash for amino acid code convertingmy %AA123= (A => 'Ala', V => 'Val', L => 'Leu', I => 'Ile',P => 'Pro', W => 'Trp', F => 'Phe', M => 'Met',G => 'Gly', S => 'Ser', T => 'Thr', Y => 'Tyr',C => 'Cys', N => 'Asn', Q => 'Gln', K => 'Lys',R => 'Arg', H => 'His', D => 'Asp', E => 'Glu',);my %AA321= reverse %AA123;my $comment_H= ">metuximab Heavy chain";my $comment_L= ">metuximab Light chain";my $heavy= <<'HEAVY';Met Asn Phe Gly Leu Ser Trp Val Phe Ile Val Ile Leu Leu Lys Gly Val Gln Ser Glu Val Lys Leu Glu Glu Ser Gly Gly Gly Leu Val Gln Pro Gly Gly Ser Met Lys Leu Ser Cys Val Ala Ser Gly Phe Thr Phe Ser Asp Ala Trp Met Asp Trp Val Arg Gln Ser Pro Glu Lys Gly Leu Glu Trp Val Ala Glu Ile Arg Ser Ala Asn Asn His Ala Pro Tyr Tyr Tyr Thr Glu Ser Val Lys Gly Arg Phe Thr Ile Ser Arg Asp Asp Ser Lys Ser Ile Ile Tyr Leu Gln Met Asn Asn Leu Arg Ala Glu Asp Thr Gly Ile Tyr Tyr Cys Thr Arg Asp Ser Thr Ala Thr His Trp Gly Gln Gly Thr Leu Val Thr Val Ser Ser Ala Ser Thr Thr Ala Pro Ser ValHEAVYmy $light= <<'LIGHT';Met Asp Ser His Thr Gln Val Phe Ile Phe Leu Leu Leu Cys Val Ser Gly Ala His Gly Ser Ile Val Met Thr Gln Thr Pro Thr Phe Leu Val Val Ser Ala Gly Asp Arg Val Thr Ile Thr Cys Lys Ala Ser Gln Ser Val Ile Asn Asp Val Ala Trp Tyr Gln Gln Lys Pro Gly Gln Ser Pro Lys Leu Leu Ile Phe Tyr Ala Ser Asn Arg Asn Thr Gly Val Pro Asp Arg Phe Thr Gly Ser Gly Tyr Gly Thr Asp Phe Thr Phe Thr Ile Ser Thr Val Gln Ala Glu Asp Leu Ala Val Tyr Phe Cys Gln Gln Asp Tyr Ser Pro Pro Phe Thr Phe Gly Ser Gly Thr Lys Leu Glu Ile Lys Arg Lys Ser Thr Ala Pro Thr Val Ser LIGHTmy @harray= split /\s/, $heavy;my @larray= split /\s/, $light;foreach my $residue(@harray){$residue= $AA321{$residue};}foreach my $residue(@larray){$residue= $AA321{$residue};}print "$comment_H\n",@harray,"\n\n";print "$comment_L\n",@larray;第三章2.提示1 - 1 # 1 minus 1 is 0, would convert to string "0", so false.0.0 # equals 0, would convert to string "0", so false."0.0" # is the string "0.0", neither "" nor "0", so true!"0" # is the string "0", so false."" # is a null string, so false." " # is a whitespace, not null string, so true."0.0" + 0 # would become the number 0 (coerced by the +), so false. \$i # is a reference to $i, so true, even if $i is false, null, 0. undef() # is a function returning the undefined value, so false.3. 办法之一#!D:\Perl\perl\bin\perl.exemy @aa= qw(C D E F G H I K L M N P Q R S T V W Y A);my @pep= qw(A A A A A A A A A);print foreach @pep, "\n";for ( my $i= 0; $i<= 8; $i++ ) {foreach my $res(@aa){$pep[$i] = $res;if ($res ne 'A'){print foreach @pep;}print "\n";}}5、（1）只见于muscle ：grep -v -f skin muscle >muscle.txt 结果见于目录下muscle.txt （2）只见于skin ：grep -v -f muscle skin >skin.txt 结果见于目录下skin.txt（3）都有：grep -Ff muscle skin >skin-muscle.txt 结果见于目录下skin-muscle.txt 6、diff -y laf update >laf-update.xls 结果见于目录下laf-update.xls。