linearm3_bak2.m

来自「线性M3网络及其示意图」· M 代码 · 共 49 行

function Z=LinearM3(X,Y)
%%GZCM3 function
global M

NC1=size(M,1);
NC2=size(M{1},1);

AllZero=0;% 0 means all M{i}{j}(5) is zero for all j in i

I=1;
for i=1:NC1
    
    for j=1:NC2
       if M{i}{j}(5)==1 
           AllZero=1;
           break;
       end
    end
    
    if j==NC2
        continue;
    end
    
    TMIN=Linear(X,Y,i,j);  %TEMP MIN
      
    for k=(j+1):NC2
        if M{i}{k}(5)==1 
            TMIN=min(TMIN,Linear(X,Y,i,k));
        end
    end
    
    
    MIN(:,:,I)=TMIN;
    I=I+1;
end


MAX=MIN(:,:,1);
if (I-1)>=2
  for i=2:(I-1)
    TMAX=MIN(:,:,i);    %TEMP MAX
    MAX=max(MAX,TMAX);
  end
end


Z=MAX;