ineqelim.m

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

%function MM=IneqElim(P)
global M BM C1 C2
% A*x>b




    
for i=1:size(C1,1)
    P=[C1(i,:);C2];
    
    for p=1:size(P,1)
        
        %%%%% MIN
        TMIN=Linear(P(p,1),P(p,2),i,1);  %TEMP MIN
      
        for k=2:size(C2,1)

                TT=Linear(P(p,1),P(p,2),i,k);
                TMIN=min(TMIN,TT);

        end
        %%%%% MIN
        
        
        for j=1:size(M{i},1)
           
                if(  TMIN== linearDist(P(p,:),M{i}{j}) )  % if M3 output = linearDist   
                    BM{i}{j}(5)=1;      % 1 will not be elim
                end 
           
           
       end
       
    end
    
end




celldisp(BM);