JP7315200B2 - Analysis device, method and program - Google Patents

Description

The present invention relates to a technique for analyzing the concentration or content of another specific substance contained in a substance, or whether or not a substance contains another substance.

例えば、米等の粉状、粒状又は不定形のバラ状の対象物の水分量や材料の特性を、体積充填率の変動に影響されず、オンラインでリアルタイムに非接触且つ非破壊で測定する技術として、例えば以下に示す特許文献１には、周波数１０ ⁵乃至１０ ¹² ［Ｈｚ］にある所定周波数の電磁波を照射する電磁波送信手段と、透過または反射した電磁波を受信する電磁波受信手段と、送受信波の位相と振幅を測定する測定器と、測定器での測定から振幅変化と位相差を検出し解析する解析装置とを備える水分量測定装置が開示されている。 This analyzer detects the amplitude change and the phase difference for a substance with a known moisture content in advance, obtains the slope of a straight line linearly approximating the relationship between the amplitude change and the phase difference, obtains a calibration curve representing the relationship between the slope and the water content, stores it in a storage means, obtains the slope of the straight line that linearly approximates a plurality of pairs of the amplitude change and the phase difference detected by irradiating the substance flowing through the processing line with an electromagnetic wave of a predetermined frequency, and detects the water content by applying the slope to the calibration curve.

This technique is premised on measurement at a frequency where linear approximation can be obtained and the linearity of the calibration curve can be maintained, resulting in the problem of poor sensitivity.

Japanese Patent No. 6253096

Therefore, one aspect of the present invention is to provide a technique for accurately measuring the concentration or content of another specific substance contained in a substance, or whether or not a substance contains another substance.

An analysis apparatus according to the present invention is an analysis apparatus for analyzing the concentration or content of another specific substance contained in a substance or whether or not a substance contains another substance, comprising: (A) means for acquiring data representing transmission characteristics, reflection characteristics or impedance characteristics at a specific frequency within a predetermined band including the original resonance frequency of the measurement system for a measurement system in which a substance causes a change in resonance frequency; and (C) means for determining the concentration or content of another specific substance corresponding to the calculated predetermined feature quantity, or whether or not the other substance is included, based on the relational data stored in the data storage unit that stores relationship data representing the relationship between the concentration or content of another specific substance and the predetermined feature quantity, or the relationship between whether or not the other substance is included and the predetermined feature quantity.

According to one aspect, it becomes possible to accurately measure the concentration or content of another specific substance contained in a certain substance, or whether or not a certain substance contains another substance.

FIG. 1 is a diagram showing an overview of the entire system. FIG. 2 is a diagram for explaining the shift of the resonance frequency due to the object to be measured. FIG. 3(a) is a diagram for explaining the amount of amplitude change, and FIG. 3(b) is a diagram for explaining the amount of phase change. FIG. 4 is a functional block diagram of the analysis device. FIG. 5 is a diagram illustrating a processing flow of the analysis device; FIG. 6 is a diagram for explaining the fitting of a quadratic function on a plane spanned by amplitude variation and phase variation. FIG. 7 is a diagram for explaining fitting of a quadratic function in the complex plane. FIG. 8 is a diagram for explaining the calibration curve. FIG. 9 is a diagram for explaining circle fitting in the complex plane. FIG. 10 is a functional block diagram of a computer. FIG. 11 (a) is a diagram showing an overview of the entire system when a microstrip line antenna is used as an example of a narrowband antenna, and (b) is a diagram showing an overview of the entire system when a patch antenna is used as an example of a narrowband antenna. FIG. 12(a) is a diagram showing an overview of the entire system when a microstrip line-type antenna is used as an example of a narrowband antenna, and FIG. 12(b) is a diagram showing an overview of the entire system when a patch-type antenna is used as an example of a narrowband antenna. FIGS. 13(a) and 13(b) are diagrams showing an example of a measurement system when no vector network analyzer is used.

In Patent Document 1 mentioned above, for example, when rice flows through a processing line with air, an electromagnetic wave of a predetermined frequency is irradiated, and the amplitude change and phase difference of the transmitted electromagnetic wave or reflected electromagnetic wave are detected at a plurality of times, and the moisture content of the rice is specified by the slope of the linear approximation on the plane spanned by the amplitude change and the phase difference.

