-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathRandTauRelimit.m
More file actions
147 lines (134 loc) · 6.64 KB
/
Copy pathRandTauRelimit.m
File metadata and controls
147 lines (134 loc) · 6.64 KB
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
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
% check if tau_i and/or R_i need to be relimited. If yes --> inform the user
tau1_limit = str2num(evalin('base','tau1_limit'));
tau2_limit = str2num(evalin('base','tau2_limit'));
tau3_limit = str2num(evalin('base','tau3_limit'));
Ri_limit = str2num(evalin('base','Ri_limit'));
if NRC == 1
if isempty(tau1_limit) == 1
if tau_1_reco >= 1000 * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau1 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
else
if tau_1_reco >= tau1_limit * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau1 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
end
RR1 = prctile (R(1,:),50); % median of R1
if isempty(Ri_limit) == 1
if RR1 >= 0.3 * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for Ri and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
else
if RR1 >= Ri_limit * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for Ri and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
end
end
if NRC == 2
if isempty(tau1_limit) == 1
if tau_1_reco >= 50 * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau1 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
else
if tau_1_reco >= tau1_limit * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau1 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
end
if isempty(tau2_limit) == 1
if tau_2_reco >= 1500 * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau2 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
else
if tau_2_reco >= tau2_limit * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau2 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
end
RR1 = prctile (R(1,:),50); % median of R1
RR2 = prctile (R(2,:),50); % median of R2
if isempty(Ri_limit) == 1
if (RR1 >= 0.3 * 0.95) || (RR2 >= 0.3 * 0.95)
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for Ri and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
else
if (RR1 >= Ri_limit * 0.95) || (RR2 >= Ri_limit * 0.95)
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for Ri and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
end
end
if NRC == 3
if isempty(tau1_limit) == 1
if tau_1_reco >= 10 * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau1 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
else
if tau_1_reco >= tau1_limit * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau1 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
end
if isempty(tau2_limit) == 1
if tau_2_reco >= 100 * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau2 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
else
if tau_2_reco >= tau2_limit * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau2 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
end
if isempty(tau3_limit) == 1
if tau_3_reco >= 1000 * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau3 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
else
if tau_3_reco >= tau3_limit * 0.95
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for tau3 and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
end
RR1 = prctile (R(1,:),50); % median of R1
RR2 = prctile (R(2,:),50); % median of R2
RR3 = prctile (R(3,:),50); % median of R3
if isempty(Ri_limit) == 1
if (RR1 >= 0.3 * 0.95) || (RR2 >= 0.3 * 0.95) || (RR3 >= 0.3 * 0.95 )
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for Ri and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
else
if (RR1 >= Ri_limit * 0.95) || (RR2 >= Ri_limit * 0.95) || (RR3 >= Ri_limit * 0.95)
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
disp(' Please increase the upper limit for Ri and rerun the codes')
disp(' ======CAUTION======CAUTION======CAUTION======CAUTION======')
end
end
end