CN106777539A - A kind of method of Automatic-searching solid propellant rocket frequency characteristic time of integration point - Google Patents
A kind of method of Automatic-searching solid propellant rocket frequency characteristic time of integration point Download PDFInfo
- Publication number
- CN106777539A CN106777539A CN201611059097.9A CN201611059097A CN106777539A CN 106777539 A CN106777539 A CN 106777539A CN 201611059097 A CN201611059097 A CN 201611059097A CN 106777539 A CN106777539 A CN 106777539A
- Authority
- CN
- China
- Prior art keywords
- time
- frequency values
- frequency
- integration
- data point
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Testing Relating To Insulation (AREA)
- Tests Of Electronic Circuits (AREA)
Abstract
The present invention proposes a kind of method of Automatic-searching solid propellant rocket frequency characteristic time of integration point, the leading zero's time of each sine wave signal is found first, further according to frequency dot sequency and each Frequency point to deserved sine wave number, take the half (being rounded downwards when not for integer) of waveform number under each frequency values, the second half section is selected as integration segment, integration segment corresponding time started and end time are frequency characteristic time of integration point, can so find all frequency characteristic time of integration points.Using this method Automatic-searching frequency characteristic time of integration point, only need to sequentially be input into all frequency values and its corresponding waveform number, and the instruction sine wave signal maximum that collects and minimum value, command signal time started and end time, all frequency characteristic time of integration points can be found according to method operation, time was fallen sharply by 2 seconds by 30 minutes, data-handling efficiency is substantially increased, labour intensity is reduced, also causes that the possibility that artificial treatment is malfunctioned due to fatigue no longer occurs.
Description
Technical field
The invention belongs to solid propellant rocket technical field of measurement and test, specially a kind of Automatic-searching solid propellant rocket
The method of frequency characteristic time of integration point.
Background technology
In solid propellant rocket test process, carry out needing to ask for jet pipe using test data when frequency characteristic is processed
The resonant frequency with engine is swung, technical method usual at present is：Obtain actuator displacement, actuator pressure difference, jet pipe pendulum
The amplitude-frequency and phase-frequency characteristic of angle correspondence input actuator command signal, draw log-frequency characteristic curve, and then obtain jet pipe pendulum
The resonant frequency of dynamic and engine.
Actuator command signal is the sine wave signal of a succession of different frequency, for the angle of cant signal of synchronous acquisition
And command signal, it is thus necessary to determine that frequency characteristic time of integration point, and be all at present to use hand in solid propellant rocket test
Dynamic method：The instruction sine wave signal figure of collection is drawn, each Frequency point later half number sine wave is found out manually as integration
Section, integration segment corresponding time started and end time are frequency characteristic time of integration point.
The content of the invention
The technical problem to be solved
, to find frequency characteristic time of integration point manually, when running into, frequency points are more for prior art, when waveform quantity is big, hand
It is time consuming that dynamic method becomes, and treatment once just takes half an hour, and error probability is also increased.The present invention is sought to using automation
Algorithm, finds all frequency characteristic time of integration points, improves efficiency.
Technical scheme
The thinking of solve problem of the present invention is：
The leading zero's time of each sine wave signal is found first, it is corresponding with each Frequency point further according to frequency dot sequency
Sine wave number is obtained, the half (being rounded downwards when not for integer) of waveform number under each frequency values is taken, the second half section is selected as product
Segmentation, integration segment corresponding time started and end time are frequency characteristic time of integration point, can so find all frequencies
Rate characteristic time of integration point.
According to above-mentioned thinking, the technical scheme is that：
A kind of method of the Automatic-searching solid propellant rocket frequency characteristic time of integration point, it is characterised in that：Bag
Include following steps：
Step 1：During solid propellant rocket testing experiment, with n points/second acquisition instructions signal, the finger for collecting
Sine wave signal Amplitude maxima is made for max, minimum value is min, the command signal time started is starttime, end time
It is endtime；And according to the command signal of input, all frequency values of determine instruction signal, and the corresponding waveform of each frequency
Number；Command signal peak frequency is mHz；
Step 2：From command signal time started to end time, sequentially find and all meet any in following two conditions
One data point of condition, and obtain the data point corresponding time：
Condition 1：The data value of data point is more than 0.4max+0.6min, and n/5m data point being close to before it
Data value be both less than 0.4max+0.6min；
Condition 2：The data value of data point is less than 0.6max+0.4min, and n/5m data point being close to after which
Data value be both greater than 0.6max+0.4min；
Step 3：The sequencing of the data point obtained according to step 2, since second data point for obtaining, calculates the
The average value of 2k time corresponding with the 2k+1 data point, k=1,2 ..., at the end of obtaining all command signal sine waves
Between；
Step 4：According to the sequencing of each frequency values of command signal, the initial time of each frequency values integration segment is determined
And the end time：
For t-th frequency values, by the corresponding waveform number of the frequency values and the waveform number phase of all frequency values before
Plus, obtain the corresponding waveform number accumulated value of t-th frequency values；
The half of the corresponding waveform number of t-th frequency values is subtracted using the corresponding waveform number accumulated value of t-th frequency values
Value, obtains t-th sine wave onset wave serial number of frequency values integration segment, and all command signals obtained according to step 3
The end time of sine wave, between determining at the beginning of t-th sine wave onset wave of frequency values integration segment；
Added up using the corresponding waveform number of t-th frequency values and be worth to the sine wave of t-th frequency values integration segment and terminate
Waveform sequence number, and all command signal sine waves obtained according to step 3 end time, determine t-th frequency values integration
The sine wave of section terminates the end time of waveform.
Beneficial effect
Using this method Automatic-searching frequency characteristic time of integration point, all frequency values and its corresponding only need to be sequentially input into
Waveform number, and collect instruction sine wave signal maximum max and minimum value min, the command signal time started and
End time, all frequency characteristic time of integration points can be at once found according to automatic algorithms operation program, the time was by 30 minutes
Once fall sharply by 2 seconds, substantially increase data-handling efficiency, reduce labour intensity, also cause artificial treatment due to fatigue
The possibility of error no longer occurs.
Additional aspect of the invention and advantage will be set forth in part in the description, and will partly become from the following description
Obtain substantially, or recognized by practice of the invention.
Brief description of the drawings
Of the invention above-mentioned and/or additional aspect and advantage will become from description of the accompanying drawings below to embodiment is combined
Substantially and be readily appreciated that, wherein：
Fig. 1：Command signal sine waveform schematic diagram；
Fig. 2：Command signal sine wave schematic diagram in embodiment；
Fig. 3：All instruction sine wave end times in embodiment；
Fig. 4：The corresponding waveform number accumulated value of each frequency values in embodiment；
Fig. 5：In embodiment at the beginning of each frequency values integration segment between and the end time.
Specific embodiment
Embodiments of the invention are described below in detail, the embodiment is exemplary, it is intended to for explaining the present invention, and
It is not considered as limiting the invention.
The purpose of the present invention be can according to collection command signal sine wave amplitude, starting and the end time, and
When the frequency dot sequency and each Frequency point of command signal sine wave find out frequency characteristic integration to deserved sine wave number automatically
Between point, principle is as follows：The command signal of collection is the sine wave signal of a succession of different frequency, and signal has certain interference, but
Its amplitude is identical, the ascent stage of a sine wave find out first more than -20% amplitude data point and last
Less than a data point for+20% amplitude, same ascent stage first is more than data point correspondence time of -20% amplitude and tight
With thereafter last less than the average value of the data point correspondence time of+20% amplitude as next sine wave signal
Initial time, as shown in Figure 1；Further according to frequency dot sequency and each Frequency point to deserved sine wave number, each frequency is taken
The half (being rounded downwards when not for integer) of the lower waveform number of value, selects the second half section as integration segment, at the beginning of integration segment is corresponding
Between and the end time be frequency characteristic time of integration point.
One ascent stage of sine wave, first data point more than -20% amplitude was so looked for：Current data point data value
More than -20% amplitude, and the several data point datas being close to before it are both less than -20% amplitude；One sine wave it is upper
Last is so looked for less than the data point of+20% amplitude to rise section：Current data point data value is less than+20% amplitude, and
And the several data point datas being close to after which are both greater than+20% amplitude.
It is of the invention to concretely comprise the following steps based on above-mentioned principle：
Step 1：During solid propellant rocket testing experiment, with n points/second acquisition instructions signal, the finger for collecting
Sine wave signal Amplitude maxima is made for max, minimum value is min, the command signal time started is starttime, end time
It is endtime；And according to the command signal of input, all frequency values of determine instruction signal, and the corresponding waveform of each frequency
Number；Command signal peak frequency is mHz；
Step 2：From command signal time started to end time, sequentially find and all meet any in following two conditions
One data point of condition, and obtain the data point corresponding time：
Condition 1：The data value of data point is more than 0.4max+0.6min, and n/5m data point being close to before it
Data value be both less than 0.4max+0.6min；
Condition 2：The data value of data point is less than 0.6max+0.4min, and n/5m data point being close to after which
Data value be both greater than 0.6max+0.4min；
Step 3：The sequencing of the data point obtained according to step 2, since second data point for obtaining, calculates the
The average value of 2k time corresponding with the 2k+1 data point, k=1,2 ..., at the end of obtaining all command signal sine waves
Between, namely at the beginning of next sine wave between；
Step 4：According to the sequencing of each frequency values of command signal, the initial time of each frequency values integration segment is determined
And the end time：
For t-th frequency values, by the corresponding waveform number of the frequency values and the waveform number phase of all frequency values before
Plus, obtain the corresponding waveform number accumulated value of t-th frequency values；
The half of the corresponding waveform number of t-th frequency values is subtracted using the corresponding waveform number accumulated value of t-th frequency values
Value (is rounded downwards) when not for integer, obtains t-th sine wave onset wave serial number of frequency values integration segment, and according to step
The end time of the 3 all command signal sine waves for obtaining, determine opening for t-th sine wave onset wave of frequency values integration segment
Time beginning；
Added up using the corresponding waveform number of t-th frequency values and be worth to the sine wave of t-th frequency values integration segment and terminate
Waveform sequence number, and all command signal sine waves obtained according to step 3 end time, determine t-th frequency values integration
The sine wave of section terminates the end time of waveform.
With 5000 points/second acquisition instructions signal in the present embodiment, command signal peak frequency is 50Hz, the instruction for collecting
Signal sine wave maximum is 0.235, and minimum value is -0.2243, and the command signal time started is 21.6s, and the end time is
32.1s, as shown in Figure 2.Command signal is included successively：0.5,1,1.592,2.387,3.183 five frequency values, it is corresponding in turn to
Waveform number be：2,2,3,3,3.
From command signal time started to end time, sequentially find and all meet any one in following two conditions
The data point of part, and obtain the data point corresponding time：
Condition 1：The data value of data point is more than -0.04058, and 20 data values of data point being close to before it
Both less than -0.04058；
Condition 2：The data value of data point is less than 0.05128, and 20 data values of data point being close to after which
Both greater than 0.05128.
According to the sequencing of data point obtained above, since second data point for obtaining, calculate 2k and
The 2k+1 average value of data point correspondence time, k=1,2 ..., the end time of all command signal sine waves is obtained,
Between at the beginning of i.e. next sine wave；As shown in Figure 3.
According to the sequencing of each frequency values of command signal, initial time and the end of each frequency values integration segment are determined
Time：
For t-th frequency values, by the corresponding waveform number of the frequency values and the waveform number phase of all frequency values before
Plus, obtain the corresponding waveform number accumulated value of t-th frequency values；Accumulated value is as shown in Figure 4.
The half of the corresponding waveform number of t-th frequency values is subtracted using the corresponding waveform number accumulated value of t-th frequency values
Value (is rounded downwards) when not for integer, obtains t-th sine wave onset wave serial number of frequency values integration segment, and according to step
The end time of the 3 all command signal sine waves for obtaining, determine opening for t-th sine wave onset wave of frequency values integration segment
Time beginning；Added up using the corresponding waveform number of t-th frequency values and be worth to the sine wave of t-th frequency values integration segment and terminate ripple
Shape serial number, and all command signal sine waves obtained according to step 3 end time, determine t-th frequency values integration segment
Sine wave terminate end time of waveform.Result is as shown in Figure 5.
Although embodiments of the invention have been shown and described above, it is to be understood that above-described embodiment is example
Property, it is impossible to limitation of the present invention is interpreted as, one of ordinary skill in the art is not departing from principle of the invention and objective
In the case of above-described embodiment can be changed within the scope of the invention, change, replace and modification.
Claims (1)
1. a kind of method of Automatic-searching solid propellant rocket frequency characteristic time of integration point, it is characterised in that：Including following
Step：
Step 1：During solid propellant rocket testing experiment, with n points/second acquisition instructions signal, the instruction for collecting is just
String ripple signal amplitude maximum is max, and minimum value is min, and the command signal time started is starttime, and the end time is
endtime；And according to the corresponding waveform of the command signal of input, all frequency values of determine instruction signal, and each frequency
Number；Command signal peak frequency is mHz；
Step 2：From command signal time started to end time, sequentially find and all meet any one in following two conditions
The data point of condition, and obtain the data point corresponding time：
Condition 1：The data value of data point be more than 0.4max+0.6min, and n/5m data point being close to before it number
0.4max+0.6min is both less than according to value；
Condition 2：The data value of data point be less than 0.6max+0.4min, and n/5m data point being close to after which number
0.6max+0.4min is both greater than according to value；
Step 3：The sequencing of the data point obtained according to step 2, since second data point for obtaining, calculates 2k
The average value of time corresponding with the 2k+1 data point, k=1,2 ..., obtain the end time of all command signal sine waves；
Step 4：According to the sequencing of each frequency values of command signal, the initial time and knot of each frequency values integration segment are determined
The beam time：
For t-th frequency values, the corresponding waveform number of the frequency values is added with the waveform number of all frequency values before, obtained
To the corresponding waveform number accumulated value of t-th frequency values；
The half value of the corresponding waveform number of t-th frequency values is subtracted using the corresponding waveform number accumulated value of t-th frequency values, is obtained
To t-th sine wave onset wave serial number of frequency values integration segment, and all command signal sine waves obtained according to step 3
End time, between determining at the beginning of t-th sine wave onset wave of frequency values integration segment；
Added up using the corresponding waveform number of t-th frequency values and be worth to the sine wave of t-th frequency values integration segment and terminate waveform
Serial number, and all command signal sine waves obtained according to step 3 end time, determine t-th frequency values integration segment
Sine wave terminates the end time of waveform.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611059097.9A CN106777539B (en) | 2016-11-28 | 2016-11-28 | A kind of method of Automatic-searching solid propellant rocket frequency characteristic time of integration point |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611059097.9A CN106777539B (en) | 2016-11-28 | 2016-11-28 | A kind of method of Automatic-searching solid propellant rocket frequency characteristic time of integration point |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106777539A true CN106777539A (en) | 2017-05-31 |
| CN106777539B CN106777539B (en) | 2019-09-03 |
Family
ID=58913131
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611059097.9A Active CN106777539B (en) | 2016-11-28 | 2016-11-28 | A kind of method of Automatic-searching solid propellant rocket frequency characteristic time of integration point |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106777539B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108168897A (en) * | 2017-12-15 | 2018-06-15 | 西安航天动力测控技术研究所 | A kind of method of solid engines experiment 1553B servo instruction automatic segmentations |
| CN110082069A (en) * | 2019-05-06 | 2019-08-02 | 吉洁阳 | A kind of engine optical fiber test apparatus |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2587321A1 (en) * | 2011-10-25 | 2013-05-01 | Siemens Aktiengesellschaft | Wave filtering using differential evolution for dynamic positioning systems |
| CN103411629A (en) * | 2013-08-15 | 2013-11-27 | 北京航天控制仪器研究所 | Rocket sled test exterior data fusion method |
| CN103455645A (en) * | 2012-05-31 | 2013-12-18 | 北京宇航系统工程研究所 | Overall-modal extraction method |
-
2016
- 2016-11-28 CN CN201611059097.9A patent/CN106777539B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2587321A1 (en) * | 2011-10-25 | 2013-05-01 | Siemens Aktiengesellschaft | Wave filtering using differential evolution for dynamic positioning systems |
| CN103455645A (en) * | 2012-05-31 | 2013-12-18 | 北京宇航系统工程研究所 | Overall-modal extraction method |
| CN103411629A (en) * | 2013-08-15 | 2013-11-27 | 北京航天控制仪器研究所 | Rocket sled test exterior data fusion method |
Non-Patent Citations (2)
| Title |
|---|
| 张敏等: "带线性隔离的微弱信号调理电路设计及应用", 《工业控制计算机》 * |
| 杨云飞等: "运载火箭摇摆发动机与全检动力学特性耦合关系研究", 《宇航学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108168897A (en) * | 2017-12-15 | 2018-06-15 | 西安航天动力测控技术研究所 | A kind of method of solid engines experiment 1553B servo instruction automatic segmentations |
| CN110082069A (en) * | 2019-05-06 | 2019-08-02 | 吉洁阳 | A kind of engine optical fiber test apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106777539B (en) | 2019-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109781412B (en) | Rolling bearing self-adaptive resonance demodulation method based on EEMD | |
| CN102680860B (en) | Automatic fault-point locating method for traveling-wave based fault location of high-voltage electric power lines | |
| CN105547698A (en) | Fault diagnosis method and apparatus for rolling bearing | |
| CN106777539B (en) | A kind of method of Automatic-searching solid propellant rocket frequency characteristic time of integration point | |
| JP2015075821A (en) | State diagnostic method and state diagnostic device | |
| WO2016026558A1 (en) | A method for assessing the condition of rotating machinery connected to an electric motor | |
| CN103424183B (en) | Method for eliminating abnormal interference on detection for mechanical vibration signals | |
| CN103584854A (en) | Extraction method of electrocardiosignal R waves | |
| CN108470018A (en) | Smoothing method and device based on the intrinsic mode functions that empirical mode decomposition decomposes | |
| He et al. | Feature extraction of acoustic signals based on complex Morlet wavelet | |
| CN110653661A (en) | Cutter state monitoring and identifying method based on signal fusion and multi-fractal spectrum algorithm | |
| CN106198013B (en) | A kind of envelope Analysis Method based on empirical mode decomposition filtering | |
| CN110530507B (en) | Edge calculation method, medium, and system for monitoring rotating device | |
| CN103310113A (en) | Universal blood glucose prediction method based on frequency band separation and data modeling | |
| CN102566555A (en) | Major appliance work state monitoring method based on pattern recognition | |
| CN106894950B (en) | A kind of power characteristic verification method based on Wind turbines drive characteristic | |
| CN111075704B (en) | Frequency hopping bandwidth detection system and intelligent algorithm for frequency-conversion compressor of air conditioner in data machine room | |
| CN114486260B (en) | Bearing fault diagnosis method based on self-adaptive variational modal decomposition | |
| CN106682791B (en) | Energy-conserving efficient refuse treatment system | |
| RU2705070C1 (en) | Method of aircraft with turbojet engine type identification in pulse-doppler radar station | |
| CN117390972B (en) | Electronic cigarette atomizer service life prediction method | |
| CN109301842A (en) | A kind of wind power plant sub-synchronous oscillation cutting method based on negative damping contribution | |
| CN110920931B (en) | Static load and dynamic load separation method for flight test of rotating part | |
| Gryllias et al. | A peak energy criterion (pe) for the selection of resonance bands in complex shifted morlet wavelet (csmw) based demodulation of defective rolling element bearings vibration response | |
| CN111523231A (en) | Subsynchronous oscillation analysis method based on EEMD and Prony method |
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 |