genpai.m

详细介绍了Mathlab的使用方法与技巧

function [l,li,r,ri,f,fi,x,xi]=genpai(gen)
gen=sort(gen);genn=length(gen);l=0;li=0;r=0;ri=0;f=0;fi=0;x=0;xi=0;
for i=1:genn
    if abs(real(gen(i)))<0.00001
       gen(i)=gen(i)-real(gen(i));
    end
    if abs(imag(gen(i)))<0.00001
       uuuu=real(gen(i));
       gen(i)=uuuu;
   end
end
%gen
for i=1:genn
   if abs(real(gen(i)))<0.00001&abs(imag(gen(i)))<0.00001           %判num的零根，sumzi0：零根在z中的下标；%lingzi0i：零根的个数；
      l(li+1)=i;li=li+1;
   else if abs(real(gen(i)))>=0.00001&abs(imag(gen(i)))<0.00001  %判num的实数根，sumzi0：实数根在z中的下标；lingzi0i：实数根的个数；
         r(ri+1)=i;ri=ri+1;
   else if abs(real(gen(i)))>=0.00001&abs(imag(gen(i)))>=0.00001  %判num的复数根，sumzi0：复数根在z中的下标；lingzi0i：复数根的个数；
         f(fi+1)=i;fi=fi+1;
      else
         x(xi+1)=i;xi=xi+1;                  %虚根
      end
   end
end 
end

if li>=1
   for i=1:li
      l(i)=0;                                           %l：零根矩阵
   end
else
   l=[];
end
if ri>=1
   for i=1:ri
      r(i)=real(gen(r(i)));                                           %l：实数根矩阵
   end
   else
   r=[];
end
if fi>=1
   for i=1:fi
      f(i)=gen(f(i));                                           %l：复数根矩阵
   end
   else
   f=[];
end
if xi>=1
   for i=1:xi
      xu(i)=imag(gen(x(i)));                                           %l：虚根矩阵
   end
   x=j*xu;
else
   x=[];
end
l=sort(l);r=sort(r);