CN1544917A - Spectrum quantitative automatic analysis method - Google Patents
Spectrum quantitative automatic analysis method Download PDFInfo
- Publication number
- CN1544917A CN1544917A CNA2003101087339A CN200310108733A CN1544917A CN 1544917 A CN1544917 A CN 1544917A CN A2003101087339 A CNA2003101087339 A CN A2003101087339A CN 200310108733 A CN200310108733 A CN 200310108733A CN 1544917 A CN1544917 A CN 1544917A
- Authority
- CN
- China
- Prior art keywords
- quantitative
- wavelength
- value
- parameter
- related coefficient
- 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
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention is an optical spectrum quantitative automatic analysis method, firstly determining an optical spectrum quantitative automatic analysis method of the four quantifying methods: non-normalization wavelength quantifying, normalization wavelength quantifying, area quantifying and speed quantifying and the quantitative parameters, then in the automatically ergodic circulating course, calculating light measurement the current values of the quantitative parameters, making linear fitting on the light measurement and measured quantity by using least squares method, recording all the quantitative parameters and corresponding related coefficients and finding the maximum related coefficient, and after the automatic ergodic course is finished, obtaining a 'related coefficient-quantitative parameter relation chart' and the suggested optimum quantitative parameter. It provides an automatic analysis method for the quantitative analysis result of optical spectrum, can evaluate quantitative measuring ability of the instrument for the special sample on given condition, gives the optimum quantitative analysis parameter and range and provides the reference for the users of spectral photometer and spectrograph.
Description
Technical field:
The present invention relates to the quantitative automatic analysis method of a kind of spectrum, can be used for evaluation appts, determine best quantitative test parameter and scope thereof automatically, realize that the sample concentration of optimum efficiency is measured the quantitative measurment ability of specific sample under specified criteria.Belong to the spectrophotometry field.
Background technology:
Quantitative test is one of important composition of spectrophotometry.A fundamental analysis content wherein is exactly the concentration that spectrum (absorbance) from sample is extrapolated sample, and its theoretical foundation is Beer law.
In real work, developed and multiple concrete grammar, as wavelength quantitatively, the normalization wavelength is quantitative, area is quantitative, speed is quantitative etc.Wavelength quantitatively is to select a wavelength points from the overall optical spectral curve of sample, and at this wavelength place, the absorbance and the concentration of sample have good linear relationship.This relation can be determined with linear fit.The normalization wavelength quantitatively selects a wavelength points that spectrum is carried out normalization earlier, selects another wavelength points then and makes linear fit.Area quantitatively then will be selected two wavelength points in the curve of spectrum, area and sample concentration that curve is enclosed between these two wavelength points have good linear relationship.Slightly different with first three methods, what speed was quantitatively utilized is the absorbance rate of change curve of noting in the example reaction process, chooses a proper time point on rate curve.At this time point place, the absorbance rate of change of sample and the concentration of sample have the good linear relation.Linear fit in these quantivative approachs often uses least square method.For convenience of explanation, the area and the absorbance rate of change of the absorbance in these methods, absorbance curve all is called light quantity, and selected wavelength points or time point are called quantitative parameter when these methods are used.Like this, wavelength quantitatively and the speed sizing technique all have only a quantitative parameter, and the normalization wavelength is quantitative and the area sizing technique need be determined two quantitative parameters.
Under complete ideal situation, quantitative parameter can be selected arbitrarily in it allows span.But during actual measurement, multiple factor, as catoptrical existence, the restriction of the precision of instrument own, the electrochemical conditions of sample etc. all can be to the light quantity and the generation of the linear relationship between the concentration influence in various degree of sample.Like this, when quantitative parameter is got some value, still have good linear relationship between light quantity and concentration, and under other values, can not obtain quantitative result accurately.Therefore, selecting suitable quantitative parameter is the steps necessary of high-quality quantitative spectrochemical analysis.
In classic method, the loaded down with trivial details manually-operated of definite needs of quantitative parameter.The analyst will attempt selecting for use different quantitative parameters repeatedly, selecting the back to calculate light quantity each time, light quantity and concentration is carried out match, investigates the quality of its linear degree.Though can utilize some existing practical experiences to help select quantitative parameter, when using wavelength quantitative, can attempt selecting a certain absorption peak wavelength of sample to carry out measurement of concetration usually, this is still a process that wastes time and energy.And because the work quantitative limitation, the number of times of trial also is limited, can't guarantee that the quantitative parameter of determining always can reach best quantitative test effect.
Summary of the invention:
The objective of the invention is at the deficiencies in the prior art and actual needs, provide a kind of spectrum quantitative automatic analysis method, can replace loaded down with trivial details manual method and definite automatically quantitative parameter, measure with the sample concentration of realizing optimum efficiency.
For realizing such purpose, the key of the inventive method is the possible value of institute that travels through quantitative parameter automatically.Under each value, automatically finishing light quantity calculates, use least square fitting light quantity and the measured linear relationship of sample, related coefficient with least square method is the quality that index is passed judgment on linear degree, provides quantitative parameter and " related coefficient～quantitative parameter graph of a relation " an of the best at last.Whole process can utilize computing machine to realize automatically.
Method of the present invention comprises following concrete steps:
1, determines quantitative analysis method of spectrum and quantitative parameter.From non-normalized wavelength quantitatively, the normalization wavelength is quantitative, area is quantitative and the quantitative four kinds of quantivative approachs of speed select a kind of quantitative test that is used for, also determined employed quantitative parameter simultaneously.Enter the cyclic process of automatic traversal then.
Method of the present invention at first provides four kinds of default quantivative approachs, comprise non-normalized wavelength quantitatively, the normalization wavelength is quantitative, area is quantitative and speed is quantitative, can select wherein a kind of on demand.When determining quantivative approach, quantitative parameter and span thereof that this method is required have also just been determined.The quantitative quantitative parameter of non-normalized wavelength is a wavelength points, and span is the wavelength coverage of the curve of spectrum.The normalization wavelength quantitatively and area quantitatively all need two wavelength points as quantitative parameter, span also is the wavelength coverage of the curve of spectrum.The quantitative quantitative parameter of speed is a time point, and span is the writing time of example reaction process.
2, enter the circulation of automatic traversal after, round-robin begins each time, in possible the span of quantitative parameter automatically value of not attempting as yet of selection as currency.Automatically sampling process is such: non-normalized wavelength quantitatively and speed quantitatively all have only a quantitative parameter, value is from the minimum value of allowed band, increase successively.The normalization wavelength quantitatively and the area quantivative approach all need two quantitative parameters, then select first wavelength earlier, first wavelength also is that the minimum value from allowed band begins value.Then, the value of second wavelength increases one by one from the currency of first wavelength.And before the possible value of second wavelength was not all attempted as yet, the value of first wavelength immobilized.
3, after the value of quantitative parameter in current circulation determined, the currency of then pressing quantitative parameter calculated light quantity.The quantitative luminosity value of non-normalized wavelength is exactly the absorbance of current wavelength correspondence.The normalization wavelength quantitatively earlier at first wavelength points place to curve of spectrum normalization, the absorbance (after the normalization) of getting second wavelength place correspondence again is as the luminosity value.Area quantitatively then will calculate the area that the curve of spectrum is enclosed between two current wavelength points, as the value of light quantity.The absorbance rate of change of the current point in time correspondence that speed is quantitative is exactly a light quantity.
4, after light quantity is calculated and finished, with least square method light quantity and sample concentration are carried out linear fit, and calculate the related coefficient of least square method.Related coefficient can reflect the quality of the linear degree of match.Current quantitative parameter and corresponding related coefficient are all noted, and the related coefficient that obtains in current related coefficient and the circulation before relatively, found out maximum related coefficient and note.
5, judge whether current quantitative parameter value has traveled through the possible value of its institute,, then continue circulation up to finishing whole ergodic process if quantitative parameter still has new value to attempt.Quantitative for non-normalized wavelength, quantitative parameter is as having got maximum wavelength points, but end loop.For the normalization wavelength quantitatively and area quantitative, if two wavelength have all been got the wavelength of maximum, expression can end loop.And in speed was quantitative, the sign of loop ends was to have attempted maximum time point.Automatically after traversal finishes,, will find the related coefficient of a maximum owing to the work in the step 4.The quantitative parameter value of this related coefficient correspondence is exactly the quantitative parameter of this kind sample required the best when using this method.Simultaneously, the value of each quantitative parameter and corresponding related coefficient have automatically recorded in step 4 also along with ergodic process.
6, display result.After having finished step 5, just obtained a series of and corresponding related coefficient of each value quantitative parameter, and the related coefficient of a maximum.The pairing quantitative parameter value of this maximum related coefficient can make the quantitative test of current sample reach optimum.Thereby, just absorbance data to be done once to analyze with this quantitative parameter, the fitting a straight line that obtains shows in the mode of mapping.Simultaneously, corresponding quantitative parameter value is mapped, obtain " related coefficient～quantitative parameter graph of a relation " with related coefficient.
On " related coefficient～quantitative parameter graph of a relation ", the analyst can further observe the influence of quantitative parameter value to quantitative analysis results, also the effect that can carry out the concentration quantitative measurment as certain test method of reference analysis or testing equipment.
The present invention is applicable to the various single component surveying works that obtain the spectral class instrument of continuous spectrum curve.Determine best quantitative parameter value for the spectrophotometric analysis personnel provide quick, an easy method, and provide auxiliary information for the measurement effect of evaluation appts.The analyst only need prepare some sample solutions that are used to calibrate and measure, and uses this method to get final product after obtaining all spectroscopic datas.The analyst can directly use the best quantitative parameter of suggestion, or serves as to select targetedly and revise with reference to the value of quantitative parameter is made with " related coefficient～quantitative parameter graph of a relation ", obtains best quantitative analysis results.
The present invention replaces manually-operated with the automation process of computing machine, has significantly reduced analysis time and workload, also is difficult for producing the situation of careless omission, and the result is more accurate, directly perceived.
Description of drawings:
Fig. 1 is the process flow diagram of the inventive method.
Fig. 2 is the protein solution spectral line that has adopted the variable concentrations that the small-sized hypervelocity ultraviolet-visual spectrometer of the inventive method measures.
Fig. 3 is " related coefficient～wavelength relationship figure " that obtains after analyzing.
Fig. 4 is the best quantitatively fitting quantitative line at wavelength place that obtains after analyzing.
Embodiment:
Below in conjunction with accompanying drawing and specific embodiment technical scheme of the present invention is further described.
Method flow of the present invention as shown in Figure 1.
As shown in Figure 1, the flow process of this method is divided into 6 steps.Wherein, step 2,3,4 and 5 has constituted the automatic ergodic process of key of the present invention.Below, to the in addition simple declaration respectively of 6 steps:
1. selection quantivative approach.A kind of by analyst's selection in 4 kinds of default methods, also just determined quantitative parameter type, quantity and the span that method is required simultaneously.
2. the automatic value of quantitative parameter.Method is quantitative parameter value in allowed band automatically.
3. light quantity is calculated.Under the current value of quantitative parameter, calculate light quantity.
4. least square fitting is sought best related coefficient.Carry out linear fit with least square method, calculate related coefficient simultaneously.Write down this related coefficient and current quantitative parameter value, and find out maximum related coefficient.
5. traversal control automatically.When the control cyclic process of traversal automatically finishes, and standard is that quantitative parameter has used the possible value of institute.
6. display result.Automatically after traversal finishes, obtained the record of a series of quantitative parameter values and corresponding related coefficient, and the related coefficient of a maximum and corresponding quantitative parameter value, also promptly best quantitative parameter value.To the mapping of quantitative parameter value, obtain intuitively " related coefficient～quantitative parameter graph of a relation " with related coefficient.With the best quantitative parameter value that method is determined absorbance data is done once to analyze, fitting a straight line also shows with the mapping mode.
Embodiment:
In the spectral analysis module of small-sized hypervelocity ultraviolet-visual spectrometer, used the present invention.The concentration that this spectrometer measurement is obtained is that the spectroscopic data of the protein solution of 1.000mol/L, 0.500mol/L, 0.250mol/L, 0.125mol/L carries out quantitative test.The spectrogram of the protein solution of four kinds of concentration is seen accompanying drawing 2.
Use method of the present invention to carry out the non-normalized wavelength of protein solution when quantitative, help to select quantitative wavelength at this spectrometer.
Implementation step is as follows:
(1) selects quantivative approach.The curve of spectrum of protein solution is seen Fig. 2.Because the spectrum of protein solution has tangible absorption peak, there is not background absorption substantially, so selected non-normalized wavelength quantivative approach.Simultaneously just determined that also quantitative parameter has only one, type is a wavelength.Its span is determined from absorbance data, is 2048 wavelength points between 190nm～820nm;
(2) the automatic value of quantitative parameter.Wavelength begins value from the 190nm of minimum, increases one by one.If the sequence number with wavelength points is represented, is exactly since 1, increase by 1 each time.
(3) light quantity is calculated.The quantitative light quantity of non-normalized wavelength is exactly an absorbance.Can directly look into the absorbance of getting current wavelength value correspondence from the spectrum of four protein solutions.
(4) linear fit.Light quantity has obtained four and with least square method these absorbances has been carried out linear fit to sample concentration corresponding to variable concentrations solution absorbency value after calculating and finishing, and calculates related coefficient.Related coefficient and current wavelength value are noted, and sought out maximum related coefficient.Related coefficient is big more, shows that the linear degree of the data that are used for match is good more, thereby the result of quantitative test more accurately and reliably;
(5) traversal control automatically.The purpose of automatic traversal control is to make quantitative parameter use the possible value of institute, and does not omit.The quantitative wavelength value of non-normalized wavelength is carried out from small to large successively, only needs to judge whether wavelength value has used 820nm.If represent, then check whether used 2048 wavelength points with the wavelength points sequence number;
(6) display result.With wavelength value of preserving in the step 4 and the mapping of related coefficient record, obtain " related coefficient～wavelength relationship figure ".The wavelength value of the maximum correlation coefficient correspondence that finds, the optimal wavelength that with non-normalized wavelength quantivative approach protein solution is carried out quantitative test exactly.Carry out a quantitative test at this wavelength place, the fitting a straight line that the obtains demonstration of also mapping.
Non-normalized quantitatively " related coefficient～wavelength relationship figure " sees Fig. 3.See that from figure 220nm～280nm and 370nm～420nm wave band has higher linear degree.Compare with spectrogram, 370nm～420nm wave band is place, place, absorption of sample peak just, and this often selects the absorption peak quantitative test for use with classic method is consistent.The best quantitatively wavelength of the suggestion that method provides is 227nm, and the fitting a straight line at this wavelength place is seen Fig. 4.Can see that on scheming the coincidence degree of data point and fitting a straight line is very good, the quantitative measurment of carrying out concentration with this relation can obtain result accurately and reliably.
By using the result as seen, the present invention can help the user to finish quantitative test to sample quickly and easily, has greatly improved work efficiency.
Claims (1)
1. automatic analysis method that spectrum is quantitative is characterized in that may further comprise the steps:
1) determines quantitative analysis method of spectrum and quantitative parameter: quantitative from default non-normalized wavelength, the normalization wavelength is quantitative, selection is a kind of quantitatively and in the quantitative four kinds of quantivative approachs of speed is used for quantitative test for area, the quantitative quantitative parameter of non-normalized wavelength is a wavelength points, span is the wavelength coverage of the curve of spectrum, the normalization wavelength quantitatively and area quantitatively all need two wavelength points as quantitative parameter, span also is the wavelength coverage of the curve of spectrum, the quantitative quantitative parameter of speed is a time point, span is the writing time of example reaction process, enters the cyclic process of automatic traversal then;
2) enter the circulation of automatic traversal after, round-robin begins each time, in the possible span of quantitative parameter, select a value of not attempting as yet as currency automatically, the value that non-normalized wavelength is quantitative and speed is quantitative is from the minimum value of allowed band, increase successively, quantitative and then selected earlier first wavelength of area quantivative approach of normalization wavelength, begin value from the minimum value of allowed band, the value of second wavelength is from the currency of first wavelength, increase one by one, before the possible value of second wavelength was not all attempted as yet, the value of first wavelength immobilized;
3) after the value of quantitative parameter in current circulation determined, the currency of pressing quantitative parameter calculates light quantity, the quantitative luminosity value of non-normalized wavelength is exactly the absorbance of current wavelength correspondence, the normalization wavelength quantitatively earlier at first wavelength points place to curve of spectrum normalization, get absorbance after the normalization of second wavelength place correspondence again as the luminosity value, area quantitatively then will calculate area that the curve of spectrum the enclosed value as light quantity between two current wavelength points, the absorbance rate of change of the current point in time correspondence that speed is quantitative is exactly a light quantity;
4) after light quantity is calculated and is finished, with least square method light quantity and sample concentration are carried out linear fit, and the related coefficient of calculating least square method, current quantitative parameter and corresponding related coefficient are all noted, and the related coefficient that obtains in current related coefficient and the circulation before relatively, finding out maximum related coefficient notes;
5) judge whether current quantitative parameter value has traveled through the possible value of its institute, if quantitative parameter still has new value to attempt, then continue circulation up to finishing whole ergodic process, quantitative for non-normalized wavelength, quantitative parameter is as having got maximum wavelength points, but end loop, for the normalization wavelength quantitatively and area quantitative, can end loop if two wavelength have all been got maximum wavelength, speed quantitatively in, the sign of loop ends is to have attempted maximum time point;
6) after traversal finishes automatically, with the corresponding related coefficient mapping of each value of a series of and quantitative parameter that obtains, obtain " related coefficient～quantitative parameter graph of a relation ", the best quantitative parameter of maximum correlation coefficient correspondence is done once to analyze to absorbance data, and the fitting a straight line that obtains shows in the mode of mapping.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200310108733 CN1235035C (en) | 2003-11-20 | 2003-11-20 | Spectrum quantitative automatic analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200310108733 CN1235035C (en) | 2003-11-20 | 2003-11-20 | Spectrum quantitative automatic analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1544917A true CN1544917A (en) | 2004-11-10 |
| CN1235035C CN1235035C (en) | 2006-01-04 |
Family
ID=34334838
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 200310108733 Expired - Fee Related CN1235035C (en) | 2003-11-20 | 2003-11-20 | Spectrum quantitative automatic analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1235035C (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101275912B (en) * | 2008-05-08 | 2010-08-25 | 中国农业大学 | Fluid food browning testing method |
| CN104568173A (en) * | 2013-10-29 | 2015-04-29 | 南开大学 | CCD wavelength calibrating method |
| CN106124420A (en) * | 2016-07-07 | 2016-11-16 | 威凯检测技术有限公司 | A kind of spectral detection instrument data analysis and information synthesis method |
| CN109187392A (en) * | 2018-09-26 | 2019-01-11 | 中南大学 | A kind of zinc liquid trace metal ion concentration prediction method based on two-zone model |
| CN110044477A (en) * | 2019-04-24 | 2019-07-23 | 中国人民解放军战略支援部队航天工程大学 | A kind of luminosity data searching method with similar spatial observation Geometrical change rule |
-
2003
- 2003-11-20 CN CN 200310108733 patent/CN1235035C/en not_active Expired - Fee Related
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101275912B (en) * | 2008-05-08 | 2010-08-25 | 中国农业大学 | Fluid food browning testing method |
| CN104568173A (en) * | 2013-10-29 | 2015-04-29 | 南开大学 | CCD wavelength calibrating method |
| CN106124420A (en) * | 2016-07-07 | 2016-11-16 | 威凯检测技术有限公司 | A kind of spectral detection instrument data analysis and information synthesis method |
| CN109187392A (en) * | 2018-09-26 | 2019-01-11 | 中南大学 | A kind of zinc liquid trace metal ion concentration prediction method based on two-zone model |
| CN110044477A (en) * | 2019-04-24 | 2019-07-23 | 中国人民解放军战略支援部队航天工程大学 | A kind of luminosity data searching method with similar spatial observation Geometrical change rule |
| CN110044477B (en) * | 2019-04-24 | 2021-05-18 | 中国人民解放军战略支援部队航天工程大学 | Photometric data searching method with similar space observation geometric change rule |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1235035C (en) | 2006-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109799203B (en) | Wide-range high-precision spectrum detection method for COD concentration in water body | |
| CN103175805B (en) | Method for determining indexes of COD and BOD5 in sewage through near infrared spectrometry | |
| CN111965140B (en) | Wavelength point recombination method based on characteristic peak | |
| CN108982406A (en) | A kind of soil nitrogen near-infrared spectral characteristic band choosing method based on algorithm fusion | |
| CN106770058A (en) | The quick special purpose device and its application method of the soil nitrate-N based on infrared spectrum | |
| CN106018331B (en) | The method for estimating stability and pretreatment optimization method of multi-channel spectral system | |
| CN103472028A (en) | Method for determining content of components in blended fiber | |
| CN108896507B (en) | Method for estimating river humification index | |
| CA2469099A1 (en) | Spectroscopic method and apparatus for total hemoglobin measurement | |
| CN1235035C (en) | Spectrum quantitative automatic analysis method | |
| CN106442474B (en) | A kind of cement slurry three ratio measurement method based on Partial Least Squares | |
| CN109001182B (en) | Raman spectrum nondestructive testing method for alcohol content in closed container | |
| CN104849234A (en) | Assay method for analyzing contents of principal components of imidacloprid based on near-infrared spectrum | |
| CN1793854A (en) | Method and system for determining sewage COD using spectrum technology | |
| CN110987845A (en) | Real-time detection method for dye concentration in single-component and double-component dye solutions | |
| CN102384786B (en) | A kind of spectrophotometer and control method thereof | |
| CN113916817B (en) | Spectrum method chromaticity online measurement method for urban living drinking water | |
| CN103398966A (en) | Method for detecting TMC concentration in organic solution by using spectrometer | |
| CN116399836A (en) | Cross-talk fluorescence spectrum decomposition method based on alternating gradient descent algorithm | |
| CN102384785B (en) | Full-wave band wavelength calibration method for spectrophotometer | |
| CN107727602A (en) | The method of middle infrared spectrum combination vector angle quantitative analysis Sucralose content | |
| CN2741042Y (en) | On-line multi-channel near-infrared spectrograph for measuring total nitrogen and acid content in soy sauce | |
| CN2854583Y (en) | System for investigating sewage chemical oxygen demand by spectral technology | |
| RU2366929C1 (en) | Method for quantative determination of methionine in water solutions | |
| CN2739619Y (en) | On-line multi-channel near infrared spectrougraph for measuring vinegar total acid content |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C19 | Lapse of patent right due to non-payment of the annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |