CN108982732A - A kind of efficient liquid phase glycosylated hemoglobin chromatographic peak area recognition methods - Google Patents

A kind of efficient liquid phase glycosylated hemoglobin chromatographic peak area recognition methods Download PDF

Info

Publication number
CN108982732A
CN108982732A CN201810966041.4A CN201810966041A CN108982732A CN 108982732 A CN108982732 A CN 108982732A CN 201810966041 A CN201810966041 A CN 201810966041A CN 108982732 A CN108982732 A CN 108982732A
Authority
CN
China
Prior art keywords
curve
gaussian curve
peak
exponential
peak value
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.)
Pending
Application number
CN201810966041.4A
Other languages
Chinese (zh)
Inventor
何刚
吴文青
俞文心
许康
董建明
叶长宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest University of Science and Technology
Original Assignee
Southwest University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southwest University of Science and Technology filed Critical Southwest University of Science and Technology
Priority to CN201810966041.4A priority Critical patent/CN108982732A/en
Publication of CN108982732A publication Critical patent/CN108982732A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating 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/02Column chromatography
    • G01N30/86Signal analysis
    • G01N30/8624Detection of slopes or peaks; baseline correction
    • G01N30/8631Peaks
    • G01N30/8634Peak quality criteria

Landscapes

  • Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

The present invention discloses a kind of efficient liquid phase glycosylated hemoglobin chromatographic peak area recognition methods, including data collection: being collected to glycosylated hemoglobin test data and removes impulse disturbances;Onset index Gaussian curve: by test data onset index Gaussian curve, and parameter in gauge index Gaussian curve;Carry out index gaussian curve approximation: the mathematical model for using the convolution of Gaussian curve and exponential decay curve as curved flow of chromatography is fitted the peak type of efficient liquid phase;Correct the index gaussian curve approximation.The present invention realizes the efficient identification of efficient liquid phase saccharification peak area, there are the abilities such as anti-impulse disturbances, Parameter adjustable, the amendment of automatic peak height in identification process, the accuracy of identification of liquid phase glycosylated hemoglobin chromatographic peak area is improved, high-precision targeting inspection result can be provided for the detection of glycosylated hemoglobin.

Description

