CN114660521A - Automatic calibration device for electric power experimental instrument - Google Patents

Abstract

The application provides electric power laboratory glassware automatic calibration device includes: the device comprises a test configuration module, an auxiliary module, a frequency synchronization module and an autonomous operation module, wherein the modules are matched with each other to detect the frequency updating characteristic of an external reference clock so as to diagnose whether the external reference clock is synchronous with a standard clock source or not, if so, the external reference clock source is determined to be accurate and reliable in frequency, and the initial configuration of a local crystal oscillator of a signal source is synchronized to the external reference clock, so that the frequency adaptive calibration is completed, and the coverage rate and the accuracy of the automatic calibration of the device are improved.

Description

Automatic calibration device for electric power experimental instrument

Technical Field

The application relates to the technical field of communication calibration, in particular to automatic calibration device for electric power experimental instruments.

Background

The statements in this section merely provide background information related to the present application and may not constitute prior art.

In the development process of mobile communication technology, in order to meet the requirement of interoperability test, auxiliary devices such as signal sources and the like are developed and used for supporting the test of protocol consistency of communication equipment.

In The communication consistency test process, frequency synchronization is a precondition, and in practical application, The interoperation test between a signal source and communication equipment comprises two schemes of cable connection and OTA (over The air) connection test, wherein The cable connection scheme can provide a reference clock by The communication equipment, then The reference clock is connected to The signal source through a cable, and The signal source locks a local crystal oscillator to The frequency of The reference clock, so that The frequency synchronization between The signal source and The communication equipment is realized, and The effective development of The subsequent interoperation test is ensured; however, for the OTA test, because the signal source and the communication device are in different directions, there are reasons such as long distance and inconvenient wiring, and it is difficult to achieve frequency synchronization between the signal source and the communication device through the cable in many cases, so that before the OTA interoperability test, the frequency calibration needs to be manually performed on the signal source device in advance, and the test can be performed after the synchronization is ensured. Obviously, this brings workload to laboratory management, and for a large-scale laboratory, due to the large variety of instruments and the small number of each instrument, the management cost is directly increased. Therefore, it is a problem to be solved in the prior art to provide an automatic calibration method for laboratory instruments, which improves the coverage rate and accuracy of automatic calibration of equipment.

Disclosure of Invention

In order to solve the above problems, the present application provides an automatic calibration apparatus for an electrical experimental instrument, which diagnoses whether an external reference clock is synchronized with a standard clock source by detecting a frequency update characteristic of the external reference clock.

The application provides electric power laboratory glassware automatic calibration device includes: the test system comprises a test configuration module, an auxiliary module, a frequency synchronization module and an autonomous operation module, wherein the modules are matched with each other to calibrate the instrument, and the steps are as follows:

step 1: the test configuration module receives the signal source configuration information, if the configuration information enables an external reference clock, a recording list k is newly built in a frequency difference recording component in a database, and the step 2 is skipped, otherwise, the step 4 is skipped;

step 2: the auxiliary module inquires a clock frequency error requirement F of a corresponding communication system from a database according to the configuration information, and then reads an initial configuration value of the crystal oscillator from the database and configures the crystal oscillator;

and step 3: the frequency synchronization module periodically calculates a crystal oscillator adjustment offset value A corresponding to the frequency synchronization of the external reference clock according to the period A, and dynamically forms a non-configuration indication signal tNoConfig; during the period that tNoConfig is 0, superposing and configuring the crystal oscillator based on the period A and the bias value Ap, and recording each time point and the bias value Ap to a PartA sublist in a list k; periodically calculating a crystal oscillator adjustment offset value Bp corresponding to the frequency synchronization of an external reference clock according to a period B during the period tNoConfig is 1, and recording each time point and the offset value Bp into a PartB sub-table in a list k; after the test is detected to be finished, writing the clock frequency error requirement F and the end time point information T into PartC in a list k, and stopping updating information to the list k;

and 4, step 4: the autonomous operation module judges whether a frequency difference recording component in the database is empty; if the configuration value is empty, reading the initial configuration value of the crystal oscillator from the database and configuring the crystal oscillator; if the frequency difference recording component is not empty, selecting a signal source frequency initial calibration optimal reference queue from each group of data of PartB in each recording list of the frequency difference recording component in the database, determining a superposition configuration value corresponding to an Nth deviation value of which the root mean square of the latest continuous N times of deviation values of PartA in the optimal reference queue is smaller than a threshold Am as an initial configuration value of the crystal oscillator, writing the initial configuration value into the database, deleting all information of the frequency difference recording component in the database, configuring the crystal oscillator according to the initial configuration value of the crystal oscillator, and outputting a corresponding signal to serve by the signal source according to the configuration information.

Preferably, in the step 3, the method for performing superposition configuration on the crystal oscillator based on the period a and the offset value Ap includes: at time T0, the crystal oscillator is initialized, the crystal oscillator configuration value is V0, and after L periods a, the crystal oscillator configuration value Vl is V0+ ∑ APl, L is 1,2 … L, and Apl is the offset value calculated for each period a.

Preferably, the frequency synchronization module comprises a reference source construction sub-module, an initial synchronization sub-module, a dynamic management module and a homology identification module.

