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

Automatic calibration device for electric power experimental instrument Download PDF

Info

Publication number
CN114660521B
CN114660521B CN202210293667.XA CN202210293667A CN114660521B CN 114660521 B CN114660521 B CN 114660521B CN 202210293667 A CN202210293667 A CN 202210293667A CN 114660521 B CN114660521 B CN 114660521B
Authority
CN
China
Prior art keywords
frequency
value
period
crystal oscillator
list
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210293667.XA
Other languages
Chinese (zh)
Other versions
CN114660521A (en
Inventor
刘爱华
连维坤
蒋珂
张建涛
刘修伟
岳秀锋
张天琦
陈泽宇
张良
赵鹏
洪云起
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Shandong Electric Power Co Linqing Power Supply Co
Original Assignee
State Grid Shandong Electric Power Co Linqing Power Supply Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Shandong Electric Power Co Linqing Power Supply Co filed Critical State Grid Shandong Electric Power Co Linqing Power Supply Co
Priority to CN202211556942.9A priority Critical patent/CN115877302A/en
Priority to CN202210293667.XA priority patent/CN114660521B/en
Publication of CN114660521A publication Critical patent/CN114660521A/en
Application granted granted Critical
Publication of CN114660521B publication Critical patent/CN114660521B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing or calibrating of apparatus covered by the other groups of this subclass

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 disclosure and may not necessarily 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 laboratory with a large scale, due to the fact that the variety of instruments is large and the number of each instrument is not small, the management cost is directly increased. Therefore, it is an object of the present invention to provide an automatic calibration method for laboratory instruments, which improves the coverage 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;
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 a crystal oscillator initial configuration value from the database and configures the crystal oscillator;
and 3, step 3: 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, 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 the Nth deviation value of which the root mean square of the latest continuous N times deviation value of the PartB in the optimal reference queue is smaller than the 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: and initializing the crystal oscillator at the time of T0, wherein the configuration value of the crystal oscillator is V0, the configuration value of the crystal oscillator is V0 plus Sigma APl, and L =1,2 \8230, L and Apl are deviation values calculated in each period A after L periods A reach the time of Tl.
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, the local crystal oscillator of the signal source is adjusted until the average value of M latest continuous 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 multiplied by Step, the tNoConfig is set to be 1, the Step 3.3 is skipped, and otherwise, the operation of the Step 3.2 is repeated;
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 the step 3.3, during the period tNoConfig is 1, the homology identification module calculates the frequency deviation DeltaB _ j _ p of the local crystal oscillator of the signal source from the external reference clock according to the period B, and writes the statistical result of the period that tNoConfig continues to be 1 each time into the group p in the PartB sublist in the list k in the database, where p represents the number of times that tNoConfig changes from 0 to 1, and j is the number of the period B during the period that tNoConfig is 1 after the p of tNoConfig changes 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 which is greater than AvgB _ w _ p _ k in the sliding window time length PeriodSel, a start-stop time point information set SetB _ down _ w _ p _ k corresponding to a continuous time period smaller than AvgB _ w _ p _ k, and then subtracting a minimum time length Min _ w _ p _ k from a maximum time length Max _ w _ p _ k in SetB _ up _ w _ p _ k and SetB _ down _ w _ p _ k to obtain a weight Q _ w _ p _ k of a p-th group w-th sliding window in a PartB sub table in a list 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 frequency characteristic values among the types, removing the types of which Q _ Type _ h is larger than a threshold Bm, sequencing the Q _ Type _ h of different types in the remaining types from large to small based on time points, calculating the time difference between each Type of time point and the maximum time point, obtaining the frequency error deteriorated value of each Type according to the mapping relation between time length and frequency difference, superposing each Type of frequency error deteriorated value and each Type of frequency error required value to obtain each Type of final frequency error value, and selecting a queue with the minimum final frequency error value as a frequency calibration optimal reference queue.
Compared with the prior art, the beneficial effects of this application do:
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 is 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 orientations or positional relationships shown in the drawings, are only terms of relationships determined for convenience in describing structural relationships of the components or elements of the present disclosure, do not refer to any components 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 record list by a frequency difference record 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 a crystal oscillator initial configuration value 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 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; 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 empty, calculating the frequency effectiveness 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 the PartB in the recording list with the highest effectiveness is less than the 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 according to the configuration information by a signal source to serve.
The method for calibrating the instrument by mutually matching the modules comprises the following steps:
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 a crystal oscillator initial configuration value from the database and configures the crystal oscillator;
and step 3: the frequency synchronization module periodically calculates a crystal oscillator adjustment bias value Ap 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 the Nth deviation value of which the root mean square of the latest continuous N times deviation value of the PartB in the optimal reference queue is smaller than the 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 fact that the configuration information enables the external reference clock means that the configuration information indicates that the outsourced reference clock is valid, which is an industry standard term.
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 the time T0, the crystal oscillator configuration value is V0, the time Tl is reached after L periods A, the crystal oscillator configuration value Vl = V0 plus sigma APl, L =1,2 \8230, L and Apl are deviation values calculated in each period A, for example, the deviation value is calculated to be Ap1 after the period A reaches the time T1, the crystal oscillator configuration value at the time T1 is V0 plus A1, the deviation value is calculated to be Ap2 after the period A reaches the time T2, the crystal oscillator configuration value at the time T2 is V0 plus Ap1 plus Ap2, and the like, wherein the values of Ap1 and Ap2 can be positive numbers, zero numbers 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 A and a reference pulse B with a period B by a reference source constructing submodule by adopting an external reference clock, wherein the time length of the period A is G times of that 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 A, 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 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.
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 which is greater than AvgB _ w _ p _ k in the sliding window time length PeriodSel, a start-stop time point information set SetB _ down _ w _ p _ k corresponding to a continuous time period smaller than AvgB _ w _ p _ k, and then subtracting a minimum time length Min _ w _ p _ k from a maximum time length Max _ w _ p _ k in SetB _ up _ w _ p _ k and SetB _ down _ w _ p _ k to obtain a weight Q _ w _ p _ k of a p-th group w-th sliding window in a PartB sub table in a list 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 subsets SetType _ h of different types according to the frequency error requirements, wherein the value of h is 0, 9, type _ Num-1;
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.
The following describes a specific embodiment of an automatic calibration device for an electrical power laboratory instrument by using a specific embodiment:
in the embodiment, the TEST includes six TESTs, that is, TEST0, TEST1, TEST2, TEST3, TEST4, and TEST5, where the first five TESTs are TESTs enabling an external reference clock, the sixth TEST is a TEST not enabling the reference clock, TEST0 and TEST1 both belong to communication system 0, the clock frequency error requirements of both TEST 2_3 and TEST3 both belong to communication system 1, the clock frequency error requirements of both TEST 2_3 and TEST4 belong to communication system 2, and the clock frequency error requirement is F4, where TEST0 TEST data is shown in table 1 in detail, in this embodiment, M is 3, K is 5, step is 10Hz, K × Step is 50hz, f0 threshold value 1 is 150Hz, as shown in table 1, period sel is 8 periods B, period a is equal to 4 time units, period B is equal to 2 × 2 time units, step B is 20% of time unit, and Step a is 20% of time unit.
First, step 1 is performed: and the TEST configuration module receives the signal source configuration information TEST0, detects that the configuration information enables an external reference clock, creates a new recording list k (at the moment, the value of k is 0) in the frequency difference recording component in the database, and skips to the step 2.
k is the list number, starting with 0, i.e., list 0, list 1,. The number 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 150 Hz) 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 a list k, and the information updating to the list k is stopped.
In the step 3, the specific implementation of dynamically forming the non-configuration period tNoConfig _ i is as follows:
first, step 3.1 is executed, the reference source building submodule uses an external reference clock to build a reference pulse a with a period a (4 time units in this embodiment, for example, 10-13 is one period and 14-17 is the next period in table 1), and a reference pulse B with a period B (2 time units in this embodiment, for example, 10-11 is one period and 12-13 is the next period in table 1), 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 periods 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 the signal source local crystal oscillator is adjusted until the average value of M continuous deviation values DeltaA _ i, which are the latest between the signal source local crystal oscillator and the external reference clock frequency, is less than K × Step, tNoConfig is set to 1, and the Step 3.3 is skipped, otherwise, the operation of Step 3.2 is repeated, where Step is a maximum frequency Step value corresponding to one adjustment granularity of the signal source local crystal oscillator, and preferably, K is not less than 2, as shown in table 1, at a time point where T is equal to 13, it is determined that K × Step is a small average value K × Step (i.e., less than 50 Hz) of M continuous deviation values DeltaA _ i (in this embodiment, the value is 3), and then, at the next time point (when instant time point T is equal to 14), tNoConfig is set to 1, and the Step 3.3 is skipped.
Then, executing Step 3.3, calculating a frequency error DeltaA _ i between the local crystal oscillator of the signal source A and the external reference clock in the current period A by the dynamic management module, if the DeltaA _ i is smaller than the frequency error (F-K multiplied by Step), not configuring the local crystal oscillator of the signal source, not updating the tNoConfig state, and repeating the operation of Step 3.3; 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) = (150 Hz-5 × 10 Hz) =100Hz, and referring to table 1, after T equals 100013, after time point, deltaA _ i equals 110+10=120hz, and is greater than (F-K × Step), that is, is greater than 100Hz, tNoConfig is set to 0, and the process jumps to Step 3.2 to perform local oscillator adjustment, and during the period from T equal to 14 to T equal to 100013, tNoConfig is set to 1, and thus, the local oscillator is not adjusted.
And then, TESTs of TEST1, TEST2, TEST3 and TEST4 are carried out, corresponding data information is generated, and 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 (different from) the data format of the list 1.
Next, 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 the PartB 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 the step 4, the optimal reference queue for initial calibration of the frequency of the signal source is selected from each group of data of the part b 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 PartB under each list, then the following operation is carried out on each group of data under each list, the frequency characteristic value Q _ p _ k of the group of data is determined, taking the data of TEST0 as an example, the specific calculation method is described as follows, firstly, only one group of data in TEST0 exists, T in the group of data corresponds to 14 to 100013 in Table 1, then the sliding window calculation module takes the period B as step, takes PeriodSel as time length, calculates the frequency characteristic value Q _ w _ p _ k in each PeriodSel time length in the group of data, referring to Table 1, then calculates average values AvgB _ w _ p _ k of 993 sliding windows of T equal to 14 to 29, T equal to 16 to 31, T equal to 9999999999998 to 100013, and W is a sliding window number, corresponds to 0 to 49993, then the information of each sliding window is counted, taking the sliding window 0 as an example, and the calculation process is as follows: t is calculated to be equal to the average AvgB _0_p _kover 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
<xnotran>, 0 PeriodSel ( T 14 29) AvgB _0_p_k SetB _ up _ w _ p _ k = { [14,17], [22,25] }, AvgB _0_p_k SetB _ up _ w _ p _ k = { [18,21], [26,29] }, Max _0_p_k Min _0_p_k, k PartB p w Q _ w _ p _ k. </xnotran>
Then, step 4.2 is executed, and the list characteristic value determining module determines the smallest value among 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 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 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 is 0, wherein h is equal to Type 1, and Type _ Num-1, in this embodiment, three types of Type _ Num are equal to 3, wherein TEST0 and TEST1 are divided into SetType _0, TEST2, TEST3 is divided into SetType _1, and TEST4 is divided into SetType _2.
Next, step 4.4 is executed, the intra-Type frequency feature value optimization module calculates a mean AvgInType _ h of the frequency feature values Q _ k in each list in the 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 in which Q _ k in the Type is smaller than AvgInType _ h as the frequency feature value Q _ Type _ h of the current Type, in this embodiment, if TEST0 is the frequency feature value of Type SetType _0, TEST2 is the frequency feature value of Type SetType _1, and TEST4 is the frequency feature value of Type SetType _2, if the frequency error requirements of the current embodiment are as follows:
SetType _0: the occurrence time point T is equal to 100013, and the frequency error requirement is 150Hz;
SetType _1: the occurrence time point T is equal to 105013, and the frequency error requirement is 300Hz;
SetType _2: the time of occurrence T is equal to 109013 and the frequency error requirement is 1000Hz.
Then, step 4.5 is executed, assuming that Q _ Type _ h is all smaller than the threshold B, so that the frequency feature value optimization module between types does not reject any Type, then sorting Q _ Type _ h of different types from large to small based on time points to obtain SetType _2>, setType \ u 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 duration and frequency difference:
assuming that the frequency changes by 10Hz every 1000 time granularities in this embodiment, converting the worsening effect of the upper time on the frequency error, at the maximum time point, i.e. the time point of SetType _2, corresponding to the time point T being equal to 109013, the final frequency error values of each type are as follows:
the final frequency error of SetType _0 is maximum (150 + (109013-100013)/1000 × 10) =240Hz;
the final frequency error of SetType _1 is maximally (300 + (109013-105013)/1000 × 10) =340Hz;
the final frequency error of SetType _2 is maximum (1000 + (109013-109013)/1000 x 10) =1000Hz;
therefore, 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 queue is smaller than the threshold A is determined as the 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 be in service.
TABLE 1 TEST0 test procedure database record information (corresponding to "frequency offset record component" List 0)
Figure GDA0003935290620000161
Figure GDA0003935290620000171
Figure GDA0003935290620000181
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 to the present application by those skilled in the art. Any modification, equivalent replacement, or improvement 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 (3)

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 established 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 Ap 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 deviation value of the PartB in the optimal reference queue is less 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;
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 A, 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;
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;
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 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, j is the number of the period B during the period of tNoConfig being 1 after the period p of tNoConfig jumps from 0 to 1, and the number 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 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 each 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 which is greater than AvgB _ w _ p _ k in the sliding window time length PeriodSel, a start-stop time point information set SetB _ down _ w _ p _ k corresponding to a continuous time period smaller than AvgB _ w _ p _ k, and then subtracting a minimum time length Min _ w _ p _ k from a maximum time length Max _ w _ p _ k in SetB _ up _ w _ p _ k and SetB _ down _ w _ p _ k to obtain a weight Q _ w _ p _ k of a p-th group w-th sliding window in a PartB sub table in a list 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;
step 4.5, a module for optimizing the frequency characteristic values among types, eliminating the types of which Q _ Type _ h is larger than a threshold Bm, sequencing the Q _ Type _ h of different types in the rest types from big to small based on time points, calculating the time difference between each Type of time point and the maximum time point, obtaining the frequency error deterioration value of each Type according to the mapping relation between time length and frequency difference, superposing each Type of frequency error deterioration value and each Type of frequency error requirement value to obtain each Type of final frequency error value, and selecting the queue with the minimum final frequency error value as the optimal reference queue for frequency calibration;
the DeltaB _ j _ p _ k is the DeltaB _ j _ p value in the list k.
2. The automatic calibration device for electric power laboratory 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 comprises the following steps:
and initializing the crystal oscillator at the time T0, wherein the crystal oscillator configuration value is V0, and the crystal oscillator configuration value Vl = V0 plus sigma APl and L =1, 2' -8230, and L and Apl are deviation values calculated in each period A after L periods A reach the time Tl.
3. The automatic calibration device for electric power laboratory 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.
CN202210293667.XA 2022-03-23 2022-03-23 Automatic calibration device for electric power experimental instrument Active CN114660521B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211556942.9A CN115877302A (en) 2022-03-23 2022-03-23 Automatic frequency synchronization method for test instrument in power test
CN202210293667.XA CN114660521B (en) 2022-03-23 2022-03-23 Automatic calibration device for electric power experimental instrument

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210293667.XA CN114660521B (en) 2022-03-23 2022-03-23 Automatic calibration device for electric power experimental instrument

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN202211556942.9A Division CN115877302A (en) 2022-03-23 2022-03-23 Automatic frequency synchronization method for test instrument in power test

