zsly.m

来自「这是本人在2005年参加全国大学生电工建模做b题的相关资料」· M 代码 · 共 37 行

function ZB=zsly(b)
lp=find(b);
t=1;
[H,L]=size(b);
lpc=length(lp);
for i=1:H
   for j=t:lpc
       if  lp(j)/L>i
           a(i)=j-t;
           t=a(i)+1;
           break;
       end
   end
end
t=1;
for i=1:L
   for j=t:lpc
       if  lp(j)/H>i
           c(i)=j-t;
           t=c(i)+1;
           break;
       end
   end
end           
    X=min(a);
    Y=min(b);
    for i=1:H
        if a(i)==X
            break;
        end
    end
    for j=1:L
        if b(j)==Y
            break;
        end
    end
    ZB=[i,j];