Preferably, in the step 3, a specific method for dynamically forming the non-configuration period tNoConfig _ i includes:

3.1, constructing a reference pulse A with a period of A and a reference pulse B with a period of B by using an external reference clock by a reference source constructing submodule, wherein the time length of the period A is G times of the time length of the period B, and the boundary of each G period B is aligned with the boundary of the period A;

3.2, the initial synchronization submodule is periodically synchronized with the reference pulse A, and the local crystal oscillator of the signal source is adjusted until the average value of M continuous deviation values DeltaA _ i which are the latest frequency of the local crystal oscillator of the signal source and the external reference clock frequency is less than KxStep, setting tNoConfig to 1, skipping to the Step 3.3, otherwise, repeating the operation of the Step 3.2;

3.3, the dynamic management module calculates the frequency error DeltaA _ i of the local crystal oscillator of the signal source A and the external reference clock in the current period, if the DeltaA _ i is smaller than the frequency error (F-K multiplied by Step), the local crystal oscillator of the signal source is not configured, the tNoConfig state is not updated, and the operation of the Step 3.3 is repeated; otherwise, configuring a local crystal oscillator of the signal source according to the DeltaA _ i, setting tNoConfig to 0, and then jumping to the step 3.2.

Preferably, Step is a maximum frequency Step value corresponding to an adjustment granularity of the local crystal oscillator of the signal source, and the value of K is greater than or equal to 2.

Preferably, in step 3.3, during the period of tNoConfig being 1, the homology identification module calculates the frequency deviation DeltaB _ j _ p between the local crystal oscillator of the signal source and the external reference clock according to the period B, and writes the statistical result of the period of tNoConfig continuing to be 1 each time into the p-th group in the PartB sublist in the list k in the database, where p represents the number of times tNoConfig jumps from 0 to 1, and j is the number of the period B during the period of tNoConfig being 1 after the p-th jump of tNoConfig from 0 to 1, and is counted from 0.

Preferably, in the step 4, the optimal reference queue for initial calibration of the signal source frequency is selected from each group of data of the PartB in each record list of the frequency difference record component in the database, and the specific method is as follows:

step 4.1, the sliding window calculation module obtains each group of data of the PartB in each list, then the following operation is carried out on each group of data in each list, and the frequency characteristic value Q _ p _ k of the group of data is determined, wherein the specific calculation method is as follows: taking the period B as a step and the period PeriodSel as a time length, calculating a frequency characteristic value Q _ w _ p _ k in each PeriodSel time length in the group of data, specifically, calculating an average value AvgB _ w _ p _ k of DeltaB _ j _ p _ k in each sliding window time length PeriodSel, calculating a start-stop time point information set SetB _ up _ w _ p _ k corresponding to a continuous time length larger than AvgB _ w _ p _ k in the sliding window time length PeriodSel, calculating a start-stop time point information set SetB _ up _ w _ p _ k corresponding to a continuous time length smaller than AvgB _ w _ p _ k, and subtracting a minimum time length Min _ w _ p _ k from a maximum time length Max _ p _ k in SetB _ up _ w _ p _ k and SetB _ down _ w _ p _ k to obtain a weight value Q _ w _ k of a parp in a parp _ p _ k list in the SetB _ w _ p _ k;

step 4.2, the list characteristic value determining module determines the smallest value in Q _ w _ p _ k corresponding to p and w values in the list k as the frequency characteristic value Q _ k of the list k;

4.3, the Type management module counts the Type number Type _ Num of the frequency error requirements according to the frequency error requirements F of PartC in each list k in the database, and then divides the list k into different types of subsets SetType _ h according to the frequency error requirements, wherein the value of h is 0,.

Step 4.4, calculating the average value AvgInType _ h of each list frequency characteristic value Q _ k in a single Type in SetType _ h by an intra-Type frequency characteristic value optimization module, and then determining the frequency characteristic value corresponding to the list with the maximum time point information T at the end of the elements with Q _ k smaller than AvgInType _ h in the Type as the frequency characteristic value Q _ Type _ h of the Type;

and 4.5, selecting a preferred module for the inter-Type frequency characteristic values, removing the types of which Q _ Type _ h is larger than the threshold Bm, sequencing the Q _ Type _ h of different types in the rest types from large to small based on time points, calculating the time difference between the time point of each Type and the maximum time point, obtaining the frequency error deterioration value of each Type according to the mapping relation between the time length and the frequency difference, superposing the frequency error deterioration value of each Type and the frequency error requirement value of each Type to obtain the final frequency error value of each Type, and selecting the queue with the minimum final frequency error value as the optimal reference queue for frequency calibration.

Compared with the prior art, the beneficial effect of this application is:

according to the method and the device, whether the external reference clock is synchronous with the standard clock source or not is diagnosed by detecting the frequency updating characteristic of the external reference clock, if so, the frequency of the external reference clock source is determined to be accurate and reliable, and the initial configuration of the local crystal oscillator of the signal source is synchronized to the external reference clock, so that the frequency self-adaptive calibration is completed, and the coverage rate and the accuracy of the automatic calibration of the equipment are improved.

