Research on the Prediction Method of Clock Tester Calibration Data Based on Radial Basis Function Neural Network

: A radial basis function (RBF) neural network-based calibration data prediction model for clock testers is proposed to address the issues of ﬁxed calibration cycles, low efﬁciency, and waste of electrical energy. This provides a new method for clock tester traceability calibration. First, analyze the mechanism of clock tester calibration parameters and the inﬂuencing factors of prediction targets. Based on the learning rules of an RBF neural network, determine the data types of training and testing sets. Second, normalize the training and testing data to avoid the adverse effects of data characteristics and distribution differences on the prediction model. Finally, based on different prediction objectives, time-driven and data-driven calibration data prediction models are constructed using RBF neural networks. Through simulation analysis, it is shown that an RBF neural network is superior to a BP neural network in predicting clock tester calibration data, and time-driven prediction accuracy is superior to data-driven prediction accuracy. Moreover, the prediction error and mean square error of both prediction models are on the order of 10 − 9 , meeting the prediction accuracy requirements.


Introduction
A clock tester is an instrument that uses an internal quartz crystal oscillator as a time reference and can display daily errors, second errors, and frequency errors.It is used to measure frequency accuracy and second pulse signals.It is composed of a waveform shaping circuit, a frequency division circuit, an electronic control circuit, a display circuit, etc.It is mainly used for measuring the second pulse signal of electric energy meters [1].The general traceability calibration interval for clock testers is one year, and the acquisition time of calibration indicators (stability and accuracy) is long (at least 8 h) [2], which poses challenges to customers and inspectors and wastes electrical energy.
At present, there are no more efficient physical measurement methods and calibration standards, and research on predictive models for clock tester calibration data is still in its infancy.A radial basis function neural network is capable of approximating any nonlinear function and can handle difficult-to-analyze laws within the system.It has good generalization ability and fast learning convergence speed.It has been successfully applied in nonlinear function approximation, time series analysis, data classification, pattern recognition, information processing, image processing, system modeling, control, and fault diagnosis [3][4][5][6][7].The prediction method based on radial basis function neural networks has been widely studied and applied in other fields [8][9][10][11][12].Lashkenari et al. [13] developed a radial basis function neural network prediction model for predicting the viscosity of Iranian crude oil that has universality and accuracy.Luo et al. [14] proposed a fast method for predicting the fatigue life of automotive wheels based on radial basis function neural networks combined with orthogonal decomposition, which has ideal accuracy.Mohammad et al. [15] established a concrete mix ratio model based on radial basis functions to estimate the compressive strength of concrete containing different amounts of fly ash at Electronics 2023, 12, 4677 2 of 14 any time in response to the fact that fly ash can enhance the mechanical properties and durability of concrete materials.Hao et al. [16] proposed a prediction model for chaotic radial basis function neural networks based on radial basis function neural networks and chaos theory, which utilizes the Earth's natural pulse electromagnetic field signals for potential strength trend prediction.Stephen et al. [17] introduced the radial basis function neural network algorithm into the path loss prediction model to solve the problems of large errors and low universality in empirical and deterministic models.Zhang et al. [18] developed an intelligent model based on radial basis function neural networks to predict the thermal conductivity of nanofluids under various conditions in response to the inherent complexity of nanofluids.This is of great significance for promoting the industrial application of nanofluids.Tao et al. [19] proposed an improved radial basis function prediction model based on a differential evolution algorithm for the complexity of the coking energy consumption process, which is of great significance for reducing energy consumption, saving costs, and improving enterprises' economic benefits.Yoon et al. [20] proposed a probabilistic motion prediction algorithm based on radial basis function neural networks for the uncertainty of vehicle motion trajectories.This algorithm can accurately calculate the likelihood of multiple target lanes and surrounding vehicle trajectories using radial basis function neural networks, which is crucial for avoiding potential risks in the future.The successful application of the prediction method based on radial basis function neural networks in other fields has verified the theoretical feasibility of the clock tester calibration data prediction model.In addition, GMDH neural networks, BP neural networks, etc., can also be used for building prediction models.Roshani et al. [21] used the GMDH neural network to build a prediction model for estimating the gas volume percentage in two-phase flow without relying on the flow pattern and improving the measurement accuracy by 2.7 times.Qin et al. [22] built a BP neural network model to predict user behavior that has high prediction efficiency and accuracy.The calibration data of the clock tester studied in this article has typical nonlinear characteristics and randomness within a certain numerical range.Therefore, an radial basis function neural network with a strong generalization ability and a good nonlinear function approximation effect is used for research.
In view of this, a radial basis function neural network is used to build a clock tester calibration data prediction model to replace traditional calibration methods.By analyzing the calibration parameters of the clock tester, time-and data-driven training functions are established to expand the prediction range of the model and improve its applicability.The clock tester calibration data prediction model provides a new method to improve the efficiency of the clock tester traceability calibration, which can effectively reduce the number of traditional calibration times, not only saving electricity but also providing theoretical support for subsequent calibration time decisions and creating conditions for condition-based traceability.

