The article analyses why colleges and universities should strengthen innovation and entrepreneurship education based on ‘mass entrepreneurship and innovation.’ First, we conduct a questionnaire survey on the status quo of college students’ innovation and entrepreneurship and use nonlinear methods to construct an evaluation model of innovation and entrepreneurship capabilities to evaluate students’ innovation and entrepreneurship capabilities quantitatively. Finally, we verify the effectiveness of the combined evaluation model through the data on the innovation and entrepreneurship of college students. The research results can provide a new idea for appraisal of college students’ innovation and entrepreneurship ability.
Keywords
 nonlinear model
 analytic hierarchy process
 colleges and universities
 innovation and entrepreneurship ability evaluation
MSC 2010
 62J12
As a distributed vibration monitoring technology, phasesensitive optical timedomain reflectometer (ΦOTDR) has the advantages of corrosion resistance, electromagnetic interference resistance, high sensitivity, wide detection range [1], and easy concealment [2] etc. ΦOTDR has been successfully used in perimeter security [3], protection of important underground infrastructures [4] and health monitoring of large buildings [5]. ΦOTDR technology is also a new way of protecting the underground cable. Because most of the cables are buried in areas of high population density, the behaviour of stealing and damaging cable will bring great inconvenience to people's life and also cause great economic losses. ΦOTDR, which is used directly on the paved communication optical fibre without additional sensing devices, does not affect the communication on fibre and also protect the fibre itself [6]. Therefore, ΦOTDR is a convenient, high costeffective cable protection way.
The primary problem of vibration monitoring is the localisation of vibration events, but the premise of localisation is to measure the distance between the vibration source and fibre. In the past, the research on ΦOTDR technology focused on the improvement of the detection range [7] as well as the sensitivity [8,9,10]. As the technology is maturing, the distance measurement to vibration signal through ΦOTDR becomes possible.
However, the distance measurement accuracy on ΦOTDR is low due to its high sensitivity along with other factors such as the inherent noises on devices, environmental noises, spatial resolution and human activities’ interference. Thus, this paper focuses on efficient distance measurement for vibration events.
This paper presents the method to determine the minimum distance between vibration source and fibre under existing optical vibration signals.
Some communication fibres are distributed in the dedicated pipe chase, nevertheless, some are directly embedded under the ground. Fibres in pipe chase are not directly contacted with the soil layers. Thus the capacity of sensing vibration is weak. In comparison, fibres under the ground have a stronger ability to sense vibration. The research of this paper is only limited to distance measurement between vibration source and fibre under the ground. In practice, evacuation is the main threat to optical fibres. Evacuation can be divided into shovel evacuation, excavator evacuation and tamping machine evacuation, etc. The experiment of this paper concentrated on shovel excavation due to its high frequency.
Vibration signals from excavation conform to the propagation law in the soil layers. The attenuation is in the exponential law [11], which means that the farther the signal is from the vibration source, the more it attenuates.
The spatial resolution of distributed fibres refers to the minimum spatial distance as measuring vibration along with the fiber distribution. One fibre can be divided into different segments. Each segment is equivalent to a sensor, and the length is considered as spatial resolution. For example, when spatial resolution is 1 m, vibration signals are distinguished on two fibre segments with an interval of 1 m.
After thorough consideration of spatial resolution and vibration attenuation law, the propagation process of distributed fibre for vibration can be described in Figure 1. The segment closest to the vibration source receives the strongest signal, vice versa. For convenient description, the fibre segment closest to the vibration source is numbered as 0—the 0th fibre or the datum fibre segment. The adjacent fibre segments are numbered as −1 and 1—the −1st fibre and the 1st fibre, and so on. Fibres are marked by using the subscript subsequence (−N. . . −1, 0, 1 . . . N). For convenient description, such a method is called MMDFS.
The vibration source closest to the fibre segment is l_{0}, the distance from each segment to the vibration source is marked as l_{i}, (i= −N, . . . 0, . . . , N). Δt_{i} refers to TDOA between the ith fibre and the 0th fibre. Δt_{0}, the TDOA of 0th fibre, is equal to 0. Δt_{i} is >0. v represents propagation velocity of vibration. d stands for spatial resolution. The relationship between Δt_{i} and l_{i} is
The goal of this paper is to substitute TDOA Δt_{i} into Eq. (4) to solve l_{0}. Since Eq. (4) is a nonlinear equation set, solving l_{0} directly is complicated. Therefore, Eq. (4) is firstly linearised, then the least square principle is used to deduce a solving equation. To ensure the estimation accuracy of Δt_{i}, the method of TDOA estimation and eliminating outliers are studied.
Distance Estimation between Vibration Source and Fibre Based on the Least Square
A variable substitution is made as follows
Accurate TDOA estimation is the key to measure the distance between vibration source and fibre. The classical method for TDOA estimation is GCC. In recent years, GCCSC has been developed to conduct correlation again based on GCC. Under the low SNR, the results from the secondary relation are much better than GCC. We conducted the simulation comparison experiments by aiming at five methods, including crosscorrelation (CC), SCOT, ML, GCCSC based on SCOT, and GCCSC based on ML. The findings of the comparative experiments indicate that the secondary relation estimation based on ML weighting has the maximum accuracy.
A signal emitting from a remote source and monitored in the presence of noise at two spatially separated receivers can be mathematically modelled as
GCC is the improvement on the basic crosscorrelation method, and it can be used under the low SNR. GCC's principle is prefiltering to eliminate noises and disturbance and then carrying out correlation operations.
According to the different weighting methods, GCC can be divided into four types, including Roth, SCOT, PHAT and ML. The ML algorithm can give the big weight to the frequency bands with the high SNR and the small weight to the bands with the small SNR, therefore it is the optimal filter statistically.
The generalised correlation between x_{1} (n) and x_{2} (n) is
The process of GCCSC is that the signals are autocorrelated and crosscorrelated firstly, and then the results of autocorrelation and crosscorrelation are used to carry out crosscorrelation operation again, so as to improve resolution power and antinoise performance.
Autocorrelation function of x_{1} (n) is
To improve the estimation accuracy, two aspects of the secondary relation are improved: one is limiting the range of delay peak; the other is using threepoint interpolation to improve the delay estimate.
The method of calculating Δt_{i} has been discussed above. However, in practice, affected by noises and interference, the estimated delay value Δt_{i} often sometimes exceeds the normal range. These points with apparent deviation are named as outliers. Too many outliers can affect the accuracy of estimation, and the removal of outliers can significantly improve the estimation accuracy.
Judging whether a data point is an outlier, it is necessary to select a suitable indicator. In this paper, the Deviation Square (DS) is used to judge whether Δt_{i} is an outlier. The definition is
According to Figure 2(a), when i is −1, 1, −6 and 6, the four DSi is relatively large. All DSi in Figure 2(a) are calculated and shown in Figure 2(b). In Figure 2(b), the four DSi marked with solid points are apparent outliers, and the rest 11 hollow points are more likely to be normal.
In Figure 2(b), the DS value of outliers is higher than normal points. The normal points gather in a specific region with smaller DS values, and the number is more. However, the outliers cluster into one or more classes in the region with larger DS values, and the number is less. Thus DS can be used to distinguish the normal points from outliers effectively.
This paper adopts the idea of eliminating outliers by clustering. The primary thought is as follows: First, all points are clustered into multiple classes according to the DS value. Second, the most likely outliers are determined by judging the number of points in the class and the distance between the classes. Third, and ared after removal of such points. Fourth, the new DS value is calculated. The above steps are repeated until the remaining data can be considered as one class. That is, no apparent outliers.
There are multiple clustering methods, including KMeans, hierarchical clustering, Clustering with Gaussian Mixture Models (GMM) and DBSCAN. However, this paper requires that only one class remain in the end, no matter how many classes there are at the beginning. KMeans and GMM cannot meet the requirement of changeable clustering number, and hierarchical clustering is not easy to determine the dividing layers. As a result, DBSCAN is used in this paper.
DBSCAN has two advantages: (1) High efficiency. It can remove multiple outliers at one time. Compared with other clustering methods, iterations are fewer. (2) The clustering number is reduced with the increase of iterations, and finally into one class.
Before using DBSCAN, it is necessary to set two parameters, one is the cluster radius
Considering all the above problems, a complete method to calculate the distance between the vibration source and the optical fibre is described below.
The experiment is divided into two parts, including simulation fibre data and practical fibre data.
Ricker wavelets were used as the analogue vibration source signal. According to the attenuation coefficient and the distance between the optical fibre and vibration source, the attenuated signal was calculated. Then the attenuated signal was mixed with zero mean Gaussian noise to obtain the simulated vibration signal received by fibre.
Ricker wavelet is the typical waveform in vibration simulation. The timedomain expression is shown as follows
Figure 3 shows that the Ricker wavelet is not only similar to real vibration signals in the time domain but also in the frequency domain. Therefore, the Ricker wavelet is the desirable waveform of simulation vibration signals.
The vibration source data was composed of five Ricker wavelets with a frequency peak of 50 Hz. To improve the data quality, Ricker wavelet amplitude was expanded five times. After attenuation, the vibration source data was mixed with Gaussian noise with zero mean and 0.05 variance. The sampling frequency was set as 10000 Hz, the duration was 3 s, the wave velocity was 10000 m/s, and the minimum distance between the optical fibre and vibration source was 1 m. There were 19 fibre segments receiving vibration signals, and their serial numbers based on MMDFS were (−9, −8. . . −1, 0, 1. . . 8, 9).
At the same time, all data was stored as a matrix, with 30,000 rows and 19 columns. The stored data is shown in Figure 4.
In Figure 4(b), each column refers to data received by one fibre segment. Among these columns, the 10th column is the datum fibre segment and has the best signal quality with an SNR of −1.05 db. But the first and 19th columns have the worst data quality with an SNR of −2.12 db. Five Ricker wavelets in signals of No. 10 column can be observed obviously, but the Ricker wavelets are extremely close to noises and challenging to distinguish in signals of No. 1 column. The simulation fibre data used in this paper is named as SimuFiberData.
The experiment used the method of GCCSC based on ML to estimate TDOA. To improve accuracy, threepoint interpolation was conducted. Table 1 shows the simulation fiber data for making a comparison on the real TDOA with five methods, including CC, SCOT, ML, GCCSC based on SCOT and GCCSC based on ML. It is worth noting that the unit of TDOA is the sampling number, instead of the second. Mean Square Error and success rate of five methods in Table 1 are shown in Table 2.
Comparison of TDOA estimation for simu fibre data.
−9  7  −2  7.3147  7.3945  7.9716  8.055 
−8  8  8  8.2805  8.8228  7.4851  7.062 
−7  5  5  4.9436  4.3144  8.02  6.071 
−6  7  −7  6.9048  6.3663  5.464  5.083 
−5  9  −6  6.0667  8.9716  5.7098  4.099 
−4  3  −3  2.7389  2.8608  3.211  3.123 
−3  5  −1  −1.1754  −1.2043  2.3093  2.162 
−2  5  −10  5.2752  5.0414  2.1822  1.236 
−1  1  1  0.7914  0.9828  0.8831  0.414 
0  0  0  0  0  0  0 
1  0  13  −0.016  0.1585  0.3109  0.414 
2  2  3  2.8405  2.7737  1.5947  1.236 
3  3  −14  3.6201  3.1296  2.7702  2.162 
4  1  9  0.9783  1.2096  1.4064  3.123 
5  4  −8  3.9291  4.1022  4.1073  4.099 
6  8  14  7.8445  7.9643  6.8807  5.083 
7  11  −13  11.2104  11.013  6.5797  6.071 
8  3  3  11.0275  3.4893  8.4726  7.062 
9  6  −10  6.4191  6.4771  7.1539  8.055 
CC, crosscorrelation.
Mean square error and success rate of five methods of Table 1.
MSE  5.943  17.0042  5.2276  5.6606  0.9272 
Success Rate  100%  42.1053%  89.4737%  94.7368%  100% 
CC, crosscorrelation.
TDOA estimation in Table 1 has two unreasonable situations: (1) TDOA value is negative. For example, the estimation value of No. −9 in SCOT is −2; (2) TDOA value dramatically exceeds the theoretical value. For example, the estimation value of No. 1 in SCOT is 13, but the actual value is 1.236. These belong to unsuccessful estimation values. Thus the success rate is also an important indicator to estimate TDOA.
The success rate of the five methods in Table 2 shows that GCCSC based on ML has the highest success rate than other methods. The mean square error estimated by GCCSC based on ML is the minimum and it has the highest success rate. For this reason, GCCSC based on ML was selected to do TDOA estimation.
To obtain satisfactory distance estimation results, not only the accurate TDOA method should be selected, but also the TDOA outliers should be removed reasonably. The entire detailed computational process of SimuFiberData is shown in Table 3.
Complete process to simu fibre data.
1  8.055  7.9716  −9  248.021  0  0 
2  7.062  7.4851  −8  28.806  0  0 
3  6.071  8.02  −7  286.574  0  0 
4  5.083  5.464  −6  1.424  −6  0 
5  4.099  5.7098  −5  153.47  0  0 
6  3.123  3.211  −4  0.612  −4  1.402 
7  2.162  2.3093  −3  0.191  −3  0.03 
8  1.236  2.1822  −2  22.245  0  0 
9  0.414  0.8831  −1  2.88  −1  1.612 
10  0  0  0  0  0  0 
11  0.414  0.3109  1  0.032  1  1.28E−01 
12  1.236  1.5947  2  2.861  2  1.199 
13  2.162  2.7702  3  10.526  3  7.343 
14  3.123  1.4064  4  124.818  0  0 
15  4.099  4.1073  5  8.236  5  7.932 
16  5.083  6.8807  6  221.649  0  0 
17  6.071  6.5797  7  3.14  7  0.673 
18  7.062  8.4726  8  70.983  0  0 
19  8.055  7.1539  9  711.649  0  0 
DSMSE  105.4509  2.2577  