Drawings

The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the application and, together with the description, serve to explain the application and are not intended to limit the application.

Figure 1 is a flow chart of a method of one embodiment of the present application,

fig. 2 is a system component diagram of an embodiment of the present application.

The specific implementation mode is as follows:

the present application will be further described with reference to the following drawings and examples.

It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments according to the present disclosure. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.

In the present disclosure, terms such as "upper", "lower", "left", "right", "front", "rear", "vertical", "horizontal", "side", "bottom", and the like indicate orientations or positional relationships based on those shown in the drawings, and are only relational terms determined for convenience in describing structural relationships of the parts or elements of the present disclosure, and do not refer to any parts or elements of the present disclosure, and are not to be construed as limiting the present disclosure.

As shown in fig. 1 to 2, the present application provides an automatic calibration device for an electrical power laboratory instrument, comprising: the system comprises a test configuration module, an auxiliary module, a frequency synchronization module and an autonomous operation module, wherein the functions of the modules are as follows:

a test configuration module: the module is responsible for receiving signal source configuration information, creating a recording list by a frequency difference recording component in a database and controlling a test flow according to the enablement of an external reference clock.

An auxiliary module: the module inquires a clock frequency error requirement F of a corresponding communication system from a database according to the configuration information, and then reads an initial configuration value of the crystal oscillator from the database and configures the crystal oscillator.

A frequency synchronization module: the module periodically calculates a crystal oscillator adjustment offset value Ap corresponding to the frequency synchronization of an external reference clock according to a period A, dynamically forms a non-configuration indicating signal tNoConfig, performs superposition configuration on the crystal oscillator based on the period A and the offset value Ap when the tNoConfig is 0, and records each time point and the offset value Ap to a PartA sub-table in a list k; periodically calculating a crystal oscillator adjustment offset value Bp corresponding to the frequency synchronization of an external reference clock according to a period B during the period tNoConfig is 1, and recording each time point and the offset value Bp into a PartB sub-table in a list k; and after the test is detected to be finished, writing the clock frequency error requirement F into PartC in the list k, and stopping updating information to the list k.

An autonomous operation module: the module is responsible for judging whether a frequency difference recording component in the database is empty, and if the frequency difference recording component in the database is empty, reading an initial configuration value of the crystal oscillator from the database and configuring the crystal oscillator; if the data is not null, calculating the frequency validity of each group of data of PartC in each recording list of the frequency difference recording component in the database, determining a superposition configuration value corresponding to the Nth deviation value of which the root mean square of the last continuous N times of deviation values of PartA in the recording list with the highest validity is less than the threshold Am as the initial configuration value of the crystal oscillator, writing the initial configuration value into the database, deleting all information of the frequency difference recording component in the database, configuring the crystal oscillator according to the initial configuration value of the crystal oscillator, and outputting a corresponding signal according to the configuration information by a signal source to serve.

The steps of calibrating the instrument by mutually matching the modules are as follows:

step 1: the test configuration module receives the signal source configuration information, if the configuration information enables an external reference clock, a recording list k is newly built in a frequency difference recording component in a database, and the step 2 is skipped, otherwise, the step 4 is skipped;

and 2, step: the auxiliary module inquires a clock frequency error requirement F of a corresponding communication system from a database according to the configuration information, and then reads an initial configuration value of the crystal oscillator from the database and configures the crystal oscillator;

and step 3: the frequency synchronization module periodically calculates a crystal oscillator adjustment offset value A corresponding to the frequency synchronization of the external reference clock according to the period A, and dynamically forms a non-configuration indication signal tNoConfig;

during the period that tNoConfig is 0, superposing and configuring the crystal oscillator based on the period A and the bias value Ap, and recording each time point and the bias value Ap to a PartA sub-table in a list k;

periodically calculating a crystal oscillator adjustment offset value Bp corresponding to the frequency synchronization of an external reference clock according to a period B during the period tNoConfig is 1, and recording each time point and the offset value Bp into a PartB sub-table in a list k; after the test is detected to be finished, writing the clock frequency error requirement F and the end time point information T into PartC in a list k, and stopping updating information to the list k;

and 4, step 4: the autonomous operation module judges whether a frequency difference recording component in the database is empty;

if the configuration value is empty, reading the initial configuration value of the crystal oscillator from the database and configuring the crystal oscillator;

if the frequency difference recording component is not empty, selecting a signal source frequency initial calibration optimal reference queue from each group of data of PartB in each recording list of the frequency difference recording component in the database, determining a superposition configuration value corresponding to an Nth deviation value of which the root mean square of the latest continuous N times deviation value of PartA in the optimal reference queue is smaller than a threshold Am as an initial configuration value of the crystal oscillator, writing the initial configuration value into the database, deleting all information of the frequency difference recording component in the database, configuring the crystal oscillator according to the initial configuration value of the crystal oscillator, and outputting a corresponding signal to serve by the signal source according to the configuration information.

The application also provides an automatic calibration method of the experimental instrument, and the specific steps are consistent with the steps 1 to 4.

In step 1, the configuration information enables the external reference clock, which is an industry standard term, to mean that the configuration information indicates that the outsourced reference clock is valid.

