【正文】 K=K+length(strfind(Sequence1{k},39。K39。))。 E=E+length(strfind(Sequence1{k},39。E39。))。 R=R+length(strfind(Sequence1{k},39。R39。))。 G=G+length(strfind(Sequence1{k},39。G39。))。 F=F+length(strfind(Sequence1{k},39。F39。))。 N=N+length(strfind(Sequence1{k},39。N39。))。 P=P+length(strfind(Sequence1{k},39。P39。))。 Y=Y+length(strfind(Sequence1{k},39。Y39。))。 I=I+length(strfind(Sequence1{k},39。I39。))。 H=H+length(strfind(Sequence1{k},39。H39。))。 C=C+length(strfind(Sequence1{k},39。C39。))。 D=D+length(strfind(Sequence1{k},39。D39。))。 sum_Aa=sum_Aa+length(Sequence1{k})。 X=X+length(strfind(Sequence1{k},39。X39。))。 B=B+length(strfind(Sequence1{k},39。B39。))。 Z+length(strfind(Sequence1{k},39。Z39。))。 U=U+length(strfind(Sequence1{k},39。U39。))。 end sum_Aa=sum_AaXBZU。 Paa(1,2)=A/sum_Aa。 Paa(2,2)=R/sum_Aa。 Paa(3,2)=D/sum_Aa。 Paa(4,2)=C/sum_Aa。 Paa(5,2)=Q/sum_Aa。 Paa(6,2)=E/sum_Aa。 Paa(7,2)=H/sum_Aa。 Paa(8,2)=I/sum_Aa。 Paa(9,2)=G/sum_Aa。 Paa(10,2)=N/sum_Aa。 Paa(11,2)=L/sum_Aa。 Paa(12,2)=K/sum_Aa。 Paa(13,2)=M/sum_Aa。 Paa(14,2)=F/sum_Aa。 Paa(15,2)=P/sum_Aa。 Paa(16,2)=S/sum_Aa。 Paa(17,2)=T/sum_Aa。 Paa(18,2)=W/sum_Aa。 Paa(19,2)=Y/sum_Aa。 Paa(20,2)=V/sum_Aa。 A=0。 W= 0。 M= 0。 S= 0。 T= 0。 V=0。 Q=0。 L=0。 K=0。 E=0。 R=0。 G=0。 F=0。 N=0。 P=0。 Y=0。 I=0。 H=0。 C=0。 D=0。 X=0。 B=0。 Z=0。 U=0。 sum_Aa=0。[seq_header2,Sequence2]=fastaread(39。文件3位置\39。)。for k=1:length(seq_header2) A=A+length(strfind(Sequence2{k},39。A39。))。 W=W+length(strfind(Sequence2{k},39。W39。))。 M=M+length(strfind(Sequence2{k},39。M39。))。 S=S+length(strfind(Sequence2{k},39。S39。))。 T= T+length(strfind(Sequence2{k},39。T39。))。 V=V+length(strfind(Sequence2{k},39。V39。))。 Q=Q+length(strfind(Sequence2{k},39。Q39。))。 L=L+length(strfind(Sequence2{k},39。L39。))。 K=K+length(strfind(Sequence2{k},39。K39。))。 E=E+length(strfind(Sequence2{k},39。E39。))。 R=R+length(strfind(Sequence2{k},39。R39。))。 G=G+length(strfind(Sequence2{k},39。G39。))。 F=F+length(strfind(Sequence2{k},39。F39。))。 N=N+length(strfind(Sequence2{k},39。N39。))。 P=P+length(strfind(Sequence2{k},39。P39。))。 Y=Y+length(strfind(Sequence2{k},39。Y39。))。 I=I+length(strfind(Sequence2{k},39。I39。))。 H=H+length(strfind(Sequence2{k},39。H39。))。 C=C+length(strfind(Sequence2{k},39。C39。))。 D=D+length(strfind(Sequence2{k},39。D39。))。 sum_Aa=sum_Aa+length(Sequence2{k})。 X=X+length(strfind(Sequence2{k},39。X39。))。 B=B+length(strfind(Sequence2{k},39。B39。))。 Z+length(strfind(Sequence2{k},39。Z39。))。 U=U+length(strfind(Sequence2{k},39。U39。))。 end sum_Aa=sum_AaXBZU。 Paa(1,3)=A/sum_Aa。 Paa(2,3)=R/sum_Aa。 Paa(3,3)=D/sum_Aa。 Paa(4,3)=C/sum_Aa。 Paa(5,3)=Q/sum_Aa。 Paa(6,3)=E/sum_Aa。 Paa(7,3)=H/sum_Aa。 Paa(8,3)=I/sum_Aa。 Paa(9,3)=G/sum_Aa。 Paa(10,3)=N/sum_Aa。 Paa(11,3)=L/sum_Aa。 Paa(12,3)=K/sum_Aa。 Paa(13,3)=M/sum_Aa。 Paa(14,3)=F/sum_Aa。 Paa(15,3)=P/sum_Aa。 Paa(16,3)=S/sum_Aa。 Paa(17,3)=T/sum_Aa。 Paa(18,3)=W/sum_Aa。 Paa(19,3)=Y/sum_Aa。 Paa(20,3)=V/sum_Aa。 bar(Paa) xlabel(39。amino acid39。)。 ylabel(39。Frequence39。)。 set(gca,39。XTickLabel39。,{39。A39。,39。W39。,39。M39。,39。S39。,39。T39。,39。V39。,39。Q39。,39。L39。,39。K39。,39。E39。,39。R39。,39。G39。,... 39。F39。,39。N39。,39。P39。,39。Y39。,39。I39。,39。H39。,39。C39。,39。D39。},39。XTick39。,1:1:20)。 set(gca,39。xlim39。,[,])。 h=legend({39。所有39。,39。無序區(qū)39。,39。有序區(qū)39。},1)。 set(h,39。location39。,39。NorthEast39。) applyhatch(gcf,39。\.x.39。)。 %hold on 23