8508.390  9370.197  

1.419  0.958 
DS, deviation square.
The 1st column in Table 3 is the serial numbering for TDOA. Successively, the 2nd column is the true TDOA; the 3rd column is the TDOA estimation value according to GCCSC based on ML. After the 3rd column, data in two columns are divided into a group. Each group represents one distance estimation. No. 1 column of each group is the serial numbering based on MMDFS. If the serial numbering is 0, it means data at this point have been eliminated from the distance estimation. No. 2 column in each group is the corresponding DS of data involved in TDOA estimation. The antepenultimate row is the mean square error of all DS values after each distance estimation, implying the effect of DS clustering. The penultimate row is the estimation value of wave velocity. The tailender row is the estimation value of the distance between the vibration source and optical fibre.
In the first estimation, wave velocity is 8508.390, the distance is 1.419, both of which belong to the reasonable value range. However, DS has the more significant mean square error and it is 105.451, indicating that the clustering result is not desirable. The TDOA value involved in estimation probably has outliers, which should be eliminated.
The primary step of eliminating outliers was to set up cluster radius
After removing outliers in the first time, the second estimation was conducted to obtain wave velocity as 9370.197 and distance as 0.958, close to the set value 1; DS mean square error was very small, that is 2.2577, indicating that the clustering result was desirable and outliers had been eliminated. Figure 5(b) shows the scatter diagram of all DS values in the No. 7 column in Table 3—the scatter diagram after removing outliers, showing that the DS values of residual points were maintained at 10 and there were no outliers. The whole calculation was over.
Fibre data were collected from the space in a college in Beijing, China. The space was an area of 800 m2 in a rectangle. The southnorth length was about 20 m and the southeast length was about 40 m. The length of fibre was 300 m, paved around the site with two circles. Its buried depth was about 10 cm under the ground.
The fibre was wireline cable. NBXS3000 from the Japanese Neubrex Company was used as the fibre collecting device. The sampling frequency was 10000 Hz. Fibre spatial resolution was 1.03 m, thus the entire fibre could be divided into 290 segments.
The excavation was carried out at 108.83 m of the fibre. A shovel was used to excavate on the ground. The excavation position was distant from the fibre about 1 m. The strongest excavation signal was at 108.83 m. Excavation signals could be distinguished at 99.59 m and 117.04 m. But it was difficult to distinguish farther excavation signals with the naked eyes. The vibration signals received from 99.59 m to 117.04 m were chosen as the analysis objects. According to MMDFS, the fibre segment at 108.83 m was the datum fibre segment with the serial numbering as 0. The serial numbering at 109.86 m, 117.04 m, 107.8 m and 99.59 m was 1, 8, −1 and −9, respectively. There were a total of 18 fibre segments, which had the sequence numbered as (−9, −8, −7, −6, −5, −4, −3, −2, −1, 0, 1, 2, 3, 4, 5, 6, 7, 8).
Practical fibre data were greatly affected by the background noises, and the SNR was low. The best SNR obtained from the strongest excavation signal at 108.83 m was 17.72 dB, the worst SNR was 9.49 dB. At 117.04 m, the best SNR was 8.45 dB; the worst SNR was −2.27 dB.
Relative to the simulation fibre data, practical fibre data were affected by noises and soil medium, resulting in greater deviation. As estimating TDOA, more unreasonable values than simulation fibre data were generated, sometimes accounting for >30% of the total data. Therefore, the main task of pretreatment was to remove unreasonable TDOA.
As shown in Figure 6, the fibre data contains four groups of apparent excavation signals circled with dotted lines. In distance estimation, four segments of data are partitioned. The length of each segment ranges from 130 to 200 sampling points. Excavation data in each group regards the No. 0 column as the benchmark to estimate TDOA with other columns and eliminate unreasonable TDOA, remaining 25 effective values.
Because the method adopted for estimating TDOA has strong noise resistance, this paper does not carry out denoising processing on optical fibre data.
The complete details of distance estimation for excavation data are recorded in Table 4.
Differing from Table 3, because true TDOA is unknown, there is no real value column in Table 4. The rest is the same as Table 4. The original number ranges from 1 to 72. After pretreatment, only 25 valid TDOA values remained. No. 2 column is the TDOA estimation value. The entire calculation conducted distance estimation five times. Outliers were eliminated four times. The setting of cluster radius
The overall process of excavation data treatment.
11  0.30  1  76.82  1  22.90  1  5.46  1  1.83  1  0.02 
12  0.40  2  73.13  2  11.87  2  0.09  2  0.92  2  6.18 
13  0.45  3  20.16  3  1.15  3  20.17  3  34.44  0  0 
14  0.61  4  1.15  4  34.72  4  103.16  0  0  0  0 
15  0.71  5  40.03  5  195.43  0  0  0  0  0  0 
16  0.82  6  246.77  6  577.32  0  0  0  0  0  0 
17  0.81  7  889.91  0  0  0  0  0  0  0  0 
18  0.93  8  1902.98  0  0  0  0  0  0  0  0 
20  0.95  −8  1878.44  0  0  0  0  0  0  0  0 
21  1.94  −7  504.32  0  0  0  0  0  0  0  0 
22  1.61  −6  63.24  −6  268.29  0  0  0  0  0  0 
23  1.56  −5  13.22  −5  24.83  −5  101.21  0  0  0  0 
24  1.05  −4  81.91  −4  0.01  −4  31.30  −4  62.44  0  0 
25  0.28  −3  0.08  −3  16.31  −3  40.42  −3  53.09  0  0 
26  0.32  −2  37.20  −2  3.69  −2  0.42  −2  2.83  −2  8.66 
27  0.25  −1  52.07  −1  14.75  −1  3.11  −1  0.87  −1  0.01 
47  0.53  1  228.79  1  77.38  1  24.15  1  11.17  1  2.14 
48  0.81  2  322.49  2  94.42  2  22.01  2  6.90  2  0.06 
49  1.36  3  388.30  3  110.21  3  24.98  3  6.87  3  0.04 
50  1.81  4  172.13  4  34.59  4  3.03  4  0.06  4  3.18 
51  2.55  5  20.72  5  9.16  5  2.03  5  2.11  5  0.13 
65  0.63  1  303.47  1  106.37  1  35.55  1  17.58  1  4.40 
66  1.13  2  520.48  0  0  0  0  0  0  0  0 
67  1.45  3  407.44  3  124.57  3  33.62  3  11.92  3  1.26 
68  1.94  4  157.81  4  40.46  4  8.47  4  1.33  4  0.00 
DSMSE  336.12  88.42  27.01  14.29  2.17  