In the step 3, the method for performing superposition configuration on the crystal oscillator based on the period a and the offset value Ap includes:

the crystal oscillator is initialized at time T0, the crystal oscillator configuration value is V0, the time T is reached after L cycles a, the crystal oscillator configuration value Vl is V0+ ∑ APl, L is 1,2 … L, Apl is the offset value calculated for each cycle a, for example, the offset value is Ap1 when the cycle a reaches time T1, the crystal oscillator configuration value at time T1 is V0+ a1, the offset value is Ap2 when the cycle a reaches time T2, the crystal oscillator configuration value at time T2 is V0+ Ap1+ Ap2, and so on, where Ap1 and Ap2 may take positive, zero, and negative numbers.

The frequency synchronization module comprises a reference source construction submodule, an initial synchronization submodule, a dynamic management module and a homologous identification module, and in the step 3, a specific method for dynamically forming a non-configuration time period tNoConfig _ i comprises the following steps:

3.1, constructing a reference pulse A with a period of A and a reference pulse B with a period of B by using an external reference clock by a reference source constructing submodule, wherein the time length of the period A is G times of the time length of the period B, and the boundary of each G period B is aligned with the boundary of the period A;

3.2, the initial synchronization submodule is periodically synchronized with the reference pulse A, and the local crystal oscillator of the signal source is adjusted until the average value of M continuous deviation values DeltaA _ i which are the latest frequency of the local crystal oscillator of the signal source and the external reference clock frequency is less than KxStep, setting tNoConfig to 1, skipping to the Step 3.3, otherwise, repeating the operation of the Step 3.2;

3.3, the dynamic management module calculates the frequency error DeltaA _ i of the local crystal oscillator of the signal source A and the external reference clock in the current period, if the DeltaA _ i is smaller than the frequency error (F-K multiplied by Step), the local crystal oscillator of the signal source is not configured, the tNoConfig state is not updated, and the operation of the Step 3.3 is repeated; otherwise, configuring a local crystal oscillator of the signal source according to the DeltaA _ i, setting tNoConfig to 0, and then jumping to the step 3.2.

In the Step 3.2, Step is a maximum frequency Step value corresponding to an adjustment granularity of the local crystal oscillator of the signal source, and preferably, the value of K is greater than or equal to 2.

In the step 3.3, during the period tNoConfig is 1, the homology identification module calculates the frequency deviation DeltaB _ j _ p between the local crystal oscillator of the signal source and the external reference clock according to the period B, and writes the statistical result of the period in which tNoConfig continues to be 1 each time into the group p in the part B sub-table in the list k in the database, wherein p represents the number of times that tNoConfig jumps from 0 to 1, and j is the number of the period B during the period in which tNoConfig is 1 after the period p of tNoConfig jumps from 0 to 1, and is counted from 0.

In the step 4, an optimal reference queue for initial calibration of the frequency of the signal source is selected from each group of data of the PartB in each record list of the frequency difference record components in the database, and the specific method is as follows:

step 4.1, the sliding window calculation module obtains each group of data of the PartB in each list, then the following operation is carried out on each group of data in each list, and the frequency characteristic value Q _ p _ k of the group of data is determined, wherein the specific calculation method is as follows: taking the period B as a step and the period PeriodSel as a time length, calculating a frequency characteristic value Q _ w _ p _ k in each PeriodSel time length in the group of data, specifically, calculating an average value AvgB _ w _ p _ k of DeltaB _ j _ p _ k in each sliding window time length PeriodSel, calculating a start-stop time point information set SetB _ up _ w _ p _ k corresponding to a continuous time length larger than AvgB _ w _ p _ k in the sliding window time length PeriodSel, calculating a start-stop time point information set SetB _ up _ w _ p _ k corresponding to a continuous time length smaller than AvgB _ w _ p _ k, and subtracting a minimum time length Min _ w _ p _ k from a maximum time length Max _ p _ k in SetB _ up _ w _ p _ k and SetB _ down _ w _ p _ k to obtain a weight value Q _ w _ k of a parp in a parp _ p _ k list in the SetB _ w _ p _ k;

step 4.2, a list characteristic value determining module determines the smallest value in Q _ w _ p _ k corresponding to p and w values in the list k as a frequency characteristic value Q _ k of the list k;

4.3, the Type management module counts the Type number Type _ Num of the frequency error requirements according to the frequency error requirements F of PartC in each list k in the database, and then divides the list k into different types of subsets SetType _ h according to the frequency error requirements, wherein the value of h is 0,.

Step 4.4, calculating the average value AvgInType _ h of the frequency characteristic values Q _ k of each list in a single Type in SetType _ h by a frequency characteristic value optimizing module in the Type, and then determining the frequency characteristic value corresponding to the list with the maximum information T at the ending time point in the elements with Q _ k smaller than AvgInType _ h in the Type as the frequency characteristic value Q _ Type _ h of the Type;