In this technology, any frequency can be selected from the frequency bands such as RF waves, microwaves, and millimeter waves, but in reality, the measurement is performed using frequencies that can obtain measurement results that can be linearly approximated, and frequencies that can maintain a linear relationship between the moisture content and the slope of the line.

Measurement under such a premise has a problem that the sensitivity is poor and the measurement accuracy is lowered.

Therefore, in the present embodiment, the measurement accuracy is increased by increasing the sensitivity of the measurement system using resonance. However, since the result of such measurement has non-linearity, analysis processing corresponding to it is performed. As a measurement system using resonance, there are cases where an interferometer is used, and cases where a narrowband antenna is used.

Note that the present embodiment is applicable not only to measuring the water content of a powdery, granular, or amorphous loose substance as in Patent Document 1, but also to measuring the concentration or content of another substance contained in a certain substance (main substance). Furthermore, it also corresponds to determining whether or not a certain substance (main substance) contains other substances (sugar, salt, alcohol, fat, carbon dioxide, urea, etc., and may be some kind of foreign substance).

An example using an interferometer is shown in FIG.

The system according to this embodiment has an interferometer 100 , a vector network analyzer 200 as a measuring instrument, and an analysis device 300 .

Interferometer 100 has coaxial cable 101 connected to the output of vector network analyzer 200 , T-branch circuit 103 , microstrip line 104 , and coaxial cable 102 connected to the input of vector network analyzer 200 .

A signal S1 output from the vector network analyzer 200 is branched into a signal S2 not passing through the microstrip line 104 and a signal S3 to the microstrip line 104 in the T branch circuit 103 . The signal S3 branched to the microstrip line 104 propagates through the microstrip line 104, is reflected by the short circuit at the other end, and returns to the T branch circuit 103 as the reflected signal S4. Then, in the T branch circuit 103 , the signals S 2 and S 3 are synthesized, and the synthesized signal S 5 is input to the vector network analyzer 200 .

The vector network analyzer 200 is a general one and measures the transmission characteristics, reflection characteristics, impedance characteristics, etc. of the device under test. Data representing these characteristics include a combination of amplitude and phase and a combination of real and imaginary parts. Both represent substantially the same thing and are convertible from one to the other.

The electrical length L _e from the branch point of the T branch circuit 103 to the short circuit of the microstrip line 104 is determined by L _e =L/{ε _r } ^0.5 from the physical length L and the effective dielectric constant ε _r of the substrate, cable, pattern, etc. The effective dielectric constant ε _r also changes depending on the object to be measured placed on the microstrip line 104 .

Based on this electrical length L _e , when the wavelength of the signal is λ/2=2*L _e , the resonance state is reached and the amplitude of the signal is minimized. This situation is shown in FIG.

In FIG. 2, the horizontal axis represents frequency and the vertical axis represents signal amplitude. As shown in FIG. 2, notches with minimum amplitude periodically appear based on the electrical length L _e . Curve A represents the characteristics when the object to be measured is not placed on the microstrip line 104 , and curve B represents the characteristics when the object to be measured is placed on the microstrip line 104 . As described above, it can be seen that the resonance frequency is shifted and the notch is shifted because the effective dielectric constant changes depending on the object to be measured.

In the present embodiment, measurement is performed using such a phenomenon. As shown in FIG. 3A, the amplitude is analyzed based on the amount of change in amplitude, which is the difference between the value of curve B and the value of curve A at a specific observation frequency f _v in a predetermined band W including the resonance frequency of interferometer 100. As shown in FIG. 3B, the phase is also analyzed based on the phase change amount, which is the difference between the value of curve C when the object to be measured is not placed on the microstrip line 104 and the value of curve D when the object to be measured is placed on the microstrip line 104 at a specific observation frequency _fv .

Note that the predetermined band W is a range in which sensitivity can be enhanced using resonance, and varies depending on the measurement system and the object to be measured. The observation frequency f _v is specified within a predetermined band W as a frequency other than the resonance frequency, for example, at which the amplitude change amount and the phase change amount are sufficiently large.

A combination of such amplitude and phase changes may be measured as a combination of real and imaginary parts.

In the measurement system shown in FIG. 1, the object to be measured is placed on the microstrip line 104 and measured. For example, if it is a liquid, it is placed in a container such as a bag. At this time, the state of the object to be measured is changed by changing the placement method or by changing the volume or amount of the object to be measured, and the measurement is performed a plurality of times.