32894.05i  21501.91i  8321.26i  8756.23  16942.8  

−5.50i  −4.884i  −6.976i  4.498  0.951  
100  40  10  10  
MinPts  5  5  5  5 
DS, deviation square; MinPts, minimum number of points.
However, this method has the possibility of failure. Especially when there are many outliers with the DS values close to each other, they are easy to aggregate into a class, resulting in the deletion of the normal points. In this case, the result has a significant deviation.
When this method is applied to the practical optical fibre data, it is affected by many factors and will produce significant errors. These factors include
Sampling affects the measurement of the time when the vibration reaches the fibre. Sampling occurs at integer multiples of the sampling period. If the time of the vibration signal arriving at the optical fibre is not an integer multiple of the sampling period, the signal arrival time is regarded as the nearest next integer time. For example, the vibration that reaches the fibre is the 5.1 sampling point theoretically, but the actual time collected is in the 6th sampling point.
It is assumed that the arrival time of the vibration to the fibre is uniformly distributed, the average arrival time obtained will be 0.5 sampling points later than the actual arrival time during sampling.
The better method to solve the sampling problem is interpolation which can greatly reduce the sampling error.
Noises will affect the accuracy of TDOA estimation. With the decrease of SNR, the accuracy of TDOA estimation also decreases sharply. Moreover, when the SNR exceeds 10 dB, the accuracy of TDOA estimation obtained is also higher.
In this paper, a fibre segment is equivalent to a sensor. However, a sensor can be thought of as a point, and a fibre segment as a line segment whose length is equal to the spatial resolution. Therefore, the spatial resolution will cause a ±d/2 error in the estimation of the closest point of the vibration source to the optical fibre. Further reduction of spatial resolution is the key to improve the accuracy of distance estimation.
Fibres in the experiment are paved in the earth ditch with a depth of 10–15 cm. In the pavement process, the fibre is not straight due to the factors such as the fluctuation of the pavement surface and the natural bending of the fibre.
The propagation speed of the vibration signal in the soil layer is related to the propagation medium. If the elasticity modulus of the propagation medium increases, the velocity will be faster. The elasticity modulus of stone is greater than that of soil, and the elasticity modulus of soil is greater than that of sand. In the experimental site, when the earth ditch was excavated, it was found that there were more stones and harder soil. Therefore, when the vibration signal propagated in the soil layer, the wave velocity in different media was different. It even happened that the nondatum fibre segment received the signal earlier than the datum fibre segment.
Safety protection for underground communication fibre is of great significance. When a threatening event to optical fibre security happens, the effective distance measurement between the threat location and the optical fibre is an important basis for further protective actions. OTDR technology is used to make use of the existing communication fibre without adding additional sensors. This is a new way to protect the security of optical fibre by using the fibre itself as the sensor.
The method of distance measurement used in this paper includes four steps: fibre marking, TDOA estimation, outliers’ removal, and distance calculation. It is proved that the method is correct and useful in theory.
However, affected by many factors, it is still not a very accurate distance estimation method in practice. To further improve the accuracy of the estimation, certain improvements, such as reducing spatial resolution and increasing SNR, need further study.
Complete process to simu fibre data.
1  8.055  7.9716  −9  248.021  0  0 
2  7.062  7.4851  −8  28.806  0  0 
3  6.071  8.02  −7  286.574  0  0 
4  5.083  5.464  −6  1.424  −6  0 
5  4.099  5.7098  −5  153.47  0  0 
6  3.123  3.211  −4  0.612  −4  1.402 
7  2.162  2.3093  −3  0.191  −3  0.03 
8  1.236  2.1822  −2  22.245  0  0 
9  0.414  0.8831  −1  2.88  −1  1.612 
10  0  0  0  0  0  0 
11  0.414  0.3109  1  0.032  1  1.28E−01 
12  1.236  1.5947  2  2.861  2  1.199 
13  2.162  2.7702  3  10.526  3  7.343 
14  3.123  1.4064  4  124.818  0  0 
15  4.099  4.1073  5  8.236  5  7.932 
16  5.083  6.8807  6  221.649  0  0 
17  6.071  6.5797  7  3.14  7  0.673 
18  7.062  8.4726  8  70.983  0  0 
19  8.055  7.1539  9  711.649  0  0 
DSMSE  105.4509  2.2577  

