% Пошук Парето-оптимальних варіантів структури вимірювальної системи

% Матриця значень експертних оцінок ефективності функціонування вимірювальної системи
% для 7 варіантів структури (рядки) та 6 критеріїв оптимальності (стовбці)  
% Оцінки за 9 бальною шкалою, 1-min; 9-max

% x - початкові оцінки 
% px - Парето-оптимальні оцінки

function [px pxvar]=paretov(x, xvar);
    n=length(x(:,1));
    i=1;
    while i<n
        for j=i+1:n
            m=min(x(i,:)-x(j,:));
            q=min(x(j,:)-x(i,:));
            if m>=0
                x(j,:)=[];
                xvar(j,:)=[];
                i=i-1;
                break
            elseif q>=0
                x(i,:)=[];
                xvar(i,:)=[];
                i=i-1;
                break
            else
            end
        end
        i=i+1;
        n=length(x(:,1));
    end
    % Парето-оптимальні варіанти
    px=x;
    pxvar=xvar;