Based on the measurement results output from the vector network analyzer 200, the analysis device 300 determines the concentration or content of other substances contained in a given substance.

FIG. 4 shows an example of the functional configuration of the analysis device 300. As shown in FIG.

The analysis device 300 has a data acquisition unit 301 , a measurement data storage unit 303 , a feature amount calculation unit 305 , a feature amount storage unit 307 , a determination unit 309 , a calibration curve data storage unit 311 and an output unit 313 .

The data acquisition section 301 acquires measurement data from the vector network analyzer 200 and stores it in the measurement data storage section 303 . The feature amount calculation unit 305 executes processing for calculating feature amounts from the measurement data stored in the measurement data storage unit 303 and stores the calculated feature amounts in the feature amount storage unit 307 . In the present embodiment, quadratic function fitting or circle fitting is performed from the measurement data for a plurality of times to calculate the coefficients of the fitting curve.

The calibration curve data storage unit 311 stores calibration curve data, for example, for each coefficient of the fitting curve. The calibration curve is also a quadratic function, for example. Based on the data stored in the calibration curve data storage unit 311, the determination unit 309 determines the content or concentration of another substance corresponding to the feature quantity stored in the feature quantity storage unit 307, or whether a certain substance contains another substance. The output unit 313 outputs the processing result of the determination unit 309 to an output device such as a display device.

Next, the processing contents of the analysis device 300 will be described with reference to FIGS. 5 to 9. FIG.

First, the data acquisition unit 301 acquires measurement data for different states of the device under test from a measuring device such as the vector network analyzer 200, and stores the data in the measurement data storage unit 303 ( FIG. 5 : step S101).

The feature amount calculation unit 305 calculates the coefficient of the fitting curve of the quadratic function as the feature amount from the data stored in the measurement data storage unit 303, and stores it in the feature amount storage unit 307 (step S103).

Although the amount of amplitude change and the amount of phase change, or the real and imaginary parts thereof, are obtained from the measuring device, these values may be calculated from the measurement data from the measuring device by the data acquisition unit 301 or the feature amount calculation unit 305 or the like.

A case in which data on the amount of amplitude change and the amount of phase change are acquired will be described with reference to FIG. In FIG. 6, the vertical axis represents the amount of amplitude change [dB], and the horizontal axis represents the amount of phase change [deg]. Here, measurement results for high-purity water and an impurity aqueous solution containing chlorine as an impurity are shown at a frequency of 300 MHz. In this way, when the measured data are plotted on the plane spanned by the amplitude variation amount and the phase variation amount, each of them is fitted to the curve of the quadratic function. In the example of FIG. 6, high-purity water is fitted to curve a, and impurity aqueous solution is fitted to curve b. As a fitting method, a well-known method such as the method of least squares may be used.

In the case of high-purity water, the change is steep near the phase change amount of −50°, and the data is missing.

The coefficients a, b and c are then calculated for each curve in the form y=ax ² +bx+c, as shown in FIG.

As for the phase, the output of the measuring instrument is generally in the range of -180° to +180°, but depending on the measurement data, it may be separated into the vicinity of -180° and the vicinity of +180°. In this case, for example, by performing folding processing to convert −170° to +190°, the data are put together on one side to facilitate fitting.

On the other hand, the case where the data of the real part and the imaginary part are acquired will be described with reference to FIG. In FIG. 7, the vertical axis represents the imaginary part (Imaginary) and the horizontal axis represents the real part (Real). Here, measurement results for high-purity water and impurity aqueous solution at a frequency of 300 MHz are shown. Thus, when the measured data are plotted in the plane spanned by the real and imaginary parts (that is, the complex plane), they are each fitted to a curve of a quadratic function. In the example of FIG. 7, high-purity water is fitted to curve c, and impurity aqueous solution is fitted to curve d.

Then, even in the case of FIG. 7, the coefficients a, b and c are calculated for each curve in the form of y=ax ² +bx+c.

In the case of real part and imaginary part data, there is an advantage that folding processing is unnecessary.

Next, based on the calibration curve data stored in the calibration curve data storage unit 311, the determination unit 309 determines the concentration or content corresponding to the feature amount data (i.e., fitting curve coefficient) stored in the feature amount storage unit 307 (step S107).

