CN104122440B - Crystal oscillator frequency drift cognitive method - Google Patents
Crystal oscillator frequency drift cognitive method Download PDFInfo
- Publication number
- CN104122440B CN104122440B CN201410354488.8A CN201410354488A CN104122440B CN 104122440 B CN104122440 B CN 104122440B CN 201410354488 A CN201410354488 A CN 201410354488A CN 104122440 B CN104122440 B CN 104122440B
- Authority
- CN
- China
- Prior art keywords
- frequency drift
- crystal oscillator
- drift value
- fitting function
- class
- 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
Links
Landscapes
- Oscillators With Electromechanical Resonators (AREA)
Abstract
The invention discloses a kind of crystal oscillator frequency drift cognitive method.Comprising: measure the crystal oscillator static phase noise when static state；Set vibration environment condition, the crystal oscillator of vibration under vibration environment condition is taken multiple measurements, obtains organizing vibration phase noise more；Carry out difference operation by often organizing vibration phase noise with static phase noise, obtain many class frequencys drift value；Every class frequency drift value is carried out Screening Treatment, to reject the frequency drift value more than predetermined threshold；Every class frequency drift value is carried out Screening Treatment, to reject the frequency drift value more than predetermined threshold；Obtain and remain each class frequency drift value, according to default response surface model and every class frequency drift value, the fitting function that previous group frequency drift value is corresponding is modified, obtains final fitting function, and determine the BURNING RATE ACCELERATION SENSITIVITY direction of final fitting function.The present invention can improve the accuracy of crystal oscillator frequency drift perception, thus excavates the vibration environment condition Influencing Mechanism to crystal oscillator electrical property.
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 particularly important reference frequency, crystal oscillator (hereinafter referred to as crystal oscillator) is usually mounted on various high-speed aircraft, and high-speed aircraft is the most all operated in extremely under rugged environment, accordingly, crystal oscillator is just highly prone to the impact of vibration environment, therefore, the crystal oscillator that high-speed aircraft is used for degree of stability in vibration environment has the harshest requirement.
In order to meet this requirement, it is necessary to the electrical property of crystal oscillator is tested, to grasp the degree of stability of crystal oscillator, it is to avoid the occurrence frequency drift when reality is applied, thus cause serious consequence.Therefore, study the vibration environment condition Influencing Mechanism to crystal oscillator electrical property, become current urgent needs.
The method measuring crystal acceleration sensitivity at present mainly has 2g to invert measurement method, vibration experiment measurement method etc., the most generally takes to survey the method the most simply averaged more and calculates BURNING RATE ACCELERATION SENSITIVITY.It is true that affected by external environment disturbance (mechanical disturbance and electromagnetic distu), test data dithering amplitude is relatively big, has directly influenced the accuracy of crystal oscillator frequency drift perception, it is difficult to excavate the vibration environment condition Influencing Mechanism to crystal oscillator electrical property.
Summary of the invention
The technical problem to be solved is: for the problem of above-mentioned existence, it is provided that a kind of crystal oscillator frequency drift cognitive method, it is possible to increase the accuracy of crystal oscillator frequency drift perception, thus excavates the vibration environment condition Influencing Mechanism to crystal oscillator electrical property.
The technical solution used in the present invention is to provide a kind of crystal oscillator frequency drift cognitive method, including: measure the crystal oscillator static phase noise when static state；Set vibration environment condition, the described crystal oscillator of vibration under described vibration environment condition is taken multiple measurements, obtains organizing vibration phase noise more；Carry out difference operation by often organizing vibration phase noise with described static phase noise, obtain many class frequencys drift value；Every class frequency drift value is carried out Screening Treatment, to reject the frequency drift value more than predetermined threshold；A class frequency drift value in many class frequencys drift value is fitted, obtains the fitting function of correspondence；Obtain and remain each class frequency drift value, according to default response surface model and every class frequency drift value, the fitting function that previous group frequency drift value is corresponding is modified, obtains 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 fitted, after the step of the fitting function obtaining correspondence, before the described step obtaining each class frequency drift value of residue, described crystal oscillator frequency drift cognitive method also includes: be smoothed described fitting function, to remove the catastrophe point on described fitting function.
Preferably, described crystal oscillator frequency drift cognitive method also includes: the final fitting function of described crystal oscillator is graphically changed expression.
In sum, owing to have employed technique scheme, the invention has the beneficial effects as follows:
1, under vibration the frequency shift (FS) Value Data obtained can be carried out real-time filtering noise reduction process；
2, processed by match choosing, it is possible to the frequency shift (FS) Value Data that automatic rejection degree of association is relatively low；
3, can the impact on crystal oscillator electrical property of the intensity of variation that various vibration environment condition is given of accurate quantification.
Accompanying drawing explanation
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 of the present invention drift cognitive method.
Fig. 2 is the three-dimensional coordinate schematic diagram of the fitting function obtained in crystal oscillator frequency of the present invention drift cognitive method.
Fig. 3 is the three-dimensional cloud atlas of the final fitting function obtained in crystal oscillator frequency of the present invention drift cognitive method.
Detailed description of the invention
All features disclosed in this specification, or disclosed all methods or during step, in addition to mutually exclusive feature and/or step, all can combine by any way.
Any feature disclosed in this specification, unless specifically stated otherwise, all can be by other equivalences or there is the alternative features of similar purpose replaced.I.e., unless specifically stated otherwise, an example during each feature is a series of equivalence or similar characteristics.
As it is shown in figure 1, be the schematic flow sheet of crystal oscillator frequency of the present invention drift cognitive method.Crystal oscillator frequency drift cognitive method comprises the following steps:
S1: measure the crystal oscillator static phase noise when static state.
Wherein, the crystal oscillator phase noise when static state should be zero in the ideal case.But in most cases, owing to crystal oscillator is affected by self-condition or external condition, even if also can there is certain phase noise when static state.Measure after static phase noise, static phase noise record is got off.
S2: set vibration environment condition, takes multiple measurements the crystal oscillator of vibration under vibration environment condition, obtains organizing vibration phase noise more.
Wherein, vibration environment condition includes temperature, electromagnetic distu amount, mechanical disturbance amount etc..Owing to when measuring, the impact of each moment crystal oscillator environmental condition vibrated is different, it is therefore desirable to take multiple measurements, measures every time and obtain one group of vibration phase noise.It is noted that repetitive measurement herein refers to that multiple measurement process, each measurement process comprise multiple pendulous frequency.
S3: carry out difference operation by often organizing vibration phase noise with static phase noise, obtain many class frequencys drift value.
Wherein, one group of vibration phase noise comprises multiple vibration phase noise, and with static phase noise, the multiple vibration phase noises in one group of vibration phase noise are carried out difference operation respectively, obtains comprising a class frequency drift value of multiple frequency drift value.
S4: every class frequency drift value is carried out Screening Treatment, to reject the frequency drift value more than predetermined threshold.
Wherein, the intermittent sudden change such as external electromagnetic disturbance can deteriorate vibration phase noise, thus deteriorates frequency drift value.Predetermined threshold can set as required, and the frequency drift value below predetermined threshold may be considered normal frequency drift value, is then the frequency drift value deteriorated more than the frequency drift value of predetermined threshold.
S5: be fitted the class frequency drift value in many class frequencys drift value, obtains the fitting function of correspondence.
Wherein, fitting function characterizes the functional relation between vibration environment condition and the crystal oscillator frequency drift value deterioration degree of current setting.
S6: obtain and remain each class frequency drift value, is modified the fitting function that previous group frequency drift value is corresponding according to default response surface model and every class frequency drift value, obtains final fitting function, and determine the BURNING RATE ACCELERATION SENSITIVITY direction of final fitting function.
Wherein, the response surface model preset is for revising further fitting function.For example, if fitting function is y=ax, when obtaining next class frequency drift value, response surface model is to deduce on y=ax, and deduction result is y=ax2+bx；When obtaining next class frequency drift value again, response surface model is at y=ax2Deducing on+bx, deduction result is y=ax2+bx+c.By that analogy, during until obtaining last class frequency drift value, deducing result is y=ax2+bx+c(x-1)2.As a rule, response surface model is the function of below the cubic term of x.Revise every time and all will obtain a fitting function, as shown in Figure 3；In Fig. 2, each bar line segment in three-dimensional system of coordinate is each bar fitting function and is tried to achieve.The direction vector of each bar line segment in Fig. 2 is the BURNING RATE ACCELERATION SENSITIVITY direction of final fitting function.
S7: the final fitting function of crystal oscillator is graphically changed expression.
Wherein, the measurement result in traditional BURNING RATE ACCELERATION SENSITIVITY direction is generally data form, does not affects the inherence of crystal oscillator frequency drift value deterioration degree from vibration environment condition and sends the vibration operating characteristic expressing crystal oscillator.By final fitting function is graphically changed table can be cheer and bright express the vibration environment condition inherent impact on crystal oscillator frequency drift value deterioration degree, the result of graph-based is as it is shown on figure 3, the curved surface in three-dimensional system of coordinate is the response surface.
In the present embodiment, after step s 5, before step S6, crystal oscillator frequency drift cognitive method also includes: be smoothed fitting function, to remove the catastrophe point on fitting function.Wherein, frequency drift value can remove the frequency drift value of deterioration after Screening Treatment, but there may be the catastrophe points such as sharp-pointed or burr on the fitting function after matching, can revise fitting function further after catastrophe point being removed.
By the way, the crystal oscillator frequency drift cognitive method of the present invention is fitted by frequency drift value is carried out Screening Treatment again, such that it is able to improve the accuracy of crystal oscillator frequency drift perception, and by the response surface model preset and every class frequency drift value, the fitting function that previous group frequency drift value is corresponding is modified such that it is able to excavate the vibration environment condition Influencing Mechanism to crystal oscillator electrical property.
The invention is not limited in aforesaid detailed description of the invention.The present invention expands to any new feature disclosed in this manual or any new combination, and the arbitrary new method that discloses or the step of process or any new combination.
Claims (3)
1. a crystal oscillator frequency drift cognitive method, it is characterised in that including:
Measure the crystal oscillator static phase noise when static state；
Set vibration environment condition, the described crystal oscillator of vibration under described vibration environment condition is repeatedly surveyed
Amount, obtains organizing vibration phase noise more；
Carry out difference operation by often organizing vibration phase noise with described static phase noise, obtain the drift of many class frequencys
Shifting value；
Every class frequency drift value is carried out Screening Treatment, to reject the frequency drift value more than predetermined threshold；
A class frequency drift value in many class frequencys drift value is fitted, obtains the fitting function of correspondence；
Obtain and remain each class frequency drift value, according to default response surface model and every class frequency drift value pair
The fitting function that previous group frequency drift value is corresponding is modified, and obtains final fitting function, and determines described
The BURNING RATE ACCELERATION SENSITIVITY direction of final fitting function.
Crystal oscillator frequency the most according to claim 1 drift cognitive method, it is characterised in that described
A class frequency drift value in many class frequencys drift value is fitted, obtains the step of the fitting function of correspondence
Afterwards, before the described step obtaining each class frequency drift value of residue, described crystal oscillator frequency drift cognitive method
Also include:
Described fitting function is smoothed, to remove the catastrophe point on described fitting function.
Crystal oscillator frequency the most according to claim 1 drift cognitive method, it is characterised in that described crystalline substance
Vibration frequency drift cognitive method also includes:
The final fitting function of described crystal oscillator is graphically changed expression.
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 CN104122440A (en) | 2014-10-29 |
| CN104122440B true 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) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104408296B (en) * | 2014-11-14 | 2017-09-22 | 中国电子科技集团公司第二十九研究所 | Crystal oscillator module regulating allocation device and crystal oscillator module environmental suitability adjusting method |
| CN107219395B (en) * | 2017-05-23 | 2019-06-18 | 郑州云海信息技术有限公司 | A kind of detection method of RTC crystal oscillator |
| CN110376446B (en) * | 2019-07-17 | 2021-06-08 | 中国电子科技集团公司第二十九研究所 | High-robustness active noise reduction method for crystal oscillator module |
| CN110808712B (en) * | 2019-10-21 | 2023-04-18 | 中国电子科技集团公司第二十九研究所 | Intelligent crystal oscillator phase noise comprehensive compensation method |
| CN114884479B (en) * | 2022-07-11 | 2022-09-09 | 成都优弗科技有限公司 | Real-time correction system for crystal oscillator vibration |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE60327633D1 (en) * | 2002-05-28 | 2009-06-25 | Capital Formation Inc | Attachment structures for crystal resonators with low acceleration sensitivity |
| US7046584B2 (en) * | 2003-07-09 | 2006-05-16 | Precision Drilling Technology Services Group Inc. | Compensated ensemble crystal oscillator for use in a well borehole system |
-
2014
- 2014-07-24 CN CN201410354488.8A patent/CN104122440B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
|---|
| Impact of Mechanical Vibrations on Laser Stability;Maxim Kuschnerov 等;《IEEE Photonics Technology Letters》;20100801;第22卷(第15期);1114-1116 * |
| 加速度对晶体振荡器的影响及补偿技术研究;纪龙蛰 等;《宇航计测技术》;20120430;第32卷(第2期);37-42 * |
| 抗振晶体振荡器相位噪声测试方法的对比研究;陈中平 等;《电子科技》;20140515;第27卷(第5期);36-38 * |
| 用响应面函数分析煤巷的稳定性;郑全九;《煤矿开采》;20060830;第11卷(第4期);7-9 * |
| 石英晶体振荡器加速度效应及敏感度的测试;唐钱 等;《测试技术学报》;20120630;第26卷(第3期);200-206 * |
| 石英晶体振荡器加速度敏感度综述;夏顺冬;《现代雷达》;19890630;13-26 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104122440A (en) | 2014-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104122440B (en) | Crystal oscillator frequency drift cognitive method | |
| CN109382822B (en) | Simulation device and simulation method for robot system | |
| US8140911B2 (en) | Dynamic software tracing | |
| US20220075919A1 (en) | Systems and methods for co-simulation | |
| RU2011102604A (en) | DEVICE AND METHOD OF NAVIGATION | |
| EP2897001A3 (en) | Multi-level/multi-threshold/multi-persistency GPS/GNSS atomic clock monitoring | |
| US11137321B2 (en) | Failure diagnosis system | |
| US10564658B2 (en) | Servo control apparatus, servo control method and computer-readable recording medium | |
| CN104702214B (en) | A kind of method of crystal oscillator frequency compensation | |
| RU2013119639A (en) | MODELING OF KARST FORMATION | |
| CN106557073A (en) | The Servocontrol device of the function with the characteristic for determining learning controller | |
| EP3244224A1 (en) | Integrated system and method for testing system timing margin | |
| WO2015145195A1 (en) | Determining location and depth of subsurface magnetic sources | |
| 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 | |
| US20080100360A1 (en) | Pulsed local clock buffer (lcb) characterization ring oscillator | |
| CN104734824B (en) | The method of Chaotic Synchronous detection weak periodic signal under strong noise background | |
| CN113064674B (en) | Method and device for expanding state machine logic, storage medium and electronic device | |
| CN104880713A (en) | GBAS integrity assessment and processing method and system | |
| US8098249B2 (en) | Apparatus and method for automatic scaling of tick marks | |
| CN108643895B (en) | Resistivity logging response calculation method and device for gradient stratum | |
| CN103401534B (en) | Based on the method for the final voltage fast prediction of RC circuit discharge and recharge | |
| JP2015170091A (en) | Processing conditions calculation apparatus, method, program, processing apparatus and measurement unit | |
| EP3226181A1 (en) | Apparatus, method, and program for supporting development of production line | |
| EP3220294A1 (en) | Calculation of electromagnetic force distribution, and device for calculating electromagnetic force distribution |
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 |