High-efficiency liquid-phase glycosylated hemoglobin chromatographic peak area identification method
Technical Field
The invention belongs to the technical field of glycosylated hemoglobin identification, and particularly relates to a high-efficiency liquid-phase glycosylated hemoglobin chromatographic peak area identification method.
Background
Diabetes is increasingly threatening the health of humans, and monitoring efforts directed to diabetes treatment are also receiving increased attention from medical workers. Among parameters reflecting diabetes treatment indexes, glycated hemoglobin is more accepted and favored by authorities and institutions due to its stability.
Glycated hemoglobin assay has been used as a medical diagnostic, originated in the last eighties of the last seventies, and mainly measures the ability of glucose molecules in blood to attach to red blood cells, as an assessment of the quality of blood glucose level control in the past 2-3 months in diabetic patients. At present, methods clinically used for measuring the HbA1c content mainly comprise immunization, ion exchange chromatography, boron affinity chromatography, capillary electrophoresis and the like. The methods comprise a separation step before the photometric quantitative determination of the HbA1c content, such as ion exchange chromatography and electrophoretic separation based on the surface charge difference between HbA1c and Hb and immune and boron affinity chromatography based on the structure difference between HbA1c and Hb, the detection process needs multiple elution and special instrument detection, the detected data cannot be effectively identified, the detection period is greatly increased, the detection precision is reduced, and therefore, the application in the aspects of early warning of diabetes and medium-long term control of the blood sugar level of a diabetic patient is limited.
Disclosure of Invention
In order to solve the problems, the invention provides a high-efficiency liquid-phase glycosylated hemoglobin chromatographic peak area identification method, which realizes the high-efficiency identification of a high-efficiency liquid-phase glycosylated peak area by introducing Gaussian curve fitting, has the capabilities of pulse interference resistance, parameter adjustability, automatic peak height correction and the like in the identification process, improves the identification precision of the liquid-phase glycosylated hemoglobin chromatographic peak area, and can provide a high-precision target detection result for the detection of the glycosylated hemoglobin.
In order to achieve the purpose, the invention adopts the technical scheme that: a method for identifying the chromatographic peak area of high-efficiency liquid-phase glycated hemoglobin comprises the following steps:
s100, data collection: collecting test data of the glycosylated hemoglobin and removing pulse interference;
s200, establishing an exponential Gaussian curve: establishing an exponential Gaussian curve through test data, and calculating parameters in the exponential Gaussian curve;
s300, performing exponential Gaussian curve fitting: fitting the peak type of the high-performance liquid phase by taking the convolution of the Gaussian curve and the exponential decay curve as a mathematical model of a chromatographic outflow curve;
and S400, correcting the exponential Gaussian curve fitting.
Further, in step S100, test data is collected, and the collected test data (x) is recorded in the form of pairsi,yi),i=1,2,N。
Further, in step S200, an exponential gaussian curve is created by recording data, and the calculation formula of the exponential gaussian curve is:
wherein S is the area under the curve; μ is data xiI is the mean of 1,2, N,σ is data xiI is the standard deviation of 1,2, N,
further, parameters in the calculated exponential gaussian curve include the area under the curve, S, the mean of the test data, μ, and the standard deviation of the test data, σ.
Further, in step S200, the calculation process of the parameters in the exponential gaussian curve includes the steps of:
s201, taking natural logarithms at two sides of an exponential Gaussian curve calculation formula, including
S202, introducing a new variable Zi, and converting the formula obtained in S201 into a quadratic polynomial function:
wherein z isi=ln yi
S203, expressing the quadratic polynomial function into a matrix form to obtain a matrix function:
s204, according to a least square method, giving a solving parameter b0,b1,b2Is expressed as
Wherein,
s205, by the relational expressionμ=σ2b1And solving the parameter values in the exponential Gaussian curve.
Further, in step S300, exponential gaussian curve fitting is performed: taking the convolution of a Gaussian curve and an exponential decay curve as a mathematical model of a chromatogram outflow curve to fit the peak type of the high-performance liquid phase, wherein the fitting calculation formula is as follows:
wherein h (x) is a function of the outflow time t and the exponential decay curve isTau is the time constant of exponential decay curve, v is convolution integral variable, and any two curve functions fX(x) And fYThe convolution operation of (y) is defined as
Furthermore, overfitting and under-fitting can occur when saccharification peak fitting is performed, probably because the model parameters are greatly influenced by the change of the peak pattern; in order to solve the problem that the fitting result of the directly used glycated peak is not accurate due to the possible subtle differences in the peak patterns in the collection process caused by the various interferences, the exponential gaussian curve fitting is modified in step S400, which includes the steps of:
calculating a peak value fitted by an exponential Gaussian curve as a saccharification peak value;
comparing the saccharification peak value with the actual peak value; when the saccharification peak value is the same as the actual peak value, no correction is carried out; when the peak value is deviated from the actual value, correcting;
the modification process comprises the following steps: the calculated glycation peak is multiplied by a factor to equalize the peak with the actual peak.
Further, the formula for determining the comparison between the saccharification peak and the actual peak is as follows:
wherein,denotes that at time t ═ PACalculating the obtained saccharification peak value;denotes that at time t ═ PAActual peak value of true glycation peak;
at the peak of saccharificationAnd the actual peak valueWhen the same, adopting the expression of the upper half part;
when the saccharification peakAnd the actual peak valueAnd the lower half expression is adopted, the coefficient is automatically corrected, and automatic peak height correction is realized.
The beneficial effects of the technical scheme are as follows:
according to the method, the high-efficiency identification of the high-efficiency liquid phase glycosylated peak area is realized by introducing Gaussian curve fitting, the method has the capabilities of pulse interference resistance, parameter adjustability, automatic peak height correction and the like in the identification process, the identification precision of the liquid phase glycosylated hemoglobin chromatographic peak area is improved, and high-precision targeted inspection data can be provided for the detection of the glycosylated hemoglobin; provides reliable basis for early warning of diabetes and medium-and-long-term control of blood sugar level of diabetic patients;
the invention improves the accuracy of peak identification by removing pulse interference; the peak height can be automatically corrected based on the corrected Gaussian function closure, so that a high-precision numerical result is obtained.
Drawings
FIG. 1 is a schematic flow chart of a method for identifying the chromatographic peak area of high performance liquid glycated hemoglobin according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described with reference to the accompanying drawings.
In this embodiment, referring to fig. 1, the present invention provides a method for identifying a chromatographic peak area of a high performance liquid glycated hemoglobin, comprising the steps of:
s100, data collection: collecting test data of the glycosylated hemoglobin and removing pulse interference;
s200, establishing an exponential Gaussian curve: establishing an exponential Gaussian curve through test data, and calculating parameters in the exponential Gaussian curve;
s300, performing exponential Gaussian curve fitting: fitting the peak type of the high-performance liquid phase by taking the convolution of the Gaussian curve and the exponential decay curve as a mathematical model of a chromatographic outflow curve;
and S400, correcting the exponential Gaussian curve fitting.
As an optimization of the above embodiment, in step S100, the test data is collected as a numberForm of pairs the test data (x) collected was recordedi,yi),i=1,2,N。
As an optimization scheme of the above embodiment, in step S200, an exponential gaussian curve is created by recording data, and the calculation formula of the exponential gaussian curve is as follows:
wherein S is the area under the curve; μ is data xiI is the mean of 1,2, N,σ is data xiI is the standard deviation of 1,2, N,
parameters in the calculated exponential gaussian curve include the area under the curve, S, the mean value of the test data, mu, and the standard deviation of the test data, sigma.
In step S200, the process of calculating parameters in the exponential gaussian curve includes the steps of:
s201, taking natural logarithms at two sides of an exponential Gaussian curve calculation formula, including
S202, introducing a new variable Zi, and converting the formula obtained in S201 into a quadratic polynomial function:
wherein z isi=ln yi
S203, expressing the quadratic polynomial function into a matrix form to obtain a matrix function:
s204, according to a least square method, giving a solving parameter b0,b1,b2Is expressed as
Wherein,
s205, by the relational expressionμ=σ2b1And solving the parameter values in the exponential Gaussian curve.
As an optimization of the above embodiment, in step S300, exponential gaussian curve fitting is performed: taking the convolution of a Gaussian curve and an exponential decay curve as a mathematical model of a chromatogram outflow curve to fit the peak type of the high-performance liquid phase, wherein the fitting calculation formula is as follows:
wherein h (x) is a function of the outflow time t and the exponential decay curve isTau is the time constant of exponential decay curve, v is convolution integral variable, and any two curve functions fX(x) And fYThe convolution operation of (y) is defined as
When saccharification peak fitting is carried out, overfitting, under-fitting and the like can occur because the model parameters are greatly influenced by the change of the peak pattern; in order to solve the problem that the fitting result of the directly used glycated peak is not accurate due to the possible subtle differences in the peak patterns in the collection process caused by the various interferences, the exponential gaussian curve fitting is modified in step S400, which includes the steps of:
calculating a peak value fitted by an exponential Gaussian curve as a saccharification peak value;
comparing the saccharification peak value with the actual peak value; when the saccharification peak value is the same as the actual peak value, no correction is carried out; when the peak value is deviated from the actual value, correcting;
the modification process comprises the following steps: the calculated glycation peak is multiplied by a factor to equalize the peak with the actual peak.
Wherein, the judgment formula for comparing the saccharification peak value with the actual peak value is as follows:
wherein,denotes that at time t ═ PACalculating the obtained saccharification peak value;denotes that at time t ═ PATrue saccharification peakA cross-peak value;
at the peak of saccharificationAnd the actual peak valueWhen the same, adopting the expression of the upper half part;
when the saccharification peakAnd the actual peak valueAnd the lower half expression is adopted, the coefficient is automatically corrected, and automatic peak height correction is realized.
The foregoing shows and describes the general principles and broad features of the present invention and advantages thereof. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (8)