For example, data as shown in FIG. 8 is stored in the calibration curve data storage unit 311 . In the example of FIG. 8, the horizontal axis represents concentration [ppm], the left vertical axis represents the value of coefficient a, and the right vertical axis represents the value of coefficient b or c/100. FIG. 8 shows calibration curve data for the amplitude variation and phase variation shown in FIG.

In this embodiment, a calibration curve fitted with a quadratic function is drawn for each coefficient. That is, a calibration curve can be obtained by measuring a sample with a known concentration using the above measurement system, calculating the coefficients of the fitting curve for the amplitude change amount and the phase change amount as described above, and fitting the relationship between the coefficients and the concentration to the curve of the quadratic function.

Since there are three coefficients, three calibration curves are obtained, but any one calibration curve with high reliability may be specified and the concentration determined from that calibration curve.

On the other hand, after specifying the concentration from the calibration curve for each of the three coefficients, the concentration may be determined by averaging or the like. A weighted average may be used.

Even if the measured data are the real part and the imaginary part, the concentration can be specified by similarly preparing the calibration curve. Also, the content is the same as that of the concentration.

When determining whether or not a certain substance contains another substance, for example, for each coefficient, the value range when the other substance is contained or the value range when the other substance is not contained is specified, and the data is stored in the calibration curve data storage unit 311. Then, for example, for all coefficients or a specific coefficient, it is sufficient to confirm whether or not they are within the value range when other substances are included, or whether they are within the value range when other substances are not included.

Depending on the situation, if the value range of each coefficient is specified for the case where substance A may contain substance B and substance C, the case where substance A only contains substance B, the case where substance A contains substance B, the case where substance A contains substance C, and the case where substance A contains substance B and substance C, it may be possible to distinguish between the above cases.

The output unit 309 outputs the processing result by the determination unit 309 to a display device or the like that is an output device (step S107). In addition, it may be output to other computers connected to the network.

In FIG. 7, fitting with a quadratic function was described, but when there is a large change in the state of the object to be measured (arrangement state, amount, density ratio with air, etc.), nonlinearity becomes strong, so circular fitting may be preferable.

FIG. 9 shows an example of circle fitting. In FIG. 9, the vertical axis represents the imaginary part, the horizontal axis represents the real part, the curve e represents the arc of the high-purity water, and the curve f represents the arc of the impurity aqueous solution. For circle fitting, calculate the coefficients A, B and C in the function x ² +y ² +Ax+By+C=0. In the case of FIG. 9, A=−0.726, B=0.830, C=0.122 for curve e and A=−1.226, B=1.448, C=0.302 for curve f.

Of course, when calculating the coefficients A, B and C of circle fitting, the calibration curve data for each of the coefficients A, B and C are prepared in the calibration curve data storage unit 311 . When determining whether or not a certain substance contains another substance, for example, for each coefficient, the value range when the other substance is contained or the value range when the other substance is not contained is specified, and the data is stored in the calibration curve data storage unit 311.

In this way, even if the state of the object to be measured changes, the sensitivity can be improved by resonance, and accurate measurement can be performed.

Although the embodiment of the present invention has been described above, the present invention is not limited to this.

For example, the functional block configuration shown in FIG. 4 is an example, and may not match the program module configuration. Also, the analysis device 300 may be realized by a plurality of devices instead of a single device. The analysis device 3 may be integrated with the measuring device.

The measurement system using the interferometer 100 is also an example, and another measurement system capable of measuring at a frequency near the resonance frequency may be adopted by utilizing the fact that the resonance frequency changes depending on the object to be measured.

A narrow band antenna as shown in FIGS. 11 and 12 may be used instead of the interferometer 100 as shown in FIG. In the example of FIGS. 11 and 12, the narrowband antenna 400 is a microstrip line antenna (FIGS. 11(a) and 12(a)) and a patch antenna (FIGS. 11(b) and 12(b)). Such a narrowband antenna 400 is based on the principle of resonance, emits electromagnetic waves in a narrow band including a specific frequency (corresponding to the resonance frequency), and receives electromagnetic waves in such a band. 11(a) and (b) show a reflection type example in which a narrowband antenna 400 is connected to one port of the vector network analyzer 200. FIG. That is, an electromagnetic wave having a frequency near the resonance frequency of narrowband antenna 400 is radiated to the object under test, a reflected wave from the object under test is received by narrowband antenna 400, and vector network analyzer 200 measures data representing reflection characteristics. On the other hand, the examples of FIGS. 12( a ) and 12 ( b ) show interference type examples in which the narrowband antenna 400 is connected to two ports of the vector network analyzer 200 . This is to measure the state of interference between the electromagnetic waves radiated from the narrowband antenna 400 to the object under test and the electromagnetic waves reflected from the object under test with the vector network analyzer 200 in the same manner as the interferometer described above.

