Abstract:
In a pulse oximeter for obtaining an oxygen saturation in a blood, a light emitter irradiates a living tissue with light beams having five different wavelengths. A light receiver receives respective light beams reflected from or transmitted through the living tissue, and converts the received light beams to electric signals. A first calculator calculates five attenuation changes of the living tissue based on fluctuations of the respective electric signals. A second calculator calculates at least four attenuation change ratios from the five attenuation changes. Each of the attenuation change ratios is defined by a ratio between any two of the five attenuation changes. A third calculator calculates the oxygen saturation based on the attenuation change ratios, while taking an oxygen saturation of arterial blood, an oxygen saturation of venous blood, a ratio between changes in arterial blood and venous blood, and a tissue term as four unknown values.

Description:
BACKGROUND OF THE INVENTION  
       [0001]     The present invention relates to a pulse oximeter which can measure an oxygen saturation of arterial blood continuously and non-invasively by utilization of variations in the volume of arterial blood by pulsation.  
         [0002]     A pulse oximeter has hitherto been widely known as an apparatus which measures an oxygen saturation of blood; more specifically, a concentration ratio of oxyhemoglobin to a sum of oxyhemoglobin and deoxyhemoglobin. Japanese Patent Publication No. 5-88609B discloses an apparatus for determining concentrations of constituents of blood with high precision without being affected by pure tissue whose thickness varies due to pulsation of blood.  
         [0003]     More specifically, this publication discloses an apparatus for determining concentrations of constituents of blood comprising: a light emitter which irradiates a living tissue with light beams of N different wavelengths; a light receiver which receives light originated from the light emitter and then reflected from or transmitted through a living tissue; a first circuit which detects an attenuation change at the living tissue based on output signals from the light receiver for each of the N wavelengths; a second circuit which calculates N-1 attenuation change ratios between different wavelengths, in accordance with detection signals pertaining to the N different wavelengths output from the first circuit; and a third circuit for calculating which calculates relative concentrations of N-1 constituents of blood, through use of arithmetic equations determined by solving simultaneous equations with N-1 unknowns in connection with concentrations of constituents of blood, on the basis of values of attenuation change ratios output from the second circuit. This calculation is based on the assumption that attenuation changes in the living tissue stem from changes in the thickness of blood and thickness of pure tissue which does not include blood.  
         [0004]     In the above apparatus, the second circuit enables calculation of a ratio between N-1 pulsation changes in different wavelengths when the second circuit receives signals output from the first circuit. The third circuit performs computation by substituting actually-measured values of ratios of pulsation changes and respective coefficient values into an equation for determining relative concentrations of N-1 constituents of blood. The equation is obtained by solving simultaneous equations pertaining to N-1 pulsation change ratios that also take into consideration the influence of pulsation of pure tissue. Accordingly, concentrations (relative concentrations) of N-1 constituents of blood can be measured with high precision and without being influenced by pulsation of pure tissue.  
         [0005]     Meanwhile, when body motion occurs during measurement of concentrations of constituents of blood by a pulse oximeter, artifacts are superposed on the transmitted light. When the artifacts stemming from such body motion arm large, difficulty is encountered in removing the artifacts by correcting pulse waveforms, or the like. Japanese Patent Publication No. 11-216133A discloses a pulse oximeter capable of performing high-precision measurement without being influenced by the artifacts stemming from body motion even when the artifacts are large.  
         [0006]     More specifically, the pulse oximeter comprises: a light emitter which irradiates a living tissue with a plurality of light beams having different wavelengths; a photoelectric converter which converts a light beam transmitted through the living tissue into an electric signal for each of the wavelengths; a first detector which detects an attenuation change at the living tissue based on fluctuations of the signals output from the photoelectric converter for each of the wavelengths; a variable filter which receives the attenuation change for each of the wavelengths and allows to pass through a component having a prescribed frequency band; a band prescriber which prescribes the frequency band for the variable filter; and a second detector which determines an oxygen saturation in accordance with an output from the variable filter.  
         [0007]     The apparatus disclosed in Japanese Patent Publication No. 5-88609B teaches that a plurality of light beams having different wavelengths to precisely measure concentrations of constituents of blood without being influenced by pure tissue whose thickness varies due to pulsation of blood. The influence of a term of the pure tissue can be eliminated by use of light beams of three different wavelengths during the measurement process, thereby enabling measurement of a relative concentration between oxyhemoglobin and deoxyhemoglobin, which are two constituents of blood.  
         [0008]     The influence of a term of the pure tissue can be eliminated by use of light beams of four different wavelengths during the measurement process, as in the case of the above, thereby enabling precise measurement of relative concentrations between oxyhemoglobin, deoxyhemoglobin and another pigment, which are three constituents of blood. Furthermore, use of light beams of five different wavelengths enables measurement of relative concentrations between four constituents of blood including carboxyhemoglobin in addition to the above three constituents of blood measured by use of the four different light wavelengths.  
         [0009]     However, this publication is silent about the influence of the artifacts stemming from body motion. Therefore, means for eliminating the influence of the artifacts is not disclosed at all.  
         [0010]     Meanwhile, Japanese Patent Publication No. 11-216133A teaches means for eliminating the influence of the artifacts stemming from body motion. However, in a pulse oximeter disclosed in this publication, changes in the thickness of arterial blood and changes in the thickness of pure tissue, which are assumed to be causes of the aromas stemming from body motion, are regarded as negligible factors. Focus is paid to changes in the thickness of venous blood as the paramount cause, and measures are taken against the changes in the thickness of venous blood. In summary, this publication teaches a method and an apparatus for performing: a two-wavelength measurement which disregards changes in tissue thickness; and a three-wavelength measurement which takes into consideration changes in tissue thickness.  
         [0011]     Accordingly, the most serious problem faced by a currently available pulse oximeter is the artifacts stemming from body motion. A primary method adopted as a countermeasure for avoiding the influence of the artifacts stemming from body motion is a statistical method. Specifically, reference is made to antecedent data and subsequent data in order to obtain a reliable measurement value at a certain point in time. Therefore, such a measurement entails a time delay and value smoothing. This is against the primary object of a pulse oximeter to detect anomalies in a patient at an early stage. That is, a function for eliminating the influence of artifacts is insufficient or involves inconveniences.  
       SUMMARY OF THE INVENTION  
       [0012]     It is therefore an object of the invention to provide a pulse oximeter capable of measuring an oxygen saturation of arterial blood with high precision by using a deterministic method which is able to eliminate the aforementioned influence of artifacts stemming from body motion, as well as eliminating the time delay and the value smoothing.  
         [0013]     In order to achieve the above object, according to the invention, there is provided a pulse oximeter for obtaining an oxygen saturation in a blood, comprising: 
        a light emitter, adapted to irradiate a living tissue with light beams having five different wavelengths;     a light receiver, adapted to receive respective light beams reflected from or transmitted through the living tissue, and to convert the received light beams to electric signals;     a first calculator, which calculates five attenuation changes of the living tissue based on fluctuations of the respective electric signals;     a second calculator, which calculates at least four attenuation change ratios from the five attenuation changes, each of the attenuation change ratios being defined by a ratio between any two of the five attenuation changes; and     a third calculator, which calculates the oxygen saturation based on the attenuation change ratios, while taking an oxygen saturation of arterial blood, an oxygen saturation of venous blood, a ratio between changes in arterial blood and venous blood, and a tissue term as four unknown values.        
 
         [0019]     With the above configuration, the influence of the pulsation when venous blood is being pulsated for some reason can be reliably eliminated, whereby an oxygen saturation of arterial blood can be measured with high precision and without involving the time delay or the value smoothing. Furthermore, even in a case where a pulse wave in a living tissue of a subject is so small that a pulse oximetry is not applicable, measurement of an oxygen saturation of arterial blood is enabled by forcibly effecting body motion against the living tissue by a vibration mechanism or the like, to thus generate artifacts on attenuation changes in light transmitted through the living tissue. Moreover, an oxygen saturation of venous blood can also be measured simultaneously with the oxygen saturation of arterial blood.  
         [0020]     According to the invention, there is also provided a method of obtaining an oxygen saturation in a blood, comprising steps of: 
        irradiating a living tissue with light beams having five different wavelengths;     receiving respective light beams reflected from or transmitted through the living tissue;     converting the received light beams to electric signals;     calculating five attenuation changes of the living tissue based on fluctuations of the respective electric signals;     calculating at least four attenuation change ratios from the five attenuation changes, each of the attenuation change ratios being defined by a ratio between any two of the five attenuation changes; and     calculating the oxygen saturation based on the attenuation change ratios, while taking an oxygen saturation of arterial blood, an oxygen saturation of venous blood, a ratio between changes in arterial blood and venous blood, and a tissue term as four unknown values.        
 
         [0027]     The venous blood is pressed by the pulsation of the arterial blood, and inverse pulsation is thus presented. This fact contributes an increase of the SpO2 value. This increase should not be ignored in a case where the SaO2 value is almost 100%. Although there is a possibility that the SpO2 value exceeds 100%, the conventional two-wavelength type pulse oxymeter ignores the amount of the SpO2 value exceeding 100%. However, the contribution from the inversed pulsation of the venous blood should be removed to obtain an accurate SpO2 value. For example, an extremely immature baby would become blind due to retinopathy if the SaO2 value is too large, whereas would become encephalopathy or die if the SaO2 value is too small. It is apparent that the detection of an accurate SpO2 value is so critical for the extremely immature baby.  
         [0028]     According to the above configurations, since the contribution from the inversed pulsation of the venous blood can be removed by the five-length measurement, so that an accurate SpO2 value can be obtained especially in the case where the SaO2 value is almost 100%. Therefore, it is advantageous for the respiration monitoring of the extremely immature baby.  
     
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0029]     The above objects and advantages of the present invention will become more apparent by describing in detail preferred exemplary embodiments thereof with reference to the accompanying drawings, wherein:  
         [0030]      FIG. 1  Is a block diagram showing an electric configuration of a pulse oximeter according to one embodiment of the present invention;  
         [0031]      FIG. 2  is a graph showing a trace on planes in a case where a subject moved his hand in a “waving” manner;  
         [0032]      FIG. 3A  is a graph showing artifact characteristics obtained by a two-wavelength calculation as a result of body motion in the case of  FIG. 2 ;  
         [0033]      FIG. 3B  is a graph showing artifact characteristics obtained by a three-wavelength calculation as a result of body motion in the case of  FIG. 2 ;  
         [0034]      FIG. 3C  is a graph showing artifact characteristics obtained by a five-wavelength calculation as a result of body motion in the case of  FIG. 2 ;  
         [0035]      FIG. 4  is a graph showing a trace on Φ-Φ planes in a case where a subject moved his hand in a “chopping” manner;  
         [0036]      FIG. 5A  is a graph showing artifact characteristics obtained by the two-wavelength calculation as a result of body motion in the case of  FIG. 4 ;  
         [0037]      FIG. 5B  is a graph showing artifact characteristics obtained by the three-wavelength calculation as a result of body motion in the case of  FIG. 4 ;  
         [0038]      FIG. 5C  is a graph showing artifact characteristics obtained by the five-wavelength calculation as a result of body motion in the case of  FIG. 4 ;  
         [0039]      FIG. 6A  is a graph showing artifact characteristics obtained by the two-wavelength calculation in a case where a subject moved his hand in the “waving” manner after holding his breath;  
         [0040]      FIG. 6B  is a graph showing artifact characteristics obtained by the three-wavelength calculation as a result of body motion in the case of  FIG. 6A ;  
         [0041]      FIG. 6C  is a graph showing artifact characteristics obtained by the five-wavelength calculation as a result of body motion in the case of  FIG. 6A ;  
         [0042]      FIG. 7A  is a graph showing artifact characteristics obtained by the two-wavelength calculation in a case where a subject moved his hand in the “chopping” manner after holding his breath;  
         [0043]      FIG. 7B  is a graph showing artifact characteristics obtained by the three-wavelength calculation as a result of body motion in the case of  FIG. 7A ; and  
         [0044]      FIG. 7C  is a graph showing artifact characteristics obtained by the five-wavelength calculation as a result of body motion in the case of  FIG. 7A . 
     
    
     DETAILED DESCRIPTION OF THE INVENTION  
       [0045]     One embodiment of the present invention will be described below in detail with reference to the accompanying drawings.  
         [0046]     In a pulse oximeter shown in  FIG. 1 , a light emitter  10  includes five light emitting elements LED  1  to LED  5  for emitting light beams of five different wavelengths. The light beams emitted from the light emitter  10  is transmitted through a living tissue  12  and received by a light receiver  14  constituted by a light receiving element PD, a current-voltage converter  16 , and an AD converter  18 . The light receiver  14  outputs transmitted light signals for the respective wavelengths.  
         [0047]     A storage  20  is formed from storage devices  20 A to  20 E which store the transmitted light signals of the respective wavelengths. A peak bottom detector  22  is configured to detect a peak value and a bottom value of changes in a transmitted light signal of one of the wavelengths, to thus output time points at which the peak value and the bottom value appear. A first calculator  24  is constituted by attenuation change calculators  24 A to  24 E which capture transmitted light signals of the respective wavelengths at the time points obtained by the peak bottom detector  22  to calculate attenuation changes ΔAi (i=1, 2, 3, 4, 5) on the basis of the transmitted light signals.  
         [0048]     A second calculator  26  is configured to calculate an oxygen saturation SpO2 on the basis of the attenuation changes ΔAi calculated by the respective attenuation change calculators  24 A to  24 E of the first calculator  24 . More specifically, the second calculator  26  has the function of an attenuation change ratio calculator which calculates attenuation change ratios between the attenuation changes Φij=ΔAi/ΔAj (e.g., Φ12=ΔA1/ΔA2, Φ34=ΔA3/ΔA4, Φ51=ΔA5/ΔA1, Φ23=ΔA2/ΔA3, Φ45=ΔA4/ΔA5) on the basis of the attenuation changes ΔAi; and the function of an oxygen saturation calculator which calculates the oxygen saturation of blood SpO 2  on the basis of the attenuation change ratios Φij by use of four unknown values consisting of: an oxygen saturation of arterial blood Sa; an oxygen saturation of venous blood Sv; a ratio V (=ΔDv/ΔDa) between changes in arterial blood ΔDa and changes in venous blood ΔDv; and a tissue term Wi.  
         [0049]     A timer  28  is configured to generate a timing signal which is in conjunction with the operation of the second calculator  26 , and controls a timing at which the respective light emitting elements LED  1  to LED  5  of the light emitter  10  are to be operated and a timing at which the transmitted light signals are to be stored in the respective storage devices  20 A to  20 E of the storage  20 .  
         [0050]     Next, operations of calculating oxygen saturation of arterial blood and venous blood in the pulse oximeter will be described.  
         [0051]     The light emitting elements LED  1  to LED  6  of the light emitter  10  are caused to illuminate sequentially and alternately at different wavelengths λ1, λ2, λ3, λ4, and λ5 in accordance with the timing signal from the timer  28 . As a result, light which has transmitted through the living tissue  12  is received at the light receiver  14 . In response to the respective wavelengths of the light emitting elements LED  1  to LED  5 , the respective transmitted light signals L 1 , L 2 , L 3 , L 4 , and L 5  are stored in the respective storage devices  20 A to  20 E of the storage  20  in accordance with the timing signal from the timer  28 . Meanwhile, the storage devices  20 A to  20 E store data pertaining to outputs (digital data) produced by the AD converter  18  of the light receiver  14  during a predetermined time period.  
         [0052]     Next, the peak bottom detector  22  detects a peak value and a bottom value of changes in the transmitted light signal L 1  stored in the storage device  20 A and corresponding to the wavelength λ1, so that time points at which the peak value and the bottom value appear (hereinafter, referred as a peak time point and a bottom time point). The peak bottom detector  22  is configured such that a peak time point and a bottom time point can be detected arbitrarily even though changes in the transmitted light signal L 1  are non-periodic; e.g., even when disturbed by body motion.  
         [0053]     Accordingly, the respective attenuation change calculators  24 A to  24 E of the first calculator  24  captures the transmitted light signals of the respective wavelengths at the peak and bottom time points which have been obtained by the peak bottom detector  22 . In this case, the peak and bottom time points are both past times, and values of transmitted light at the time points are represented as L 1 (t-t 0 ), L 2 (t-t 0 ), L 3 (t-t 0 ), L 4 (t-t 0 ), and L 5 (t-t 0 ) (t 0  denotes a time in the past). The respective attenuation change calculators  24 A to  24 E calculate the attenuation changes ΔAi corresponding to the peak and bottom time points with the following approximations:  
         Δ   ⁢           ⁢   Ai     =         Li   ⁡     (   tp   )       -     Li   ⁡     (   tb   )         Li         
 
 where, i=1, 2, 3, 4, 5; tp denotes the peak time point; and tb denotes the bottom time point. 
 
         [0055]     From the attenuation changes ΔAi obtained by the respective attenuation change calculators  24 A to  24 E of the first calculator  24 , ratios between the attenuation changes Φij are obtained from the following equations: 
 