Principles of Radial Basis Function Neural Networks
Radial basis function neural networks, abbreviated as RBF neural networks, is a neural network structure proposed by J. Moody and C. Darken in the late 1980s.An RBF neural network is a high-performance feedforward network with the ability to achieve optimal approximation and overcome local minimum problems [7].
According to the different number of radial basis functions, RBF neural networks can be divided into regularized RBF neural networks and generalized RBF neural networks.When using a regularized RBF neural network structure, the number of hidden layer nodes is the number of samples, and the data center of the basis function is the sample itself.Parameter design only needs to consider the extension constant and the weight of the output nodes.The network topology is shown in Figure 1 [5].where X is the input data vector and Cn is the center of the data sample.The third layer is the output layer, where d1, d2, …, dn are the output data, and their activation function is a linear function.
When using a generalized RBF neural network structure, the learning algorithm of the RBF neural network should address issues such as how to determine the number of hidden layer nodes in the network, how to determine the data centers and extension constants of each radial basis function, and how to correct the output weights.The topology structure of the RBF neural network is shown in Figure 2 [23].

Input layer
Hidden layer Output layer Topological structure of generalized RBF neural networks.
In Figure 2, the number of hidden layer nodes in the RBF neural network is smaller than the number of samples.This is because when the number of samples is large, the computational complexity of the network will be large, and the larger the matrix, the greater the likelihood of its ill condition.Therefore, the generalized RBF neural network reduces the likelihood of a matrix ill condition by reducing the number of hidden layer nodes.Figure 2 φ0 is the offset value.
The paper builds a clock tester calibration data prediction model based on the RBF neural network, which has a small amount of data.Therefore, a regularized RBF neural network is used.The Gaussian function is chosen as the radial basis function, as shown in Formula (1).wherein r is the distance from the data to the center of the sample, δ the width of the center of the basis function.

Problem Description
Using a clock tester as the measurement object, analyze the measured data and construct a prediction model.
The calibration of the clock tester mainly involves measuring its stability and accuracy, and the main calibration methods [2] are as follows: The regularized RBF neural network is usually a three-layer feedforward network.The first layer is the input layer, where x 1 , x 2 , . .., x n are the input data.The second layer is the hidden layer, where the node basis function ϕ(X, C n ) is the radial basis function, where X is the input data vector and C n is the center of the data sample.The third layer is the output layer, where d 1 , d 2 , . .., d n are the output data, and their activation function is a linear function.
When using a generalized RBF neural network structure, the learning algorithm of the RBF neural network should address issues such as how to determine the number of hidden layer nodes in the network, how to determine the data centers and extension constants of each radial basis function, and how to correct the output weights.The topology structure of the RBF neural network is shown in Figure 2 [23].
itself.Parameter design only needs to consider the extension constant and the weight of the output nodes.The network topology is shown in Figure 1 [5].

Input layer
Hidden layer Output layer Topological structure of regularized RBF neural networks.
The regularized RBF neural network is usually a three-layer feedforward network.The first layer is the input layer, where x1, x2, …, xn are the input data.The second layer is the hidden layer, where the node basis function ( , ) where X is the input data vector and Cn is the center of the data sample.The third layer is the output layer, where d1, d2, …, dn are the output data, and their activation function is a linear function.
When using a generalized RBF neural network structure, the learning algorithm of the RBF neural network should address issues such as how to determine the number of hidden layer nodes in the network, how to determine the data centers and extension constants of each radial basis function, and how to correct the output weights.The topology structure of the RBF neural network is shown in Figure 2 [23].