Publications (2)

Publication Number Publication Date
CN114660521A CN114660521A (en) 2022-06-24
CN114660521B true CN114660521B (en) 2023-01-24

Family

ID=82031328

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202211556942.9A Pending CN115877302A (en) 2022-03-23 2022-03-23 Automatic frequency synchronization method for test instrument in power test
CN202210293667.XA Active CN114660521B (en) 2022-03-23 2022-03-23 Automatic calibration device for electric power experimental instrument

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN202211556942.9A Pending CN115877302A (en) 2022-03-23 2022-03-23 Automatic frequency synchronization method for test instrument in power test

Country Status (1)

Country Link
CN (2) CN115877302A (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002139585A (en) * 2000-10-31 2002-05-17 Fujitsu Ltd Correcting method and device for clock
CN101039145A (en) * 2007-03-30 2007-09-19 华为技术有限公司 Method and apparatus for realizing clock
CN104158616A (en) * 2014-07-28 2014-11-19 京信通信系统(中国)有限公司 Clock synchronization method and device in optical fiber transmission system
CN105187061A (en) * 2015-08-28 2015-12-23 京信通信系统(中国)有限公司 Crystal oscillator control method and device thereof
CN108400861A (en) * 2018-02-11 2018-08-14 京信通信系统(中国)有限公司 Crystal oscillator original frequency method of adjustment and device and clock system
CN109884877A (en) * 2019-01-15 2019-06-14 华南理工大学 A kind of high-precision GPS sync identification system and method
CN113014350A (en) * 2021-03-02 2021-06-22 北京机电工程研究所 PMC interface-based time synchronization method between simulation devices
CN113552794A (en) * 2021-06-24 2021-10-26 南方电网科学研究院有限责任公司 Automatic calibration device and method for clock signal in power chip

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002139585A (en) * 2000-10-31 2002-05-17 Fujitsu Ltd Correcting method and device for clock
CN101039145A (en) * 2007-03-30 2007-09-19 华为技术有限公司 Method and apparatus for realizing clock
CN104158616A (en) * 2014-07-28 2014-11-19 京信通信系统(中国)有限公司 Clock synchronization method and device in optical fiber transmission system
CN105187061A (en) * 2015-08-28 2015-12-23 京信通信系统(中国)有限公司 Crystal oscillator control method and device thereof
CN108400861A (en) * 2018-02-11 2018-08-14 京信通信系统(中国)有限公司 Crystal oscillator original frequency method of adjustment and device and clock system
CN109884877A (en) * 2019-01-15 2019-06-14 华南理工大学 A kind of high-precision GPS sync identification system and method
CN113014350A (en) * 2021-03-02 2021-06-22 北京机电工程研究所 PMC interface-based time synchronization method between simulation devices
CN113552794A (en) * 2021-06-24 2021-10-26 南方电网科学研究院有限责任公司 Automatic calibration device and method for clock signal in power chip

Also Published As

Publication number Publication date
CN114660521A (en) 2022-06-24
CN115877302A (en) 2023-03-31

Similar Documents

Publication Publication Date Title
JP5331656B2 (en) Power supply
CN101304265B (en) Method and apparatus for data reception
CN114660521B (en) Automatic calibration device for electric power experimental instrument
JP2023123601A (en) Systems, devices, and methods of cyclic coulometry
US20040125874A1 (en) Method of determining an eye diagram of a digital signal
CN107273322B (en) Parallel data output method and device
CN103698639B (en) Crystal oscillator waveform parameter automatic measurement system and method
CN110542722B (en) Fault detection method and device for ultrasonic probe
CN112134593B (en) Zero crossing point signal output and power line data sending method and equipment
CN103368786B (en) The method of testing of CAN bus data and device
CN1199049C (en) Regular calibration method and device for integrated circuit testing unit
US20190222437A1 (en) Switching device and determination method
CN115118618A (en) Intelligent gateway performance test method and system
US11404024B2 (en) Signal adjusting device, signal adjusting method, and test equipment
CN217766678U (en) High-low temperature test system of touch unit
CN112539523B (en) Control method and control device for air conditioner communication, communication system and storage medium
CN116755021B (en) Calibrating device and calibrating method for direct-current high-voltage divider
US10965441B1 (en) Frame trigger recreation method and frame trigger recreator
CN113315443B (en) Three-electric-system resonance point identification method, motor controller and storage medium
CN114337892B (en) Multi-equipment clock synchronization method based on low-voltage distribution network
CN113125961B (en) Battery detection method and device
CN112129993B (en) Zero crossing signal output and power line data transmission method and device
EP4293920A1 (en) Transceiver
CN117811449A (en) Motor driving method, motor model generating method and electronic device
CN106940432B (en) Method and system for testing instrument connectivity

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant