]> AND Private Git Repository - LiCO.git/blob - PeCO-EO/reponse.tex
Logo AND Algorithmique Numérique Distribuée

Private GIT Repository
ff5a80ce947dbd71d43ca51b9a4cbae0a0cafe71
[LiCO.git] / PeCO-EO / reponse.tex
1 \documentclass[14]{article}
2
3 \usepackage{color}
4 \usepackage{times}
5 \usepackage{titlesec}
6 \usepackage{pifont}
7 %\usepackage[T1]{fontenc}
8 %\usepackage[latin1]{inputenc}
9
10 \renewcommand{\labelenumii}{\labelenumi\arabic{enumii}}
11 %\titleformat*{\section}{\Large\bfseries}
12
13 %\title{Response to the reviewers of \bf "Perimeter-based Coverage Optimization to Improve Lifetime in Wireless Sensor Networks"}
14 %\author{Ali Kadhum Idrees, Karine Deschinkela, Michel Salomon and Raphael Couturier}
15
16 \begin{document}
17
18 \begin{flushright}
19 \today
20 \end{flushright}%
21
22 \vspace{-0.5cm}\hspace{-2cm}FEMTO-ST Institute, UMR 6714
23
24 \hspace{-2cm}University of Franche-Comt\'e
25
26 \hspace{-2cm}IUT Belfort-Montbéliard, BP 527, 90016 Belfort Cedex, France.
27
28 \bigskip
29
30 \begin{center}
31 Detailed changes and addressed issues in the revision of the article
32
33 "Perimeter-based Coverage Optimization to Improve Lifetime in Wireless Sensor Networks"\\
34
35 by Ali Kadhum Idrees, Karine Deschinkel Michel Salomon, Rapahel Couturier
36
37 \bigskip
38 \end{center}
39 Dear Editor and Reviewers,
40
41 First of all, we would like to thank you very much for your kind help to improve our article named:
42 "Perimeter-based Coverage Optimization to Improve Lifetime in Wireless Sensor Networks". We highly appreciate the detailed valuable comments of the reviewers on our article. The suggestions are quite helpful for us and we incorporate them in the revised article. We are happy to submit to you a revised version that considers most of your remarks and suggestions of improvement to improve the quality of our article.
43
44 As below, we would like to clarify some of the points raised by the reviewers and we hope the reviewers and the editors will be satisfied with our responses to the comments and the revision for the original manuscript. \\
45
46
47 %Journal: Engineering Optimization
48 %Reviewer's Comment to the Author Manuscript id GENO-2015-0094
49 %Title: \bf "Perimeter-based Coverage Optimization to Improve Lifetime in Wireless Sensor Networks"
50 %Authors: Ali Kadhum Idrees, Karine Deschinkela, Michel Salomon and Raphael Couturier
51
52
53 \section*{Response to Reviewer No. 1 Comments}
54
55 This paper proposes a scheduling technique for WSN to maximize coverage and network lifetime. The novelty of this paper is the integration of an existing perimeter coverage measure with an existing integer linear programming model. Here are few comments:\\
56
57 \noindent {\bf 1.} The paper makes use of the existing integer optimization model to govern the state of each sensor node within the WSN to maximize coverage and network lifetime. This formulation of the  coverage problem is different from the literature in the sense that they use the perimeter coverage measures to optimize coverage as opposed to the targets/points coverage. The methodology uses existing methods and the original contribution lies only in the application of these methods for the coverage scheduling problem.\\
58
59 \textcolor{blue}{\textbf{\textsc{Answer:} To the best of our knowledge, no integer linear programming based on perimeter coverage has been already proposed in the literature. As specified in the paper, in section 4, it is inspired from one model developed for brachytherapy treatment planning for optimizing dose distribution. In this model the deviation between an actual dose distribution and a required dose distribution in each organ is minimized. In WSN the deviations between the actual level of coverage and the required level are minimized. Outside this parallel between these two applications the mathematical formulation is completly different.                            }}\\
60
61
62 \noindent {\bf 2.} The theory seems mathematically sound. However, the assumption made on the selection criteria for the leader seems too vague.  \\
63
64 \textcolor{blue}{\textbf{\textsc{Answer:}  The selection criteria for the leader inside each subregion is explained in page 8, lines 50-51. After information exchange among the sensor nodes in the subregion, each node will have all required information to decide if it is a leader or not. The decision is based on selecting the sensor node that has a larger number of one-hop neighbors. If this value is the same for many sensors, the node that has the largest remaining energy will be selected as a leader. If there exists sensors with  the same number of neighbors and the same value for the remaining energy, the sensor node that has the largest  index will be selected as a leader. }}\\
65 %{\bf In fact, we gave a high priority to the number of neighbors to reduce the communication energy consumption - PAS CLAIR  }}.\\
66
67
68 \noindent {\bf 3.} The communication and information sharing required to cooperate and make these
69 decisions was not discussed.  \\
70
71 \textcolor{blue}{\textbf{\textsc{Answer:}  The communication and information sharing required to cooperate and make these decisions was discussed in page 8, lines 48-49. Position coordinates, remaining energy, sensor node ID and number of one-hop neighbors are exchanged.}}\\
72
73
74
75 \noindent {\bf 4.} The definitions of the undercoverage and overcoverage variables are not clear. I suggest
76 adding some information about these values, since without it, you cannot understand how M and V are computed for the optimization problem.  \\
77
78 \textcolor{blue}{\textbf{\textsc{Answer:} The perimeter of each sensor may be cut in parts called coverage intervals (CI). The level of coverage of one CI is defined as  the number of active sensors neighbours covering this part of the perimeter. If a given level of coverage $l$ is required  for one sensor, the sensor is said to be undercovered (respectively overcovered) if the level of coverage of one of its CI is less (respectively greater) than $l$. In other terms, we define undercoverage and overcoverage through the use of variables $M_{i}^{j}$ and $V_{i}^{j}$ for one sensor $j$ and its coverage interval $i$. If the sensor $j$ is undercovered, there exists at least one of its CI (say $i$) for which the number of active sensors (denoted by $l^{i}$) covering this part of the perimeter is less than $l$ and in this case : $M_{i}^{j}=l-l^{i}$, $V_{i}^{j}=0$. In the contrary, if the sensor $j$ is overcovered, there exists at least one of its CI (say $i$) for which the number of active sensors (denoted by $l^{i}$) covering this part of the perimeter is greater than $l$ and in this case : $M_{i}^{j}=0$, $V_{i}^{j}=l^{i}-l$.                       }}\\
79
80
81
82 \noindent {\bf 5.} Can you mathematically justify how you chose the values of alpha and beta? This is not
83 very clear. I would suggest possibly adding more results showing how the algorithm performs with different alphas and betas.  \\
84
85 \textcolor{blue}{\textbf{\textsc{Answer:} The choice of alpha and beta should be made according to the needs of the application. Alpha should be enough large to prevent undercoverage and so to reach the highest possible coverage ratio. Beta should be enough large to prevent overcoverage and so to activate a minimum number of sensors. The values of $\alpha_{i}^{j}$ can be identical for all coverage intervals $i$ of one sensor $j$ in order to express that the perimeter of each sensor should be uniformly covered, but  $\alpha_{i}^{j}$ values can be differenciated between sensors to force some regions to be better covered than others. The choice of $\beta \gg \alpha$  prevents the overcoverage, and so limit the activation of a large number of sensors, but as $\alpha$ is  low, some areas may be poorly covered. This explains the results obtained for {\it Lifetime50} with $\beta \gg \alpha$: a large number of periods with low coverage ratio. With $\alpha \gg \beta$, we priviligie the coverage even if some areas may be overcovered, so high coverage ratio is reached, but a large number of sensors are activated to achieve this goal. Therefore network lifetime is reduced. The choice $\alpha=0.6$ and $\beta=0.4$ seems to achieve the best compromise between lifetime and coverage ratio.                }}\\
86
87
88
89 \noindent {\bf 6.} The authors have performed a thorough review of existing coverage methodologies.
90 However, the clarity in the literature review is a little off. Some of the descriptions of the method
91 s used are very vague and do not bring out their key contributions. Some references are not consistent and I suggest using the journals template to adjust them for overall consistency.   \\
92
93 \textcolor{blue}{\textbf{\textsc{Answer:}                               }}\\
94
95
96
97 \noindent {\bf 7.} The methodology is implemented in OMNeT++ (network simulator) and tested against 2 existing algorithms and a previously developed method by the authors. The simulation results are thorough and show that the proposed method improves the coverage and network lifetime compared with the 3 existing methods. The results are similar to previous work done by their team.   \\
98
99 \textcolor{blue}{\textbf{\textsc{Answer:} Although the study conducted in this paper reuses the same protocol presented in our previous work, we focus in this paper on the mathematical optimization model developed to schedule nodes activities.  We deliberately chose to keep the same performance indicators to compare the results obtained with this new formulation with other existing algorithms.                            }}\\
100
101
102 \noindent {\bf 8.}  Since this paper is attacking the coverage problem, I would like to see more information on the amount of coverage the algorithm is achieving. It seems that there is a tradeoff in this algorithm that allows the network to increase its lifetime but does not improve the coverage ratio. This may be an issue if this approach is used in an application that requires high coverage ratio.   \\
103
104 \textcolor{blue}{\textbf{\textsc{Answer:}  Your remark is interesting. Indeed, figures 8(a) and (b) highlight this result. PeCO methods allows to achieve a coverage ratio greater than $50\%$ for many more periods than the others three methods, but for applications requiring an high level of coverage (greater than $95\%$), DilCO method is more efficient.                                 }}\\
105
106 %%%%%%%%%%%%%%%%%%%%%%  ENGLISH and GRAMMER %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
107
108 \noindent\textcolor{black}{\textbf{\Large English and Grammar:}} \\
109
110 \noindent {\ding{90} The first paragraph of every section is not indented.  }  \\
111
112 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed. The first paragraph of every section is indented in the new version. }}\\
113
114
115 \noindent {\ding{90} You seem to be writing in the first person. I suggest rewriting sentences that include “we” “our” or “I” in the third person. (There are too many instances to list them all. They are easily found using the find tool.)  }  \\
116
117 \textcolor{blue}{\textbf{\textsc{Answer:}  It is very common to find sentences with "we" and "our" in scientific papers to explain the work made by the authors. Nevertheless we agree with the reviewer and we reformulated some sentences in the paper to avoid too many uses of the first person. }}\\
118
119
120 \noindent {\ding{90}  Run-on sentence: Page 2 lines 43-48}  \\
121
122 \textcolor{blue}{\textbf{\textsc{Answer:} We rewrote this sentence in two separated sentences.     }}\\
123
124
125
126 \noindent {\ding{90} Add an “and” after the comma on page 3 line 34.}  \\
127
128 \textcolor{blue}{\textbf{\textsc{Answer:}   Right, fixed  }}\\
129
130
131 \noindent {\ding{90}  “model as” instead of “Than” on page 10 line 12.}  \\
132
133 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}\\
134
135
136 \noindent {\ding{90} “no longer” instead of “no more” on page 10 line 31.}  \\
137
138 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}\\
139
140
141 \noindent {\ding{90} “in the active state” add the on page 10 line 34. }  \\
142
143 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}\\
144
145
146
147 \noindent { \ding{90} Lots of English and grammar mistakes. I recommend rereading the paper line by line and adjusting the sentences that do not make sense.} \\
148
149 \textcolor{blue}{\textbf{\textsc{Answer:} ?????? relecture par Ingrid     }}.\\
150
151
152
153
154 \section*{Response to Reviewer No. 2 Comments}
155 The paper entitled "Perimeter-based Coverage Optimization to Improve Lifetime in Wireless Sensor Networks", by Ali Kadhum Idrees, Karine Deschinkela, Michel Salomon and Raphael Couturier proposes a new protocol for Wireless Sensor Networks called PeCO (Perimeter-based Coverage Optimization protocol) that aims at optimizing the use of energy by conjointly exploiting a spatial and temporal subdivision. The protocol is based on solving a Mixed Integer Linear Program at each leader node, and at each iteration of the protocol. The results obtained by PeCO are compared with three other competitors.\\
156
157
158 \noindent\textcolor{black}{\textbf{MAJOR COMMENTS:}} \\
159
160 \noindent {\bf 1.} The protocol framework is not described in details. In particular, the spatial and temporal subdivision (page 2, line 11) that is at the core of PeCO, is not described nor justified in much detail. How to implement an efficient spatial subdivision? On page 10, line 11, the number of subdivisions is said to be equal to 16, but the clustering algorithm used is not mentioned. Is this number dependent of the size of the sensing area? Of the number of sensors? Of the sensing range? The proposed protocol cannot be adopted by practitioners if such an important step is not documented. Temporal subdivision suffers from the same lack of description and justification: why should time intervals have the same duration? If they have the same duration, how should this common duration should be chosen?    \\
161
162 \textcolor{blue}{\textbf{\textsc{Answer:} Spatial and temporal choices of subdivision are not the topics of the paper. In the study, we assume that the deployment of sensors is almost uniformly over the region. So we only need to fix a regular division of the region into subregions to make the problem tractable. The subdivision is made such that the number of hops between any pairs of sensors inside a subregion is less than or equal to 3. Concerning the choice of the sensing period duration, it is correlated with the types of applications, with the amount of initial energy in sensors batteries and also with the duration of the exchange phase. All applications do not have the same requirements of Quality of Service. Here information exchange is executed every hour but the length of the sensing period could be reduced and adapted dynamically. On the one hand a small sensing period would allow to be more reliable but would have higher communication costs. On the other hand the choice of a long duration may cause problems in case of nodes failure during the sensing period.}}\\
163
164
165 \noindent {\bf 2.}Page 9, Section 4, is the Perimeter-based coverage problem NP-hard? This question is important for justifying the use of a Mixed Integer Linear Programming model.   \\
166
167 \textcolor{blue}{\textbf{\textsc{Answer:} The perimeter scheduling coverage problem is NP-hard in general, it  has been proved 
168  in the paper entitled "Perimeter Coverage Scheduling in Wireless Sensor Networks Using Sensors with a Single Continuous Cover Range" from  Ka-Shun Hung and King-Shan Lui (EURASIP Journal on Wireless Communications and Networking 2010, 2010:926075  doi:10.1155/2010/926075). In this paper, authors study the coverage of the perimeter of a large object requiring to be monitored. In our study, the large object to be monitored is the sensor itself (or more precisely its sensing area).                                        }}\\
169
170
171 \noindent {\bf 3.} Page 9, the major problem with the present paper is, in my opinion, the objective function of the Mixed Integer Linear Program (2). It is not described in the paper, and looks like an attempt to address a multiobjective problem (like minimizing overcoverage and undercoverage). However, using a weighted sum is well known not to be an efficient way to address biobjective problems. The introduction of various performance metrics in Section 5.1 also suggests that the authors have not decided exactly which objective function to use, and compare their protocols against competitors without mentioning the exact purpose of each of them. If the performance metrics list given in Section 5.1 is exhaustive, then the authors should mention at the beginning of the paper what are the aims of the protocol, and explain how the protocol is built to optimize these objectives.  \\
172
173 \textcolor{blue}{\textbf{\textsc{Answer:}  As far as we know, representing the objective function as a weighted sum of criteria to be minimized in case of multicriteria optimization is a classical method.                                        }}\\
174
175
176 \noindent {\bf 4.}Page 11 Section 5.2, the sensor nodes are said to be based on Atmels AVR ATmega103L microcontroller. If I am not mistaken, these devices have 128 KBytes of memory, and I didn't find any clue that they can run an operating system like Linux. This point is of primary importance for the proposed protocol, since GLPK (a C API) is supposed to be executed by the cluster leader. In addition to that, GLPK requires a non negligible amount of memory to run properly, and the Atmels AVR ATmega103L microcontroller might be insufficient for that purpose. The authors are urged to provide references of previous works showing that these technical constraints are not preventing their protocol to be implemented on the aforementioned microcontroller. Then, on page 13, in Section "5.2.3 Energy Consumption", the estimation of $E_p^{com}$ for the considered microcontroller seems quite challenging and should be carefully documented. Indeed, this is a key point in providing a fair comparison of PeCO with its competitors.    \\
177
178 \textcolor{blue}{\textbf{\textsc{Answer:}
179 To implement PeCO on real sensors nodes with limited memories capacities, we can act on :
180 \begin{itemize}
181 \item the solver : GLPK is memory consuming for the resolution of integer programming (IP) compared with other commercial solvers like CPLEX\texttrademark. Commercial solvers generally outperform open source solvers  (See "Analysis of commercial and free and open source
182 solvers for linear optimization problems" by B. Meindl and M. Templ from Vienna University of Technology). Memory use depends on the number of variables and number of  constraints.  For linear programs (LP), a reasonable estimate of memory use with CPLEX\texttrademark is to allow one megabyte per thousand constraints. For integer programs, no simple formula exists since memory use depends so heavily on the size of the branch and bound tree. But, the estimate for linear programs still provides a lower bound. In our case, the characteristics of the integer programming (2) are the following:\\
183 \begin{itemize}
184 \item number of variables : $S* (2*I+1)$
185 \item number of constraints : $2* I *S$
186 \item number of non-zero coefficients : $2* I *S * B$
187 \item number of parameters (in the objective function) : $2* I *S$
188 \end{itemize}
189 where $S$ denotes the number of sensors in the subregion, $I$ the average number of cover intervals per sensor, $B$ the average number of sensors involved in a cover interval. The following table gives the memory used with GLPK to solve the integer program (column 3) and its LP-relaxation (column 4) for different problem sizes. The sixth column gives an estimate of the memory used with CPLEX to solves the LP-relaxation according to the number of constraints.
190 \\
191 \begin{tabular}{|c|c|c|c|c|c|}
192 \hline
193 Total number of nodes& S & I & GLPK IP & GLPK LP & number of nodes in the &CPLEX\\
194 of nodes &&&&relaxation &branch-and-bound tree &\\
195 \hline
196 100 & 6.25& 5&0.2 Mb & 0.2 Mb &1 &  64 Kb\\
197 \hline
198 200 & 12.5& 11&1.7 Mb & 1.6 Mb &1 & 281 Kb\\
199 \hline
200 300 &18.5 & 17&3.6 Mb & 3.5 Mb & 3 &644 Kb\\
201 \hline
202 \end{tabular}
203 It is noteworthy that the difference of memory used with GLPK between the resolution of the IP and its LP-relaxation is very weak (not more than 0.1 Mb). The size of the branch and bound tree dos not exceed 3 nodes. This result leads one to believe the memory used with CPLEX for solving the IP would be very close to that for the LP-relaxation, that is to say around 100 Kb for a subregion containing $S=10$ sensors. Moreover the IP seems to have some specifities that encourage us to develop our own solver (coefficents matrix is very sparse) or to use an existing heuristic to find good approximate solution ().
204
205
206
207 \item the subdivision of the region of interest. To make the resolution of integer programming tractable by a leader sensor, we need to limit the number of nodes in each subregion (the number of variables and constraints of the integer programming is directly depending on the number of nodes and neigbors). It is therefore necessary to adapt the subdvision  according to the number of sensors deployed in the area and their sensing range (impact on the number of  cover intervals).  
208 \item heuristic
209
210 \end{itemize}}}\\
211
212
213
214
215 \noindent\textcolor{black}{\textbf{MINOR COMMENTS:}} \\
216
217
218 \noindent {\ding{90} Page 12, lines 7-15, the authors mention that DiLCO protocol is close to PeCO. This should be mentioned earlier in the paper, ideally in Section 2 (Related Literature), along with the detailed description of DESK and GAF, the competitors of the proposed protocol, PeCO.  }  \\
219
220 \textcolor{blue}{\textbf{\textsc{Answer:}       }}.\\
221
222
223
224 \noindent {\ding{90} Page 2, line 20, "An optimal scheduling" should be replaced with "An optimal schedule" }  \\
225
226 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}.\\
227
228
229
230 \noindent {\ding{90}  Page 4, we first read (line 23) "we assume that each sensor node can directly transmit its measurements to a mobile sink", then on line 30, "We also assume that the communication range Rc satisfies $Rc >=2Rs$. In fact, Zhang and Hou (2005) proved that if the transmission range fulfills the previous hypothesis, the complete coverage of a convex area implies connectivity among active nodes.". These two assumptions seems redundant. }  \\
231
232 \textcolor{blue}{\textbf{\textsc{Answer:} Yes, you are right and we removed sentences about the sink. Indeed we consider multi-hops communication.}}.\\
233
234
235
236 \noindent {\ding{90}   Page 4, line 37, a definition for k-covered is missing (the sentence is an equivalence property).}  \\
237
238 \textcolor{blue}{\textbf{\textsc{Answer:} Right. A network area is said to be $k$-covered if every point in the area covered by at least k sensors. We added this definition in the paper}}.\\
239
240
241
242
243 \noindent {\ding{90}  Page 5, lines 34 and 37, replace [0, $2\pi$] with [0, $2\pi$) }  \\
244
245 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}.\\
246
247
248
249
250 \noindent {\ding{90}  Page 5, line 36 and 43, replace "figure 2" with "Figure 2" }  \\
251
252 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}.\\
253
254
255
256
257 \noindent {\ding{90}  Page 5, line 50, replace "section 4" with "Section 4" }  \\
258
259 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}.\\
260
261
262 \noindent {\ding{90}   Page 5, line 51, replace "figure 3" with "Figure 3"}  \\
263
264 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}.\\
265
266
267 \noindent {\ding{90}  Page 7, line 20 "regular homogeneous subregions" is too vague. }  \\
268
269 \textcolor{blue}{\textbf{\textsc{Answer:} As mentioned in the previous remark, the spatial subdivision was not clearly explained in the paper. We added a discussion about this question in the article. Thank you for highlighting it. A FAIRE          }}.\\
270
271
272 \noindent {\ding{90}   Page 7, line 24, replace "figure 4" with "Figure 4"}  \\
273
274 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}.\\
275
276
277 \noindent {\ding{90}  Page 7, line 47, replace "Five status" with "Five statuses" }  \\
278
279 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}.\\
280
281
282
283 \noindent {\ding{90}  Page 9, the constraints of the Mixed Integer Linear Program (2) are not numbered. There are two inequalities for overcoverage and undercoverage that are used to define Mij and Vij. Why not using replacing these inequalities by equalities? }  \\
284
285 \textcolor{blue}{\textbf{\textsc{Answer:}  For minimizing the objective function, $M_{i}^{j}$ and $V_{i}^{j}$ should be set to the smallest possible value such that the inequalities are satisfied. It is explained in the answer 4 for the reviewer 1. But, at optimality, constraints are not necessary satisfied with equality. For instance, if a sensor $j$ is overcovered, there exists at least one of its coverage interval (say $i$) for which the number of active sensors (denoted by $l^{i}$) covering this part of the perimeter is greater than $l$. In this case, $M_{i}^{j}=0$, $V_{i}^{j}=l^{i}-l$, the corresponding inequality $\sum_{k \in A} ( a^j_{ik} ~ X_{k}) + M^j_i  \geq l$ is a strict inequality since $\sum_{k \in A} ( a^j_{ik} ~ X_{k})=l^{i} > l$. }}\\
286
287 \noindent {\ding{90}  Page 10, line 50, "or if the network is no more connected". In order to assess this, the communication range should be known, but it is not given in Table 2. }  \\
288
289 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed}}.\\
290
291
292 \noindent {\ding{90}  Page 10, line 53, the "Coverage ratio" definition is provided for a given period p? Then in the formula on top of page 11, N is set to 51 times 26, why? Is it somehow related to the sensing area having size 50 times 25? }  \\
293
294 \textcolor{blue}{\textbf{\textsc{Answer:} Yes, the "Coverage ratio" definition is provided for a given period p. N is set to 51 times 26 = 1326 grid points because  we discretized the sensing field as a regular grid, a point on the contour and a point every meter. Yes, it is related to the sensing area having size 50 meters times 25 meters. }}\\
295
296
297 \noindent {\ding{90}  Page 11, line 17 in the formula of ASR, |S| should be replaced with J (where J is defined page 4 line 16). }  \\
298
299 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}.\\
300
301
302 \noindent {\ding{90}  Page 13, line 41 and 43, replace "figure 8" with "Figure 8" }  \\
303
304 \textcolor{blue}{\textbf{\textsc{Answer:} Right, fixed }}.\\
305
306
307
308 We are very grateful to the reviewers who, by their recommendations, allowed us to improve the quality of our article.
309
310 \begin{flushright}
311 Best regards\\
312 The authors
313 \end{flushright} 
314  
315
316
317 \end{document}