Objective-Extraction-via-Fuzzy-Clustering-in-Evolutionary-Many-Objective-Optimization/MF.m at master · wangyirui/Objective-Extraction-via-Fuzzy-Clustering-in-Evolutionary-Many-Objective-Optimization · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
function MF(testName,problemDim,objectiveDim)

global mop weightStatus U Cluster;

switch testName
    case 'MF1'
        mop = MF1(mop,problemDim,objectiveDim);
    case 'MF2'
        mop = MF2(mop,problemDim,objectiveDim);
    otherwise
        error('Undefined test problem name');
end

%% MF1 Problem
    function p = MF1(p, ~, ~)
        p.name               = 'MF1';
        p.odim               = 3;
        p.pdim               = 5;
        p.varMin            = 0;
        p.varMax           = 1;
        p.func               =@evaluate;
        function [testRawCost, testCost, testConstrains]= evaluate(x)
            testConstrains = zeros(p.odim - 1,1);
            a = x(1) * x(2);
            b = (1 - x(2)) * x(1);
            c = 1- x(1);
            if a == max(a, max(b,c))
                tmp = a;
            else
                tmp = 0.5 * (a + min(b, c));
            end
            f(1) = tmp + sum( (x(3:p.pdim) - 0.5) .^2 );
            a = x(1) * (1 - x(2));
            b = 1 - x(1);
            c = x(1) * x(2);
            if a == max(a, max(b,c))
                tmp = a;
            else
                tmp = 0.5 * (a + min(b, c));
            end
            f(2) = tmp + sum( (x(3:p.pdim) - 0.5) .^2 );
            a = 1 - x(1);
            b = x(1) * x(2);
            c = x(1) * (1 - x(2));
            if a == max(a, max(b,c))
                tmp = a;
            else
                tmp = 0.5 * (a + min(b, c));
            end
            f(3) = tmp + sum( (x(3:p.pdim) - 0.5) .^2 );
            testRawCost = f';

            % Extra Operation
            if strcmp(weightStatus, 'Initialize')
                testCost = testRawCost;
            else
                tmp1 = sum (f .* U(1, :));
                tmp2 = sum (f .* U(2, :));
                f(1) = tmp1;
                f(2) = tmp2;
                if Cluster == 3
                    tmp3 = sum (f .* U(3, :));
                    f(3) = tmp3;
                    testCost = (f(1:3))';
                else
                    testCost = (f(1:2))';
                end
            end
        end
    end

%% MF2 Problem
    function p = MF2(p, ~, ~)
        p.name               = 'MF2';
        p.odim               = 3;
        p.pdim               = 5;
        p.varMin            = 0;
        p.varMax           = 1;
        p.func               =@evaluate;
        function [testRawCost, testCost, testConstrains]= evaluate(x)
            testConstrains = zeros(p.odim - 1,1);
            a = x(1) * x(2);
            b = (1 - x(2)) * x(1);
            c = 1- x(1);
            if a == min(a, min(b,c))
                tmp = 0;
            else
                tmp = ( a * (a + b + c) ) / ( a + max(b,c));
            end
            f(1) = tmp + sum( (x(3:p.pdim) - 0.5) .^2 );
            a = x(1) * (1 - x(2));
            b = 1 - x(1);
            c = x(1) * x(2);
            if a == min(a, min(b,c))
                tmp = 0;
            else
                tmp = ( a * (a + b + c) ) / ( a + max(b,c));
            end
            f(2) = tmp + sum( (x(3:p.pdim) - 0.5) .^2 );
            a = 1 - x(1);
            b = x(1) * x(2);
            c = x(1) * (1 - x(2));
            if a == min(a, min(b,c))
                tmp = 0;
            else
                tmp = ( a * (a + b + c) ) / ( a + max(b,c));
            end
            f(3) = tmp + sum( (x(3:p.pdim) - 0.5) .^2 );
            testRawCost = f';

            % Extra Operation
            if strcmp(weightStatus, 'Initialize')
                testCost = testRawCost;
            else
                tmp1 = sum (f .* U(1, :));
                tmp2 = sum (f .* U(2, :));
                f(1) = tmp1;
                f(2) = tmp2;
                if Cluster == 3
                    tmp3 = sum (f .* U(3, :));
                    f(3) = tmp3;
                    testCost = (f(1:3))';
                else
                    testCost = (f(1:2))';
                end
            end
        end
    end

end