In addition, by changing the volume or amount of the object to be measured, the state of the object to be measured is changed and the measurement is performed a plurality of times, but the object to be measured placed in a container such as a bag may be measured for each container while it is being moved by a conveyer such as a belt conveyor. Also, by devising the arrangement of the interferometer 100 or the narrowband antenna 400, it is possible to measure an object flowing through a pipe or the like at a plurality of timings. In this way, online analysis becomes possible.

The vector network analyzer 200 is a device that outputs a signal and detects the characteristics of the signal after it has been transmitted through a transmission line or an object or the reflected signal from them as amplitude and phase or real and imaginary parts. In the above example, only the example using the vector network analyzer 200 was shown, but as shown in FIG. In FIG. 13(a), the receiver 1 receives a signal from a signal source, and the receiver 2 receives a reflected signal from a sensor or a circuit (resonant circuit, etc.), and the amplitude change amount and the phase change amount are calculated from those outputs, for example, by the analysis device 300. Similarly, in FIG. 13B, the receiver 1 receives the signal from the signal source, and the receiver 2 receives the transmission signal, which is the output from the sensor or the circuit.

10, a memory 2501, a CPU (Central Processing Unit) 2503, a hard disk drive (HDD) 2505, a display control unit 2507 connected to a display device 2509, a drive device 2513 for a removable disk 2511, an input device 2515, and a communication control unit 2517 for connecting to a network are connected to a bus. 2519 are connected. Note that the HDD may be a storage device such as a solid state drive (SSD). An operating system (OS) and an application program for performing processing in the embodiments of the present invention are stored in the HDD 2505 and read from the HDD 2505 to the memory 2501 when executed by the CPU 2503 . The CPU 2503 controls the display control unit 2507, the communication control unit 2517, and the drive device 2513 according to the processing content of the application program to perform predetermined operations. In addition, data in the middle of processing is mainly stored in the memory 2501, but may be stored in the HDD 2505 as well. In the embodiment of the present invention, an application program for carrying out the processing described above is stored and distributed in computer-readable removable disk 2511 and installed from drive device 2513 to HDD 2505 . It may be installed in the HDD 2505 via a network such as the Internet and the communication control unit 2517 . Such a computer apparatus implements the various functions described above through organic cooperation between hardware such as the CPU 2503 and memory 2501 described above and programs such as the OS and application programs.

It should be noted that data used by executing the above-described processing is stored in a storage device such as the memory 2501 or the HDD 2505 regardless of whether it is in the middle of processing or processing results.

The embodiments described above are summarized as follows.

The analysis apparatus according to the present embodiment is an analysis apparatus that analyzes the concentration or content of another specific substance contained in a certain substance, or whether or not a certain substance contains other substances (either singular or plural), comprising: (A) means for acquiring data representing transmission characteristics, reflection characteristics, or impedance characteristics at a specific frequency within a predetermined band including the original resonance frequency of the measurement system for a measurement system in which the resonance frequency is changed by a certain substance, for a plurality of different states of the substance; (C) means for determining the concentration or content of another specific substance corresponding to the calculated predetermined feature quantity, or whether or not another substance is included, based on the relational data stored in the data storage section storing the relationship data representing the relationship between the concentration or content of another specific substance and the predetermined feature quantity, or the relationship between the presence of another substance and the predetermined feature quantity.

By improving the sensitivity through resonance in this way, it becomes possible to accurately identify the concentration or content of another specific substance contained in a substance, or whether or not a substance contains another substance.

The above-described calculation means may calculate a fitting curve on a plane spanned by the amplitude change amount and the phase change amount from the data of the amplitude change amount and the phase change amount specified by the data representing the transmission characteristics, the reflection characteristics, or the impedance characteristics for a plurality of times, and calculate the coefficient of the fitting curve as a predetermined feature amount. This makes it possible to appropriately extract the feature quantity of a certain substance. In this case, quadratic function fitting is preferred.