Input layer
Hidden layer Output layer . Topological structure of generalized RBF neural networks.
In Figure 2, the number of hidden layer nodes in the RBF neural network is smaller than the number of samples.This is because when the number of samples is large, the computational complexity of the network will be large, and the larger the matrix, the greater the likelihood of its ill condition.Therefore, the generalized RBF neural network reduces the likelihood of a matrix ill condition by reducing the number of hidden layer nodes.Figure 2 φ0 is the offset value.
The paper builds a clock tester calibration data prediction model based on the RBF neural network, which has a small amount of data.Therefore, a regularized RBF neural network is used.The Gaussian function is chosen as the radial basis function, as shown in Formula (1).
wherein r is the distance from the data to the center of the sample, δ the width of the center of the basis function.

Problem Description
Using a clock tester as the measurement object, analyze the measured data and construct a prediction model.
The calibration of the clock tester mainly involves measuring its stability and accuracy, and the main calibration methods [2] are as follows: In Figure 2, the number of hidden layer nodes in the RBF neural network is smaller than the number of samples.This is because when the number of samples is large, the computational complexity of the network will be large, and the larger the matrix, the greater the likelihood of its ill condition.Therefore, the generalized RBF neural network reduces the likelihood of a matrix ill condition by reducing the number of hidden layer nodes.Figure 2 ϕ 0 is the offset value.
The paper builds a clock tester calibration data prediction model based on the RBF neural network, which has a small amount of data.Therefore, a regularized RBF neural network is used.The Gaussian function is chosen as the radial basis function, as shown in Formula (1).
wherein r is the distance from the data to the center of the sample, δ the width of the center of the basis function.

Problem Description
Using a clock tester as the measurement object, analyze the measured data and construct a prediction model.
The calibration of the clock tester mainly involves measuring its stability and accuracy, and the main calibration methods [2] are as follows: Calibration basis: JJF1662-2017 "Clock Tester" calibration specification.Calibration instrument: GPS-locked rubidium atomic frequency standard, universal counter.
Stability calibration method: After starting the machine, measure eight times using 500 kHz as the standard, with a measurement interval of 1 h.We randomly sample three times during each measurement, take the arithmetic mean of them as the measurement value, and calculate the stability according to Formula (2).
(2) wherein S represents the calculated stability, x i represents the measured data, and x i max and x i min are divided into the maximum and minimum values in the measured data.Accuracy calibration method: After starting the machine, randomly measure 10 times, take three samples of data each time, take their arithmetic mean as the measurement value, and calculate the accuracy according to Formula (3).
wherein A is the calculated accuracy, and x 0 is the reference standard.Taking a clock tester as an example, collect data measured by the clock tester.Build a prediction model based on the RBF neural network algorithm, train the collected data, and predict the accuracy and stability of the clock tester after meeting the requirements.

Factors Affecting Calibration Parameters
According to Section 2.2, the calibration parameters of the clock tester are accuracy and stability.Through analysis of the calibration process, it can be seen that the measurement parameters are mainly affected by the calibration instrument.Therefore, in the measurement process of the clock tester, it is necessary to consider the calibration uncertainty [2].According to Formulas (2) and (3), the uncertainty of measurement results mainly comes from measurement repeatability and accuracy of reference standards.Therefore, the measurement data are shown in Formula (4), wherein x s , x 0 , ∆x represent the general counter display data (i.e., measurement data), output data standard values, and data measurement errors.
(1) Measurement Repeatability Use a universal counter to repeatedly sample the output data of the clock tester n times.According to the Bessel formula, calculate the measurement uncertainty component as shown in Formula (5).
wherein u 1 is the calculated uncertainty component, x i is the collected data, and x is the average value of the collected data.
(2) Accuracy of Reference Standards According to uniform distribution, if the reference standard accuracy is c b , then k = √ 3, and the uncertainty component is calculated as shown in Formula (6).
Electronics 2023, 12, 4677 5 of 14 (3) Uncertainty Calculation Due to the independence of the introduced measurement uncertainty components, the combined measurement uncertainty is shown in Formula (7).
The expanded uncertainty of measurement is shown in Formula ( 8), wherein k is the expansion coefficient, generally taken as k = 2.In the actual calibration work, if we estimate that the measured value follows the normal distribution, when k = 2, U is the half-width of the confidence interval with a confidence probability of approximately 95%, representing that 95% of the measured value is within the range of ±U of the actual value.

