السلام عليكـم …
مسـاء الخير
لو سمحتو أنـا اريد شـرح الكود هذا
وشكراً لكــم اذا تكرمتــو
for jj=1:lengthofseq:50*lengthofseq
the_test_seq=the_all_new_seq(:,jj:jj+(lengthofseq-1))
the_corrected_seq(1,:)=the_test_seq(1,:);
gf=2;
while gf<=numseq
the_line_seq=the_test_seq(gf,:)
[COLOR=red]the_ref_seq=The_original(gf,:);[/color]
[COLOR=red]gapsloc=findstr(the_line_seq,'_');[/color]
[COLOR=red]the_line_seq(gapsloc)=[];[/color]
x=strcmp(the_ref_seq,the_line_seq);
if x==1
disp(' it is okey')
the_line_seq=the_test_seq(gf,:);
else
disp('it must be corrected')
the_line_seq=the_test_seq(gf,:);
[the_line_seq]=correct(the_line_seq,the_ref_seq,numofgap)
end
the_corrected_seq(gf,:)= the_line_seq;
gf=gf+1;
end
the_corrected_seq
score_of_the_corected_seq=calculatescore(the_corrected_seq,1,-1,-2)
the_all_new_corected_seq(:,jj:jj+(lengthofseq-1))=the_corrected_seq;
end
the_all_new_corected_seq
وخاصه الكود المكتوب بالون الاحمر , وشكراً