1. A method for identifying the chromatographic peak area of high-efficiency liquid-phase glycated hemoglobin is characterized by comprising the following steps:
s100, data collection: collecting test data of the glycosylated hemoglobin and removing pulse interference;
s200, establishing an exponential Gaussian curve: establishing an exponential Gaussian curve through test data, and calculating parameters in the exponential Gaussian curve;
s300, performing exponential Gaussian curve fitting: fitting the peak type of the high-performance liquid phase by taking the convolution of the Gaussian curve and the exponential decay curve as a mathematical model of a chromatographic outflow curve;
and S400, correcting the exponential Gaussian curve fitting.
2. The method of claim 1, wherein the test data is collected and the collected test data (x) is recorded in the form of a plurality of pairs in step S100i,yi),i=1,2,N。
3. The method of claim 2, wherein in step S200, an exponential gaussian curve is created by recording data, wherein the exponential gaussian curve is calculated by the formula:
wherein S is the area under the curve; μ is data xiI is the mean of 1,2, N,σ is data xiI is the standard deviation of 1,2, N,
4. the method of claim 3, wherein the parameters of the exponential Gaussian curve include an area S under the curve, a mean μ of the test data, and a standard deviation σ of the test data.
5. The method of claim 4, wherein the step of calculating parameters of the exponential Gaussian curve in step S200 comprises the steps of:
s201, taking natural logarithms at two sides of an exponential Gaussian curve calculation formula, including
S202, introducing a new variable Zi, and converting the formula obtained in S201 into a quadratic polynomial function:
wherein z isi=lnyi
S203, expressing the quadratic polynomial function into a matrix form to obtain a matrix function:
s204, according to a least square method, giving a solving parameter b0,b1,b2Is expressed as
Wherein,
s205, by the relational expressionμ=σ2b1And solving the parameter values in the exponential Gaussian curve.
6. The method of claim 5, wherein in step S300, an exponential Gaussian curve fitting is performed: taking the convolution of a Gaussian curve and an exponential decay curve as a mathematical model of a chromatogram outflow curve to fit the peak type of the high-performance liquid phase, wherein the fitting calculation formula is as follows:
wherein h (x) is a function of the outflow time t and the exponential decay curve isTau is the time constant of exponential decay curve, v is convolution integral variable, and any two curve functions fX(x) And fYThe convolution operation of (y) is defined as
7. The method of claim 6, wherein the step S400 of modifying the exponential Gaussian curve fit comprises the steps of:
calculating a peak value fitted by an exponential Gaussian curve as a saccharification peak value;
comparing the saccharification peak value with the actual peak value; when the saccharification peak value is the same as the actual peak value, no correction is carried out; when the peak value is deviated from the actual value, correcting;
the modification process comprises the following steps: the calculated glycation peak is multiplied by a factor to equalize the peak with the actual peak.
8. The method of claim 7, wherein the formula for comparing the glycated peak with the actual peak is:
wherein,denotes that at time t ═ PACalculating the obtained saccharification peak value;denotes that at time t ═ PAActual peak value of true glycation peak;
at the peak of saccharificationAnd the actual peak valueWhen the same, adopting the expression of the upper half part;
when the saccharification peakAnd the actual peak valueAnd the lower half expression is adopted to automatically correct the coefficient.
CN201810966041.4A 2018-08-23 2018-08-23 A kind of efficient liquid phase glycosylated hemoglobin chromatographic peak area recognition methods Pending CN108982732A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810966041.4A CN108982732A (en) 2018-08-23 2018-08-23 A kind of efficient liquid phase glycosylated hemoglobin chromatographic peak area recognition methods

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810966041.4A CN108982732A (en) 2018-08-23 2018-08-23 A kind of efficient liquid phase glycosylated hemoglobin chromatographic peak area recognition methods

