我正在执行以下代码,以根据其中一行(指定行 q)中所有字符的值过滤从一个文件中提取的四行组,并将过滤后的行写入新文件。
fin = fopen(fname_in,'r') ;
fout = fopen(fname_out,'w') ;
y = 0 ;
w = 0 ;
fprintf(1,'Filtering') ;
while ~feof (fin)
w=w+1 ;
if ~mod(w,50000)
fprintf(1,'.');
end
t = fgets(fin) ;
s = fgets(fin) ;
p = fgets(fin) ;
q = fgets(fin) ;
if all(q(1:len) >= phred+phred_code)
y = y + 1 ;
fwrite(fout,[t s([1:len end]) p q([1:len end])]) ;
end
end
fclose(fin) ;
fclose(fout) ;
代码可以运行,但是需要很长时间才能运行。
还有什么其他方法可以加速代码?
谢谢