8508.390  9370.197  

1.419  0.958 
Mean square error and success rate of five methods of Table 1.
MSE  5.943  17.0042  5.2276  5.6606  0.9272 
Success Rate  100%  42.1053%  89.4737%  94.7368%  100% 
The overall process of excavation data treatment.
11  0.30  1  76.82  1  22.90  1  5.46  1  1.83  1  0.02 
12  0.40  2  73.13  2  11.87  2  0.09  2  0.92  2  6.18 
13  0.45  3  20.16  3  1.15  3  20.17  3  34.44  0  0 
14  0.61  4  1.15  4  34.72  4  103.16  0  0  0  0 
15  0.71  5  40.03  5  195.43  0  0  0  0  0  0 
16  0.82  6  246.77  6  577.32  0  0  0  0  0  0 
17  0.81  7  889.91  0  0  0  0  0  0  0  0 
18  0.93  8  1902.98  0  0  0  0  0  0  0  0 
20  0.95  −8  1878.44  0  0  0  0  0  0  0  0 
21  1.94  −7  504.32  0  0  0  0  0  0  0  0 
22  1.61  −6  63.24  −6  268.29  0  0  0  0  0  0 
23  1.56  −5  13.22  −5  24.83  −5  101.21  0  0  0  0 
24  1.05  −4  81.91  −4  0.01  −4  31.30  −4  62.44  0  0 
25  0.28  −3  0.08  −3  16.31  −3  40.42  −3  53.09  0  0 
26  0.32  −2  37.20  −2  3.69  −2  0.42  −2  2.83  −2  8.66 
27  0.25  −1  52.07  −1  14.75  −1  3.11  −1  0.87  −1  0.01 
47  0.53  1  228.79  1  77.38  1  24.15  1  11.17  1  2.14 
48  0.81  2  322.49  2  94.42  2  22.01  2  6.90  2  0.06 
49  1.36  3  388.30  3  110.21  3  24.98  3  6.87  3  0.04 
50  1.81  4  172.13  4  34.59  4  3.03  4  0.06  4  3.18 
51  2.55  5  20.72  5  9.16  5  2.03  5  2.11  5  0.13 
65  0.63  1  303.47  1  106.37  1  35.55  1  17.58  1  4.40 
66  1.13  2  520.48  0  0  0  0  0  0  0  0 
67  1.45  3  407.44  3  124.57  3  33.62  3  11.92  3  1.26 
68  1.94  4  157.81  4  40.46  4  8.47  4  1.33  4  0.00 
DSMSE  336.12  88.42  27.01  14.29  2.17  