Publications (1)

Publication Number Publication Date
CN108982732A true CN108982732A (en) 2018-12-11

Family

ID=64548046

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810966041.4A Pending CN108982732A (en) 2018-08-23 2018-08-23 A kind of efficient liquid phase glycosylated hemoglobin chromatographic peak area recognition methods

Country Status (1)

Country Link
CN (1) CN108982732A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109682907A (en) * 2019-01-14 2019-04-26 华中科技大学 A method of it is demanded perfection the first dimension chromatographic peak in two-dimensional chromatography by two modulation peaks
CN109917061A (en) * 2019-03-11 2019-06-21 嘉兴迈维代谢生物科技有限公司 A kind of method of quick point in wide target metabolism group
CN112444589A (en) * 2020-12-04 2021-03-05 深圳普门科技股份有限公司 Chromatographic peak detection method, device, computer equipment and storage medium
CN113419020A (en) * 2021-06-30 2021-09-21 成都师范学院 Glycated hemoglobin overlapping peak recognition method, apparatus, system, device, and medium
CN114002367A (en) * 2021-10-15 2022-02-01 无锡博慧斯生物医药科技有限公司 Algorithm of chromatographic peak area of high-efficiency liquid-phase glycosylated hemoglobin
CN114431856A (en) * 2022-01-28 2022-05-06 上海乾康医疗设备股份有限公司 Neural feedback rehabilitation training system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10170492A (en) * 1996-12-09 1998-06-26 Jasco Corp Apparatus for measuring time-developed spectrum
CN1712955A (en) * 2004-06-25 2005-12-28 中国科学院大连化学物理研究所 Precisive measurement for parameter of chromatography spike and area of overlapped peak
CN103389352A (en) * 2013-07-25 2013-11-13 上海申瑞继保电气有限公司 Method for smoothing chromatographic data of oil gas
CN105067732A (en) * 2015-07-07 2015-11-18 天津大学 Overlapped chromatographic peak quantitative method based on second-order Gaussian function fitting
CN106442830A (en) * 2016-09-29 2017-02-22 广州供电局有限公司 Method and system for detecting alarm value of gas content of transformer oil
CN106770861A (en) * 2016-11-21 2017-05-31 浙江大学 The evaluation method of oil-filled transformer on-line monitoring availability of data

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10170492A (en) * 1996-12-09 1998-06-26 Jasco Corp Apparatus for measuring time-developed spectrum
CN1712955A (en) * 2004-06-25 2005-12-28 中国科学院大连化学物理研究所 Precisive measurement for parameter of chromatography spike and area of overlapped peak
CN103389352A (en) * 2013-07-25 2013-11-13 上海申瑞继保电气有限公司 Method for smoothing chromatographic data of oil gas
CN105067732A (en) * 2015-07-07 2015-11-18 天津大学 Overlapped chromatographic peak quantitative method based on second-order Gaussian function fitting
CN106442830A (en) * 2016-09-29 2017-02-22 广州供电局有限公司 Method and system for detecting alarm value of gas content of transformer oil
CN106770861A (en) * 2016-11-21 2017-05-31 浙江大学 The evaluation method of oil-filled transformer on-line monitoring availability of data

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"简单但不平凡"(C914620529): "《CSDN博客》", 25 December 2015 *
AGILENT TECHNOLOGIES, INC.: "《Technical Note》", 31 December 2005 *
涂佳丽 等: "基于液相色谱的糖化血红蛋白浓度的检测", 《微计算机信息》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109682907A (en) * 2019-01-14 2019-04-26 华中科技大学 A method of it is demanded perfection the first dimension chromatographic peak in two-dimensional chromatography by two modulation peaks
CN109917061A (en) * 2019-03-11 2019-06-21 嘉兴迈维代谢生物科技有限公司 A kind of method of quick point in wide target metabolism group
CN112444589A (en) * 2020-12-04 2021-03-05 深圳普门科技股份有限公司 Chromatographic peak detection method, device, computer equipment and storage medium
CN112444589B (en) * 2020-12-04 2021-10-08 深圳普门科技股份有限公司 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
CN113419020A (en) * 2021-06-30 2021-09-21 成都师范学院 Glycated hemoglobin overlapping peak recognition method, apparatus, system, device, and medium
CN114002367A (en) * 2021-10-15 2022-02-01 无锡博慧斯生物医药科技有限公司 Algorithm of chromatographic peak area of high-efficiency liquid-phase glycosylated hemoglobin
CN114002367B (en) * 2021-10-15 2023-05-09 无锡博慧斯生物医药科技有限公司 Algorithm for high performance liquid glycosylated hemoglobin chromatographic peak area
CN114431856A (en) * 2022-01-28 2022-05-06 上海乾康医疗设备股份有限公司 Neural feedback rehabilitation training system

