数学建模 DNA序列分类模型（终稿）(7)

if pa<0.1& pa>0.01 pa=pa*10;

zhishu_a=zhishu_a+1; elseif pa<0.01 &pa>1e-3 pa=pa*100;

zhishu_a=zhishu_a+2; end

if pb<0.1 & pb>.01 pb=pb*10;

zhishu_b=zhishu_b+1; elseif pb<.01 pb=pb*100;

zhishu_b=zhishu_b+2; end

pa=round(pa*100)/100;%保留两位有效数字 pb=round(pb*100)/100; end

ee_a=[ee_a zhishu_a]; %记录数据 ee_b=[ee_b zhishu_b]; PA=[PA pa];

PB=[PB pb]; end

ee_ab=ee_a-ee_b;

det_ab=PA-PB.*10.^(ee_ab); %b_type =find(ee_a-ee_b>10); %a_type =find(ee_b-ee_a>10); aa=find(det_ab>0) bb=find(det_ab<0)