## 3 Theory of Proposed Method
### Filtering Outliers by Weight Estimation  
#### 3.1.1 Weight Estimation  
The main idea of weight estimation is that a small weight is given to the measurement with a large residual, which in turn will reduce the influence of the outlier. Therefore, it is crucial to establish the relationship between weights and residuals. However, the residuals of the real model cannot be obtained directly, so the weights are also difficult to be determined. In this work, the AE source is located by the LLS method, and each measurement is treated with equal weights firstly. Then, the residuals of the fitting model can be obtained and used to replace the residuals of the true model:  
\[\Delta_{i}=b_{i}-A_{i}\theta,\quad i=1,2,\cdots,n \tag{6}\]  
where, \(\Delta_{i}\) is the \(i\)th element of a residual vector of the fitting model. After obtaining the residuals, the weight estimation is developed by establishing the relationship between weights and residuals. The simplest method is to exploit the reciprocal value of the residuals as the weights directly. However, when a residual tends to be zero, the corresponding measurement is assigned with an unreasonably high weight, resulting in overfitting. Therefore, the weight should have a maximum upper limit to avoid overfitting, which can be realized by setting a lower limit of \(\lambda_{m}\). Second, residuals exceeding the threshold \(\lambda_{L}\) are considered outliers, and the relative measurements are also regarded as outliers and filtered by setting the weight to zero. Therefore, the weight estimation can be expressed as follows [52]:  
\[w_{i}=\left\{ \begin{array}{l}1/\lambda_{m},\,if\,|\Delta_{i}|\leq\lambda_{m }\\ 1/\Delta_{i},\,if\,\lambda_{m}<\mid\Delta_{i}\mid<\lambda_{L}\\ 0,\quad\quad\quad else\end{array} \right. \tag{7}\]  
To determine the specific parameters of the weight estimation, it is assumed that residuals obey the normal distribution. Therefore, the threshold \(\lambda_{m}\) and \(\lambda_{L}\) in (7) can be defined by:  
\[\lambda_{m}=k_{m}\sigma \tag{8}\]  
\[\lambda_{L}=k_{L}\sigma \tag{9}\]  
where \(k_{m}\) and \(k_{L}\) are the correlation factors, and \(\sigma\) is the estimated value of the standard deviation of residuals, with expression as follows:  
\[\sigma=\sqrt{\frac{1}{n-m}\frac{\sum\limits_{i=1}^{n}w_{i}(b_{i}-A_{i}\theta) ^{2}}{\frac{1}{n}\sum\limits_{i=1}^{n}w_{i}}} \tag{10}\]  
where \(m\) is the number of the unknown parameters. The denominator \(\frac{1}{n}\sum\limits_{i=1}^{n}w_{i}\) is used to standardize the weights. Figure 1 shows the regional division of residuals according to the normal distribution. The parameters \(k_{m}\) and \(k_{L}\) divide the complete residuals into three regions, including minor residuals, medium residuals, and outlier residuals. Besides, the corresponding weights transformed from these residuals are divided into three parts, including an equal weight part, a variable weight part, and a zero weight part, as shown in Figure 2. The first part has the equal and highest weights out of the consideration that at least half of the measurement data are accurate. Therefore, \(k_{m}=0.765\) is selected, so that 50% of the normal residuals fall in the equal weight area. To avoid the unbalance of weighting, \(\lambda_{m}\) should have a lowest limitation equal to \(0.05\text{max}(|\Delta_{i}|)\) when there are many minor residuals. The weights of the second part varies with the relative residuals according to Equation (7). The third part is the zero weight region. Any observations from the \(n\) residuals that lie inside the probability band (\(1\bigtriangleup\)\(1>k_{L}\sigma\)) can be considered outliers. Besides, the corresponding measurements are regarded as outlier measurements, which are filtered on the following calculation by setting the weights to zeros. Figure 1: Residual division according to the normal distribution. In this figure, PDF is the normal probability density function and CDF means the cumulative distribution function.