-
Notifications
You must be signed in to change notification settings - Fork 0
/
map1.m
88 lines (86 loc) · 2.25 KB
/
map1.m
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
function [sel1,sel2,mm] = map1(x,y)
%%%%%%% E1
switch x
case 1
mm = 1;
DD0 = [-1];
DD1 = [1];
DD0 = y *DD0;
DD1 = y*DD1;
sel1=mm*DD0;
sel2=mm*DD1;
%%%%%%%%%%% E04
case 2
mm = 4;
DD0 = [-1,1,0,0;0,-1,1,0;0,0,-1,1;0,0,0,-1];
DD1 = zeros(4,4);
DD1(4,1) = 1;
DD0 = y *DD0;
DD1 = y*DD1;
sel1=mm*DD0;
sel2=mm*DD1;
% mm = 3;
% DD0 = [-1,1,0;0,-1,1;0,0,-1];
% DD1 = zeros(3,3)
% DD0 = y *DD0;
% DD1 = y*DD1;
% sel1=mm*DD0;
% sel2=mm*DD1;
%%%%%%%% Hyperexponential
case 3
mm = 2;
DD0 = (1/5.263157894736842)*[-10,0;0,-1];
DD1 = (1/5.263157894736842)*[9,1;0.9,0.1];
DD0 = y*DD0;
DD1 = y*DD1;
sel1=DD0;
sel2=DD1;
%%%%%%%%%%%%% MAP with Negative correlation
case 4
mm=5;
DD0 = [-2.25,2.25,0,0,0;0,-2.25,2.25,0,0;0,0,-2.25,2.25,0;0,0,0,-2.25,0;0,0,0,0,-4.5];
DD1 = [0,0,0,0,0;0,0,0,0,0;0,0,0,0,0;0.0225,0,0,0,2.2275;4.455,0,0,0,0.045];
DD0 = y*DD0;
DD1 = y*DD1;
sel1=DD0;
sel2=DD1;
% %%%%%%%%%%%%% MAP with Negative correlation
% case 5
% mm = 3;
% DD0 = (1/3.811765)*[-2,2,0;0,-81,0;0,0,-81];
% DD1 = (1/3.811765)*[0,0,0;25.25,0,55.75;55.75,0,25.25];
% % DD0 = (1/1.995570321151716)*[-2,2,0;0,-2,0;0,0,-450.50];
% % DD1 = (1/1.995570321151716)*[0,0,0;0.02,0,1.98;445.995,0,4.505];
% DD0 = y *DD0;
% DD1 = y*DD1;
% sel1=DD0;
% sel2=DD1;
% %%%%%%%%%%%%% MAP with Positive correlation
% case 6
% mm = 3;
% DD0 = (1/3.811765)*[-2,2,0;0,-81,0;0,0,-81];
% DD1 = (1/3.811765)*[0,0,0;55.25,0,25.75;25.75,0,55.25];
% % DD0 = (1/1.995570321151716)*[-2,2,0;0,-2,0;0,0,-450.50];
% % DD1 = (1/1.995570321151716)*[0,0,0;1.98,0,0.02;4.505,0,445.995];
% DD0 = y *DD0;
% DD1 = y*DD1;
% sel1=DD0;
% sel2=DD1;
%
% case 7
% mm = 3;
% DD0 = (3/25)*[-15,0,0;0,-5,0;0,0,-10];
% DD1 = (3/25)*[4.5,4.5,6;1.5,1.5,2;3,3,4];
% DD0 = y*DD0;
% DD1 = y*DD1;
% sel1=DD0;
% sel2=DD1;
case 5
mm=5;
DD0 = [-2.25,2.25,0,0,0;0,-2.25,2.25,0,0;0,0,-2.25,2.25,0;0,0,0,-2.25,0;0,0,0,0,-4.5];
DD1 = [0,0,0,0,0;0,0,0,0,0;0,0,0,0,0;2.2275,0,0,0,0.0225;0.045,0,0,0,4.455];
DD0 = y*DD0;
DD1 = y*DD1;
sel1=DD0;
sel2=DD1;
end