CN116136518A - Chromatograph - Google Patents
Chromatograph Download PDFInfo
- Publication number
- CN116136518A CN116136518A CN202310427488.5A CN202310427488A CN116136518A CN 116136518 A CN116136518 A CN 116136518A CN 202310427488 A CN202310427488 A CN 202310427488A CN 116136518 A CN116136518 A CN 116136518A
- Authority
- CN
- China
- Prior art keywords
- chromatogram
- chromatographic
- wavelet
- function
- denoising
- 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
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000001914 filtration Methods 0.000 claims abstract description 20
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000001514 detection method Methods 0.000 claims abstract description 12
- 238000005070 sampling Methods 0.000 claims description 29
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 4
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 238000004587 chromatography analysis Methods 0.000 claims 1
- 230000002159 abnormal effect Effects 0.000 abstract description 8
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000004445 quantitative analysis Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 239000012080 ambient air Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004817 gas chromatography Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000001819 mass spectrum Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8603—Signal analysis with integration or differentiation
- G01N30/8617—Filtering, e.g. Fourier filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8624—Detection of slopes or peaks; baseline correction
- G01N30/8627—Slopes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8624—Detection of slopes or peaks; baseline correction
- G01N30/8631—Peaks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8603—Signal analysis with integration or differentiation
- G01N2030/862—Other mathematical operations for data preprocessing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
- Y02A50/30—Against vector-borne diseases, e.g. mosquito-borne, fly-borne, tick-borne or waterborne diseases whose impact is exacerbated by climate change
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Treatment Of Liquids With Adsorbents In General (AREA)
Abstract
The invention relates to a chromatograph, which comprises a sample inlet, a chromatographic column, a detector and a signal processing module, wherein the signal processing module is used for processing an electric signal detected by the detector to obtain a chromatogram; the signal decomposition and filtering module is used for performing empirical mode decomposition on chromatographic signals of a chromatogram of a target sample to be detected to obtain an eigenmode function and a high-frequency signal of the target sample to be detected, and performing Bayes wavelet filtering on the high-frequency signal to obtain a wavelet reconstructed signal; the method is also used for carrying out empirical mode decomposition on chromatographic signals of a chromatogram of the standard sample to obtain a standard intrinsic base mode function; the method is also used for fitting the eigenmode function of the target sample to be measured by utilizing the standard eigenmode function to obtain a fitted eigenmode function; and the signal reconstruction module is used for carrying out signal reconstruction on the wavelet reconstructed signal and the eigen-mode function after fitting to obtain a denoising chromatogram. The invention is beneficial to eliminating the influence of abnormal peaks and improving the detection precision.
Description
Technical Field
The invention belongs to the technical field of detection and analysis, and particularly relates to a chromatograph.
Background
The chromatograph is used for separating components of the substance to be detected, so that chromatograms of the corresponding substance at different positions are obtained, and quantitative analysis is carried out based on chromatographic peak areas or heights. In recent years, with the development of chromatographic technology and matching detectors, the accuracy and sensitivity of detection are continuously improved, and the method is widely applied to detection, research and analysis in various industries.
The determination of the peak type and the base line of the chromatographic peak is the key of accurate quantitative analysis, the traditional peak area calculation generally determines the rising start point of the chromatographic peak according to the set threshold value and the times higher than the threshold value, and simultaneously determines the falling end point of the chromatographic peak according to the times lower than the corresponding threshold value, thereby determining the base line and obtaining the corresponding peak type; wherein the setting of the threshold is determined based on empirical values. However, for low-concentration sample analysis and liquid-mass or gas-mass series systems (often with higher noise), the conventional method has deviation in determining the base line and the peak, so that the error of calculating the peak area to realize quantitative calculation is larger, and the later stage needs to rely on operators with abundant experience to perform manual calculation, however, the processing habit of the operators still causes larger calculation difference, the accuracy of the quantitative result cannot be ensured, and the detection efficiency is reduced.
In addition, background components of the detected substances under different working conditions are complex, and as the situation that a mixed peak or abnormal noise is superimposed on a detected chromatographic peak is caused due to overlong use time or improper maintenance of an instrument, the base line on two sides of the shoulder peak or the chromatographic peak is fluctuated too much, and the base line and the peak form are determined by the traditional method, so that the calculation deviation of the peak area is caused, and the accuracy of quantitative analysis is affected.
Disclosure of Invention
Based on the above-mentioned drawbacks and deficiencies of the prior art, it is an object of the present invention to at least solve one or more of the above-mentioned problems of the prior art, in other words, to provide a chromatograph satisfying one or more of the aforementioned needs.
In order to achieve the aim of the invention, the invention adopts the following technical scheme:
a chromatograph, comprising:
the sample inlet is used for carrying out sample injection on the detection sample; the sample to be detected is a standard sample or a target sample to be detected;
a chromatographic column for separating components of the test sample;
a detector for detecting the separated components and outputting an electrical signal;
the signal processing module is used for processing the electric signals to obtain a chromatogram;
the signal decomposition and filtering module is used for performing empirical mode decomposition on chromatographic signals of a chromatogram of a target sample to be detected to obtain an eigenmode function and a high-frequency signal of the target sample to be detected, and performing Bayes wavelet filtering on the high-frequency signal to obtain a wavelet reconstructed signal; the method is also used for carrying out empirical mode decomposition on chromatographic signals of a chromatogram of the standard sample to obtain a standard intrinsic base mode function; the method is also used for fitting the eigenmode function of the target sample to be measured by utilizing the standard eigenmode function to obtain a fitted eigenmode function;
the signal reconstruction module is used for carrying out signal reconstruction on the wavelet reconstructed signals and the fitted eigenmode functions to obtain a denoising chromatogram;
and the calculation module is used for calculating and obtaining the concentration of each component in the target sample to be detected according to the denoising chromatogram.
Preferably, the chromatograph further comprises:
the interval dividing module is used for dividing the chromatogram of the standard sample into the target characteristic peaksMA plurality of intervals, each interval including at least one target characteristic peak;
the standard eigen-mode function corresponding to each interval is:
wherein ,is the firstmThe gaussian function corresponding to each interval,a m is the firstmThe function coefficients of the gaussian function corresponding to the respective intervals,T 1 、T 2 respectively the firstmThe time corresponding to the start and end of each interval,tin order to be able to take time,m∈[1,M],Mis a positive integer.
Preferably, the bayesian wavelet filtering of the high-frequency signal includes:
(1) Performing multi-level wavelet decomposition on the high-frequency signals corresponding to each section; the method comprises the steps of denoising a high-frequency signal obtained by wavelet decomposition of each stage by using a Bayes threshold as a wavelet threshold, and performing next-stage wavelet decomposition on the obtained denoised high-frequency signal;
(2) And carrying out wavelet reconstruction on the low-frequency signal and the de-dried high-frequency signal obtained by decomposing the wavelet of each stage to obtain a wavelet reconstructed signal.
Preferably, the number of wavelet decomposition stepsZThe value is an integer between 3 and 8.
As a preferred embodiment, the firstzThe Bayes threshold corresponding to denoising after the level wavelet decomposition is as follows:
wherein ,respectively the firstmCorresponding to the first intervalzHigh-frequency signal and low-frequency signal obtained by level wavelet decomposition, < >>Respectively the firstmCorresponding to the first intervalz+A high-frequency signal and a low-frequency signal obtained by 1-level wavelet decomposition,th1 and th2 are respectively a preset lower threshold value and an upper threshold value, < + >>Is->Is used for the normal deviation of the (c) in the (c),μas the weighting coefficient(s),signas a function of the sign of the symbol,medianas a function of the median value,z∈[1,Z]。
preferably, the method comprisesth1 and th2 respectively takes the value of and />,μThe value is 0.5.
As a preferred scheme, the process of calculating the concentration of each component in the target sample to be measured according to the denoising chromatogram comprises the following steps:
based on the denoising chromatogram, a preset threshold value is utilizedth3. Slope ofkLine widthL w Determining chromatographic peaks of the components;
calculating the peak area according to the chromatographic peak;
the concentration of the component was calculated from the peak area.
Preferably, the process of determining the chromatographic peak of the component from the denoising chromatogram comprises:
denoising the chromatogram according to the sampling frequencyIs converted from time to the number of sampling pointsiThe corresponding chromatogram is +.>;
iStarting from zero when it appearsAnd->The slope formed by the subsequent sampling points is greater than the slopekWhen in use, then->As the left starting point of the chromatographic peak, the corresponding coordinate of the chromatographic peak in the denoising chromatogram is%L 1 ,FL 1 );
Along withiContinue to increase when it occursAnd->The slope formed by the subsequent sampling points is smaller than the slopekWhen the negative value of (2) is->As the right end point of the chromatographic peak, the corresponding coordinate of the chromatographic peak in the denoising chromatogram is%R 1 ,FR 1 );
JudgingR 1 -L 1 Whether or not to be not smaller than line widthL w The method comprises the steps of carrying out a first treatment on the surface of the If yes, then the first chromatographic peak is determinedPerforming least square fitting by using a preset number of sampling points before and after the left starting point and a preset number of sampling points before and after the right ending point to obtain a baseline function of a first chromatographic peak;
after the first chromatographic peak is determined, followingiContinuing to increase, determination of the remaining chromatographic peaks and so on, yielding left onset, right endpoint and baseline functions for all chromatographic peaks.
Preferably, the peak area is calculated from the chromatographic peaks as:
wherein ,L n 、R n respectively the firstnThe left starting point and the right ending point of each chromatographic peak are positioned on the abscissa of the denoising chromatogram,is the firstnThe baseline function of the individual chromatographic peaks,n∈[1,N],Nis the total number of chromatographic peaks;
the concentration of the components was calculated from the peak areas as:
wherein ,F(P,T,V) As to pressurePTemperature (temperature)TAnd sampling flow rateVIs used to determine the fitting coefficients of (a),βis a correction coefficient.
Preferably, the detector is a FID detector or an MS detector.
Compared with the prior art, the invention has the beneficial effects that:
(1) According to the chromatograph, the standard eigenmode function obtained by empirical mode decomposition of chromatographic signals of a chromatogram of a standard sample is utilized to fit the eigenmode function of a target sample to be measured to obtain a fitted eigenmode function, and the signal reconstruction is carried out on the wavelet reconstructed signals obtained by carrying out Bayesian wavelet filtering on the fitted eigenmode function and the target sample to be measured, so that the influence of abnormal peaks (such as shoulder peaks, pits and the like) can be eliminated;
(2) The chromatograph provided by the invention utilizes the improved Bayes threshold to carry out multistage wavelet decomposition, and has the advantages of eliminating shoulder peaks, being good in noise reduction effect, high in instantaneity, easy to realize and the like;
(3) The chromatograph of the invention determines the chromatographic peak of each component to be detected by using a threshold value, slope and line width three-element judgment method, which is beneficial to eliminating abnormal peaks and eliminating the influence of shoulder peaks caused by noise.
Drawings
FIG. 1 is a block diagram of a chromatograph according to embodiment 1 of the present invention;
FIG. 2 is a flow chart of chromatographic signal filtering and reconstruction of the chromatograph of example 1 of the present invention;
FIG. 3 is a graph showing the comparison of a denoising chromatogram obtained by wavelet filtering with an improved Bayesian threshold with an original chromatogram and a chromatogram obtained by wavelet filtering with a generic Bayesian threshold according to embodiment 1 of the present invention;
fig. 4 is a comparison of the denoising chromatogram of example 1 of the present invention with a comparison chromatogram obtained by processing using the existing denoising method.
Detailed Description
In order to more clearly illustrate the embodiments of the present invention, specific embodiments of the present invention will be described below with reference to the accompanying drawings. It is evident that the drawings in the following description are only examples of the invention, from which other drawings and other embodiments can be obtained by a person skilled in the art without inventive effort.
Example 1:
as shown in fig. 1, the chromatograph of the present embodiment includes a sample inlet, a chromatographic column, a detector, a signal processing module, a signal decomposing and filtering module, a signal reconstructing module, a section dividing module and a calculating module. The sample inlet, the chromatographic column, the detector and the signal processing module are all basic components of a conventional chromatograph, and the specific structure of the sample inlet, the chromatographic column, the detector and the signal processing module is not repeated herein.
The sample inlet of the embodiment is used for sampling the detection sample; the sample to be detected is a standard sample or a target sample to be detected, and the standard sample is formed according to all target components to be detected in the target sample to be detected.
The chromatographic column of the embodiment is used for separating components of a detection sample, then entering a detector, and detecting the separated components to obtain corresponding electric signals; the electric signal is input into a signal processing module, and the electric signal is processed to obtain a chromatogram. Wherein the detector is an FID detector.
The chromatograph of the embodiment firstly carries out chromatographic detection on a standard sample to obtain a chromatogram of the standard sample; and then carrying out chromatographic detection on the target sample to be detected to obtain a chromatogram of the target sample to be detected.
The signal decomposition and filtering module of the embodiment is used for performing empirical mode decomposition on chromatographic signals of a chromatogram of a standard sample to obtain a standard intrinsic base mode function.
Specifically, the interval division module of the present embodiment is configured to divide the chromatogram of the standard sample into the target characteristic peaksMEach interval including at least one target characteristic peak. The abscissa of the chromatogram is time, the time can be converted into the number of sampling points according to the sampling frequency, the number of the sampling points in each interval is unequal, and the determination is specifically performed according to the actual distribution condition of the target characteristic peak.
Correspondingly, the standard eigen-basis mode function corresponding to each interval is:
wherein ,is the firstmThe gaussian function corresponding to each interval,a m is the firstmThe function coefficients of the gaussian function corresponding to the respective intervals,T 1 、T 2 respectively the firstmThe time corresponding to the start and end of each interval,tin order to be able to take time,m∈[1,M],Mis a positive integer.
As shown in fig. 2, the signal decomposition and filtering module of the present embodiment is further configured to perform empirical mode decomposition on a chromatographic signal of a chromatogram of a target sample to be measured, obtain an eigenmode function and a high-frequency signal of the target sample to be measured, and perform bayesian wavelet filtering on the high-frequency signal to obtain a signal after wavelet reconstruction.
Specifically, the interval division module of the present embodiment also corresponds to the chromatogram of the target sample to be measured to the chromatogram of the standard sampleMDividing each interval in one-to-one correspondence to obtain chromatographic signals of each interval. Chromatographic signals for the respective intervals +.>Empirical mode decomposition is carried out to obtain: />,/>Is the firstmChromatographic signal of individual interval->Corresponding eigenmode function,/->Is the firstmChromatographic signal of individual interval->Corresponding high frequency signals.
In addition, the signal decomposition and filtering module of the present embodiment is also used for utilizing standard eigen-basis mode functionsEigenmode function of the target sample to be tested>Fitting to obtain a fitted eigenmode function, and performing subsequent signal reconstruction by using the fitted eigenmode function to facilitate eliminating abnormal peaks (such as shoulder peak, pit, etc.)) Is a function of (a) and (b).
The process of bayesian wavelet filtering of the high-frequency signal according to the embodiment comprises the following steps:
(1) Performing multi-level wavelet decomposition on the high-frequency signals corresponding to each section; the method comprises the steps of denoising a high-frequency signal obtained by wavelet decomposition of each stage by using a Bayes threshold as a wavelet threshold, and performing next-stage wavelet decomposition on the obtained denoised high-frequency signal; the number of wavelet decomposition steps of the present embodimentZThe value is an integer between 3 and 8, which can be determined according to practical application conditions, and the embodiment usesZTake the value 3 as an example, as shown in fig. 2, namely, a threshold value 1, a threshold value 2 and a threshold value 3. In addition, the specific process of wavelet decomposition may refer to the prior art, and is not described herein. The present embodiment improves the wavelet threshold used for denoising, specifically, the firstzThe Bayes threshold corresponding to denoising after the level wavelet decomposition is as follows:
wherein ,respectively the firstmCorresponding to the first intervalzHigh-frequency signal and low-frequency signal obtained by level wavelet decomposition, < >>Respectively the firstmCorresponding to the first intervalz+A high-frequency signal and a low-frequency signal obtained by 1-level wavelet decomposition,th1 and th2 are respectively a preset lower threshold value and an upper threshold value, < + >>Is->Is used for the normal deviation of the (c) in the (c),μas the weighting coefficient(s),signas a function of the sign of the symbol,medianas a function of the median value,z∈[1,Z]the method comprises the steps of carrying out a first treatment on the surface of the As an example, the present embodimentth1 and th2 is respectively given the value-> and />,μThe value is 0.5;
(2) And carrying out wavelet reconstruction on the low-frequency signal and the de-dried high-frequency signal obtained by decomposing the wavelet of each stage to obtain a wavelet reconstructed signal.
The signal reconstruction module of the embodiment is used for carrying out signal reconstruction on the signal after wavelet reconstruction and the eigen-mode function after fitting to obtain a denoising chromatogram. The Bayes wavelet filter of the embodiment is adopted for denoising, and has the advantages of eliminating shoulder peaks, being good in denoising effect, high in instantaneity, easy to realize and the like.
As shown in fig. 3, in the low-range sample (ambient air) measurement, the bayesian wavelet filtering process of the embodiment has obvious suppression on noise, and the core is that the selection of the bayesian wavelet threshold is completely different from the suppression degree of the noise, and the improved bayesian threshold of the embodiment is adoptedSubstitute general Bayesian threshold->Fluctuation and abnormal peaks caused by noise can be effectively suppressed, and the accuracy of baseline judgment is effectively improved, so that the accuracy and the degree of automation of peak area calculation are improved.
In addition, as shown in fig. 4, with respect to pits caused by noise superimposed on a chromatographic peak, pits cannot be eliminated by the denoising method disclosed in the document of publication No. CN115700544a in the prior art; the pit can be effectively eliminated by adopting a Gaussian fundamental mode (namely a standard intrinsic fundamental mode function), and the detection precision is improved.
The calculation module of the embodiment is used for calculating and obtaining the concentration of each component in the target sample to be detected according to the denoising chromatogram. Specifically, the process of calculating the concentration of each component in the target sample to be measured according to the denoising chromatogram comprises the following steps: based on the denoising chromatogram, a preset threshold value is utilizedth3. Slope ofkLine widthL w Determining chromatographic peaks of the components; calculating the peak area according to the chromatographic peak; the concentration of the component was calculated from the peak area. The present embodiment utilizesth3. Slope ofkLine widthL w The three-element judgment method determines the chromatographic peak of each component to be detected, obtains the starting point and the end point of the chromatographic peak, eliminates the abnormal peak and eliminates the influence of the shoulder peak caused by noise.
Specifically, the process of determining the chromatographic peak of a component from a denoising chromatogram includes:
denoising the chromatogram according to the sampling frequencyIs converted from time to the number of sampling pointsiThe corresponding chromatogram is +.>;
iStarting from zero when it appearsAnd->The slope formed by the six later sampling points (the number of the sampling points is not limited to six and can be determined between four and eight according to actual requirements) is larger than the slopekWhen in use, then->As the left starting point of the chromatographic peak, the corresponding coordinate of the chromatographic peak in the denoising chromatogram is%L 1 ,FL 1 );
Along withiContinue to increase when it occursAnd->The slope formed by the six later sampling points (the number of the sampling points is not limited to six and can be determined between four and eight according to actual requirements) is smaller than the slopekNegative value (i.e.)k) When in use, then->As the right end point of the chromatographic peak, the corresponding coordinate of the chromatographic peak in the denoising chromatogram is%R 1 ,FR 1 );
JudgingR 1 -L 1 Whether or not to be not smaller than line widthL w The method comprises the steps of carrying out a first treatment on the surface of the If yes, determining a first chromatographic peak, and performing least square fitting by using three sampling points (the number of the sampling points is not limited to three and can be determined according to actual requirements) before and after the left starting point and three sampling points (the number of the sampling points is not limited to three and can be determined according to actual requirements) before and after the right ending point to obtain a baseline function of the first chromatographic peakThe method comprises the steps of carrying out a first treatment on the surface of the If not, the peak is a mixed peak or an abnormal peak, the peak record and calculation are not carried out, and the chromatographic peak is continuously searched;
after the first chromatographic peak is determined, followingiContinuing to increase, determination of the remaining chromatographic peaks and so on, yielding left onset, right endpoint and baseline functions for all chromatographic peaks.
The peak area calculated from the chromatographic peaks in this example is:
wherein ,L n 、R n respectively the firstnThe left starting point and the right ending point of each chromatographic peak are positioned on the abscissa of the denoising chromatogram,is the firstnThe baseline function of the individual chromatographic peaks,n∈[1,N],Nis the total number of chromatographic peaks.
The concentration of the components calculated from the peak area in this example was:
wherein ,F(P,T,V) As to pressurePTemperature (temperature)TAnd sampling flow rateVIs used to determine the fitting coefficients of (a),βfor the determination of the correction coefficients, the fitting coefficients and the correction coefficients, reference is made in particular to the prior art.
The chromatograph of the embodiment carries out peak searching and baseline determining on the denoising chromatogram obtained by the improved Bayesian wavelet filtering treatment, reduces the misjudgment probability, ensures more accurate calculation of peak area and concentration, and meets the automatic and real-time online monitoring requirements.
Example 2:
the chromatograph of this embodiment is different from embodiment 1 in that:
the detector adopts an MS detector, is suitable for gas chromatography and mass spectrum combination, and meets the requirements of different applications;
other structures can be referred to embodiment 1.
The foregoing is only illustrative of the preferred embodiments and principles of the present invention, and changes in specific embodiments will occur to those skilled in the art upon consideration of the teachings provided herein, and such changes are intended to be included within the scope of the invention as defined by the claims.
Claims (10)
1. A chromatograph, comprising:
the sample inlet is used for carrying out sample injection on the detection sample; the sample to be detected is a standard sample or a target sample to be detected;
a chromatographic column for separating components of the test sample;
a detector for detecting the separated components and outputting an electrical signal;
the signal processing module is used for processing the electric signals to obtain a chromatogram;
the signal decomposition and filtering module is used for performing empirical mode decomposition on chromatographic signals of a chromatogram of a target sample to be detected to obtain an eigenmode function and a high-frequency signal of the target sample to be detected, and performing Bayes wavelet filtering on the high-frequency signal to obtain a wavelet reconstructed signal; the method is also used for carrying out empirical mode decomposition on chromatographic signals of a chromatogram of the standard sample to obtain a standard intrinsic base mode function; the method is also used for fitting the eigenmode function of the target sample to be measured by utilizing the standard eigenmode function to obtain a fitted eigenmode function;
the signal reconstruction module is used for carrying out signal reconstruction on the wavelet reconstructed signals and the fitted eigenmode functions to obtain a denoising chromatogram;
and the calculation module is used for calculating and obtaining the concentration of each component in the target sample to be detected according to the denoising chromatogram.
2. The chromatograph of claim 1, further comprising:
the interval dividing module is used for dividing the chromatogram of the standard sample into the target characteristic peaksMA plurality of intervals, each interval including at least one target characteristic peak;
the standard eigen-mode function corresponding to each interval is:
wherein ,is the firstmThe gaussian function corresponding to each interval,a m is the firstmThe function coefficients of the gaussian function corresponding to the respective intervals,T 1 、T 2 respectively the firstmThe time corresponding to the start and end of each interval,tin order to be able to take time,m∈[1,M],Mis a positive integer.
3. The chromatograph of claim 2, wherein the bayesian wavelet filtering of the high frequency signal comprises:
(1) Performing multi-level wavelet decomposition on the high-frequency signals corresponding to each section; the method comprises the steps of denoising a high-frequency signal obtained by wavelet decomposition of each stage by using a Bayes threshold as a wavelet threshold, and performing next-stage wavelet decomposition on the obtained denoised high-frequency signal;
(2) And carrying out wavelet reconstruction on the low-frequency signal and the de-dried high-frequency signal obtained by decomposing the wavelet of each stage to obtain a wavelet reconstructed signal.
4. A chromatograph according to claim 3, characterized in that the number of wavelet decomposition stepsZThe value is an integer between 3 and 8.
5. The chromatograph of claim 4, wherein the firstzThe Bayes threshold corresponding to denoising after the level wavelet decomposition is as follows:
wherein ,respectively the firstmCorresponding to the first intervalzHigh-frequency signal and low-frequency signal obtained by level wavelet decomposition, < >>Respectively the firstmCorresponding to the first intervalz+A high-frequency signal and a low-frequency signal obtained by 1-level wavelet decomposition,th1 and th2 are respectively a preset lower limit threshold value and an upper limitThreshold value (S)>Is->Is used for the normal deviation of the (c) in the (c),μas the weighting coefficient(s),signas a function of the sign of the symbol,medianas a function of the median value,z∈[1,Z]。
7. The chromatograph according to claim 5, wherein the process of calculating the concentration of each component in the target sample to be measured from the denoising chromatogram comprises:
based on the denoising chromatogram, a preset threshold value is utilizedth3. Slope ofkLine widthL w Determining chromatographic peaks of the components;
calculating the peak area according to the chromatographic peak;
the concentration of the component was calculated from the peak area.
8. The chromatograph of claim 7, wherein the process of determining the chromatographic peak of the component from the de-noised chromatogram comprises:
denoising the chromatogram according to the sampling frequencyIs converted from time to the number of sampling pointsiThe corresponding chromatogram is +.>;
iStarting from zero when it appearsAnd->The slope formed by the subsequent sampling points is greater than the slopekWhen in use, then->As the left starting point of the chromatographic peak, the corresponding coordinate of the chromatographic peak in the denoising chromatogram is%L 1 ,FL 1 );
Along withiContinue to increase when it occursAnd->The slope formed by the subsequent sampling points is smaller than the slopekWhen the negative value of (2) is->As the right end point of the chromatographic peak, the corresponding coordinate of the chromatographic peak in the denoising chromatogram is%R 1 ,FR 1 );
JudgingR 1 -L 1 Whether or not to be not smaller than line widthL w The method comprises the steps of carrying out a first treatment on the surface of the If yes, determining a first chromatographic peak, and performing least square fitting by using a preset number of sampling points before and after a left starting point and a preset number of sampling points before and after a right ending point to obtain a baseline function of the first chromatographic peak;
after the first chromatographic peak is determined, followingiContinuing to increase, determination of the remaining chromatographic peaks and so on, yielding left onset, right endpoint and baseline functions for all chromatographic peaks.
9. The chromatograph of claim 8, wherein the peak area is calculated from the chromatographic peaks as:
wherein ,L n 、R n respectively the firstnThe left starting point and the right ending point of each chromatographic peak are positioned on the abscissa of the denoising chromatogram,is the firstnThe baseline function of the individual chromatographic peaks,n∈[1,N],Nis the total number of chromatographic peaks;
the concentration of the components was calculated from the peak areas as:
wherein ,F(P,T,V) As to pressurePTemperature (temperature)TAnd sampling flow rateVIs used to determine the fitting coefficients of (a),βis a correction coefficient.
10. The chromatograph according to any of claims 1-9, wherein said detector is a FID detector or an MS detector.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310427488.5A CN116136518B (en) | 2023-04-20 | 2023-04-20 | Chromatograph |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310427488.5A CN116136518B (en) | 2023-04-20 | 2023-04-20 | Chromatograph |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116136518A true CN116136518A (en) | 2023-05-19 |
CN116136518B CN116136518B (en) | 2023-08-01 |
Family
ID=86333758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310427488.5A Active CN116136518B (en) | 2023-04-20 | 2023-04-20 | Chromatograph |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116136518B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116399996A (en) * | 2023-06-02 | 2023-07-07 | 海能未来技术集团股份有限公司 | Organic element analysis method and device based on thermal conductivity detection and electronic equipment |
Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0694696A (en) * | 1992-09-17 | 1994-04-08 | Hitachi Ltd | Method for analyzing chromatogram and chromatographic device |
JP2006177980A (en) * | 2006-03-27 | 2006-07-06 | Hitachi Ltd | Chromatographic data processing device, chromatographic data processing method, and chromatographic analyzer |
CN101131380A (en) * | 2006-08-21 | 2008-02-27 | 安捷伦科技有限公司 | Method and apparatus for detecting trace substance using chromatographic instrument |
JP2009008582A (en) * | 2007-06-29 | 2009-01-15 | Shimadzu Corp | Chromatogram data processor |
US20100161238A1 (en) * | 2007-05-31 | 2010-06-24 | Politecnico Di Milano | Method and System for Detecting Peptide Peaks in HPLC-MS Signals |
US20110161053A1 (en) * | 2009-12-31 | 2011-06-30 | Industrial Technology Research Institute | Slope-based fast intrinsic mode functions decomposition method and apparatus |
US20120179389A1 (en) * | 2009-08-20 | 2012-07-12 | Spectrosense Ltd. | Gas Chromatographic Analysis Method and System |
CN103902844A (en) * | 2014-04-24 | 2014-07-02 | 国家电网公司 | Transformer vibration signal de-noising method based on EEMD kurtosis threshold value |
US20140257712A1 (en) * | 2011-09-05 | 2014-09-11 | Shimadzu Corporation | System and method for processing chromatogram data |
CN104297504A (en) * | 2014-10-22 | 2015-01-21 | 上海申腾信息技术有限公司 | Automatic gas chromatographic control system |
US20150308992A1 (en) * | 2012-12-07 | 2015-10-29 | Shimadzu Corporation | Data processing system and data processing method for chromatograph |
US20170336370A1 (en) * | 2014-09-03 | 2017-11-23 | Shimadzu Corporation | Chromatogram data processing method and chromatogram data processing apparatus |
CN108898155A (en) * | 2018-05-18 | 2018-11-27 | 浙江工业大学 | A kind of SAR image target recognition method of Threshold Denoising combination convolutional neural networks |
CN109085282A (en) * | 2018-06-22 | 2018-12-25 | 东南大学 | A kind of chromatographic peaks analytic method based on wavelet transformation and Random Forest model |
US20210048417A1 (en) * | 2018-03-14 | 2021-02-18 | Hitachi High-Tech Corporation | Chromatography mass spectrometry and chromatography mass spectrometer |
CN112444589A (en) * | 2020-12-04 | 2021-03-05 | 深圳普门科技股份有限公司 | Chromatographic peak detection method, device, computer equipment and storage medium |
CN112730712A (en) * | 2020-12-31 | 2021-04-30 | 杭州谱育科技发展有限公司 | Method for improving LC-MS data signal-to-noise ratio |
WO2021106356A1 (en) * | 2019-11-27 | 2021-06-03 | アルプスアルパイン株式会社 | Chromatogram data processing device, chromatogram data processing method, chromatogram data processing program, and storage medium |
CN114280212A (en) * | 2022-01-10 | 2022-04-05 | 华谱科仪(北京)科技有限公司 | Chromatographic detection correction method, storage medium and electronic equipment |
CN114358502A (en) * | 2021-12-08 | 2022-04-15 | 北京工商大学 | Edible oil risk assessment method based on optimized ELM |
CN114814036A (en) * | 2022-05-09 | 2022-07-29 | 上海谱锐赛思生物技术有限公司 | Method for measuring concentration of azilsartan and amlodipine in blood plasma |
CN115586280A (en) * | 2022-10-31 | 2023-01-10 | 国电南瑞科技股份有限公司 | Chromatographic peak identification method based on self-adaptive threshold |
CN115700544A (en) * | 2022-10-25 | 2023-02-07 | 山东科技大学 | Chromatographic signal denoising method combining empirical mode decomposition and wavelet soft threshold |
-
2023
- 2023-04-20 CN CN202310427488.5A patent/CN116136518B/en active Active
Patent Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0694696A (en) * | 1992-09-17 | 1994-04-08 | Hitachi Ltd | Method for analyzing chromatogram and chromatographic device |
JP2006177980A (en) * | 2006-03-27 | 2006-07-06 | Hitachi Ltd | Chromatographic data processing device, chromatographic data processing method, and chromatographic analyzer |
CN101131380A (en) * | 2006-08-21 | 2008-02-27 | 安捷伦科技有限公司 | Method and apparatus for detecting trace substance using chromatographic instrument |
US20100161238A1 (en) * | 2007-05-31 | 2010-06-24 | Politecnico Di Milano | Method and System for Detecting Peptide Peaks in HPLC-MS Signals |
JP2009008582A (en) * | 2007-06-29 | 2009-01-15 | Shimadzu Corp | Chromatogram data processor |
US20120179389A1 (en) * | 2009-08-20 | 2012-07-12 | Spectrosense Ltd. | Gas Chromatographic Analysis Method and System |
US20110161053A1 (en) * | 2009-12-31 | 2011-06-30 | Industrial Technology Research Institute | Slope-based fast intrinsic mode functions decomposition method and apparatus |
US20140257712A1 (en) * | 2011-09-05 | 2014-09-11 | Shimadzu Corporation | System and method for processing chromatogram data |
US20150308992A1 (en) * | 2012-12-07 | 2015-10-29 | Shimadzu Corporation | Data processing system and data processing method for chromatograph |
CN103902844A (en) * | 2014-04-24 | 2014-07-02 | 国家电网公司 | Transformer vibration signal de-noising method based on EEMD kurtosis threshold value |
US20170336370A1 (en) * | 2014-09-03 | 2017-11-23 | Shimadzu Corporation | Chromatogram data processing method and chromatogram data processing apparatus |
CN104297504A (en) * | 2014-10-22 | 2015-01-21 | 上海申腾信息技术有限公司 | Automatic gas chromatographic control system |
US20210048417A1 (en) * | 2018-03-14 | 2021-02-18 | Hitachi High-Tech Corporation | Chromatography mass spectrometry and chromatography mass spectrometer |
CN108898155A (en) * | 2018-05-18 | 2018-11-27 | 浙江工业大学 | A kind of SAR image target recognition method of Threshold Denoising combination convolutional neural networks |
CN109085282A (en) * | 2018-06-22 | 2018-12-25 | 东南大学 | A kind of chromatographic peaks analytic method based on wavelet transformation and Random Forest model |
WO2021106356A1 (en) * | 2019-11-27 | 2021-06-03 | アルプスアルパイン株式会社 | Chromatogram data processing device, chromatogram data processing method, chromatogram data processing program, and storage medium |
CN112444589A (en) * | 2020-12-04 | 2021-03-05 | 深圳普门科技股份有限公司 | Chromatographic peak detection method, device, computer equipment and storage medium |
WO2022116673A1 (en) * | 2020-12-04 | 2022-06-09 | 深圳普门科技股份有限公司 | Chromatographic peak detection method and apparatus, computer device, and storage medium |
CN112730712A (en) * | 2020-12-31 | 2021-04-30 | 杭州谱育科技发展有限公司 | Method for improving LC-MS data signal-to-noise ratio |
CN114358502A (en) * | 2021-12-08 | 2022-04-15 | 北京工商大学 | Edible oil risk assessment method based on optimized ELM |
CN114280212A (en) * | 2022-01-10 | 2022-04-05 | 华谱科仪(北京)科技有限公司 | Chromatographic detection correction method, storage medium and electronic equipment |
CN114814036A (en) * | 2022-05-09 | 2022-07-29 | 上海谱锐赛思生物技术有限公司 | Method for measuring concentration of azilsartan and amlodipine in blood plasma |
CN115700544A (en) * | 2022-10-25 | 2023-02-07 | 山东科技大学 | Chromatographic signal denoising method combining empirical mode decomposition and wavelet soft threshold |
CN115586280A (en) * | 2022-10-31 | 2023-01-10 | 国电南瑞科技股份有限公司 | Chromatographic peak identification method based on self-adaptive threshold |
Non-Patent Citations (4)
Title |
---|
MUHAMMAD FAROOQ,WAHAB;THOMAS C,O\'HAVER: "Wavelet transforms in separation science for denoising and peak overlap detection.", J SEP SCI, vol. 43, no. 9, pages 1998 - 2010, XP055881254, DOI: 10.1002/jssc.202000013 * |
MUHAMMAD FAROOQ,WAHAB;THOMAS C,O\'HAVER: "Wavelet transforms in separation science for denoising and peak overlap detection.", JOURNAL OF SEPARATION SCIENCE * |
熊智新,胡上序,路文初: "基于小波变换的重叠色谱峰检测方法研究", 浙江大学学报(工学版), no. 08 * |
石松;刘升;: "高斯拟合结合小波降噪处理IMS峰信号", 淮北师范大学学报(自然科学版), no. 02 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116399996A (en) * | 2023-06-02 | 2023-07-07 | 海能未来技术集团股份有限公司 | Organic element analysis method and device based on thermal conductivity detection and electronic equipment |
Also Published As
Publication number | Publication date |
---|---|
CN116136518B (en) | 2023-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116136518B (en) | Chromatograph | |
CN105223482B (en) | The wavelet decomposition two-value denoising method of partial-discharge ultrahigh-frequency signal waveform | |
CN108918446B (en) | Ultra-low concentration sulfur dioxide ultraviolet difference feature extraction algorithm | |
CN113588847B (en) | Biological metabonomics data processing method, analysis method, device and application | |
CN112345759B (en) | Method for detecting fluorescence intensity peak | |
CN113325277A (en) | Partial discharge processing method | |
CN117030654B (en) | Method for measuring concentration of trace-level sulfur dioxide in air | |
CN111089856A (en) | Post-processing method for extracting Raman spectrum weak signal | |
CN116432051A (en) | Evaluation method for weak difference signals of near infrared spectrum data | |
CN115420726A (en) | Method for rapidly identifying target object by using reconstructed SERS spectrum | |
CN113686811A (en) | Spectral data processing method based on double sensors | |
CN105954223A (en) | Method for improving prediction accuracy of gasoline properties | |
CN114614825B (en) | Low-cost high-speed pulse signal data sampling and peak value detection method | |
CN114550843B (en) | Prediction model of monosaccharide composition and content in traditional Chinese medicine polysaccharide, construction method and application thereof | |
CN111210876A (en) | Disturbed metabolic pathway determination method and system | |
Ji et al. | Pure ion chromatogram extraction via optimal k-means clustering | |
CN110632191A (en) | Transformer chromatographic peak qualitative method and system based on decision tree algorithm | |
CN110553988B (en) | NO low-concentration signal extraction method based on superposition state idea | |
CN110308109B (en) | Method for distinguishing perchlorate content range in tea based on near infrared spectrum | |
CN113408383A (en) | Audible noise invalid data judgment method based on bounded beta (g, h) distribution and MWKPCA | |
CN118010649B (en) | Pollution detection method for food | |
CN112505779B (en) | Method for removing collected footprints based on feature decomposition | |
CN117825601B (en) | Method for measuring sulfur dioxide in food | |
CN117454095B (en) | Bridge dynamic deflection data noise reduction method | |
CN109917458B (en) | Oil gas analysis method and device |
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 | ||
CP03 | Change of name, title or address |
Address after: 310053 Room 301, 3 / F, building 1, 22 Zhiren street, Puyan street, Binjiang District, Hangzhou City, Zhejiang Province Patentee after: Hangzhou Zetian Chunlai Technology Co.,Ltd. Country or region after: China Address before: 310053 Room 301, 3 / F, building 1, 22 Zhiren street, Puyan street, Binjiang District, Hangzhou City, Zhejiang Province Patentee before: Hangzhou Zetian Chunlai Technology Co.,Ltd. Country or region before: China |