CN104122440A - Method for sensing crystal oscillator frequency shift - Google Patents
Method for sensing crystal oscillator frequency shift Download PDFInfo
- Publication number
- CN104122440A CN104122440A CN201410354488.8A CN201410354488A CN104122440A CN 104122440 A CN104122440 A CN 104122440A CN 201410354488 A CN201410354488 A CN 201410354488A CN 104122440 A CN104122440 A CN 104122440A
- Authority
- CN
- China
- Prior art keywords
- crystal oscillator
- frequency drift
- fitting function
- drift value
- frequency shift
- 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.)
- Granted
Links
Abstract
The invention discloses a method for sensing crystal oscillator frequency shift. The method comprises the steps of measuring static phase noise of a crystal oscillator in a static state; setting vibration environment conditions, performing multiple metering on the crystal oscillator vibrating under the vibration environment conditions, obtaining a plurality of groups of vibration phase noise; performing difference calculation on each group of vibration phase noise and the static phase noise, obtaining a plurality of groups of frequency shift values; performing screening processing on each group of frequency shift values so as to remove frequency shift values larger than a preset threshold; obtaining each group of remaining frequency shift values, correcting fitting functions corresponding to a last group of frequency shift values according to a preset response surface model and each group of frequency shift values, obtaining a final fitting function, and determining the acceleration sensitivity direction of the final fitting function. The method can improve accuracy of crystal oscillator frequency shift sensing, and accordingly an influence mechanism influencing crystal oscillator electrical property under the vibration environment conditions is excavated.
Description
Technical field
The present invention relates to electronic technology field, especially relate to a kind of crystal oscillator frequency drift cognitive method.
Background technology
As very important reference frequency, crystal oscillator (hereinafter to be referred as crystal oscillator) is usually mounted on various high-speed aircrafts, and high-speed aircraft is all operated in very under rugged environment conventionally, accordingly, crystal oscillator is just very easily subject to the impact of vibration environment, therefore, to the crystal oscillator of its use, the degree of stability in vibration environment has very harsh requirement to high-speed aircraft.
In order to meet this requirement, just need to test the electrical property of crystal oscillator, to grasp the degree of stability of crystal oscillator, avoid occurrence frequency drift in the time of practical application, thereby cause serious consequence.Therefore, the mechanism that affects of research vibration environment condition on crystal oscillator electrical property, becomes current active demand.
The method of measuring at present crystal acceleration sensitivity mainly contains 2g reversion mensuration, vibration experiment mensuration etc., generally takes in the industry the methods of simply averaging several times of surveying to calculate BURNING RATE ACCELERATION SENSITIVITY more.In fact, be subject to the impact of external environment disturbance (mechanical disturbance and electromagnetic distu), test data jitter amplitude is larger, has directly had influence on the accuracy of crystal oscillator frequency drift perception, is difficult to excavate the affect mechanism of vibration environment condition on crystal oscillator electrical property.
Summary of the invention
Technical matters to be solved by this invention is: for the problem of above-mentioned existence, provide a kind of crystal oscillator frequency drift cognitive method, can improve the accuracy of crystal oscillator frequency drift perception, thereby excavate the affect mechanism of vibration environment condition on crystal oscillator electrical property.
The technical solution used in the present invention is to provide a kind of crystal oscillator frequency drift cognitive method, comprising: measure the static phase noise of crystal oscillator in the time of static state; Set vibration environment condition, the described crystal oscillator vibrating under described vibration environment condition is taken multiple measurements, obtain many group vibration phase noises; Every group of vibration phase noise and described static phase noise are carried out to difference computing, obtain many class frequencys drift value; Every class frequency drift value is carried out to Screening Treatment, to reject the frequency drift value that is greater than predetermined threshold; A class frequency drift value in many class frequencys drift value is carried out to matching, obtain corresponding fitting function; Obtain the each class frequency drift value of residue, according to default response surface model and every class frequency drift value, fitting function corresponding to last class frequency drift value revised, obtain final fitting function, and determine the BURNING RATE ACCELERATION SENSITIVITY direction of described final fitting function.
Preferably, described, a class frequency drift value in many class frequencys drift value is carried out to matching, after obtaining the step of corresponding fitting function, before the described step of obtaining the each class frequency drift value of residue, described crystal oscillator frequency drift cognitive method also comprises: described fitting function is carried out to smoothing processing, to remove the catastrophe point on described fitting function.
Preferably, described crystal oscillator frequency drift cognitive method also comprises: by the final fitting function of described crystal oscillator with graph-based.
In sum, owing to having adopted technique scheme, the invention has the beneficial effects as follows:
1, can, under vibration environment, carry out real-time filtering noise reduction process to the frequency shift (FS) Value Data obtaining;
2, by Screening Treatment, can automatic rejection and the lower frequency shift (FS) Value Data of the degree of correlation;
3, the impact of the intensity of variation that provides various vibration environment conditions that can accurate quantification on crystal oscillator electrical property.
Brief description of the drawings
Examples of the present invention will be described by way of reference to the accompanying drawings, wherein:
Fig. 1 is the schematic flow sheet of crystal oscillator frequency drift cognitive method of the present invention.
Fig. 2 is the three-dimensional coordinate schematic diagram of the fitting function that obtains in crystal oscillator frequency of the present invention drift cognitive method.
Fig. 3 is the three-dimensional cloud atlas of the final fitting function that obtains in crystal oscillator frequency of the present invention drift cognitive method.
Embodiment
Disclosed all features in this instructions, or step in disclosed all methods or process, except mutually exclusive feature and/or step, all can combine by any way.
Disclosed arbitrary feature in this instructions, unless narration especially all can be replaced by other equivalences or the alternative features with similar object.,, unless narration especially, each feature is an example in a series of equivalences or similar characteristics.
As shown in Figure 1, be the schematic flow sheet of crystal oscillator frequency drift cognitive method of the present invention.Crystal oscillator frequency drift cognitive method comprises the following steps:
S1: measure the static phase noise of crystal oscillator in the time of static state.
Wherein, the phase noise of crystal oscillator in the time of static state should be zero in the ideal case.But in most cases, because being subject to self-condition or external condition, crystal oscillator affects, even if also can there is certain phase noise in the time of static state.Measure after static phase noise, static phase noise files is got off.
S2: set vibration environment condition, the crystal oscillator vibrating under vibration environment condition is taken multiple measurements, obtain many group vibration phase noises.
Wherein, vibration environment condition comprises temperature, electromagnetic distu amount, mechanical disturbance amount etc.Due in the time measuring, each moment crystal oscillator is subject to the difference that affects of vibration environment condition, therefore needs to take multiple measurements, and measures one group of vibration phase noise at every turn.It is noted that repeatedly measurement herein refers to multiple measuring processes, each measuring process comprises multiple measurement number of times.
S3: every group of vibration phase noise and static phase noise are carried out to difference computing, obtain many class frequencys drift value.
Wherein, one group of vibration phase noise comprises multiple vibration phase noises, and the multiple vibration phase noises in one group of vibration phase noise are carried out to difference computing respectively at static phase noise, obtains a class frequency drift value that comprises multiple frequency drift values.
S4: every class frequency drift value is carried out to Screening Treatment, to reject the frequency drift value that is greater than predetermined threshold.
Wherein, the intermittent sudden change such as external electromagnetic disturbance can worsen vibration phase noise, thereby worsens frequency drift value.Predetermined threshold can be set as required, and the frequency drift value below predetermined threshold can be thought normal frequency drift value, and the frequency drift value that is greater than predetermined threshold is the frequency drift value worsening.
S5: the class frequency drift value in many class frequencys drift value is carried out to matching, obtain corresponding fitting function.
Wherein, fitting function has characterized the vibration environment condition of current setting and crystal oscillator frequency drift value worsens the functional relation between Chengdu.
S6: obtain the each class frequency drift value of residue, according to default response surface model and every class frequency drift value, fitting function corresponding to last class frequency drift value revised, obtain final fitting function, and determine the BURNING RATE ACCELERATION SENSITIVITY direction of final fitting function.
Wherein, default response surface model is for further revising fitting function.For instance, if fitting function is y=ax, in the time obtaining next class frequency drift value, response surface model is to deduce on y=ax, and deduction result is y=ax
2+ bx; While obtaining next class frequency drift value, response surface model is at y=ax again
2on+bx, deduce, deduction result is y=ax
2+ bx+c.By that analogy, until while obtaining last class frequency drift value, deduction result is y=ax
2+ bx+c (x-1)
2.As a rule, the function below the cubic term that response surface model is x.Each correction all will obtain a fitting function, and as shown in Figure 2, each line segment in figure in three-dimensional system of coordinate is each fitting function.Be the BURNING RATE ACCELERATION SENSITIVITY direction of determining final fitting function according to the direction vector of each line segment in figure.
S7: by the final fitting function of crystal oscillator with graph-based.
Wherein, the measurement result of traditional BURNING RATE ACCELERATION SENSITIVITY direction is generally data form, from vibration environment condition, the inherence of crystal oscillator frequency drift value deterioration degree is not affected and to send the vibration operating characteristic of expressing crystal oscillator.By by final fitting function with graphically show can be cheer and bright express the inherence impact of vibration environment condition on crystal oscillator frequency drift value deterioration degree, as shown in Figure 3, the curved surface in three-dimensional system of coordinate is response surface to the result of graph-based.
In the present embodiment, after step S5, before step S6, crystal oscillator frequency drift cognitive method also comprises: fitting function is carried out to smoothing processing, to remove the catastrophe point on fitting function.Wherein, frequency drift value can be removed the frequency drift value of deterioration after Screening Treatment, but may have the catastrophe points such as sharp-pointed or burr on the fitting function after matching, after catastrophe point is removed, can further revise fitting function.
By the way, crystal oscillator frequency of the present invention drift cognitive method carries out Screening Treatment by frequency drift value just and carries out matching again, thereby can improve the accuracy of crystal oscillator frequency drift perception, and by default response surface model and each class frequency drift value, fitting function corresponding to last class frequency drift value revised, thereby can excavate the affect mechanism of vibration environment condition on crystal oscillator electrical property.
The present invention is not limited to aforesaid embodiment.The present invention expands to any new feature or any new combination disclosing in this manual, and the arbitrary new method disclosing or step or any new combination of process.
Claims (3)
1. a crystal oscillator frequency drift cognitive method, is characterized in that, comprising:
Measure the static phase noise of crystal oscillator in the time of static state;
Set vibration environment condition, the described crystal oscillator vibrating under described vibration environment condition is taken multiple measurements, obtain many group vibration phase noises;
Every group of vibration phase noise and described static phase noise are carried out to difference computing, obtain many class frequencys drift value;
Every class frequency drift value is carried out to Screening Treatment, to reject the frequency drift value that is greater than predetermined threshold;
A class frequency drift value in many class frequencys drift value is carried out to matching, obtain corresponding fitting function;
Obtain the each class frequency drift value of residue, according to default response surface model and every class frequency drift value, fitting function corresponding to last class frequency drift value revised, obtain final fitting function, and determine the BURNING RATE ACCELERATION SENSITIVITY direction of described final fitting function.
2. crystal oscillator frequency drift cognitive method according to claim 1, it is characterized in that, described, a class frequency drift value in many class frequencys drift value is carried out to matching, after obtaining the step of corresponding fitting function, before the described step of obtaining the each class frequency drift value of residue, described crystal oscillator frequency drift cognitive method also comprises:
Described fitting function is carried out to smoothing processing, to remove the catastrophe point on described fitting function.
3. crystal oscillator frequency drift cognitive method according to claim 1, is characterized in that, described crystal oscillator frequency drift cognitive method also comprises:
By the final fitting function of described crystal oscillator with graph-based.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410354488.8A CN104122440B (en) | 2014-07-24 | 2014-07-24 | Crystal oscillator frequency drift cognitive method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410354488.8A CN104122440B (en) | 2014-07-24 | 2014-07-24 | Crystal oscillator frequency drift cognitive method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104122440A true CN104122440A (en) | 2014-10-29 |
CN104122440B CN104122440B (en) | 2016-09-14 |
Family
ID=51767930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410354488.8A Active CN104122440B (en) | 2014-07-24 | 2014-07-24 | Crystal oscillator frequency drift cognitive method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104122440B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104408296A (en) * | 2014-11-14 | 2015-03-11 | 中国电子科技集团公司第二十九研究所 | Crystal oscillator module adjusting configuration device and crystal oscillator module environmental adaptability adjusting method |
CN107219395A (en) * | 2017-05-23 | 2017-09-29 | 郑州云海信息技术有限公司 | A kind of detection method of RTC crystal oscillators |
CN110376446A (en) * | 2019-07-17 | 2019-10-25 | 中国电子科技集团公司第二十九研究所 | Crystal oscillator module high robust active denoising method |
CN110808712A (en) * | 2019-10-21 | 2020-02-18 | 中国电子科技集团公司第二十九研究所 | Intelligent crystal oscillator phase noise comprehensive compensation method |
CN114884479A (en) * | 2022-07-11 | 2022-08-09 | 成都优弗科技有限公司 | Real-time correction system for crystal oscillator vibration |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003100876A2 (en) * | 2002-05-28 | 2003-12-04 | Vectron International | Low acceleration sensitivity mounting structures for crystal resonators |
US20050007887A1 (en) * | 2003-07-09 | 2005-01-13 | Martin Sorrells | Compensated ensemble crystal oscillator for use in a well borehole system |
CN103293376A (en) * | 2013-05-31 | 2013-09-11 | 江汉大学 | Frequency stability measuring method and device |
CN103760446A (en) * | 2014-01-24 | 2014-04-30 | 吴建堂 | Crystal oscillator performance detector |
-
2014
- 2014-07-24 CN CN201410354488.8A patent/CN104122440B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003100876A2 (en) * | 2002-05-28 | 2003-12-04 | Vectron International | Low acceleration sensitivity mounting structures for crystal resonators |
US20050007887A1 (en) * | 2003-07-09 | 2005-01-13 | Martin Sorrells | Compensated ensemble crystal oscillator for use in a well borehole system |
CN103293376A (en) * | 2013-05-31 | 2013-09-11 | 江汉大学 | Frequency stability measuring method and device |
CN103760446A (en) * | 2014-01-24 | 2014-04-30 | 吴建堂 | Crystal oscillator performance detector |
Non-Patent Citations (6)
Title |
---|
MAXIM KUSCHNEROV 等: "Impact of Mechanical Vibrations on Laser Stability", 《IEEE PHOTONICS TECHNOLOGY LETTERS》 * |
唐钱 等: "石英晶体振荡器加速度效应及敏感度的测试", 《测试技术学报》 * |
夏顺冬: "石英晶体振荡器加速度敏感度综述", 《现代雷达》 * |
纪龙蛰 等: "加速度对晶体振荡器的影响及补偿技术研究", 《宇航计测技术》 * |
郑全九: "用响应面函数分析煤巷的稳定性", 《煤矿开采》 * |
陈中平 等: "抗振晶体振荡器相位噪声测试方法的对比研究", 《电子科技》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104408296A (en) * | 2014-11-14 | 2015-03-11 | 中国电子科技集团公司第二十九研究所 | Crystal oscillator module adjusting configuration device and crystal oscillator module environmental adaptability adjusting method |
CN104408296B (en) * | 2014-11-14 | 2017-09-22 | 中国电子科技集团公司第二十九研究所 | Crystal oscillator module regulating allocation device and crystal oscillator module environmental suitability adjusting method |
CN107219395A (en) * | 2017-05-23 | 2017-09-29 | 郑州云海信息技术有限公司 | A kind of detection method of RTC crystal oscillators |
CN107219395B (en) * | 2017-05-23 | 2019-06-18 | 郑州云海信息技术有限公司 | A kind of detection method of RTC crystal oscillator |
CN110376446A (en) * | 2019-07-17 | 2019-10-25 | 中国电子科技集团公司第二十九研究所 | Crystal oscillator module high robust active denoising method |
CN110376446B (en) * | 2019-07-17 | 2021-06-08 | 中国电子科技集团公司第二十九研究所 | High-robustness active noise reduction method for crystal oscillator module |
CN110808712A (en) * | 2019-10-21 | 2020-02-18 | 中国电子科技集团公司第二十九研究所 | Intelligent crystal oscillator phase noise comprehensive compensation method |
CN110808712B (en) * | 2019-10-21 | 2023-04-18 | 中国电子科技集团公司第二十九研究所 | Intelligent crystal oscillator phase noise comprehensive compensation method |
CN114884479A (en) * | 2022-07-11 | 2022-08-09 | 成都优弗科技有限公司 | Real-time correction system for crystal oscillator vibration |
CN114884479B (en) * | 2022-07-11 | 2022-09-09 | 成都优弗科技有限公司 | Real-time correction system for crystal oscillator vibration |
Also Published As
Publication number | Publication date |
---|---|
CN104122440B (en) | 2016-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104122440A (en) | Method for sensing crystal oscillator frequency shift | |
MY156075A (en) | Full wavefield inversion using time varying filters | |
CN103322931A (en) | System and method for measuring gap and offset of point cloud | |
CN106959105B (en) | Method for calibrating compass and device | |
JP2016517582A5 (en) | ||
CN109115257B (en) | Method, device, equipment and storage medium for correcting sensor characteristic curve | |
CN104716904B (en) | A kind of method of crystal oscillator frequency compensation | |
MX2019000245A (en) | System and method for collecting operational vibration data for a mining machine. | |
CN104702214B (en) | A kind of method of crystal oscillator frequency compensation | |
CN108345953A (en) | Temperature predicting method and device | |
CN106557073A (en) | The Servocontrol device of the function with the characteristic for determining learning controller | |
CN105674943B (en) | A kind of non-linear whole Deformation Prediction method of general multiple spot | |
AU2017264981B2 (en) | Calculation method for compressed air-flow rate, calculation device thereof, and storage medium | |
RU2014127119A (en) | AIRCRAFT CONTAINING MEASURING SENSOR AND METHOD FOR DETERMINING THE FLIGHT PARAMETERS OF SUCH AIRCRAFT | |
US20210190272A1 (en) | Piping diagnostic device, piping diagnostic method, discriminant-problem optimizing device, discriminant-problem optimizing method, criteria distribution generating device, criteria distribution generating method, and storage medium | |
KR20200033727A (en) | Data processing method, data processing device, data processing system, and computer readable recording medium with data processing program thereon | |
CN105571610B (en) | A kind of method and device of setting pedometer parameter | |
EP2667206A1 (en) | Transient recovery voltage measuring device, transient recovery voltage measuring method, and transient recovery voltage measuring program | |
KR102093287B1 (en) | Method for measuring indirectly tool wear of CNC machine | |
CN102865882A (en) | Method for determining starting time of fiber-optic gyroscope | |
CN115790665A (en) | Gyro error compensation method and device, electronic equipment and storage medium | |
CN108988832B (en) | Method for detecting a delay associated with an electronic device and corresponding electronic device | |
US9746961B2 (en) | Background signal processing system and background signal processing method | |
JP2016031568A (en) | Abnormality diagnosis apparatus, abnormality diagnosis method, and abnormality diagnosis program | |
KR20150015200A (en) | Device for measuring frequency of power based on zero-crossing and method for measuring frequency of power based on zero-crossing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |