duixiao.m

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

function [z,p]=duixiao(sumz,nzn,sump,npn)
sumz=sort(sumz);sump=sort(sump);iz=0;ip=0;nz=0;np=0;
if nzn>0&npn>0
for i=1:nzn
   for j=1:npn
      %c=[i,j]
      if abs(real(sumz(i))-real(sump(j)))<0.0001&abs(imag(sumz(i))-imag(sump(j)))<0.0001
         if j~=1
            %j=j
             if i~=1
               %i=i
                 if sumz(i)==sumz(i-1)&sump(j)==sump(j-1)
                    nz(iz+1)=i;np(ip+1)=j;iz=iz+1;ip=ip+1;   %nz：z中相等的根在z中的下标；iz：z中相等根的个数；
                                                           %np：p中相等的根在p中的下标；ip：p中相等根的个数；
                    break
                 elseif sumz(i)==sumz(i-1)&sump(j)~=sump(j-1)
                    break
                 end
           end
             if sump(j)==sump(j-1)                         %&sumz(i)==sumz(i-1)
              break
             end
         end
         if i~=1&sumz(i)==sumz(i-1)
             break
         end
         nz(iz+1)=i;np(ip+1)=j;iz=iz+1;ip=ip+1;   %nz：z中相等的根在z中的下标；iz：z中相等根的个数；
         break   
      end                                         %np：p中相等的根在p中的下标；ip：p中相等根的个数；
   end
end
else
   z=sumz;p=sump;
end
a=[iz,ip];z=sumz;p=sump;
if iz>=1 
   for i=1:iz
      for j=nz(iz-i+1):nzn-1
         z(j)=z(j+1);
      end
      for j=np(iz-i+1):npn-1
         p(j)=p(j+1);
      end
   end
%end
z=z(1:nzn-iz);p=p(1:npn-iz);
else 
   z=sumz;p=sump;
end