Data Analysis and Preprocessing
The data in the paper came from the measurement results of five clock testers of the same model produced by the same manufacturer and year by a certain unit for 4 consecutive years.Part of the data is shown in Table 1.Taking the clock tester as the research object, the order of error is small.Therefore, it is necessary to preprocess the measured data, remove common features, and amplify the error, as shown in Formula (9).
wherein x dea is the processed data, c d is the error amplification coefficient, x mea is the measurement data, and x com is the common feature data.
In order to avoid the adverse effects of data features and distribution differences on the model, it is necessary to normalize and denormalize the training samples.The normalization process is shown in Formula (10), and the denormalization process is shown in Formula (11).
x norm = x dea − x min x max − x min (10) x dea = x norm (x max − x min ) + x min (11) wherein x norm is the normalized data, Electronics 2023, 12, 4677 6 of 14

Model Construction
From Section 2.1, it can be seen that the RBF neural network has a multi-layer forward structure with one hidden layer.According to the different number of hidden layer nodes, it can be divided into regularized RBF neural networks and generalized RBF neural networks.The paper selects the regularized RBF neural network structure.
Let the N-dimensional space have P data points X p , p = 1, 2, . . .P, and the corresponding target values in the output space are d p , p = 1, 2, . . .P.
Find a nonlinear mapping function F(X) that satisfies the condition of Formula (12).
Select P basis functions, each corresponding to one data, and the form of each basis function is shown in Formula (13), wherein the basis function ϕ is a nonlinear function, and X p is ϕ the center of the data, and the independent variable is the distance X − Xp between the data X and the center X p .F(X) based on radial basis functions is defined as a linear combination of radial basis functions, as shown in Formula (14).
wherein w p is the corresponding weight value.In the regularized RBF neural network, X p is the sample itself, and X is known.ϕ is the Gaussian function; therefore, as long as we find w p , we can obtain F(X).
Bring the data into F(X) to obtain a linear system of Equation (15).
In the formula, ϕ is an interpolation matrix, and the basis function chosen in the paper is the Gaussian function.Therefore, ϕ is reversible.
Then, the nonlinear mapping function F(X) is solved.The prediction method model is shown in Figure 3.
Write Formula (17) in vector form as shown in Formula (18).

W d ϕ = (18)
In the formula, φ is an interpolation matrix, and the basis function chosen in the paper is the Gaussian function.Therefore, φ is reversible.
Then, the nonlinear mapping function F(X) is solved.The prediction method model is shown in Figure 3.The specific steps of the calibration data prediction method are as follows: Step 1: Obtain the measurement data of five clock testers of the same model produced by the same manufacturer and year for a certain unit for 4 consecutive years.
Step 2: Classify the obtained data according to the measurement method, preprocess the data according to the calibration specifications of the clock tester, and obtain stabilityand accuracy-related data.
Step 3: Decommon the data related to measurement stability and accuracy and amplify the error to obtain the feature data.
Step 4: Normalize the feature data.
Step 5: Build time-driven prediction models and data-driven prediction models and classify feature data.Step 6: Train the two models based on the prediction target to obtain corresponding prediction models and achieve calibration data prediction.
Step 7: Verify the effectiveness of the RBF neural network prediction method using test samples.