Similar Documents

Publication Publication Date Title
CN108982732A (en) A kind of efficient liquid phase glycosylated hemoglobin chromatographic peak area recognition methods
Lin et al. Effects of hemoglobin C, D, E, and S traits on measurements of HbA1c by six methods
WO1997020496A9 (en) Medical diagnostic analysis system
CN1954207A (en) Methods for performing hematocrit adjustment in glucose assays and devices for same
EP3221693B1 (en) A method and apparatus for performing liquid chromatography purification
Marinova et al. Multicenter evaluation of hemoglobin A1c assay on capillary electrophoresis
US20100036791A1 (en) Examination value predicting device using electrophoresis waveform, prediction method, and predicting program
CN110444294B (en) Auxiliary analysis method and equipment for prostate cancer based on perception neural network
JP2017519226A (en) Hemolysis detection method and system
CN109030803B (en) Biochemical detection quality control method
CN112083175A (en) Method for correcting biochemical index by hemolytic index measured by full-automatic dry biochemical analyzer
CN115691722A (en) Quality control method, apparatus, device, medium and program product for medical data detection
Ucar et al. Estimation of biological variation and reference change value of glycated hemoglobin (HbA1c) when two analytical methods are used
CN113470753A (en) Primary central nervous system lymphoma prognosis model establishment method based on albumin and ECOG-PS and application
CN108268752A (en) A kind of chromosome abnormality detection device
CN115101171A (en) Patient information analysis management system for hemodialysis
Nelson et al. Simplified micro-scale procedure for preparing samples for theophylline determination by liquid chromatography.
CN103760159B (en) A kind of method and system of Bacteria Identification and Analysis of Drug Susceptibility
CN112816425B (en) Method for optimizing whole blood sample detection flow by utilizing HGB calibration capability
JP7081385B2 (en) How to calculate the similarity of chromatograms
US5384239A (en) Method for analyzing the glycation of hemoglobin
Roselli et al. Analysis of noninvasive macromolecular transport measurements in the lung
CN109001354B (en) Wave crest identification method and device, chromatographic analyzer and storage medium
CN110841152A (en) Heatable pressure device for blood transfusion of hematology department
Juhola The effect of digital lowpass filters on the maximum velocity of saccadic eye movements

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20181211

RJ01 Rejection of invention patent application after publication