and 4.5, selecting a preferred module for the inter-Type frequency characteristic values, removing the types of which Q _ Type _ h is larger than the threshold Bm, sequencing the Q _ Type _ h of different types in the rest types from large to small based on time points, calculating the time difference between the time point of each Type and the maximum time point, obtaining the frequency error deterioration value of each Type according to the mapping relation between the time length and the frequency difference, superposing the frequency error deterioration value of each Type and the frequency error requirement value of each Type to obtain the final frequency error value of each Type, and selecting the queue with the minimum final frequency error value as the optimal reference queue for frequency calibration.

The following describes a specific implementation of an automatic calibration device for an electrical experimental instrument by using specific examples:

in the embodiment, six TESTs are included, namely TEST0, TEST1, TEST2, TEST3, TEST4, and TEST5, wherein the first five TESTs are TESTs enabling an external reference clock, the sixth TEST is a TEST disabling the reference clock, TEST0 and TEST1 both belong to the TEST of communication system 0, and both the clock frequency error requirements are F0_1, TEST2 and TEST3 both belong to the TEST of communication system 1, both the clock frequency error requirements are F2_3 and TEST4 belong to the TEST of communication system 2, and the clock frequency error requirement is F4, wherein the TEST data of TEST0 is shown in table 1 in detail, in the embodiment, M takes a value of 3, K takes a value of 5, Step takes a value of 10Hz, K × Step takes a value of 50Hz, F0_1 takes a value of 150Hz, as shown in table 1, Step sel takes a value of PeriodSel, cycle B is equal to 4 cycles A, cycle B is equal to 2, a unit cycle, K takes a value of 2 unit time duration A takes a value, and K takes a value of 2 unit time duration is equal to 2 unit A, the threshold B is 20% of the duration of the period B.

First, step 1 is performed: the TEST configuration module receives the signal source configuration information TEST0, detects that the configuration information enables the external reference clock, creates a new recording list k (k takes a value of 0 at this time) in the frequency offset recording component in the database, and jumps to step 2.

k is the list number, starting with 0, i.e., list 0, list 1,. the numbering is incremented in sequence, embodiments begin with list 0, list 0 corresponds to TEST0, and the created list number is incremented by 1 each time a TEST is performed.

Then, step 2 is executed: the auxiliary module inquires a clock frequency error requirement F (corresponding to F0_1, namely 150Hz) of a corresponding communication system from the database according to the configuration information, and then reads an initial configuration value of the crystal oscillator from the database and configures the crystal oscillator.

Then, step 3 is executed: the frequency synchronization module periodically calculates a crystal oscillator adjustment offset value Ap corresponding to the frequency synchronization of the external reference clock according to the period A, dynamically forms a non-configuration indicating signal tNoConfig, performs superposition configuration on the crystal oscillator based on the period A and the offset value Ap when the tNoConfig is 0, and records each time point and the offset value Ap to a PartA sublist in a list k; periodically calculating a crystal oscillator adjustment offset value Bp corresponding to the frequency synchronization of an external reference clock according to a period B during the period tNoConfig is 1, and recording each time point and the offset value Bp into a PartB sub-table in a list k; after the test is detected to be finished, the clock frequency error requirement F and the finishing time point information T are written into PartC in the list k, and the information updating to the list k is stopped.

In step 3, the specific method for dynamically forming the non-configuration period tNoConfig _ i includes:

first, step 3.1 is executed, the reference source building submodule uses the external reference clock to build a reference pulse a with a period a (4 time units in this embodiment, for example, 10-13 in table 1 is one period, and 14-17 is the next period), and a reference pulse B with a period B (2 time units in this embodiment, for example, 10-11 in table 1 is one period, and 12-13 is the next period), where the duration of the period a is G times (i.e., 2 times) the duration of the period B, and the boundary of each G period B is aligned with the boundary of the period a.

Then, Step 3.2 is executed, the initial synchronization submodule is periodically synchronized with the reference pulse a, and adjusts the local crystal oscillator of the signal source, until the average value of the latest continuous M deviation values DeltaA _ i of the local crystal oscillator of the signal source and the frequency of the external reference clock is smaller than K × Step, then tNoConfig is set to 1, and Step 3.3 is skipped, otherwise, the operation of Step 3.2 is repeated, Step is a maximum frequency Step value corresponding to one adjustment granularity of the local crystal oscillator of the signal source, preferably, K is not smaller than 2, as shown in table 1, at a time point when T is equal to 13, a small area K × Step (i.e., smaller than 50Hz) of the average value of the continuous M (in this embodiment, 3) deviation values DeltaA _ i is determined, and then at the next time point (i.e., when T is equal to 14), tNoConfig is set to 1, and Step 3.3 is skipped.

Then, Step 3.3 is executed, the dynamic management module calculates the frequency error DeltaA _ i between the local crystal oscillator of the signal source A and the external reference clock in the current period, if the DeltaA _ i is smaller than the frequency error (F-K multiplied by Step), the local crystal oscillator of the signal source is not configured, the tNoConfig state is not updated, and the operation of Step 3.3 is repeated; otherwise, configuring a local crystal oscillator of the signal source according to the DeltaA _ i, setting tNoConfig to 0, and then jumping to the step 3.2; in this embodiment, (F-K × Step) ═ F0_1-K × Step (150Hz-5 × 10Hz) ═ 100Hz, refer to table 1, after the time point T equals 100013, tNoConfig is set to 0 and shifted to Step 3.2 to perform local oscillator adjustment since DeltaA _ i equals 110+10 Hz, which is greater than (F-K × Step), i.e., greater than 100Hz, and in the period from T, etc. 14 to T equals 100013, tNoConfig is set to 1, and therefore, local oscillator adjustment is not performed.