Simulation Analysis
To verify the effectiveness of the proposed theory, an RBF neural network model was constructed using MATLAB (https://www.mathworks.com).Modeling can be divided into time-driven models and data-driven models based on the different prediction objectives and input data.The input data of the time-driven model is time and collection frequency, and the output is corresponding measurement data, which can predict the stability and accuracy of the measurement data of the same type of device in the training sample within the year.The data-driven model inputs measurement data for a certain year and outputs measurement data for the corresponding next year, which can predict the stability and accuracy of the measurement data device for the next year.
Based on the physical characteristics of the calibration device, the system calibration uncertainty U rel = 2 × 10 −8 can be obtained from Formulas ( 5)-(8).

Time-Driven Model
Using a time-driven model, the data from devices 20163320, 20163340, 20163366, and 20163355 were used as training samples, and the data from device 20163352 were used as test samples for simulation.The simulation results were compared with the BP neural network simulation results, as shown in Figures 4 and 5 and Tables 2 and 3.
Using a time-driven model and using data from four devices as training samples, the stability fitting curves for the years 2019 to 2022 were obtained, as shown in Figure 4a. Figure 4b shows the comparison between the actual measurement results and predicted results of the fifth device from 2019 to 2022 under stability measurement conditions.The stability of the corresponding year calculated is shown in Table 2.The simulation results show that the predicted results are basically consistent with the measured values.Among them, the RBF neural network has a better-fitting performance compared to the BP neural network, and the predicted stability is closer to the actual value.
Step3: Decommon the data related to measurement stability and accuracy and amplify the error to obtain the feature data.Step4: Normalize the feature data.
Step5: Build time-driven prediction models and data-driven prediction models and classify feature data.Step6: Train the two models based on the prediction target to obtain corresponding prediction models and achieve calibration data prediction.Step7: Verify the effectiveness of the RBF neural network prediction method using test samples.

Simulation Analysis
To verify the effectiveness of the proposed theory, an RBF neural network model was constructed using MATLAB (https://www.mathworks.com).Modeling can be divided into time-driven models and data-driven models based on the different prediction objectives and input data.The input data of the time-driven model is time and collection frequency, and the output is corresponding measurement data, which can predict the stability and accuracy of the measurement data of the same type of device in the training sample within the year.The data-driven model inputs measurement data for a certain year and outputs measurement data for the corresponding next year, which can predict the stability and accuracy of the measurement data device for the next year.
Based on the physical characteristics of the calibration device, the system calibration uncertainty Urel = 2 × 10 −8 can be obtained from Formulas ( 5)-( 8).

Time-Driven Model
Using a time-driven model, the data from devices 20163320, 20163340, 20163366, and 20163355 were used as training samples, and the data from device 20163352 were used as test samples for simulation.The simulation results were compared with the BP neural network simulation results, as shown in Figures 4 and 5 and Tables 2 and 3.        Using a time-driven model and using data from four devices as training samples, the accuracy fitting curves for the years 2019 to 2022 were obtained, as shown in Figure 5a. Figure 5b shows the comparison between the actual measurement results and predicted results of the fifth device from 2019 to 2022 under accuracy measurement conditions.The accuracy of the corresponding year calculated is shown in Table 3.The simulation results show that the predicted results are basically consistent with the measured values.Among them, the RBF neural network has a better-fitting performance compared to the BP neural network, and the predicted accuracy is closer to the actual value.

Data-Driven Model
Using a data-driven model, data from five devices from 2019 to 2021 were used as training samples, and data from five devices from 2022 were used as test samples for training and predictive simulation.The results are shown in Figures 6 and 7 and Tables 4 and 5.
Figure 6a shows the fitting curve of the stability training sample after using a datadriven model, where the abscissa s11-1 represents the measurement data of the first device in the first hour of the first year (2019), and other abscissa annotations are similar.Figure 6b shows the comparison between predicted data and measured data, where the abscissa 1-1 represents the measurement data of the first device in the first hour of the fourth year (2022), and other abscissa annotations are similar.Calculate the corresponding stability of the measured and predicted data according to Formula (2), as shown in Table 4.The simulation results show that the predicted stability is basically consistent with the measured stability, with the RBF neural network having better fitting performance compared to the BP neural network, and the predicted stability is closer to the actual value.
The fitting curve of the training samples for accuracy after using a data-driven model is shown in Figure 7a, and the comparison between predicted data and measured data is shown in Figure 7b.Calculate the data according to Formula (3) to obtain the corresponding accuracy, as shown in Table 5.The simulation results show that the prediction accuracy is basically consistent with the actual measurement accuracy, with RBF neural network having a better-fitting performance compared to BP neural network, and the prediction accuracy is closer to the actual value.Using a time-driven model and using data from four devices as training samples, the stability fitting curves for the years 2019 to 2022 were obtained, as shown in Figure 4a. Figure 4b shows the comparison between the actual measurement results and predicted results of the fifth device from 2019 to 2022 under stability measurement conditions.The stability of the corresponding year calculated is shown in Table 2.The simulation results show that the predicted results are basically consistent with the measured values.Among them, the RBF neural network has a better-fitting performance compared to the BP neural network, and the predicted stability is closer to the actual value.
Using a time-driven model and using data from four devices as training samples, the accuracy fitting curves for the years 2019 to 2022 were obtained, as shown in Figure 5a. Figure 5b shows the comparison between the actual measurement results and predicted results of the fifth device from 2019 to 2022 under accuracy measurement conditions.The accuracy of the corresponding year calculated is shown in Table 3.The simulation results show that the predicted results are basically consistent with the measured values.Among them, the RBF neural network has a better-fitting performance compared to the BP neural network, and the predicted accuracy is closer to the actual value.

Data-Driven Model
Using a data-driven model, data from five devices from 2019 to 2021 were used as training samples, and data from five devices from 2022 were used as test samples for training and predictive simulation.The results are shown in Figures 6 and 7 and Tables 4 and 5.The fitting curve of the training samples for accuracy after using a data-driven model is shown in Figure 7a, and the comparison between predicted data and measured data is shown in Figure 7b.Calculate the data according to Formula (3) to obtain the corresponding accuracy, as shown in Table 5.The simulation results show that the prediction accuracy is basically consistent with the actual measurement accuracy, with RBF neural network having a better-fitting performance compared to BP neural network, and the prediction accuracy is closer to the actual value.

Error Analysis
Based on the predicted and measured stability/accuracy values, calculate the prediction error to obtain the stability and accuracy errors under time-driven and data-driven conditions, as shown in Figures 8 and 9.

Error Analysis
Based on the predicted and measured stability/accuracy values, calculate the prediction error to obtain the stability and accuracy errors under time-driven and data-driven conditions, as shown in Figures 8 and 9.

Error Analysis
Based on the predicted and measured stability/accuracy values, calculate the predic tion error to obtain the stability and accuracy errors under time-driven and data-driven conditions, as shown in Figures 8 and 9.  From Figures 8 and 9, it can be seen that both the time-driven model and the data driven model have predicted data errors of stability and accuracy on the order of 10 −9 .Th mean squared deviations of stability and accuracy prediction errors for RBF neural net works driven by time are 1.077 × 10 −9 and 1.729 × 10 −10 , respectively.The mean squared deviations of stability and accuracy prediction errors for BP neural networks driven by time are 2.269 × 10 −9 and 2.666 × 10 −9 , respectively.The mean squared deviations of stability and accuracy prediction errors for RBF neural networks driven by data are 0.650 × 10 −1 and 2.077 × 10 −9 , respectively.The mean squared deviations of stability and accuracy Error value of accuracy predicted using RBF neural network Error value of stability prediction using RBF neural network Error value of accuracy predicted using BP neural network Error value of stability prediction using BP neural network

Error Analysis
Based on the predicted and measured stability/accuracy values, calculate the predic tion error to obtain the stability and accuracy errors under time-driven and data-driven conditions, as shown in Figures 8 and 9.  From Figures 8 and 9, it can be seen that both the time-driven model and the data driven model have predicted data errors of stability and accuracy on the order of 10 −9 .Th mean squared deviations of stability and accuracy prediction errors for RBF neural net works driven by time are 1.077 × 10 −9 and 1.729 × 10 −10 , respectively.The mean squared deviations of stability and accuracy prediction errors for BP neural networks driven by time are 2.269 × 10 −9 and 2.666 × 10 −9 , respectively.The mean squared deviations of stability and accuracy prediction errors for RBF neural networks driven by data are 0.650 × 10 −1 and 2.077 × 10 −9 , respectively.The mean squared deviations of stability and accuracy Error value of accuracy predicted using RBF neural network Error value of stability prediction using RBF neural network Error value of accuracy predicted using BP neural network Error value of stability prediction using BP neural network From Figures 8 and 9, it can be seen that both the time-driven model and the datadriven model have predicted data errors of stability and accuracy on the order of 10 −9 .The mean squared deviations of stability and accuracy prediction errors for RBF neural networks driven by time are 1.077 × 10 −9 and 1.729 × 10 −10 , respectively.The mean squared deviations of stability and accuracy prediction errors for BP neural networks driven by time are 2.269 × 10 −9 and 2.666 × 10 −9 , respectively.The mean squared deviations of stability and accuracy prediction errors for RBF neural networks driven by data are 0.650 × 10 −10 and 2.077 × 10 −9 , respectively.The mean squared deviations of stability and accuracy prediction errors for BP neural networks driven by data are 5.374 × 10 −9 and 2.112 × 10 −9 , respectively.From the data, it can be seen that the mean square deviation of stability and accuracy prediction errors for RBF and BP neural networks is much smaller than the uncertainty of system calibration.The predicted results meet the requirements.The prediction error and mean square error of an RBF neural network are both smaller than those of a BP neural network.Therefore, in the prediction method of clock tester calibration data, an RBF neural network is superior to a BP neural network.The mean square error of the time-driven model prediction data is smaller than that of the data-driven model prediction data.Therefore, when predicting the stability and accuracy of equipment, when there are measurement data of other devices of the same model in the prediction year, a time-driven model can be selected.When there are no predicted annual data, a data-driven model can be selected.

Conclusions
In this study, we aimed to build a calibration data prediction model for clock testers using RBF neural networks and explore new ways to improve clock tester traceability calibration.We analyzed the mechanism of calibration parameter generation, determined the type of model training data, and, through normalization processing, avoided the impact of training data differences on the model, thereby improving the model's prediction accuracy.Based on the different prediction objectives and training sets, we established timedriven prediction models and data-driven prediction models to predict the stability and accuracy of the same type of device in the measurement data year, as well as the stability and accuracy of the measurement data equipment in the next year.Through simulation analysis, the prediction model built using an RBF neural network has more advantages in data fitting and calibration data prediction compared to the prediction model built using a BP neural network.The prediction accuracy of time-driven prediction models is better than that of data-driven prediction models.The prediction errors of both time-driven and data-driven prediction models are in the order of 10 −9 , which is much smaller than the uncertainty of system calibration and meets the prediction accuracy requirements.When the clock tester meets the prediction requirements and the equipment usage is standardized during the calibration year, the research method in the paper can be considered to replace traditional calibration methods to improve calibration efficiency and lay the foundation for the next step of tracing according to the situation.

Figure 3 .
Figure 3. Predictive model for calibration data.Figure 3. Predictive model for calibration data.

Figure 3 .
Figure 3. Predictive model for calibration data.Figure 3. Predictive model for calibration data.

Figure 4 .
Figure 4. Fitting curve and predictive data for stability of time-driven models.(a) Stability fitting curve; (b) Stability prediction results.

Figure 5 .
Figure 5. Fitting curve and predictive data for accuracy of time-driven model.(a) Accuracy fitting curve; (b) Accuracy prediction results.

Figure 6 .
Figure 6.Fitting curve and predictive data for stability of data-driven models.(a) Stability fitting curve; (b) Stability prediction results.

Figure 7 .
Figure 7. Fitting curve and predictive data for accuracy of data-driven model.(a) Accuracy fitting curve; (b) Accuracy prediction results.

Figure 8 .
Figure 8. Stability and accuracy prediction errors using time-driven models.

Figure 9 .
Figure 9. Stability and accuracy prediction errors using data-driven models.

Figure 7 .
Figure 7. Fitting curve and predictive data for accuracy of data-driven model.(a) Accuracy fitting curve; (b) Accuracy prediction results.

Figure 8 .
Figure 8. Stability and accuracy prediction errors using time-driven models.

Figure 9 .
Figure 9. Stability and accuracy prediction errors using data-driven models.

Figure 8 .
Figure 8. Stability and accuracy prediction errors using time-driven models.

Figure 8 .
Figure 8. Stability and accuracy prediction errors using time-driven models.

Figure 9 .
Figure 9. Stability and accuracy prediction errors using data-driven models.

Figure 9 .
Figure 9. Stability and accuracy prediction errors using data-driven models.

Table 2 .
Comparison between the actual value and the predicted value of stability.

Table 2 .
Comparison between the actual value and the predicted value of stability.

Table 3 .
Comparison between the actual value and the predicted value of accuracy.

Table 3 .
Comparison between the actual value and the predicted value of accuracy.

Table 4 .
Comparison between the actual value and the predicted value of stability (2022).

Table 5 .
Comparison between the actual value and the predicted value of accuracy (2022).