32894.05i  21501.91i  8321.26i  8756.23  16942.8  

−5.50i  −4.884i  −6.976i  4.498  0.951  
100  40  10  10  
MinPts  5  5  5  5 
Comparison of TDOA estimation for simu fibre data.
−9  7  −2  7.3147  7.3945  7.9716  8.055 
−8  8  8  8.2805  8.8228  7.4851  7.062 
−7  5  5  4.9436  4.3144  8.02  6.071 
−6  7  −7  6.9048  6.3663  5.464  5.083 
−5  9  −6  6.0667  8.9716  5.7098  4.099 
−4  3  −3  2.7389  2.8608  3.211  3.123 
−3  5  −1  −1.1754  −1.2043  2.3093  2.162 
−2  5  −10  5.2752  5.0414  2.1822  1.236 
−1  1  1  0.7914  0.9828  0.8831  0.414 
0  0  0  0  0  0  0 
1  0  13  −0.016  0.1585  0.3109  0.414 
2  2  3  2.8405  2.7737  1.5947  1.236 
3  3  −14  3.6201  3.1296  2.7702  2.162 
4  1  9  0.9783  1.2096  1.4064  3.123 
5  4  −8  3.9291  4.1022  4.1073  4.099 
6  8  14  7.8445  7.9643  6.8807  5.083 
7  11  −13  11.2104  11.013  6.5797  6.071 
8  3  3  11.0275  3.4893  8.4726  7.062 
9  6  −10  6.4191  6.4771  7.1539  8.055 
Law of interest rate changes in financial markets based on the differential equation model of liquidity Basalt fibre continuous reinforcement composite pavement reinforcement design based on finite element model Industrial transfer and regional economy coordination based on multiple regression model Response model for the psychological education of college students based on nonlinear finite element equations Satisfactory consistency judgement and inconsistency adjustment of linguistic judgement matrix Analysis of the relationship between industrial agglomeration and regional economic growth based on the multiobjective optimisation model Constraint effect of enterprise productivity based on constrained form variational computing The impact of urban expansion in Beijing and Metropolitan Area urban heat Island from 1999 to 2019 Ultrasonic wave promoting ice melt in ice storage tank based on polynomial fitting calculation model Regarding new wave distributions of the nonlinear integropartial Ito differential and fifthorder integrable equations Badminton players’ trajectory under numerical calculation method Innovations to Attribute Reduction of Covering Decision System Based on Conditional Information Entropy Nonlinear Differential Equations in the Teaching Model of Educational Informatisation The evaluation of college students’ innovation and entrepreneurship ability based on nonlinear model Smart Communities to Reduce Earthquake Damage: A Case Study in Xinheyuan, China Institutional investor company social responsibility report and company performance Mathematical analysis of China's birth rate and research on the urgency of deepening the reform of art education Firstprinciples calculations of magnetic and mechanical properties of Febased nanocrystalline alloy Fe_{80}Si_{10}Nb_{6}B_{2}Cu_{2} Has the belt and road initiative boosted the resident consumption in cities along the domestic route? – evidence from credit card consumption Attitude control for the rigid spacecraft with the improved extended state observer Cognitive Computational Model Using Machine Learning Algorithm in Artificial Intelligence Environment Research on tourism income index based on ordinary differential mathematical equation Application of HigherOrder Ordinary Differential Equation Model in Financial Investment Stock Price Forecast Sports Science Teaching of Athletics Based on Nonlinear Mathematical Equation Informatisation of educational reform based on fractional differential equations Research on the control of quantitative economic management variables under the numerical method based on stochastic ordinary differential equations Network monitoring and processing accuracy of big data acquisition based on mathematical model of fractional differential equation System dynamics model of output of ball mill Sensitivity Analysis of the Waterproof Performance of Elastic Rubber Gasket in Shield Tunnel Design of Morlet wavelet neural network to solve the nonlinear influenza disease system Motion about equilibrium points in the JupiterEuropa system with oblateness Badminton players’ trajectory under numerical calculation method Optimal preview repetitive control for impulsefree continuoustime descriptor systems Development of main functional modules for MVB and its application in rail transit Study on the impact of forest fire prevention policy on the health of forest resources Value Creation of Real Estate Company Spinoff Property Service Company Listing Selection by differential mortality rates Digital model creation and image meticulous processing based on variational partial differential equation The modelling and implementation of the virtual 3D animation scene based on the geometric centreofmass algorithm The policy efficiency evaluation of the Beijing–Tianjin–Hebei regional government guidance fund based on the entropy method The transfer of stylised artistic images in eye movement experiments based on fuzzy differential equations Research on behavioural differences in the processing of tenant listing information: An eyemovement experiment A review of the treatment techniques of VOC Some classes of complete permutation polynomials in the form of ( x ^{pm} −x +δ )^{s} +ax ^{pm} +bx overF _{p2m}Deformation and stress theory of surrounding rock of shallow circular tunnel based on complex variable function method The consistency method of linguistic information and other four preference information in group decisionmaking Research on the willingness of Forest Land’s Management Rights transfer under the Beijing Forestry Development A mathematical model of the fractional differential method for structural design dynamics simulation of lower limb force movement step structure based on Sanda movement Numerical calculation and study of differential equations of muscle movement velocity based on martial articulation body ligament tension Study on the maximum value of flight distance based on the fractional differential equation for calculating the best path of shot put Sports intensity and energy consumption based on fractional linear regression equation Translog function in government development of lowcarbon economy Analysis of the properties of matrix rank and the relationship between matrix rank and matrix operations Research on the Psychological Distribution Delay of Artificial Neural Network Based on the Analysis of Differential Equation by Inequality Expansion and Contraction Method Study on Establishment and Improvement Strategy of Aviation Equipment Research on Financial Risk Early Warning of Listed Companies Based on Stochastic Effect Mode The Model of Sugar Metabolism and Exercise Energy Expenditure Based on Fractional Linear Regression Equation Constructing Artistic Surface Modeling Design Based on Nonlinear Overlimit Interpolation Equation Numerical Simulation Analysis Mathematics of Fluid Mechanics for Semiconductor Circuit Breaker Characteristics of Mathematical Statistics Model of Student Emotion in College Physical Education Human Body Movement Coupling Model in Physical Education Class in the Educational Mathematical Equation of Reasonable Exercise Course The contribution of structural equation model analysis to higher education agglomeration and innovation and entrepreneurship Study on the evolutionary game theory of the psychological choice for online purchase of fresh produce under replicator dynamics formula The influence of X fuzzy mathematics method in basketball tactics scoring Mathematical statistics algorithm in the bending performance test of corroded reinforced concrete beams under fatigue load Nonlinear strategic human resource management based on organisational mathematical model Back propagation mathematical model for stock price prediction Evolutionary game research on the psychological choice of online shopping of fresh agricultural products based on dynamic simulation model Differential equation model of financial market stability based on big data Multiattribute decisionmaking methods based on normal random variables in supply chain risk management Linear fractional differential equations in bank resource allocation and financial risk management model Construction and reform of art design teaching mode under the background of the integration of nonlinear equations and the internet Spatial–temporal graph neural network based on node attention A contrastive study on the production of double vowels in Mandarin Financial accounting measurement model based on numerical analysis of rigid normal differential equation and rigid generalised functional equation Research of cascade averaging control in hydraulic equilibrium regulation of heating pipe network Mathematical analysis of civil litigation and empirical research of corporate governance Health monitoring of Bridges based on multifractal theory College students’ innovation and entrepreneurship ability based on nonlinear model Health status diagnosis of the bridges based on multifractal detrend fluctuation analysis Mathematical simulation analysis of optimal testing of shot puter's throwing path Performance evaluation of college laboratories based on fusion of decision tree and BP neural network Application and risk assessment of the energy performance contracting model in energy conservation of public buildings The term structure of economic management rate under the parameter analysis of the estimation model based on common differential equation Sensitivity analysis of design parameters of envelope enclosure performance in the dryhot and drycold areas The Spatial Form of Digital Nonlinear Landscape Architecture Design Based on Computer Big Data The improvement of museum information flow based on paste functional mapping method The art design of industrialised manufacturing furniture products based on the simulation of mathematical curves TOPSIS missile target selection method supported by the posterior probability of target recognition Research on Evaluation of Intercultural Competence of Civil Aviation College Students Based on Language Operator The incentive contract of subject librarians in university library under the nonlinear task importance Modelling and Simulation of Collaborative Innovation System in Colleges and Universities Based on Interpreted Structural Equation Model Small amplitude periodic solution of Hopf Bifurcation Theorem for fractional differential equations of balance point in group competitive martial arts The Optimal Solution of Feature Decomposition Based on the Mathematical Model of Nonlinear Landscape Garden Features Composite mechanical performance of prefabricated concrete based on hysteresis curve equation Higher education innovation and reform model based on hierarchical probit Application of Fuzzy Mathematics Calculation in Quantitative Evaluation of Students’ Performance of Basketball Jump Shot The teaching of sports science of track and fieldbased on nonlinear mathematical equations Visual error correction of continuous aerobics action images based on graph difference function Ecological balance model of effective utilization of agricultural water resources based on fractional differential equations Application of Higher Order Ordinary Differential Equation Model in Financial Investment Stock Price Forecast Precision algorithms in secondorder fractional differential equations Application of Forced Modulation Function Mathematical Model in the Characteristic Research of Reflective Intensity Fibre Sensors Fractional differential equations in National Sports Training in Colleges and Universities Radioactive source search problem and optimisation model based on metaheuristic algorithm Visualized calculation of regional power grid power data based on multiple linear regression equation Application of mathematical probabilistic statistical model of base – FFCA financial data processing Leastsquares method and deep learning in the identification and analysis of nameplates of power equipment Research on a method of completeness index based on complex model Distribution network monitoring and management system based on intelligent recognition and judgement Fake online review recognition algorithm and optimisation research based on deep learning Research on the sustainable development and renewal of Macao inner harbour under the background of digitisation Support design of main retracement passage in fully mechanised coal mining face based on numerical simulation Study on the crushing mechanism and parameters of the twoflow crusher Topological optimisation technology of gravity dam section structure based on ANSYS partial differential equation operation Interaction design of financial insurance products under the Era of AIoT Modeling the pathway of breast cancer in the Middle East Corporate social responsibility fulfilment, productmarket competition and debt risk: Evidence from China ARMA analysis of the green innovation technology of core enterprises under the ecosystem – Time series data Reconstruction of multimodal aesthetic critical discourse analysis framework Image design and interaction technology based on Fourier inverse transform What does students’ experience of eportfolios suggest Research on China interregional industrial transformation slowdown and influencing factors of industrial transformation based on numerical simulation The medical health venture capital network community structure, information dissemination and the cognitive proximity The optimal model of employment and entrepreneurship models in colleges and universities based on probability theory and statistics A generative design method of building layout generated by path Analysis of the causes of the influence of the industrial economy on the social economy based on multiple linear regression equation Research of neural network for weld penetration control Analysing the action techniques of basketball players’ shooting training using calculus method Engineering project management based on multiple regression equation and building information modelling technology Research on predictive control of students’ performance in PE classes based on the mathematical model of multiple linear regression equation Beam control method for multiarray antennas based on improved genetic algorithm The influence of X fuzzy mathematical method on basketball tactics scoring Mathematical model of back propagation for stock price forecasting Application of regression function model based on panel data in bank resource allocation financial risk management Application of Logical Regression Function Model in Credit Business of Commercial Banks Research on aerobics training posture motion capture based on mathematical similarity matching statistical analysis Application of SobolevVolterra projection and finite element numerical analysis of integral differential equations in modern art design Research on motion capture of dance training pose based on statistical analysis of mathematical similarity matching Application of Btheory for numerical method of functional differential equations in the analysis of fair value in financial accounting Research on the influence of fuzzy mathematics simulation model in the development of Wushu market Study on audiovisual family restoration of children with mental disorders based on the mathematical model of fuzzy comprehensive evaluation of differential equation Differenceindifferences test for micro effect of technological finance cooperation pilot in China Application of multiattribute decisionmaking methods based on normal random variables in supply chain risk management Exploration on the collaborative relationship between government, industry, and university from the perspective of collaborative innovation The impact of financial repression on manufacturing upgrade based on fractional Fourier transform and probability AtanKA New SVM Kernel for Classification Validity and reliability analysis of the Chinese version of planned happenstance career inventory based on mathematical statistics Visual positioning system for marine industrial robot assembly based on complex variable function Application of LaneEmden differential equation numerical method in fair value analysis of financial accounting Regression function model in risk management of bank resource allocation Application of numerical method of functional differential equations in fair value of financial accounting Optimal solution of fractional differential equations in solving the relief of college students’ mental obstacles Risk contagion in financial markets based on copula model Calculating university education model based on finite element fractional differential equations and macrocontrol analysis Automatic parameter selection ZVD shaping algorithm for crane vibration suppression based on particle swarm optimisation Educational research on mathematics differential equation to simulate the model of children's mental health prevention and control system Analysis of enterprise management technology and innovation based on multilinear regression model Verifying the validity of the whole person model of mental health education activities in colleges based on differential equation RETRACTION NOTE Calculation of tourism development income index based on finite element ordinary differential mathematical equation Adoption of deep learning Markov model combined with copula function in portfolio risk measurement Radar system simulation and nonGaussian mathematical model under virtual reality technology Comparison of compression estimations under the penalty functions of different violent crimes on campus through deep learning and linear spatial autoregressive models Research and application of constructing football training linear programming based on multiple linear regression equation Research on management evaluation of enterprise sales cash flow percentage method based on the application of quadratic linear regression equations Mathematical simulation analysis of optimal detection of shotputters’ best path Determination of the minimum distance between vibration source and fibre under existing optical vibration signals: a study Mathematical modelling of enterprise financial risk assessment based on risk conduction model Nonlinear differential equations based on the BSM model in the pricing of derivatives in financial markets Mathematical simulation experiment based on optimisation of heat treatment process of aluminium alloy materials Mathematical model of transforming image elements to structured data based on BP neural network Educational reform informatisation based on fractional differential equation 3D Mathematical Modelling Technology in Visual Rehearsal System of Sports Dance MCM of Student’s Physical Health Based on Mathematical Cone Sports health quantification method and system implementation based on multiple thermal physiology simulation Research on visual optimization design of machine–machine interface for mechanical industrial equipment based on nonlinear partial equations Informationisation of teaching model for track and field education based on finite element higherorder fractional differential equation Information technology of preschool education reform of fine arts based on fractional differential equation Information Teaching Model of Preschool Art Education in Colleges and Universities Based on Finite Element HigherOrder Fractional Differential Equation Application of artificial intelligence algorithm in mathematical modelling and solving College Students’ Mental Health Climbing Consumption Model Based on Nonlinear Differential Equations Communication architecture of power monitoring system based on incidence matrix model Differential equation to verify the validity of the model of the wholeperson mental health education activity in Universities Optimisation of Modelling of Finite Element Differential Equations with Modern Art Design Theory Analysis and Prediction of College Students’ Mental Health Based on Kmeans Clustering Algorithm Mathematical function data model analysis and synthesis system based on shortterm human movement Human gait modelling and tracking based on motion functionalisation Analysis and synthesis of function data of human movement Energysaving technology of BIM green buildings using fractional differential equation Study on the training model of football movement trajectory drop point based on fractional differential equation Financial Accounting Measurement Model Based on Numerical Analysis of Rigid Normal Differential Equation and Rigid Functional Equation User online consumption behaviour based on fractional differential equation Differential equation model of financial market stability based on Internet big data Multiattribute Decision Method Based on Normal Random Variable in Economic Management Risk Control Children’s cognitive function and mental health based on finite element nonlinear mathematical model Dichotomy model based on the finite element differential equation in the educational informatisation teaching reform model Nonlinear Dissipative System Mathematical Equations in the Multiregression Model of Informationbased Teaching Stock price analysis based on the research of multiple linear regression macroeconomic variables Fractional Linear Regression Equation in Agricultural Disaster Assessment Model Based on Geographic Information System Analysis Technology