CN106370347B - A kind of Space Pyrotechnics Devices product ignition test pressure sensor drift self-adapting correction method - Google Patents
A kind of Space Pyrotechnics Devices product ignition test pressure sensor drift self-adapting correction method Download PDFInfo
- Publication number
- CN106370347B CN106370347B CN201610814157.7A CN201610814157A CN106370347B CN 106370347 B CN106370347 B CN 106370347B CN 201610814157 A CN201610814157 A CN 201610814157A CN 106370347 B CN106370347 B CN 106370347B
- Authority
- CN
- China
- Prior art keywords
- curve
- test
- pressure sensor
- matrix
- pressure
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L27/00—Testing or calibrating of apparatus for measuring fluid pressure
- G01L27/002—Calibrating, i.e. establishing true relation between transducer output value and value to be measured, zeroing, linearising or span error determination
Abstract
It is under fire the invention discloses a kind of Space Pyrotechnics Devices product and tests pressure sensor drift self-adapting correction method, belong to Space Pyrotechnics Devices product ignition testing field.The following steps are included: Space Pyrotechnics Devices product are under fire when test, strain pressure transducer acquires the pressure signal in confined explosion space, sends computer for pressure signal by Acquisition Instrument.Computer first obtains some historical measurement datas from database, and each historical data takes partial parameters feature to form a sample matrix, obtains the main characteristic component of sample matrix using PCA algorithm, this Partial Feature can characterize the drift features of strain pressure transducer.Projection matrix is used to correct the data newly measured by the projection matrix for calculating PCA dimensionality reduction, is realized and is utilized correction of the historical data to sensor drift.New measured value, which is updated to database, obtains new sample matrix simultaneously, carries out PCA dimensionality reduction to new samples matrix and obtains new projection matrix, realizes the adaptively correcting of pressure sensor drift.
Description
Technical field
The invention belongs to Space Pyrotechnics Devices product ignition testing fields, are related to a kind of based on PCA algorithm and historical measurement data point
Analysis obtains pressure sensor drift features method in test of being under fire, and is then corrected drift according to measured data.
Background technique
Priming system is a kind of element and device for referring to ignite, ignite or provide power source, it is weapons and ammunitions, artificial
One of the most key element of the igniting elements such as satellite, spaceship, rocket and civil blast, powder power device.Firer
All too busy to get away experiment and test are developed and are produced in enormous quantities in product design.Up to the present, the design of priming system is largely
Experience is still relied on, in order to evaluate the superiority and inferiority of priming system performance, in addition to necessary quality inspection, finally also to lean on experiment and test
As a result it determines.By experiment and test, obtain various test parameters and data, with this determine it is designed, develop and production
Whether properties of product reach defined indices.
P-t curve is an index critically important in fire-working article test experiment, can be calculated according to P-t curve and start to build pressure
Time, surge pressure, surge pressure time, maximum climbing speed, quasi- static pressure, shape product gradation parameter, to infer product matter
The superiority and inferiority of amount.Strain Method is a kind of method of more commonly used measurement P-t curve, and this method uses strain pressure transducer
Air pressure in the closed confined space is measured to change with time situation.This method has system frequency response slower, cheap, is easy mark
Fixed feature.Compared with slow, the longer igniter of pressure duration or for passing fire tool, strain is selected to those pressure rise times
Method measuring system, which measures pressure time curve, has good economy.
Strain pressure transducer is the device of Strain Method most critical, and reliability is to determine test experiments success or failure
It is crucial.In order to enhance the reliability of sensor, common method is to be examined before test experiments using pressure calibrating device at present
The linearity of sensor just carries out ignition test using the sensor after sensor linearity is up to standard.
But there are still some problems for pressure sensor during the test:
(1) transducer calibration is carried out under a metastable environment, and true Exploding test environment then pole
It is severe, and with the increase of testing time, the fragment for generation of exploding can cause the deformation of strain gauge transducer surface membrane, thus
Cause sensor linearity to decline, influences the accuracy of measurement parameter;
(2) a large amount of heat can be generated when Exploding test, steeply rises the temperature of confined space, so as to cause strain-type
Temperature drift when sensor measurement, this phenomenon can not be improved with calibration.
PCA (Principal Component Analysis) is a kind of common data analysing method.PCA passes through linear
Initial data is transformed to the expression of one group of each dimension linear independence by transformation, can be used for extracting the main feature component of data.
Summary of the invention
In view of this, the purpose of the present invention is to provide a kind of ignition test pressure sensor drift of Space Pyrotechnics Devices product is adaptive
Answer bearing calibration.This method analyzes the main feature of historical measurement data using PCA, can be with by counting a large amount of test sample
It was found that main feature is not a fixed value, and it is distributed across in a region, the drift of exactly strain pressure transducer causes
This uncertainty of main feature.When measuring new data, it can use the main feature distribution region that historical data obtains and come school
Positive sensor drift bring error, improves the accuracy and reliability of ignition test experiments data.
Space Pyrotechnics Devices product, which are under fire, tests pressure sensor drift self-adapting correction method, it is comprised the steps of:
Step 1: obtaining the P-t curve of historical test data from database, then analyzes P-t curve and obtains each characterization
The characteristic parameter of firer's quality, such as: starting to build the pressure time, surge pressure, surge pressure time, maximum climbing speed, standard are quiet
The frequency spectrum maximum value for the P-t curve zero center normalization instantaneous amplitude that traditional characteristics parameter and the present invention such as pressure use, P-t
The difference of twice of envelope mean square of envelope of curve square value, the variance of P-t curve instantaneous frequency, P-t curve sizes wave system number
Variance, the new features parameter such as the ratio between fourth order cumulant of P-t curve signal.N number of sample is extracted from historical measurement data, often
A sample includes M feature, is organized into the sample matrix X0 of a N*M;
Step 2: dimension-reduction treatment is carried out to the sample matrix using PCA pivot analysis algorithm, analysis, which is found, accounts for all features
The main feature of those of energy 90%, this Partial Feature can characterize the drift features of strain pressure transducer, and find out dimensionality reduction
Transition matrix P;
Step 3: when test of being under fire starts, using real-time in strain pressure transducer acquisition confined explosion space
Pressure signal is filtered amplification to the signal of pressure sensor output, is then believed pressure after filter and amplification by Acquisition Instrument
It number is input to fusion center (computer), fusion center records the pressure of this test of being under fire change situation at any time, obtains one
Then P-t curve analyzes P-t curve and obtains M characteristic parameter of characterization firer's quality, be organized into measuring test data
Vector X is corrected the data vector X newly measured using the projection matrix P obtained from step 2, obtains corrected value Xcorr=X-
XPP′;
Step 4: by corrected value XcorrUpdate sample matrix X0, again to X0New projection matrix is obtained using PCA algorithm
P0, by seeking P-P0Two norms compare P and P0Difference then use P if more than θ and compared with given threshold θ0Instead of P, no person
It is constant to retain P.Continue test if necessary and then return to step 3, no person exits.
Advantageous effects of the invention are as follows: the present invention can make full use of historical measurement data to correct because of strain-type pressure
Error caused by force snesor drifts about, the measurement data after this method corrects are substantially better than the measurement data of not correction up,
The uncertainty as caused by strain pressure transducer temperature drift is efficiently reduced, the precision of fire-working article test parameter is provided
And reliability.
Detailed description of the invention
In order to keep the purpose of the present invention, technical scheme and beneficial effects clearer, the present invention provides following attached drawing and carries out
Illustrate:
Fig. 1 is a kind of Space Pyrotechnics Devices product ignition test pressure sensor drift self-adapting correction method of the present invention
System construction drawing.
Fig. 2 is a kind of Space Pyrotechnics Devices product ignition test pressure sensor drift self-adapting correction method of the present invention
System flow chart.
Specific embodiment
Below in conjunction with attached drawing, Space Pyrotechnics Devices product test pressure sensor drift correction method of the invention is carried out detailed
Description.
Fig. 1 be a kind of Space Pyrotechnics Devices product of the present invention be under fire that test pressure sensor drifts about self-adapting correction method be
System structure chart, comprises the following modules: strain pressure transducer, Acquisition Instrument, fusion center (computer), is shown dynamic strain indicator
Show module, memory module, database.Wherein, fusion center is connected with Acquisition Instrument, display module, memory module.Utilize Acquisition Instrument
The pressure information after filter and amplification is sent fusion center by equipment, and fusion center obtains historical data by memory module,
And test data is corrected in conjunction with PCA-CC algorithm, the mistake caused by capable of improving strain pressure transducer drift
Difference effectively improves the precision and reliability of priming system ignition test experiments data.
Fig. 2 be a kind of Space Pyrotechnics Devices product of the present invention be under fire that test pressure sensor drifts about self-adapting correction method be
System flow chart, the specific implementation steps are as follows:
(1) N (N > 100) a homologous series fire-working article test data sample is obtained from memory module, each sample includes M
Parameter attribute is organized into matrix form, obtains the sample matrix X0 of a N*M;
(2) PCA dimensionality reduction is carried out to sample matrix X0.Sample matrix is subjected to centralization first, then calculates sample matrix
Covariance matrix C diagonalization is obtained diagonal matrix Λ by the covariance matrix C of X0, and it is strongest several that energy is chosen in diagonal matrix Λ
A characteristic value makes the energy of these characteristic values account for 90% or more of total characteristic value energy, it is assumed that k characteristic value is had chosen altogether,
Wherein k < M.The corresponding feature vector of this k characteristic value forms an eigenvectors matrixThis new feature vector
Matrix P is the projection matrix of PCA dimensionality reduction;
(3) ignition test experiments start, and strain pressure transducer measures the real-time pressure in confined explosion space, due to
The signal of pressure sensor output is very faint, is filtered amplification by output signal of the dynamic strain indicator to sensor;
(4) output signal of dynamic strain indicator is uploaded to fusion center (computer) by Acquisition Instrument, fusion center note
The case where this lower pressure for being under fire experiment of record changes over time, forms a P-t curve;
(5) the P-t curve of analysis this time ignition experiment, finds out M parameter attribute, such as: starting to build pressure time, peak value pressure
Power, in the P-t curve zero that traditional characteristics parameter and the present invention such as surge pressure time, maximum climbing speed, quasi- static pressure use
The heart normalizes the frequency spectrum maximum value of instantaneous amplitude, the difference of twice of envelope mean square of P-t envelope of curve square value, P-t curve
The variance of instantaneous frequency, the variance of P-t curve sizes wave system number, the new features ginseng such as the ratio between fourth order cumulant of P-t curve signal
Number.This M parameter attribute is organized into a vector X, it may be assumed that
(6) the data X newly measured is corrected by the projection matrix P obtained from step (2), obtains Xcorr=X-XPP '.
XcorrThe strain pressure transducer drift correction data as realized by PCA-CC algorithm combination historical data;
(7) correction data X is savedcorrTo database and display is exported, by corrected value XcorrUpdate sample matrix X0, weight
Newly to X0New projection matrix P is obtained using PCA algorithm0, by seeking P-P0Two norms compare P and P0Difference, and with setting threshold
Value θ compares, and then replaces P with P0 if more than θ, it is constant that no person retains P.Continue test if necessary and then return to step (3), no person moves back
Out.
Finally, it is stated that the above selection for being under fire test priming system parameter attribute is only to illustrate technical side of the invention
Case rather than limit, although be under fire test priming system parameter attribute selection the present invention is described in detail, this
It field, such as can be from P-t song it will be appreciated by the skilled person that various changes can be made to it in the form and details
Line obtains frequency domain character (wavelet analysis), and frequency domain-energy feature (Hilbert-Huang transform) etc. is wanted without departing from right of the present invention
Seek book limited range.
Claims (4)
- The test pressure sensor drift self-adapting correction method 1. a kind of Space Pyrotechnics Devices product are under fire, it is characterised in that: including following Module: module one: strain pressure transducer;Module two: dynamic strain indicator: module three: Acquisition Instrument;Module four: fusion center; Module five: display module;Module six: memory module;Module seven: database;Including following data processing step:(1) step 1: obtaining P-t curve of historical test data from database, then analyzes P-t curve and obtains each characterization The characteristic parameter of firer's quality, including start to build pressure time, surge pressure, surge pressure time, maximum climbing speed and standard Static pressure traditional characteristic parameter, and the P-t curve zero center used normalize the frequency spectrum maximum value of instantaneous amplitude, P-t curve The difference of twice of envelope mean square of envelope square value, the variance of P-t curve instantaneous frequency, P-t curve sizes wave system number Variance, the ratio between the fourth order cumulant new feature parameter of P-t curve signal;N number of sample is extracted from historical measurement data, each Sample includes M feature, is organized into the sample matrix X of a N*M0;(2) step 2: using PCA algorithm to sample matrix X0Dimension-reduction treatment is carried out, analysis, which is found, accounts for all characteristic energies 90% Those main features, this Partial Feature can characterize the drift features of strain pressure transducer, and find out the projection matrix P of dimensionality reduction;(3) step 3: when test of being under fire starts, the pressure value that pressure sensor acquires is uploaded into computer in real time, is recorded Pressure obtains P-t curve, and analysis P-t curve obtains M characteristic parameter of characterization firer's quality, is organized into X is measured, so The data vector X newly measured is corrected by projection matrix P afterwards, obtains corrected value Xcorr=X-XPP':(4) step 4: corrected value Xcorr is updated to sample matrix X0, again to X0New projection matrix is obtained using PCA algorithm P0, by seeking P-P0Two norms compare P and P0Difference, and compared with given threshold θ, if more than 0 with P0Instead of P, no person It is constant to retain P;Continue test if necessary and then return to step 3, no person exits.
- The test pressure sensor drift self-adapting correction method 2. a kind of Space Pyrotechnics Devices product according to claim 1 are under fire, It is characterized by: being carried out on the basis of the traditional priming system characteristic parameter obtained from P-t tracing analysis in the step 1 The expansion and fusion of feature, frequency spectrum maximum value, curve including obtaining curve zero center normalization instantaneous amplitude from P-t curve The difference of twice of envelope mean square of envelope square value, the ratio between the fourth order cumulant new feature parameter of curve signal.
- The test pressure sensor drift self-adapting correction method 3. a kind of Space Pyrotechnics Devices product according to claim 1 are under fire, It is characterized by: PCA-CC algorithm is applied to Space Pyrotechnics Devices product and is under fire in test, analysis was tested in the step 3 The drift features of strain pressure transducer in journey, the pressure sensor linear fit calibration relatively only tested in the past before more can Grasp the dynamic property variation of pressure sensor, PCA-CC combination historical data to the drift of pressure sensor in test process into Row correction, improves the reliability of Space Pyrotechnics Devices product test parameter.
- The test pressure sensor drift self-adapting correction method 4. a kind of Space Pyrotechnics Devices product according to claim 1 are under fire, It is characterized by: in the step 4 new measured value being updated to sample database obtain new sample matrix in time, to new sample This matrix carries out PCA dimensionality reduction and obtains new projection matrix P0, compares P and P0 difference, updates and throws when difference reaches given threshold Shadow matrix, to realize the adaptively correcting of pressure sensor drift.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610814157.7A CN106370347B (en) | 2016-09-09 | 2016-09-09 | A kind of Space Pyrotechnics Devices product ignition test pressure sensor drift self-adapting correction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610814157.7A CN106370347B (en) | 2016-09-09 | 2016-09-09 | A kind of Space Pyrotechnics Devices product ignition test pressure sensor drift self-adapting correction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106370347A CN106370347A (en) | 2017-02-01 |
CN106370347B true CN106370347B (en) | 2019-03-19 |
Family
ID=57900274
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610814157.7A Expired - Fee Related CN106370347B (en) | 2016-09-09 | 2016-09-09 | A kind of Space Pyrotechnics Devices product ignition test pressure sensor drift self-adapting correction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106370347B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109387124A (en) * | 2018-08-23 | 2019-02-26 | 邢立平 | A kind of machinery priming system stab sensitivity test method |
CN110907061B (en) * | 2019-12-05 | 2021-10-08 | 成都理工大学 | Nuclear waste barrel temporary storage warehouse heat source distribution monitoring simulation device and monitoring method |
CN111934797B (en) * | 2020-07-09 | 2022-06-24 | 广东工业大学 | Collaborative spectrum sensing method based on covariance eigenvalue and mean shift clustering |
CN113686486B (en) * | 2021-08-19 | 2024-01-16 | 上海卫星工程研究所 | Method and system for calibrating whole star of satellite propulsion system pressure sensor |
CN114136537A (en) * | 2021-11-04 | 2022-03-04 | 歌尔微电子股份有限公司 | Pressure sensor |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201339A (en) * | 2006-12-13 | 2008-06-18 | 天津科技大学 | Apparatus and method for monitoring resistance spot welding quality |
CN202018302U (en) * | 2011-03-09 | 2011-10-26 | 中国人民解放军总装备部军械技术研究所 | Initiating explosive device action pressure-temperature test instrument |
CN102959378A (en) * | 2010-06-28 | 2013-03-06 | 通用电气公司 | Temperature independent pressure sensor and associated methods thereof |
CN103559401A (en) * | 2013-11-08 | 2014-02-05 | 渤海大学 | Failure monitoring method based on semi-supervised principal component analysis |
CN203732198U (en) * | 2013-12-20 | 2014-07-23 | 重庆川仪自动化股份有限公司 | Pressure transmitter |
CN104776947A (en) * | 2015-04-03 | 2015-07-15 | 袁川来 | Compensation method for temperature drift of robot wrist force sensor |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101103904B1 (en) * | 2009-10-01 | 2012-01-12 | 한전케이피에스 주식회사 | Pressure Transmitter Auto Calibration |
-
2016
- 2016-09-09 CN CN201610814157.7A patent/CN106370347B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201339A (en) * | 2006-12-13 | 2008-06-18 | 天津科技大学 | Apparatus and method for monitoring resistance spot welding quality |
CN102959378A (en) * | 2010-06-28 | 2013-03-06 | 通用电气公司 | Temperature independent pressure sensor and associated methods thereof |
CN202018302U (en) * | 2011-03-09 | 2011-10-26 | 中国人民解放军总装备部军械技术研究所 | Initiating explosive device action pressure-temperature test instrument |
CN103559401A (en) * | 2013-11-08 | 2014-02-05 | 渤海大学 | Failure monitoring method based on semi-supervised principal component analysis |
CN203732198U (en) * | 2013-12-20 | 2014-07-23 | 重庆川仪自动化股份有限公司 | Pressure transmitter |
CN104776947A (en) * | 2015-04-03 | 2015-07-15 | 袁川来 | Compensation method for temperature drift of robot wrist force sensor |
Non-Patent Citations (1)
Title |
---|
压力传感器温度补偿的一种新方法;孙艳梅、刘树东;《光通信研究》;20110210(第1期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106370347A (en) | 2017-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106370347B (en) | A kind of Space Pyrotechnics Devices product ignition test pressure sensor drift self-adapting correction method | |
CN105913079B (en) | Electronic nose isomeric data recognition methods based on the study of the target domain migration limit | |
CN107490397B (en) | High-accuracy self-adaptation filters the quick Peak Search Method of FBG spectrum | |
CN109685334A (en) | A kind of new hydrological model simulation evaluation method based on Multiscale Theory | |
CN109993398A (en) | A kind of air-to-ground guided missile Accuracy Asse ssment method based on Error Spectrum | |
CN112329142B (en) | Method for identifying uncertain dynamic loads of wing structure based on support vector regression | |
CN104458478B (en) | System and method for analysis in modulated thermogravimetry | |
CN109738392A (en) | Compressed sensing reconstructing method and system towards oxygen concentration in TDLAS on-line checking bottle | |
CN116989972A (en) | Method for evaluating uncertainty of heat measurement test data based on thermoelectric simulation network test | |
Oz et al. | Electrochemical impedance spectroscopy of supercapacitors: a novel analysis approach using evolutionary programming | |
CN103256999A (en) | Distributed type optical fiber temperature measuring method | |
Xu et al. | Dynamic compensation of piezoresistive pressure sensor based on sparse domain | |
Martinović et al. | Solar wind electron parameters determination on Wind spacecraft using quasi‐thermal noise spectroscopy | |
CN107203677B (en) | Electronic system multi-degradation process research method | |
US10001579B2 (en) | High frequency electromagnetic induction | |
Chen et al. | The implementation of the data validation process in a gas turbine performance monitoring system | |
Li et al. | Investigation of filtering algorithm for noise reduction in displacement sensing signal | |
CN116644658B (en) | Gas discharge temperature inversion method based on physical information neural network | |
Klevtsov | Using the Method of Normalized Amplitude for Assessing the Quality of the Calibration Tests of the Pressure Sensor | |
CN109145417A (en) | A method of Indentation strain method Stress calculation function is directly determined based on material mechanical performance | |
Peng et al. | A novel pressure sensor calibration system based on a neural network | |
CN111751037B (en) | Electric measurement explosion test data compression method | |
Tuttle et al. | Evaluation of XHVRB for capturing transition to detonation as measured by embedded gauges | |
GUO et al. | Cross-correlation analysis method on explosion seismic waves data | |
CN110389037A (en) | A kind of diesel engine intake system for measuring quantity |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190319 Termination date: 20190909 |