And then TESTs of TEST1, TEST2, TEST3 and TEST4 are carried out, corresponding data information is generated, list 1, list 2, list 3 and list 4 data are formed, the data generation process is the same as the principle, and the data formats of the list 1, the list 2, the list 3 and the list 4 are different from the data format of the list 1 (different columns).

Next, the TEST of TEST5 is performed, and since the configuration of TEST5 does not enable the external reference clock, a jump is made directly to step 4.

And step 4 is executed, the autonomous operation module judges that the frequency difference recording component in the database is not empty, a signal source frequency initial calibration optimal reference queue is selected from each group of data of PartB in each recording list of the frequency difference recording component in the database, a superposition configuration value corresponding to an Nth deviation value of which the root mean square of the last continuous N times of deviation values of PartA in the optimal reference queue is smaller than a threshold Am is determined as an initial configuration value of the crystal oscillator and is written into the database, all information of the frequency difference recording component in the database is deleted, the crystal oscillator is configured according to the initial configuration value of the crystal oscillator, and then a corresponding signal is output according to the configuration information of the signal source to serve.

In step 4, the optimal reference queue for initial calibration of the signal source frequency is selected from each group of data of the PartB in each record list of the frequency difference record component in the database, and the specific method is as follows:

step 4.1, the sliding window calculation module obtains each group of data of the PartB under each list, then performs the following operations on each group of data under each list, determines the frequency characteristic value Q _ p _ k of the group of data, taking the data of TEST0 as an example, and explains a specific calculation method as follows, firstly, only one group of data in TEST0 corresponds to T in table 1 to be 14 to 100013, then the sliding window calculation module calculates the frequency characteristic value Q _ w _ p _ k in each period of period B as stepping, and taking period of period dsel as duration, referring to table 1, then calculates average AvgB _ w _ p _ k of 993 sliding windows of T being 14 to 29, T being 16 to 31, and T being 99998 to 100013 in turn, w being a sliding window number, corresponding to 0 to 49993, and then counts the information of each sliding window, taking sliding window 0 as an example, and calculates the following steps: t is calculated to be equal to the average AvgB _0_ p _ k over the 8 time granularities 14 to 29,

AvgB_0_p_k＝(DeltaB_0_p_k+DeltaB_1_p_k+DeltaB_2_p_k+DeltaB_3_p_k+DeltaB_4_p_k+DeltaB_5_p_k+DeltaB_6_p_k+DeltaB_7_p_k)/8

＝((52+10)+(52+9)+(51-8)+(51-9)+(50+11)+(50+11)+(52-9)+(52-11))/8

＝51.75

then, the start-stop time point information set SetB _ up _ w _ p _ k corresponding to the continuous time period greater than AvgB _0_ p _ k in the time period PeriodSel (corresponding to T being equal to 14 to 29) of the sliding window 0 (corresponding to T being equal to 14 to 29) { [14,17], [22,25] }, the start-stop time point information set SetB _ up _ w _ p _ k corresponding to the continuous time period less than AvgB _0_ p _ k { [18,21], [26,29] }, and then the minimum time period Min _0_ p _ k is subtracted from the maximum time period Max _0_ p _ k to obtain the weight Q _ w _ p _ k of the w-th sliding window of the p-th group in the PartB sublist in the list k.

Then, step 4.2 is executed, and the list characteristic value determining module determines the smallest value of Q _ w _ p _ k corresponding to p and w values in the list k as the frequency characteristic value Q _ k of the list k (if a plurality of Q _ w _ p _ k values are the same, the largest time point is taken).

Then, step 4.3 is executed, the Type management module counts the number of frequency error requirement types Type _ Num according to the frequency error requirement F of the PartC in each list k in the database, and then divides the list k into different types of subsets SetType _ h according to the frequency error requirement, wherein h takes values of 0, 9, Type _ Num-1, in this embodiment, three types of Type _ Num are equal to 3, wherein TEST0, TEST1 are divided into SetType _0, TEST2, TEST3 are divided into SetType _1, and TEST4 is divided into SetType _ 2.

Then, step 4.4 is executed, the intra-Type frequency feature value optimization module calculates a mean AvgInType _ h of the respective list frequency feature values Q _ k in a single Type in SetType _ h, and then determines the frequency feature value corresponding to the list with the largest ending time point information T in the element with Q _ k smaller than AvgInType _ h in the Type as the frequency feature value Q _ Type _ h of the Type, in this embodiment, if the selected TEST0 is the frequency feature value of Type SetType _0, the TEST2 is the frequency feature value of Type SetType _1, and the TEST4 is the frequency feature value of Type SetType _2, if the frequency error requirements of the respective types in this embodiment are as follows:

SetType _ 0: the occurrence time point T is equal to 100013, and the frequency error requirement is 150 Hz;