Further, the calculation means described above may calculate a fitting curve on a plane spanned by the real and imaginary parts from the data of the real and imaginary parts specified by the data representing the transmission characteristics, the reflection characteristics, or the impedance characteristics for a plurality of times, and calculate the coefficient of the fitting curve as a predetermined feature amount. By focusing on the real and imaginary parts instead of the phase and amplitude, it is possible to avoid the range conversion associated with the phase. Again, a quadratic function fitting may be performed, but in other cases it may be preferable to perform a circle fitting.

A program for causing a computer to execute the processing of the analysis apparatus described above can be created, and the program is stored in various storage media.

100 interferometer 200 vector network analyzer 300 analysis device 301 data acquisition unit 303 measurement data storage unit 305 feature amount calculation unit 307 feature amount storage unit
309 determination unit 311 calibration curve data storage unit 313 output unit

Claims (7)

An analysis device for analyzing the concentration or content of another specific substance contained in a substance, or whether or not the substance contains another substance,
Means for acquiring data representing transmission characteristics, reflection characteristics, or impedance characteristics at a specific frequency within a predetermined band including the original resonance frequency of the measurement system for a plurality of different states of the substance while the substance exists within the measurement range of the measurement system;
a calculation means for calculating a predetermined feature amount from the data representing the transmission characteristics, reflection characteristics, or impedance characteristics acquired a plurality of times;
Means for determining the concentration or content of the other specific substance corresponding to the calculated predetermined feature quantity, or whether the other substance is included, based on the relational data stored in the data storage unit storing relational data representing the relationship between the concentration or content of the other specific substance and the predetermined feature quantity or the relationship between whether the other substance is included and the predetermined feature quantity;
An analysis device having The calculation means is
2. The analysis device according to claim 1, wherein a fitting curve on a plane spanned by the amplitude change amount and the phase change amount is calculated from the data of the amplitude change amount and the phase change amount specified by the data representing the transmission characteristics, the reflection characteristics, or the impedance characteristics for a plurality of times, and the coefficient of the fitting curve is calculated as the predetermined feature amount. The calculation means is
2. The analysis device according to claim 1, wherein a fitting curve on a plane spanned by the real and imaginary parts is calculated from the data of the real and imaginary parts specified by the data representing the transmission characteristics, the reflection characteristics, or the impedance characteristics for a plurality of times, and the coefficient of the fitting curve is calculated as the predetermined feature quantity. The analysis device according to claim 2, wherein the fitting curve is a quadratic function fitting curve. The analysis device according to claim 3, wherein the fitting curve is a quadratic function fitting curve or a circle fitting curve. An analysis method for analyzing the concentration or content of another specific substance contained in a substance, or whether or not the substance contains another substance,
the computer
obtaining data representing transmission characteristics, reflection characteristics, or impedance characteristics at a specific frequency within a predetermined band including the original resonance frequency of the measurement system for a plurality of different states of the substance while the substance exists within the measurement range of the measurement system ;
a step of calculating a predetermined feature amount from the data representing the transmission characteristics, reflection characteristics, or impedance characteristics acquired a plurality of times;
determining the concentration or content of the other specific substance corresponding to the calculated predetermined feature quantity, or whether or not the other substance is included, based on the relational data stored in the data storage unit storing relationship data representing the relationship between the concentration or content of the other specific substance and the predetermined feature quantity or the relationship between whether the other substance is included and the predetermined feature quantity;
The analysis method to perform. An analysis program for analyzing the concentration or content of another specific substance contained in a certain substance, or whether or not the certain substance contains another substance,
to the computer,
obtaining data representing transmission characteristics, reflection characteristics, or impedance characteristics at a specific frequency within a predetermined band including the original resonance frequency of the measurement system for a plurality of different states of the substance while the substance exists within the measurement range of the measurement system ;
a step of calculating a predetermined feature amount from the data representing the transmission characteristics, reflection characteristics, or impedance characteristics acquired a plurality of times;
determining the concentration or content of the other specific substance corresponding to the calculated predetermined feature quantity, or whether or not the other substance is included, based on the relational data stored in the data storage unit storing relationship data representing the relationship between the concentration or content of the other specific substance and the predetermined feature quantity or the relationship between whether the other substance is included and the predetermined feature quantity;
Analysis program for executing

Images