Φ12=Δ A 1/Δ A 2 
 
Φ34=Δ A 3/Δ A 4 
 
Φ51=Δ A 5/Δ A 1 
 
Φ23=Δ A 2/Δ A 3 
 
Φ45=Δ A 4/Δ A 5 
 
 In this case, the respective attenuation changes are represented by the following theoretical formulae:  
             Φ12   =         Δ   ⁢           ⁢   A1       Δ   ⁢           ⁢   A2       =       Eba1   +   Ebv1V   +   W       Eba2   +   Ebv2V   +   B2   +   W                 (   1   )               Φ34   =         Δ   ⁢           ⁢   A3       Δ   ⁢           ⁢   A4       =       Eba3   +   Ebv3V   +   W       Eba4   +   Ebv4V   +   W                 (   2   )               Φ51   =         Δ   ⁢           ⁢   A5       Δ   ⁢           ⁢   A1       =       Eba5   +   Ebv5V   +   W       Eba1   +   Ebv1V   +   W                 (   3   )               Φ23   =         Δ   ⁢           ⁢   A2       Δ   ⁢           ⁢   A3       =       Eba2   +   Ebv2V   +   B2   +   W       Eba3   +   Ebv3V   +   W                 (   4   )                 Φ45   =         Δ   ⁢           ⁢   A4     ΔA5     =       Eba4   +   Ebv4V   +   W       Eba5   +   Ebv5V   +   W           ⁢     
     ⁢     where   ,     
     ⁢           Ebai   =         [     SaEoi   +       (     1   -   Sa     )     ⁢   Eri       ]     ·     [     SaEoi   +       (     1   -   Sa     )     ⁢   Eri     +   F     ]                     Ebvi   =         [     SvEoi   +       (     1   -   Sv     )     ⁢   Eri       ]     ·     [     SvEoi   +       (     1   -   Sv     )     ⁢   Eri     +   F     ]                           (   5   )             
 
 Sa denotes an oxygen saturation of arterial blood; Sv denotes an oxygen saturation of venous blood; Eo denotes an extinction coefficient of oxyhemoglobin; Er denotes an extinction coefficient of deoxyhemoglobin; F denotes a scattering coefficient; and V denotes a ratio of the amplitude of a change in arterial blood ΔDv to the amplitude of a change in venous blood ΔDa, that is, ΔDv/ΔDa. 
 
         [0058]     W denotes changes in living tissue. Comparison of actual data and theoretical formulae on a certain probe has revealed that W has a slight wavelength dependency. In view of the dependency, “W+B2” is used instead of W for the second wavelength.  
         [0059]     When W is eliminated from the equations (1) to (5), the following equations hold:  
                         Φ12   ⁢     (     Eba2   +   B2   +   Ebv2V     )       -               (     Eba1   +   Ebv1V     )           θ12     =               Φ34   ⁢     (     Eba4   +   Ebv4V     )       -               (     Eba3   +   Ebv3V     )           θ34             (   6   )                           Φ34   ⁢     (     Eba4   +   Ebv4V     )       -               (     Eba3   +   Ebv3V     )           θ34     =               Φ51   ⁢     (     Eba1   +   Ebv1V     )       -               (     Eba5   +   Ebv5V     )           θ51             (   7   )                           Φ51   ⁢     (     Eba1   +   Ebv1V     )       -               (     Eba5   +   Ebv5V     )           θ51     =               Φ23   ⁢     (     Eba3   +   Ebv3V     )       -               (     Eba2   +   B2   +   Ebv2V     )           θ23             (   8   )                           Φ23   ⁢     (     Eba3   +   Ebv3V     )       -               (     Eba2   +   B2   +   Ebv2V     )           θ23     =               Φ45   ⁢     (     Eba5   +   Ebv5V     )       -               (     Eba4   +   Ebv4V     )           θ45             (   9   )             
 
 where, θij=(1−Φij). 
 
         [0061]     Subsequently, when V is eliminated from the above equations, the following nonlinear simultaneous equations with two unknowns (10) and (11) can be obtained.  
                               Φ34   ⁢           ⁢   Eba4     -   Eba3     )     /   θ34     -                 [       Φ12   ⁡     (     Eba2   +   B2     )       -   Eba1     ]     /   θ12                       (       Φ12   ⁢           ⁢   Ebv2     -   Ebv1     )     /   θ12     -                 (       Φ34   ⁢           ⁢   Ebv4     -   Ebv3     )     /   θ34             =                 (       Φ51   ⁢           ⁢   Eba1     -   Eba5     )     /   θ51     -                 (       Φ34   ⁢           ⁢   Eba4     -   Eba3     )     /   θ34                       (       Φ34   ⁢           ⁢   Ebv4     -   Ebv3     )     /   θ34     -                 (       Φ51   ⁢           ⁢   Ebv1     -   Ebv5     )     /   θ51                     (   10   )                         [       Φ23   ⁢           ⁢   Eba3     -                       (     Eba2   +   B2     )     ]     /   θ23     -                 (       Φ51   ⁢           ⁢   Eba1     -   Eba5     )     /   θ51                       (       Φ51   ⁢           ⁢   Ebv1     -   Ebv5     )     /   θ51     -                 (     Φ23Ebv3   -   Ebv2     )     /   θ23             =             (       Φ45   ⁢           ⁢   Eba5     -                     Eba4   )     /   θ45     -     [       Φ23   ⁢           ⁢   Eba3     -                       (     Eba2   +   B2     )     ]     /   θ23                       (       Φ23   ⁢           ⁢   Ebv3     -   Ebv2     )     /   θ23     -                 (       Φ45   ⁢           ⁢   Ebv5     -   Ebv4     )     /   θ45                     (   11   )             
 
 Here, the following approximations are made: 
        Eba1=(A13Eba3+B13);     Eba2=(A23Eba3+B23);     Eba4=(A43Eba3+B43);     Eba5=(A53Eba3+B53);     Ebv1=(A13Ebv3+B13);     Ebv2=(A23Ebv3+B23);     Ebv4=(A43Ebv3+B43);     Ebv5=(A53Ebv3+B53); 
 
 When these approximations are substituted into equations (10) and (11), equations (10) and (11) are transformed into simultaneous linear equations with two unknowns for Eba3. When Ebv3 is eliminated from the equation, the following linear equation for Eba3 is obtained:  
             Eba3   =           (     NL   -   JP     )     /     (     NI   -   JM     )       -       F   ⁡     (     B   +   D     )       /     (     FA   -   BE     )               (     AF   +   ED     )     /     (     FA   -   BE     )       -       (     ML   -   IF     )     /     (     NI   -   JM     )                   (   12   )             
 
 where, 
    A: (Φ34A43−1)/θ34−(Φ12A23−A13)θ12;     B: Φ34B43/θ34−[Φ12(B23+B)−B13]/θ12;     D: (Φ12B23−B13)/θ12−Φ34B43/θ34;     E: (Φ51A13−A53)/θ51−(Φ34A43−1)/θ34;     F: (Φ51B13−B53)/θ51−Φ34B43/θ34;     I: (Φ23−A23)/θ23−(Φ51A13−A53)/θ51;     J: −(B23+B)/θ23−(Φ51B13 −B53)/θ51;     L: (Φ51B13−B53)/θ51+B23/θ23;     M: (Φ45A53−A43)/θ45−(Φ23−A23)/θ23;     N: (Φ45B53−B43)/θ45+(B23+B)/θ23; and     P: −B23/θ23−(Φ45B53−B43)/θ45. 
 
 Furthermore, details of the above definitions are as follows: 
    A13: [Eb1(0.8)−Eb1(1)]/[Eb3(0.8)−Eb3(1)];     A23: [Eb2(0.8)−Eb2(1)]/[Eb3(0.8)−Eb3(1)];     A43: [Eb4(0.8)−Eb4(1)]/[Eb3(0.8)−Eb3(1)];     A53: [Eb5(0.8)−Eb5(1)]/[Eb3(0.8)−Eb3(1)];     B13: −A13Eb3(1)+Eb1(1);     B23: −A23Eb3(1)+Eb2(1);     B43: −A43Eb3(1)+Eb4(1); and     B53: −A53Eb3(1)+Eb5(1). 
 
 Meanwhile, Ebi(0.8) and Ebi(1) denote constants; Ebi(0.8) is a value of Ebi where an oxygen saturation is 0.8, and Ebi(1) is a value of Ebi where an oxygen saturation is 1. 
       
 
         [0094]     Next, Eba3 is converted into an oxygen saturation of arterial blood by the following equation:  
             Sa   =         -   B     -         B   3     -     4   ⁢           ⁢   AC             2   ⁢   A               (   13   )             
 
 where, 
        A: (Eo3−Er3) 2       B: (Eo3−Er3)(2Er3+F)     C: Er3(F+Er3)−Eba3 2  
 
 Application of the above equations provides a pulse oximeter which eliminates the influence of artifacts stemming from body motion so as to obtain an oxygen saturation of arterial blood with high precision on the basis of transmitted light of five wavelengths through living tissue. 
       
 
         [0100]     In the first calculator  24 , the respective attenuation change calculators  24 A to  24 E execute calculating operations of the attenuation changes ΔAi on the basis of the aforementioned calculating operations.  
         [0101]     Further, calculating operations of Eba3 and Sa are executed in the second calculator  26  based on the above equations (12) and (13).  
         [0102]     An example operation of the pulse oximeter configured as described above will be described below.  
         [0103]     The following five wavelengths were used: λ1=805 nm; λ2=875 nm; λ3=660 nm; λ4=700 nm; and λ5=730 nm. Meanwhile, with regard to a single point of original data, λ2 and λ3 were used for a two-wavelength calculation, and λ1, λ2, and λ3 were used for a three-wavelength calculation. A subject shook his/her hand from the wrist strongly as body motion while wearing a probe attached to the tip of his/her finger, whereby artifact characteristics of the motion were determined respectively.  
         [0104]     The subject moved his/her hand in a “waving” manner. Specifically, the subject lies on a bed in a supine position such that his/her arm extends downward over the edge of the bed. The subject shakes the hand in the front and back motion while the palm is directed backward (i.e., waving).  FIG. 2  shows a trace on Φ-Φ planes as results of the body motion. These changes show artifacts mainly stemming from changes of tissue.  FIG. 3A  shows artifact characteristics obtained by a two-wavelength calculation which indicates large artifacts.  FIG. 3B  shows artifact characteristics obtained by a three-wavelength calculation which indicates that most of the artifacts are eliminated.  FIG. 3C  shows artifact characteristics obtained by a five-wavelength calculation which indicates that artifacts are further reduced.  
         [0105]     Further, the subject moved his/her hand in a “chopping” manner. Specifically, the subject lies on a bed in a supine position such that his/her arm extends parallel to the side of the body on the bed. The subject shakes the hand in the up and down motion while the palm is directed inward (i.e., chopping).  FIG. 4  shows a trace on Φ-Φ planes as results of the body motion.  FIG. 5A  shows artifact characteristics obtained by the two-wavelength calculation which indicates large artifacts.  FIG. 5B  shows artifact characteristics obtained by the three-wavelength calculation which indicates slight improvement in artifacts.  FIG. 5C  shows artifact characteristics obtained by the five-wavelength calculation which indicates that most of the artifacts are eliminated.  
         [0106]      FIGS. 6A  to  6 C show artifact characteristics obtained in a case where the subject started moving his/her hand in the “waving” manner, after about 30 seconds from the start of measurement in a state where oxygen saturation of arterial blood had been lowered by holding his breath.  FIG. 6A  show artifact characteristics obtained by the two-wavelength calculation which indicates relatively large artifacts.  FIG. 6B  shows artifact characteristics obtained by the three-wavelength calculation which indicates improvements with respect to artifacts.  FIG. 6C  shows artifact characteristics obtained by the five-wavelength calculation which indicates further improvements with respect to artifacts.  
         [0107]     Meanwhile, in this case, attenuation changes ΔAi have a wide range of magnitude during periods of the artifacts. Therefore; attenuation changes of a predetermined level or lower have been removed in calculations. Such processing is simple and has been applied to numerous pulse oximeters. Accordingly, detailed descriptions of the processing will be omitted.  
         [0108]      FIGS. 7A  to  7 C show artifact characteristics obtained in a case where the subject started moving his/her hand in the “chopping” manner, after about 30 seconds from the start of measurement in a state where oxygen saturation of arterial blood had been lowered by holding his breath.  FIG. 7A  shows artifact characteristics obtained by the two-wavelength calculation which indicates relatively large artifacts.  FIG. 7B  shows artifact characteristics obtained by the three-wavelength calculation which Indicates little improvement in artifacts.  FIG. 7C  shows artifact characteristics obtained by the five-wavelength calculation which indicates that most of the artifacts are eliminated.  
         [0109]     Even in a case where a pulse wave in a living tissue of a subject is so small that a pulse oximetry is not applicable, according to the invention, measurement of an oxygen saturation of arterial blood is enabled by forcibly causing body motion against the living tissue by a vibration mechanism or the like, to thus generate artifacts on attenuation changes in a light signal transmitted through the living tissue.  
         [0110]     In addition, simultaneously with the oxygen saturation of arterial blood, an oxygen saturation of venous blood can also be measured.  
         [0111]     The preferred embodiment of the invention has been described; however, it should be understood that the present invention is not limited to the embodiment, and may be modified in various manners without departing from the scope of the invention.