SetType _ 1: the occurrence time point T is equal to 105013, and the frequency error requirement is 300 Hz;

SetType _ 2: the time of occurrence T equals 109013 and the frequency error requirement is 1000 Hz.

Then, step 4.5 is executed, assuming that Q _ Type _ h is all smaller than the threshold B, so that the frequency characteristic value optimization module between types does not reject any Type, then sorting different types of Q _ Type _ h from large to small based on time points to obtain SetType _2> SetType _1> SetType _0, calculating to obtain that the time difference between SetType _1 and SetType _2 is 5000 time granularities, and the time difference between SetType _0 and SetType _2 is 9000 time granularities, according to the mapping relationship between time length and frequency difference:

assuming that every 1000 time granularities in this embodiment and the frequency changes by 10Hz, the maximum time point, i.e. the time point of SetType _2, corresponds to the time point T equal to 109013, and the final frequency error values of each type are as follows:

the final frequency error of SetType _0 is maximum (150+ (109013-;

the final frequency error of SetType _1 is at most (300+ (109013-;

the final frequency error of SetType _2 is at most (1000+ (109013-;

therefore, the SetType _0 with the minimum value is selected as a queue with the minimum error, a superposition configuration value corresponding to the Nth deviation value of which the root mean square of the last continuous N deviation values of the PartA of the queue is smaller than the threshold A is determined as an initial configuration value of the crystal oscillator and is written into a database, all information of the frequency difference recording assembly in the database is deleted, the crystal oscillator is configured according to the initial configuration value of the crystal oscillator, and then a signal source outputs a corresponding signal according to the configuration information to serve.

TABLE 1 TEST0 TEST procedures database record information (corresponding to the frequency offset record component List 0)

It can be seen from the above embodiments that, by using the method of the present invention, through detecting the frequency updating characteristic of the external reference clock, it is diagnosed whether the external reference clock is synchronized with the standard clock source, if so, the frequency of the external reference clock source is determined to be accurate and reliable, and the initial configuration of the local crystal oscillator of the signal source is synchronized to the external reference clock, thereby completing the frequency adaptive calibration and improving the coverage rate and accuracy of the automatic calibration of the device.

The above description is only a preferred embodiment of the present application and is not intended to limit the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Although the embodiments of the present application have been described with reference to the accompanying drawings, it is not intended to limit the scope of the present application, and it should be understood by those skilled in the art that various modifications and variations can be made without inventive effort by those skilled in the art.

Claims (7)

1. Electric power laboratory glassware automatic calibration device, its characterized in that includes: the test system comprises a test configuration module, an auxiliary module, a frequency synchronization module and an autonomous operation module, wherein the modules are matched with each other to calibrate the instrument, and the steps are as follows:

step 1: the test configuration module receives the signal source configuration information, if the configuration information enables an external reference clock, a recording list k is newly built in a frequency difference recording component in a database, and the step 2 is skipped, otherwise, the step 4 is skipped;

step 2: the auxiliary module inquires a clock frequency error requirement F of a corresponding communication system from a database according to the configuration information, and then reads an initial configuration value of the crystal oscillator from the database and configures the crystal oscillator;

and step 3: the frequency synchronization module periodically calculates a crystal oscillator adjustment offset value A corresponding to the frequency synchronization of the external reference clock according to the period A, and dynamically forms a non-configuration indication signal tNoConfig;

during the period that tNoConfig is 0, superposing and configuring the crystal oscillator based on the period A and the bias value Ap, and recording each time point and the bias value Ap to a PartA sublist in a list k;

periodically calculating a crystal oscillator adjustment offset value Bp corresponding to the frequency synchronization of an external reference clock according to a period B during the period tNoConfig is 1, and recording each time point and the offset value Bp into a PartB sub-table in a list k; after the test is detected to be finished, writing a clock frequency error requirement F and end time point information T into PartC in a list k, and stopping updating information to the list k;

and 4, step 4: the autonomous operation module judges whether a frequency difference recording component in the database is empty;

if the configuration value is empty, reading the initial configuration value of the crystal oscillator from the database and configuring the crystal oscillator;

if the frequency difference recording component is not empty, selecting a signal source frequency initial calibration optimal reference queue from each group of data of PartB in each recording list of the frequency difference recording component in the database, determining a superposition configuration value corresponding to an Nth deviation value of which the root mean square of the latest continuous N times deviation value of PartA in the optimal reference queue is smaller than a threshold Am as an initial configuration value of the crystal oscillator, writing the initial configuration value into the database, deleting all information of the frequency difference recording component in the database, configuring the crystal oscillator according to the initial configuration value of the crystal oscillator, and outputting a corresponding signal to serve by the signal source according to the configuration information.

2. The automatic calibration device for electric power experimental instrument according to claim 1, characterized in that:

in the step 3, the method for performing superposition configuration on the crystal oscillator based on the period a and the offset value Ap includes:

at time T0, the crystal oscillator is initialized, the crystal oscillator configuration value is V0, and after L cycles a to time Tl, the crystal oscillator configuration value Vl is V0+ ∑ APl, L is 1,2 … L, and Apl is the offset value calculated for each cycle a.

3. The automatic calibration device for electric power experimental instrument according to claim 2, wherein:

the frequency synchronization module comprises a reference source construction submodule, an initial synchronization submodule, a dynamic management module and a homologous identification module.

4. The automatic calibration device for electric power experimental instrument according to claim 3, characterized in that:

in step 3, the specific method for dynamically forming the non-configuration period tNoConfig _ i includes:

3.1, constructing a reference pulse A with a period of A and a reference pulse B with a period of B by using an external reference clock by a reference source constructing submodule, wherein the time length of the period A is G times of the time length of the period B, and the boundary of each G period B is aligned with the boundary of the period A;

3.2, the initial synchronization submodule is periodically synchronized with the reference pulse A, and the local crystal oscillator of the signal source is adjusted until the average value of M continuous deviation values DeltaA _ i which are the latest frequency of the local crystal oscillator of the signal source and the external reference clock frequency is less than KxStep, setting tNoConfig to 1, skipping to the Step 3.3, otherwise, repeating the operation of the Step 3.2;

3.3, the dynamic management module calculates the frequency error DeltaA _ i of the local crystal oscillator of the signal source A and the external reference clock in the current period, if the DeltaA _ i is smaller than the frequency error (F-K multiplied by Step), the local crystal oscillator of the signal source is not configured, the tNoConfig state is not updated, and the operation of the Step 3.3 is repeated; otherwise, configuring a local crystal oscillator of the signal source according to the DeltaA _ i, setting tNoConfig to 0, and then jumping to the step 3.2.

5. The automatic calibration device for electric power experimental instrument according to claim 4, characterized in that:

step is a maximum frequency Step value corresponding to the adjustment granularity of the local crystal oscillator of the signal source, and the value of K is more than or equal to 2.

6. The automatic calibration device for electric power experimental instrument according to claim 4, characterized in that:

in the step 3.3, during the period of tNoConfig being 1, the homology identification module calculates the frequency deviation DeltaB _ j _ p between the local crystal oscillator of the signal source and the external reference clock according to the period B, and writes the statistical result of the period of tNoConfig continuing to be 1 each time into the p-th group in the PartB sublist in the list k in the database, where p represents the number of times that tNoConfig jumps from 0 to 1, and j is the number of the period B during the period of tNoConfig being 1 after the p-th time that tNoConfig jumps from 0 to 1, and counted from 0.

7. The automatic calibration device for electric power experimental instrument according to claim 6, characterized in that:

in the step 4, an optimal reference queue for initial calibration of the frequency of the signal source is selected from each group of data of the PartB in each record list of the frequency difference record components in the database, and the specific method is as follows:

step 4.1, the sliding window calculation module obtains each group of data of the PartB in each list, then the following operation is carried out on each group of data in each list, and the frequency characteristic value Q _ p _ k of the group of data is determined, wherein the specific calculation method is as follows: taking the period B as a step and the period PeriodSel as a time length, calculating a frequency characteristic value Q _ w _ p _ k in each PeriodSel time length in the group of data, specifically, calculating an average value AvgB _ w _ p _ k of DeltaB _ j _ p _ k in each sliding window time length PeriodSel, calculating a start-stop time point information set SetB _ up _ w _ p _ k corresponding to a continuous time length larger than AvgB _ w _ p _ k in the sliding window time length PeriodSel, calculating a start-stop time point information set SetB _ up _ w _ p _ k corresponding to a continuous time length smaller than AvgB _ w _ p _ k, and subtracting a minimum time length Min _ w _ p _ k from a maximum time length Max _ p _ k in SetB _ up _ w _ p _ k and SetB _ down _ w _ p _ k to obtain a weight value Q _ w _ k of a parp in a parp _ p _ k list in the SetB _ w _ p _ k;

step 4.2, the list characteristic value determining module determines the smallest value in Q _ w _ p _ k corresponding to p and w values in the list k as the frequency characteristic value Q _ k of the list k;

4.3, the Type management module counts the Type number Type _ Num of the frequency error requirements according to the frequency error requirements F of PartC in each list k in the database, and then divides the list k into different types of subsets SetType _ h according to the frequency error requirements, wherein the value of h is 0,.

Step 4.4, calculating the average value AvgInType _ h of each list frequency characteristic value Q _ k in a single Type in SetType _ h by an intra-Type frequency characteristic value optimization module, and then determining the frequency characteristic value corresponding to the list with the maximum time point information T at the end of the elements with Q _ k smaller than AvgInType _ h in the Type as the frequency characteristic value Q _ Type _ h of the Type;

and 4.5, selecting a preferred module for the inter-Type frequency characteristic values, removing the types of which Q _ Type _ h is larger than the threshold Bm, sequencing the Q _ Type _ h of different types in the rest types from large to small based on time points, calculating the time difference between the time point of each Type and the maximum time point, obtaining the frequency error deterioration value of each Type according to the mapping relation between the time length and the frequency difference, superposing the frequency error deterioration value of each Type and the frequency error requirement value of each Type to obtain the final frequency error value of each Type, and selecting the queue with the minimum final frequency error value as the optimal reference queue for frequency calibration.

Images