Title:

Kind
Code:

A1

Abstract:

The present invention involves method and apparatus for analyzing two measured signals that are modeled as containing primary and secondary portions. Coefficients relate the two signals according to a model defined in accordance with the present invention. In one embodiment, the present invention involves utilizing a transformation which evaluates a plurality of possible signal coefficients in order to find appropriate coefficients. Alternatively, the present invention involves using statistical functions or Fourier transform and windowing techniques to determine the coefficients relating to two measured signals. Use of this invention is described in particular detail with respect to blood oximetry measurements.

Inventors:

Diab, Mohamed K. (Laguna Niguel, CA, US)

Kiani-azarbayjany, Esmaiel (Laguna Niguel, CA, US)

Elfadel, Ibrahim M. (Laguna Niguel, CA, US)

Mccarthy, Rex J. (Mission Viejo, CA, US)

Weber, Walter M. (Los Angeles, CA, US)

Smith, Robert A. (Corona, CA, US)

Kiani-azarbayjany, Esmaiel (Laguna Niguel, CA, US)

Elfadel, Ibrahim M. (Laguna Niguel, CA, US)

Mccarthy, Rex J. (Mission Viejo, CA, US)

Weber, Walter M. (Los Angeles, CA, US)

Smith, Robert A. (Corona, CA, US)

Application Number:

10/062859

Publication Date:

09/12/2002

Filing Date:

01/30/2002

Export Citation:

Assignee:

DIAB MOHAMED K.

KIANI-AZARBAYJANY ESMAIEL

ELFADEL IBRAHIM M.

MCCARTHY REX J.

WEBER WALTER M.

SMITH ROBERT A.

KIANI-AZARBAYJANY ESMAIEL

ELFADEL IBRAHIM M.

MCCARTHY REX J.

WEBER WALTER M.

SMITH ROBERT A.

Primary Class:

Other Classes:

600/336, 600/509, 600/529

International Classes:

View Patent Images:

Related US Applications:

20090227856 | Electrocardiogram sensor | September, 2009 | Russell et al. |

20090312652 | ELECTRONIC MANOMETER FOR APPROPRIATELY ADJUSTING INTERNAL PRESSURE OF CUFF AND METHOD FOR CONTROLLING THE SAME | December, 2009 | Yamakoshi et al. |

20030127090 | Active pump bronchial implant devices and methods of use thereof | July, 2003 | Gifford et al. |

20090275799 | AXIAL VISUALIZATION SYSTEMS | November, 2009 | Saadat et al. |

20070100247 | Combined wrist blood pressure and ecg monitor | May, 2007 | Platt et al. |

20040153006 | Intracorporeal devices with ionomeric polymer sleeves | August, 2004 | Vrba et al. |

20060211943 | Ultrasonic blade with terminal end balance features | September, 2006 | Beaupre |

20080045796 | HISTOFLUORESCENT STAIN COMPOSITION FOR ENDOSCOPY | February, 2008 | Yamamoto et al. |

20080146942 | Catheter Position Tracking Methods Using Fluoroscopy and Rotational Sensors | June, 2008 | Dala-krishna |

20100030031 | SWING PRISM ENDOSCOPE | February, 2010 | Goldfarb et al. |

20040254428 | Intermuscular guide for retractor insertion and method of use | December, 2004 | Ritland |

Primary Examiner:

WINAKUR, ERIC FRANK

Attorney, Agent or Firm:

KNOBBE MARTENS OLSON & BEAR LLP (620 NEWPORT CENTER DRIVE, NEWPORT BEACH, CA, 92660, US)

Claims:

1. A method of determining blood oxygen saturation comprising: sensing physiological signals resulting from the attenuation of light of at least first and second wavelengths by body tissue carrying pulsing blood; determining at least two values corresponding to oxygen saturation based upon at least two alternative methods of using the physiological signals; and determining a resulting value for oxygen saturation from the at least two values corresponding to oxygen saturation.

2. The method of claim 1, wherein the step of determining a resulting value comprises selecting from the at least two values at least one value that is a maximum among the at least two values.

3. The method of claim 1, wherein the step of determining a resulting value comprises averaging at least some of the at least two values.

4. The method of claim 1, wherein one of the alternative methods comprises at least one calculation in the frequency domain.

5. The method of claim 4, wherein the calculation in the frequency domain comprises performing a Fourier Transform on the physiological signals.

6. The method of claim 1, wherein at least one of the at least two alternative methods comprises a calculation based on a ratio of a normalized representation of the physiological signal resulting from the first wavelength to a normalized representation of the physiological signal resulting from the second wavelength.

7. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on the physiological signals after they have been effected by a recursive polyphase bandpass filter.

8. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on the physiological signals after they have been effected by an adaptive implementation of a recursive polyphase bandpass filter.

9. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on the physiological signals after they have been effected by a bank of filters.

10. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on the physiological signals after they have been effected by a sum of squares analysis.

11. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on a scan of possible saturation values.

12. The method of claim 11, wherein the calculation based on a scan of possible saturation values comprises a discrete saturation transform.

13. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on determining values for saturation that minimize the correlation between a signal portion and a noise portion of at least one of the physiological signals.

14. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on the physiological signals after they have been effected by a Kalman filter.

15. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on the physiological signals after they have been effected by a neural network.

16. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on the physiological signals after they have been effected with spectral estimation techniques.

17. The method of claim 6, wherein at least one of the at least two alternative methods comprises selecting at least one of the at least two values based on characteristics of the physiological signals indicative of the quality of the physiological signals.

18. The method of claim 6, wherein step of determining comprises averaging the resulting value over time, said averaging dependent upon characteristics of the physiological signals indicative of the quality of the physiological signal.

19. The method of claim 18, wherein the averaging is based on confidence in the quality of the physiological signals.

20. The method of claim 19, wherein the confidence is determined by analyzing whether there is significant motion noise present in the physiological signals.

21. The method of claim 6, wherein at least one of the at least two alternative methods comprises a calculation based on the physiological signals after they have been effected by an adaptive algorithm.

22. The method of claim 21, wherein at least one of the at least two alternative methods comprises a calculation based upon a scan of values potentially indicative of said physiological parameter.

23. A method of determining pulse rate comprising: sensing physiological signals resulting from the attenuation of light of at least first and second wavelengths by body tissue carrying pulsing blood; determining at least two values corresponding to pulse rate based upon at least two alternative methods of processing the physiological signals; and determining a resulting value for pulse rate from the at least two values corresponding to pulse rate.

24. The method of claim 23, wherein the step of determining comprises selecting at least one of the at least two values based on a determination of confidence in the accuracy of physiological signals.

25. The method of claim 23, wherein determining a resulting value comprises averaging the at least two values.

26. The method of claim 25, wherein said step of averaging comprises averaging over a time window, wherein said window is increased for potential of said physiological parameter having a lower confidence of accuracy and decreased for potential values of said physiological parameter having a higher confidence of accuracy.

Description:

[0001] This application is a continuation of application Ser. No. 09/195,791, filed Nov. 17, 1998, which is a continuation of application Ser. No. 08/859,837, filed May 16, 1997 (now U.S. Pat. No. 6,157,850), which is a continuation of application Ser. No. 08/320,154, filed Oct. 7, 1994 (now U.S. Pat. No. 5,632,272 issued May 27, 1997), which is a continuation-in-part of application Ser. No. 08/132,812, filed Oct. 6, 1993 (now U.S. Pat. No. 5,490,505 issued Feb. 13, 1996), which is a continuation-in-part of application Ser. No. 08/249,690 filed May 26, 1994 (now U.S. Pat. No. 5,482,036), which is a continuation of application Ser. No. 07/666,060, filed Mar. 7, 1991 (now abandoned).

[0002] 1. Field of the Invention

[0003] The present invention relates to the field of signal processing. More specifically, the present invention relates to the processing of measured signals, containing a primary signal portion and a secondary signal portion, for the removal or derivation of either the primary or secondary signal portion when little is known about either of these components. More particularly, the present invention relates to modeling the measured signals in a novel way which facilitates minimizing the correlation between the primary signal portion and the secondary signal portion in order to produce a primary and/or secondary signal. The present invention is especially useful for physiological monitoring systems including blood oxygen saturation systems.

[0004] 2. Description of the Related Art

[0005] Signal processors are typically employed to remove or derive either the primary or secondary signal portion from a composite measured signal including a primary signal portion and a secondary signal portion. For example, a composite signal may contain noise and desirable portions. If the secondary signal portion occupies a different frequency spectrum than the primary signal portion, then conventional filtering techniques such as low pass, band pass, and high pass filtering are available to remove or derive either the primary or the secondary signal portion from the total signal. Fixed single or multiple notch filters could also be employed if the primary and/or secondary signal portion(s) exist at a fixed frequency(s).

[0006] It is often the case that an overlap in frequency spectrum between the primary and secondary signal portions exists. Complicating matters further, the statistical properties of one or both of the primary and secondary signal portions change with time. In such cases, conventional filtering techniques are ineffective in extracting either the primary or secondary signal. If, however, a description of either the primary or secondary signal portion can be derived, correlation canceling, such as adaptive noise canceling, can be employed to remove either the primary or secondary signal portion of the signal isolating the other portion. In other words, given sufficient information about one of the signal portions, that signal portion can be extracted.

[0007] Conventional correlation cancelers, such as adaptive noise cancelers, dynamically change their transfer function to adapt to and remove portions of a composite signal. However, correlation cancelers require either a secondary reference or a primary reference which correlates to either the secondary signal portion only or the primary signal portion only. For instance, for a measured signal containing noise and desirable signal, the noise can be removed with a correlation canceler if a noise reference is available. This is often the case. Although the amplitude of the reference signals are not necessarily the same as the amplitude of the corresponding primary or secondary signal portions, they have a frequency spectrum which is similar to that of the primary or secondary signal portions.

[0008] In many cases, nothing or very little is known about the secondary and/or primary signal portions. One area where measured signals comprising a primary signal portion and a secondary signal portion about which no information can easily be determined is physiological monitoring. Physiological monitoring generally involves measured signals derived from a physiological system, such as the human body. Measurements which are typically taken with physiological monitoring systems include electrocardiographs, blood pressure, blood gas saturation (such as oxygen saturation), capnographs, other blood constituent monitoring, heart rate, respiration rate, electro-encephalograph (EEG) and depth of anesthesia, for example. Other types of measurements include those which measure the pressure and quantity of a substance within the body such as cardiac output, venous oxygen saturation, arterial oxygen saturation, bilirubin, total hemoglobin, breathalyzer testing, drug testing, cholesterol testing, glucose testing, extra vasation, and carbon dioxide testing, protein testing, carbon monoxide testing, and other in-vivo measurements, for example. Complications arising in these measurements are often due to motion of the patient, both external and internal (muscle movement, vessel movement, and probe movement, for example), during the measurement process.

[0009] Many types of physiological measurements can be made by using the known properties of energy attenuation as a selected form of energy passes through a medium.

[0010] A blood gas monitor is one example of a physiological monitoring system which is based upon the measurement of energy attenuated by biological tissues or substances. Blood gas monitors transmit light into the test medium and measure the attenuation of the light as a function of time. The output signal of a blood gas monitor which is sensitive to the arterial blood flow contains a component which is a waveform representative of the patient's arterial pulse. This type of signal, which contains a component related to the patient's pulse, is called a plethysmographic wave, and is shown in

[0011] Typically, a digit such as a finger, an ear lobe, or other portion of the body where blood flows close to the skin, is employed as the medium through which light energy is transmitted for blood gas attenuation measurements. The finger comprises skin, fat, bone, muscle, etc., shown schematically in

[0012] An example of a more realistic measured waveform S is shown in

[0013] A pulse oximeter is a type of blood gas monitor which non-invasively measures the arterial saturation of oxygen in the blood. The pumping of the heart forces freshly oxygenated blood into the arteries causing greater energy attenuation. As well understood in the art, the arterial saturation of oxygenated blood may be determined from the depth of the valleys relative to the peaks of two plethysmographic waveforms measured at separate wavelengths. Patient movement introduces motion artifacts to the composite signal as illustrated in the plethysmographic waveform illustrated in

[0014] This invention provides improvements upon the methods and apparatus disclosed in U.S. patent application Ser. No. 08/132,812, filed Oct. 6, 1993, entitled Signal Processing Apparatus, which earlier application has been assigned to the assignee of the instant application. The present invention involves several different embodiments using the novel signal model in accordance with the present invention to isolate either a primary signal portion or a secondary signal portion of a composite measured signal. In one embodiment, a signal processor acquires a first measured signal and a second measured signal that is correlated to the first measured signal. The first signal comprises a first primary signal portion and a first secondary signal portion. The second signal comprises a second primary signal portion and a second secondary signal portion. The signals may be acquired by propagating energy through a medium and measuring an attenuated signal after transmission or reflection. Alternatively, the signals may be acquired by measuring energy generated by the medium.

[0015] In one embodiment, the first and second measured signals are processed to generate a secondary reference which does not contain the primary signal portions from either of the first or second measured signals. This secondary reference is correlated to the secondary signal portion of each of the first and second measured signals. The secondary reference is used to remove the secondary portion of each of the first and second measured signals via a correlation canceler, such as an adaptive noise canceler. The correlation canceler is a device which takes a first and second input and removes from the first input all signal components which are correlated to the second input. Any unit which performs or nearly performs this function is herein considered to be a correlation canceler.

[0016] An adaptive correlation canceler can be described by analogy to a dynamic multiple notch filter which dynamically changes its transfer function in response to a reference signal and the measured signals to remove frequencies from the measured signals that are also present in the reference signal. Thus, a typical adaptive correlation canceler receives the signal from which it is desired to remove a component and receives a reference signal of the undesired portion. The output of the correlation canceler is a good approximation to the desired signal with the undesired component removed.

[0017] Alternatively, the first and second measured signals may be processed to generate a primary reference which does not contain the secondary signal portions from either of the first or second measured signals. The primary reference may then be used to remove the primary portion of each of the first and second measured signals via a correlation canceler. The output of the correlation canceler is a good approximation to the secondary signal with the primary signal removed and may be used for subsequent processing in the same instrument or an auxiliary instrument. In this capacity, the approximation to the secondary signal may be used as a reference signal for input to a second correlation canceler together with either the first or second measured signals for computation of, respectively, either the first or second primary signal portions.

[0018] Physiological monitors can benefit from signal processors of the present invention. Often in physiological measurements a first signal comprising a first primary portion and a first secondary portion and a second signal comprising a second primary portion and a second secondary portion are acquired. The signals may be acquired by propagating energy through a patient's body (or a material which is derived from the body, such as breath, blood, or tissue, for example) or inside a vessel and measuring an attenuated signal after transmission or reflection. Alternatively, the signal may be acquired by measuring energy generated by a patient's body, such as in electrocardiography. The signals are processed via the signal processor of the present invention to acquire either a secondary reference or a primary reference which is input to a correlation canceler, such as an adaptive noise canceler.

[0019] One physiological monitoring apparatus which benefits from the present invention is a monitoring system which determines a signal which is representative of the arterial pulse, called a plethysmographic wave. This signal can be used in blood pressure calculations, blood constituent measurements, etc. A specific example of such a use is in pulse oximetry. Pulse oximetry involves determining the saturation of oxygen in the blood. In this configuration, the primary portion of the signal is the arterial blood contribution to attenuation of energy as it passes through a portion of the body where blood flows close to the skin. The pumping of the heart causes blood flow to increase and decrease in the arteries in a periodic fashion, causing periodic attenuation wherein the periodic waveform is the plethysmographic waveform representative of the arterial pulse. The secondary portion is noise. In accordance with the present invention, the measured signals are modeled such that this secondary portion of the signal is related to the venous blood contribution to attenuation of energy as it passes through the body. The secondary portion also includes artifacts due to patient movement which causes the venous blood to flow in an unpredictable manner, causing unpredictable attenuation and corrupting the otherwise periodic plethysmographic waveform. Respiration also causes the secondary or noise portion to vary, although typically at a lower frequency than the patients pulse rate. Accordingly, the measured signal which forms a plethysmographic waveform is modeled in accordance with the present invention such that the primary portion of the signal is representative of arterial blood contribution to attenuation and the secondary portion is due to several other parameters.

[0020] A physiological monitor particularly adapted to pulse oximetry oxygen saturation measurement comprises two light emitting diodes (LED's) which emit light at different wavelengths to produce first and second signals. A detector registers the attenuation of the two different energy signals after each passes through an absorptive media, for example a digit such as a finger, or an earlobe. The attenuated signals generally comprise both primary (arterial attenuator) and secondary (noise) signal portions. A static filtering system, such as a bandpass filter, removes a portion of the secondary signal which is outside of a known bandwidth of interest, leaving an erratic or random secondary signal portion, often caused by motion and often difficult to remove, along with the primary signal portion.

[0021] A processor in accordance with one embodiment of the present invention removes the primary signal portions from the measured signals yielding a secondary reference which is a combination of the remaining secondary signal portions. The secondary reference is correlated to both of the secondary signal portions. The secondary reference and at least one of the measured signals are input to a correlation canceler, such as an adaptive noise canceler, which removes the random or erratic portion of the secondary signal. This yields a good approximation to a primary plethysmographic signal as measured at one of the measured signal wavelengths. As is known in the art, quantitative measurements of the amount of oxygenated arterial blood in the body can be determined from the plethysmographic signal in a variety of ways.

[0022] The processor of the present invention may also remove the secondary signal portions from the measured signals yielding a primary reference which is a combination of the remaining primary signal portions. The primary reference is correlated to both of the primary signal portions. The primary reference and at least one of the measured signals are input to a correlation canceler which removes the primary portions of the measured signals. This yields a good approximation to the secondary signal at one of the measured signal wavelengths. This signal may be useful for removing secondary signals from an auxiliary instrument as well as determining venous blood oxygen saturation.

[0023] In accordance with the signal model of the present invention, the two measured signals each having primary and secondary signal portions can be related by coefficients. By relating the two equations with respect to coefficients defined in accordance with the present invention, the coefficients provide information about the arterial oxygen saturation and about the noise (the venous oxygen saturation and other parameters). In accordance with this aspect of the present invention, the coefficients can be determined by minimizing the correlation between the primary and secondary signal portions as defined in the model. Accordingly, the signal model of the present invention can be utilized in many ways in order to obtain information about the measured signals as will be further apparent in the detailed description of the preferred embodiments.

[0024] One aspect of the present invention is a method for use in a signal processor in a signal processor for processing at least two measured signals S_{1 }_{2 }_{1 }_{2 }

_{1}_{1}_{1}

_{2}_{2}_{2}

[0025] where s_{1 }_{2}_{1 }_{2 }

_{1}_{a}_{2 }_{1}_{v}_{2}

[0026] and where r_{a }_{v }

[0027] The method comprises a number of steps. A value of coefficient r_{a }_{1 }_{1}_{a }

[0028] In one embodiment, the clean signal is displayed on a display. In another embodiment, wherein the first and second signals are physiological signals, the method further comprises the step of processing the clean signal to determine a physiological parameter from the first or second measured signals. In one embodiment, the parameter is arterial oxygen saturation. In another embodiment, the parameter is an ECG signal. In yet another embodiment, wherein the first portion of the measured signals is indicative of a heart plethysmograph, the method further comprises the step of calculating the pulse rate.

[0029] Another aspect of the present invention involves a physiological monitor. The monitor has a first input configured to receive a first measured signal S_{1 }_{1}_{1}_{2 }_{2 }_{2}_{1 }_{2 }

_{1}_{1}_{1}

_{2}_{2}_{2}

[0030] where s_{1 }_{2}_{1 }_{2 }

_{1}_{a}_{2 }_{1}_{v}_{2}

[0031] and where r_{a }_{v }

[0032] The monitor further has a scan reference processor, the scan reference processor responds to a plurality of possible values for r_{a }_{a }

[0033] In one embodiment, the plurality of possible values correspond to a plurality of possible values for a selected blood constituent. In one embodiment the, the selected blood constituent is arterial blood oxygen saturation. In another embodiment, the selected blood constituent is venous blood oxygen saturation. In yet another embodiment, the selected blood constituent is carbon monoxide.

[0034] Another aspect of the present invention involves a physiological monitor. The monitor has a first input configured to receive a first measured signal S_{1 }_{1}_{1}_{2 }_{2 }_{2}_{1 }_{2 }

_{1}_{1}_{1}

_{2}_{2}_{2}

[0035] where s_{1 }_{2}_{1 }_{2 }

_{1}_{a}_{2 }

[0036] and

_{1}_{v}_{s}

[0037] and where r_{a }_{v }

[0038] A transform module is responsive to the first and the second measured signals and responsive to a plurality of possible values for r_{a }_{a }_{a }

[0039]

[0040]

[0041]

[0042]

[0043]

[0044]

[0045]

[0046]

[0047]

[0048]

[0049]

[0050]

[0051] _{1}_{2}_{n}_{3}_{a }_{7}_{v}

[0052] _{1}_{2}_{n}_{3}_{a }_{7}_{v}

[0053]

[0054]

[0055]

[0056]

[0057]

[0058]

[0059]

[0060]

[0061]

[0062]

[0063]

[0064]

[0065]

[0066]

[0067]

[0068]

[0069]

[0070]

[0071]

[0072]

[0073]

[0074]

[0075] FIGS.

[0076] _{λa}_{λa}

[0077] _{IR}_{λb}_{λb}

[0078]

[0079] _{λa}_{λa}_{λa}

[0080] _{λb}_{λb}_{λb}

[0081] _{1}_{2 }_{3}

[0082] The present invention involves a system which utilizes first and second measured signals that each contain a primary signal portion and a secondary signal portion. In other words, given a first and second composite signals S_{1}_{1}_{1}_{2}_{2}_{2}

[0083] The system of the present invention is particularly useful where the primary and/or secondary signal portion n(t) may contain one or more of a constant portion, a predictable portion, an erratic portion, a random portion, etc. The primary signal approximation s″(t) or secondary signal approximation n″(t) is derived by removing as many of the secondary signal portions n(t) or primary signal portions s(t) from the composite signal S(t) as possible. The remaining signal forms either the primary signal approximation s″(t) or secondary signal approximation n″(t), respectively. The constant portion and predictable portion of the secondary signal n(t) are easily removed with traditional filtering techniques, such as simple subtraction, low pass, band pass, and high pass filtering. The erratic portion is more difficult to remove due to its unpredictable nature. If something is known about the erratic signal, even statistically, it could be removed, at least partially, from the measured signal via traditional filtering techniques. However, often no information is known about the erratic portion of the secondary signal n(t). In this case, traditional filtering techniques are usually insufficient.

[0084] In order to remove the secondary signal n(t), a signal model in accordance with the present invention is defined as follows for the first and second measured signals S_{1 }_{2}

[0085] where s_{1 }_{1 }_{2 }_{2 }_{1 }_{2 }_{a }_{v }

[0086] In accordance with one aspect of the present invention, this signal model is used in combination with a correlation canceler, such as an adaptive noise canceler, to remove or derive the erratic portion of the measured signals.

[0087] Generally, a correlation canceler has two signal inputs and one output. One of the inputs is either the secondary reference n′(t) or the primary reference s′(t) which are correlated, respectively, to the secondary signal portions n(t) and the primary signal portions s(t) present in the composite signal S(t). The other input is for the composite signal S(t). Ideally, the output of the correlation canceler s″(t) or n″(t) corresponds, respectively, to the primary signal s(t) or the secondary signal n(t) portions only. Often, the most difficult task in the application of correlation cancelers is determining the reference signals n′(t) and s′(t) which are correlated to the secondary n(t) and primary s(t) portions, respectively, of the measured signal S(t) since, as discussed above, these portions are quite difficult to isolate from the measured signal S(t). In the signal processor of the present invention, either a secondary reference n′(t) or a primary reference s′(t) is determined from two composite signals measured simultaneously, or nearly simultaneously, at two different wavelengths, λa and λb.

[0088] A block diagram of a generic monitor incorporating a signal processor according to the present invention, and a correlation canceler is shown in _{λa}_{λb}_{λa}_{λa}_{λa}_{λb}_{λa}_{λb}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0089] The signals S_{λa}_{λb}_{λb}_{a}_{λa}_{λb}_{v}_{λa}_{λb}_{λb}_{λa}_{a }_{v }_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{a}_{λb}_{λa}_{λb}_{λa}_{s}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0090] In one embodiment, an adaptive noise canceler _{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λa}_{λa}_{λb}_{λb}_{λb}

[0091] The adaptive noise canceler _{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λa}_{λa}_{λb}_{λb}_{λb}

[0092] The adaptive noise canceler _{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0093]

[0094] The adaptive noise canceler _{λa}_{λb}_{λa}_{λb}_{λ}_{λ}_{λa}_{λb}_{λa}_{λb}_{λ}_{λa}_{λb}_{λa}_{λa}_{λa}_{λa}_{λb}_{λb}_{λb}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λ}_{λa}_{λb}_{λa}_{λa}_{λa}_{λa}_{λb}_{λb}_{λb}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0095] One algorithm which may be used for the adjustment of the transfer function of the internal processor

[0096] Adaptive processors _{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0097] Generalized Determination of Primary and Secondary Reference Signals

[0098] An explanation which describes how the reference signals n′(t) and s′(t) may be determined follows. A first signal is measured at, for example, a wavelength λa, by a detector yielding a signal S_{λa}

_{λa}_{λa}_{λa}

[0099] where s_{λa}_{λa}

[0100] A similar measurement is taken simultaneously, or nearly simultaneously, at a different wavelength, λb, yielding:

_{λb}_{λb}_{λb}

[0101] Note that as long as the measurements, S_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0102] To obtain the reference signals n′(t) and s′(t), the measured signals S_{λa}_{λb}_{a }_{v}_{λa}_{λb}_{λa}_{λb}

_{λa}_{a}_{λb}

_{λa}_{v}_{λb}

[0103] In accordance with the inventive signal model of the present invention, these proportionality relationships can be satisfied in many measurements, including but not limited to absorption measurements and physiological measurements. Additionally, in accordance with the signal model of the present invention, in most measurements, the proportionality constants r_{a }_{v }

_{λa}_{a}_{λb}

_{λa}_{v}_{λb}

[0104] Multiplying equation (2) by r_{a }_{λa}_{λb}

_{λa}_{a}_{λb}_{λa}_{a}_{λb}

[0105] a non-zero signal which is correlated to each secondary signal portion n_{λa}_{λb}

[0106] Multiplying equation (2) by r_{v }_{λa}_{λb}

_{λa}_{v}_{λb}_{λa}_{v}_{λb}

[0107] a non-zero signal which is correlated to each of the primary signal portions s_{λa}_{λb}

[0108] Example of Determination of Primary and Secondary Reference Signals in an Absorptive System

[0109] Correlation canceling is particularly useful in a large number of measurements generally described as absorption measurements. An example of an absorption type monitor which can advantageously employ correlation canceling, such as adaptive noise canceling, based upon a reference n′(t) or s′(t) determined by a processor of the present invention is one which determines the concentration of an energy absorbing constituent within an absorbing material when the material is subject to change. Such changes can be caused by forces about which information is desired or primary, or alternatively, by random or erratic secondary forces such as a mechanical force on the material. Random or erratic interference, such as motion, generates secondary components in the measured signal. These secondary components can be removed or derived by the correlation canceler if a suitable secondary reference n′(t) or primary reference s′(t) is known.

[0110] A schematic N constituent absorbing material comprising a container _{1}_{2}_{3}_{N}_{1 }_{N }

[0111] Initially transforming the signal by taking the natural logarithm of both sides and manipulating terms, the signal is transformed such that the signal components are combined by addition rather than multiplication, i.e.:

[0112] where I_{0 }_{i,λa }^{th }_{i}^{th }^{th }_{i}^{th }_{i}_{1 }_{N }_{1}_{N}_{i}_{i}

[0113] When the material is not subject to any forces which cause change in the thicknesses of the layers, the optical path length of each layer, x_{i}_{i}_{i}

[0114] The correlation canceler may selectively remove from the composite signal, measured after being transmitted through or reflected from the absorbing material, either the secondary or the primary signal components caused by forces which perturb or change the material differently from the forces which perturbed or changed the material to cause respectively, either the primary or secondary signal component. For the purposes of illustration, it will be assumed that the portion of the measured signal which is deemed to be the primary signal s_{λa}_{5}_{5}_{5}_{5}_{5 }_{1 }_{4 }_{6 }_{N}_{5 }_{5 }_{5}_{1 }_{4 }_{6 }_{N }_{λa}_{5}_{5 }_{1 }_{4 }_{6 }_{N}

[0115] It is often the case that the total perturbation affecting the layers associated with the secondary signal components is caused by random or erratic forces. This causes the thickness of layers to change erratically and the optical path length of each layer, x_{i}_{λa}_{λa}_{λa}_{5 }_{5}_{λa}_{λa}_{5}_{5}

[0116] The correlation canceler utilizes either the secondary reference n′(t) or the primary reference s′(t) determined from two substantially simultaneously measured signals S_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa(t)}_{5,λa}_{5}_{5}_{λa}

[0117]
_{λb(t)}_{5,λb}_{5}_{5}_{λb}

[0118] Further transformations of the signals are the proportionality relationships in accordance with the signal model of the present invention defining r_{a }_{v}

_{5,λa}_{a}_{5,λb}

_{λa}_{v}_{λb}

[0119] where

_{λa}_{a}_{λb}

_{5,λa}_{v}_{5,λb}

[0120] It is often the case that both equations (12) and (13) can be simultaneously satisfied. Multiplying equation (11) by r_{a }

_{λa}_{a}_{λb}_{λa}_{a}_{λb}

[0121]

[0122] Multiplying equation (11) by r_{v }

_{λa}_{v}_{λb}_{λa}_{v}_{λb}

_{5}_{5}_{5,λa}_{v}_{5}_{5}_{5,λb}

_{5}_{5}_{5,λa}_{v}_{5,λb}

[0123] A sample of either the secondary reference n′(t) or the primary reference s′(t), and a sample of either measured signal S_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{5,λa}_{5}_{5}_{λb}_{5,λb}_{5}_{5}_{λa}_{λa}_{λb}_{λb}_{5}_{λa}_{λb}

_{5}_{λa}_{5,λa}_{5}

[0124] or

_{5}_{λb}_{5,λb}_{5}

[0125] As discussed previously, the absorption coefficients are constant at each wavelength λa and λb and the thickness of the primary signal component, x_{5}_{5}_{5}

[0126] Determination of Concentration or Saturation in a Volume Containing More Than one Constituent

[0127] Referring to _{5 }_{6 }_{5,6}_{5}_{6}_{5 }_{6 }_{5 }_{6, }_{5}_{6}

[0128] Often it is desirable to find the concentration or the saturation, i.e., a percent concentration, of one constituent within a given thickness which contains more than one constituent and is subject to unique forces. A determination of the concentration or the saturation of a constituent within a given volume may be made with any number of constituents in the volume subject to the same total forces and therefore under the same perturbation or change. To determine the saturation of one constituent in a volume comprising many constituents, as many measured signals as there are constituents which absorb incident light energy are necessary. It will be understood that constituents which do not absorb light energy are not consequential in the determination of saturation. To determine the concentration, as many signals as there are constituents which absorb incident light energy are necessary as well as information about the sum of concentrations.

[0129] It is often the case that a thickness under unique motion contains only two constituents. For example, it may be desirable to know the concentration or saturation of A_{5 }_{5 }_{6}_{λa}_{λb}_{5 }_{6 }_{5 }_{6 }_{5 }_{5 }_{6 }_{5}_{6}_{λa}_{λb}

_{λa}_{5,λa}_{5}_{5,6}_{6,λa}_{6}_{5,6}_{λa}

_{λa}_{λa}

_{λb}_{5,λb}_{5}_{5,6}_{6,λb}_{6}_{5,6}_{λb}

_{λb}_{λb}

[0130] It is also often the case that there may be two or more thicknesses within a medium each containing the same two constituents but each experiencing a separate motion as in _{5 }_{5 }_{6 }_{3 }_{3 }_{4}_{3 }_{4 }_{5 }_{6 }_{λa}_{λb}_{5 }_{6 }_{λa}_{λb}_{3 }_{4}_{3 }_{4}_{3,λa}_{5,λa}_{3,λb}_{5,λb}_{4,λa}_{6,λa }_{4,λb}_{6,λb}_{3 }_{4 }_{5 }_{6 }

_{λa}_{5,λa}_{5}_{6,λa}_{6}_{5,6}

[0131]
_{λa}_{5,λa}_{3}_{6,λa}_{4}_{3,4}_{λa}

_{λb}_{5,λb}_{5}_{6,λb}_{6}_{5,6}

[0132]
_{λb}_{5,λb}_{3}_{6,λb}_{4}_{3,4}_{λb}

[0133] where signals n_{λa}_{λb}_{λa}_{λb}

[0134] Any signal portions whether primary or secondary, outside of a known bandwidth of interest, including the constant undesired secondary signal portion resulting from the generally constant absorption of the constituents when not under perturbation, should be removed to determine an approximation to either the primary signal or the secondary signal within the bandwidth of interest. This is easily accomplished by traditional band pass filtering techniques. As in the previous example, it is often the case that the total perturbation or change affecting the layers associated with the secondary signal components is caused by random or erratic forces, causing the thickness of each layer, or the optical path length of each layer, x_{i}_{λa}_{λa}_{λa}_{5 }_{6 }_{5 }_{6}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0135] Determination of Primary and Secondary Reference Signals for Saturation Measurements

[0136] One method for determining reference signals s′(t) or n′(t) from the measured signals S_{λa}_{λb}_{5 }_{5 }_{6 }_{3 }_{3 }_{4 }

_{5}_{5}_{5}_{6}

_{3}_{3}_{3}_{4}

_{5}_{6}_{5}^{−1}

_{3}_{4}_{3}^{−1}

[0137] are substantially constant over many samples of the measured signals S_{λa }_{λb}

[0138] The constant saturation assumption is equivalent to assuming that:

_{5}_{6}_{1}

_{3}_{4}_{2}

[0139] since the only other term in equations (23a) and (23b) is a constant, namely the numeral 1.

[0140] Using this assumption, the proportionality constants r_{a }_{v }_{λa}_{λb}

[0141]
_{λa}_{λb}

[0142]
_{λa}_{λb}_{4}

_{λa}_{v}_{λb}

[0143] In accordance with the present invention, it is often the case that both equations (26) and (30) can be simultaneously satisfied to determine the proportionality constants r_{a }_{v}_{5,λa}_{6,λa}_{5,λb}_{6,λb }_{5}_{6}_{3}_{4}_{a }_{v }_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{a }_{v}_{λa}_{λb}

[0144] Multiplying equation (19) by r_{a }

_{λa}_{a}_{λa}_{a}_{λb}

[0145] Multiplying equation (19) by r_{v }

_{λa}_{v}_{λb}_{λa}_{v}_{λb}

[0146] When using the constant saturation method in patient monitoring, initial proportionality coefficients can be determined as further explained below. It is not necessary for the patient to remain motionless even for an initialization period. With values for the proportionality coefficients r_{a }_{v }

[0147] Determination of Signal Coefficients for Primary and Secondary Reference Signals Using the Constant Saturation Method

[0148] In accordance with one aspect of the present invention, the reference processor _{λb}_{λb}_{λb}_{1}_{2}_{n }_{λa}_{λa}_{λa}

_{λa}_{λb}

[0149] for r=r_{1}_{2}_{n }

[0150] In order to determine either the primary reference s′(t) or the secondary reference n′(t) from the above plurality of reference signals of equation (32), signal coefficients r_{a }_{v }_{1}_{2}_{n}_{a }_{v }_{λa}_{λb}_{λa}_{λb}

_{λa}_{a}_{λb}

_{λa}_{v}_{λb}

_{a}_{λa}_{a}_{λb}

_{v}_{λa}_{v}_{λb}

[0151] In other words, coefficients r_{a }_{v }_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{1}_{2}_{n }_{a}_{λa}_{λb}_{v}_{λa}_{λb}

[0152] One approach to determine the signal coefficients r_{a }_{v }_{1}_{2}_{n }_{λa}_{λb}_{1}_{2}_{n}_{1}_{2}_{n}_{a }_{v }_{1}_{2}_{n }_{a }_{v }_{1}_{2}_{n }_{a }_{v}

[0153] Use of a plurality of coefficients in the processor of the present invention in conjunction with a correlation canceler _{a }_{v }

[0154] With properties (1), (2) and (3) it is easy to demonstrate that the energy or power output of a correlation canceler with a first input which corresponds to one of the measured signals S_{λa}_{λb}_{1}_{2}_{n}_{a }_{v }_{λa}_{1}_{2}_{n}_{λa}_{j}

_{λa}_{λa}_{j}_{λb}_{λa}_{j}_{λb}

[0155] where j=1, 2, . . . n and we have used the expressions

_{λa}_{λb}_{λa}

_{λa}_{λa}

_{λb}_{λb}_{λb}

[0156] The use of property (3) allows one to expand equation (35) into two terms

_{λa}_{λa}_{λa}_{λb}_{λa}_{λb}

_{λa}_{λa}_{λb}_{λa}_{λb}

[0157] so that upon use of properties (1) and (2) the correlation canceler output is given by

_{λa}_{j}_{λa}_{j}_{a}_{λa}_{j}_{v}

[0158] where δ(x) is the unit impulse function

[0159] The time variable, t, of the correlation canceler output C(S_{λa}_{j}

_{λa}_{j}^{2}_{λa}_{j}

_{j}_{a}^{2}_{λa}_{j}_{v}^{2}_{λa}

[0160] It should be understood that one could, equally well, have chosen the measured signal S_{λb}_{1}_{2}_{n}

_{λb}_{j}^{2}_{λb}

_{j}_{a}^{2}_{λb}_{j}_{v}^{2}_{λb}

[0161] It should also be understood that in practical situations the use of discrete time measurement signals may be employed as well as continuous time measurement signals. A system which performs a discrete transform (e.g., a saturation transform in the present example) in accordance with the present invention is described with reference to FIGS.

[0162] where t_{i }^{th }_{0 }_{n }

[0163] The energy functions given above, and shown in _{λa}_{λb}_{1}_{2}_{n}_{j }_{λa}_{λb}_{λa}_{λb}_{j}_{a }_{v}

[0164] It should be understood that there may be instances in time when either the primary signal portions s_{λa}_{λb}_{λa}_{λb}

[0165] Since there may be more than one signal coefficient value which provides maximum correlation canceler energy or power output, an ambiguity may arise. It may not be immediately obvious which signal coefficient together with the reference function R′(r, t) provides either the primary or secondary reference. In such cases, it is necessary to consider the constraints of the physical system at hand. For example, in pulse oximetry, it is known that arterial blood, whose signature is the primary plethysmographic wave, has greater oxygen saturation than venous blood, whose signature is the secondary erratic or random signal. Consequently, in pulse oximetry, the ratio of the primary signals due to arterial pulsation r_{a}_{λa}_{λb}_{v}_{λa}_{λb}

[0166] It should also be understood that in practical implementations of the plurality of reference signals and cross correlator technique, the ideal features listed as properties (1), (2) and (3) above will not be precisely satisfied but will be approximations thereof. Therefore, in practical implementations of this embodiment of the present invention, the correlation canceler energy curves depicted in

[0167] It should also be understood that it is possible to have more than two signal coefficient values which produce maximum energy or power output from a correlation canceler. This situation arises when the measured signals each contain more than two components each of which are related by a ratio as follows:

[0168] Thus, reference signal techniques together with a correlation cancellation, such as an adaptive noise canceler, can be employed to decompose a signal into two or more signal components each of which is related by a ratio.

[0169] Preferred Correlation Canceler Using a Joint Process Estimator Implementation

[0170] Once either the secondary reference n′(t) or the primary reference s′(t) is determined by the processor of the present invention, the correlation canceler can be implemented in either hardware or software. The preferred implementation of a correlation canceler is that of an adaptive noise canceler using a joint process estimator.

[0171] The least mean squares (LMS) implementation of the internal processor

[0172] The function of the joint process estimator is to remove either the secondary signal portions n_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0173] The joint process estimator _{λa}_{λb}_{λa}_{λa}_{λa}_{λb}_{λb}_{λb}

[0174] The joint process estimator _{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}_{λa}

[0175] The joint process estimator ^{th}^{th}^{th}_{λa}_{λa}^{th}

[0176] The least-squares lattice predictor _{f,m}_{b,m}_{m}_{m}

[0177] For each set of samples, i.e. one sample of the reference signal n′(t) or s′(t) derived substantially simultaneously with one sample of the measured signal S_{λa}_{0}_{0}_{0}_{0}_{f,1}_{1}_{0}_{b,1}_{0}_{1}_{m−1}_{m−1}_{m}_{m}

[0178] The backward prediction error b_{m}_{m,λa}_{0}_{0,λa}_{λa}_{1,λa}_{1,λa}_{1,λa}_{1}_{1,λa}_{1,λa}_{2,λa}_{2,λa}_{λa}_{λa}

[0179] The same processes are repeated in the least-squares lattice predictor _{m,λa}_{λa}_{λa}_{m}_{m}_{m,λa}_{f,m}_{b,m}_{m,λa}_{m}_{m}_{m,λa}_{f,m}_{b,m}_{m,λa}

[0180] Intermediate variables include a weighted sum of the forward prediction error squares ℑ_{m}_{m}_{m}_{m}_{m,λa}_{m}

[0181] where λ without a wavelength identifier, a or b, is a constant multiplicative value unrelated to wavelength and is typically less than or equal to one, i.e., λ≦1. The weighted sum of the backward prediction errors β_{m}

[0182] where, again, λ without a wavelength identifier, a or b, is a constant multiplicative value unrelated to wavelength and is typically less than or equal to one, i.e., λ≦1. These weighted sum intermediate error signals can be manipulated such that they are more easily solved for, as described in Chapter 9, §9.3 of the Haykin book referenced above and defined hereinafter in equations (59) and (60).

[0183] Description of the Joint Process Estimator

[0184] The operation of the joint process estimator _{m−1}_{m−1}_{m−1}_{m,λa}_{0,λa}

_{m−1}

_{m−1}

_{m−1}

_{m,λa}

_{0,λa}_{λa}

[0185] After initialization, a simultaneous sample of the measured signal S_{λa}_{λb}_{0}_{0}_{0}_{0}_{0}

_{0}_{0}

_{0}_{0}_{0}^{2}

_{0}

[0186] if a secondary reference n′(t) is used or according to:

_{0}_{0}

_{0}_{0}_{0}^{2}

^{65 }_{0}

[0187] if a primary reference s′(t) is used where, again, λ without a wavelength identifier, a or b, is a constant multiplicative value unrelated to wavelength.

[0188] Forward reflection coefficient Γ_{f,m}_{b,m}_{m,λa}_{f,1}_{b,1}_{1,λa}_{m−1}_{f,m}_{b,m}_{m}_{m}_{f,m}_{b,m}_{m}_{m,λa}_{m,λa}_{m+1λa}

_{m−1}_{m−1}_{m−1}_{m−1}_{m−1}

_{f,m}_{m−1}_{m−1}

_{b,m}_{m−1}_{m−1}

_{m}_{m−1}_{f,m}_{m−1}

_{m}_{m−1}_{b,m}_{m−1}

_{m}_{m−1}_{m−1}^{2}_{m−1}

_{m}_{m−1}_{m−1}^{2}_{m−1}

_{m}_{m−1}_{m−1}^{2}_{m−1}

_{m,λa}_{m,λa}_{m}_{m,λa}_{m}

_{m,λa}_{m,λa}_{m}

_{m+1,λa}_{m,λa}_{m}_{m}

[0189] where a (*) denotes a complex conjugate.

[0190] These equations cause the error signals f_{m}_{m}_{m,λa}_{m−1}

[0191] After a good approximation to either the primary signal s_{λa}_{λa}_{λa}_{f,m}_{b,m}_{m,λa}_{λa}_{λa}_{λa}

[0192] In a more numerically stable and preferred embodiment of the above described joint process estimator, a normalized joint process estimator is used. This version of the joint process estimator normalizes several variables of the above-described joint process estimator such that the normalized variables fall between −1 and 1. The derivation of the normalized joint process estimator is motivated in the Haykin text as problem 12 on page 640 by redefining the variables defined according to the following conditions:

[0193] This transformation allows the conversion of Equations (54)-(64) to the following normalized equations:

_{m−1}_{m−1}_{m−1}^{2}^{½}_{m−1}^{2}^{½+}_{m−1}_{m−1}

[0194]
_{m}_{m−1}^{2}_{m−1}

_{m}_{m−1}_{m−1}^{2}

[0195]
_{m+1,λa}_{m,.λa}_{m}_{m}

[0196] Initialization of Normalized Joint Process Estimator

[0197] Let N(t) be defined as the reference noise input at time index n and U(t) be defined as combined signal plus noise input at time index t the following equations apply (see Haykin, p. 619):

[0198] 1. To initialize the algorithm, at time t=0 set

_{m−1}

_{m−1}^{−6}

_{0}

[0199] 2. At each instant t≧1, generate the various zeroth-order variables as follows:

_{0}

_{0}_{0}^{2}

[0200] 3. For regression filtering, initialize the algorithm by setting at time index t=0

_{m}

[0201] 4. At each instant t≧1, generate the zeroth-order variable

_{0}

[0202] Accordingly, a normalized joint process estimator can be used for a more stable system.

[0203] In yet another embodiment, the correlation cancellation is performed with a QRD algorithm as shown diagrammatically in

[0204] The following equations adapted from the Haykin book correspond to the QRD-LSL diagram of

[0205] Computations

[0206] a. Predictions: For time t=1, 2, . . . , and prediction order m=1, 2, . . . , M, where M is the final prediction order, compute:

_{m−1}_{m−1}_{b,m−1}^{2 }_{f,m}_{b,m−1}_{f,m−1}^{·}_{b,m−1}^{½}^{·}_{f,m−1}

_{f,m−1}^{·}_{b,m−1}^{½}_{f,m−1}^{·}_{b,m−1}_{f,m−1}

_{m}^{½}_{b,m−1}_{m−1}^{½}

_{m−1}_{m−1}_{f,m−1}^{2}

[0207]
_{b,m}_{f,m−1}_{b,m−1}_{f,m−1}^{·}^{½}_{b,m−1}^{·}

_{b,m−1}^{19 }_{f,m−1}^{½}_{b,m−1}^{·}_{f,m−1}_{b,m−1}

[0208] b. Filtering: For order m=0, 1, . . . , M−1; and time t=1, 2, . . . , compute

_{m}_{m}_{b,m}^{2 }_{m+1}_{b,m}_{m}_{b,m}^{·}^{½}_{m}^{·(}

_{m}^{·}_{b,m}^{½}_{m}^{·}_{b,m}_{m}

[0209] 5. Initialization

[0210] a. Auxiliary parameter initialization: for order m=1, 2, . . . , M, set

_{f,m−1}_{b,m−1}

_{m}

[0211] b. Soft constraint initialization: For order m=0, 1, . . . , M, set

_{m}

_{m}

[0212] where δ is a small positive constant.

[0213] c. Data initialization: For t=1, 2, . . . , compute

_{f,0}_{b,0}

_{0}

_{0}

[0214] where μ(t) is the input and d(t) is the desired response at time t.

[0215] Flowchart of Joint Process Estimator

[0216] In a signal processor, such as a physiological monitor incorporating a reference processor of the present invention to determine a reference n′(t) or s′(t) for input to a correlation canceler, a joint process estimator

[0217] A flow chart of a subroutine to estimate the primary signal portion s_{λa}_{λa}_{λa }

[0218] A one-time initialization is performed when the physiological monitor is powered-on, as indicated by an “INITIALIZE NOISE CANCELER” action block

[0219] Next, a set of simultaneous samples of the composite measured signals S_{λa}_{λb}^{−1}_{0}_{1}

[0220] Then, using the set of measured signal samples S_{λa}_{λb}

[0221] A zero-stage order update is performed next as indicated in a “ZERO-STAGE UPDATE” action block _{0}_{0}_{m}_{m}

[0222] Next, a loop counter, m, is initialized as indicated in a “m=0” action block

[0223] Within the loop, the forward and backward reflection coefficient Γ_{f,m}_{b,m}

[0224] The calculation of regression filter register _{m,λa}_{m}_{m}

[0225] A new set of samples of the two measured signals S_{λa}_{λb}_{λa}_{λb}_{λa}_{λa}_{λb}_{λb}

[0226] A corresponding flowchart for the QRD algorithm of

[0227] Calculation of Saturation from Correlation Canceler Output

[0228] Physiological monitors may use the approximation of the primary signals s″_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0229] A joint process estimator _{λa}_{λb}_{m}

[0230] The second regression filter

_{m,λb}_{·m,λb}_{m}_{m,λb}_{m }

_{0,λb}

[0231] The second regression filter _{m+1,λa}

_{m+1,λb}_{m,λb}_{m,λb}_{m}

_{0,λb}_{λb}

[0232] The second regression filter has a regression coefficient κ_{m,λb}

_{m,λb}_{m,λb}_{m}

[0233] These values are used in conjunction with those intermediate variable values, signal values, register and register values defined in equations (46) through (64). These signals are calculated in an order defined by placing the additional signals immediately adjacent a similar signal for the wavelength λa.

[0234] For the constant saturation method, S_{λb}_{λb}_{λb}

[0235] The addition of the second regression filter ^{th }^{th }^{th }

[0236] An alternative diagram for the joint process estimator of

[0237] Calculation of Saturation

[0238] Once good approximations to the primary signal portions s″_{λa}_{λb}_{λa}_{λb}_{5 }_{5 }_{6}

_{λa}_{5,λa}_{5}_{5,6}_{6,λa}_{6}_{5,6}

_{λb}_{5,λb}_{5}_{5,6}_{6,λb}_{6}_{5,6}

[0239] Equations (70) and (71) are equivalent to two equations having three unknowns, namely c_{5}_{6}_{5,6}_{1 }_{2 }_{5 }_{5 }_{6 }_{3 }_{3 }_{4 }_{1 }_{2}

_{λa}_{1}_{5,λa}_{5}_{5,6}_{1}_{6,λa}_{6}_{5,6}_{1}

_{λb}_{1}_{5,λb}_{5}_{5,6}_{1}_{6,λb}_{6}_{5,6}_{1}

_{λa}_{2}_{5,λa}_{5}_{5,6}_{2}_{6,λa}_{6}_{5,6}_{2}

_{λb}_{2}_{5,λb}_{5}_{5,6}_{2}_{6,λb}_{6}_{5,6}_{2}

[0240] Then, difference signals may be determined which relate the signals of equations (72) through (75), i.e.:

_{λa}_{λa}_{1}_{λa}_{2}_{5,λa}_{5}_{6,λa}_{6}

[0241] and

_{λb}_{λb}_{1}_{λb}_{2}_{5,λb}_{5}_{,λb}_{6}

[0242] where Δx=x_{5,6}_{1}_{5,6}_{2}_{1}_{2}

_{5}_{5}_{6}

[0243] It will be understood that the Δx term drops out from the saturation calculation because of the division. Thus, knowledge of the thickness of the primary constituents is not required to calculate saturation.

[0244] Pulse Oximetry Measurements

[0245] A specific example of a physiological monitor utilizing a processor of the present invention to determine a secondary reference n′(t) for input to a correlation canceler that removes erratic motion-induced secondary signal portions is a pulse oximeter. Pulse oximetry may also be performed utilizing a processor of the present invention to determine a primary signal reference s′(t) which may be used for display purposes or for input to a correlation canceler to derive information about patient movement and venous blood oxygen saturation.

[0246] A pulse oximeter typically causes energy to propagate through a medium where blood flows close to the surface for example, an ear lobe, or a digit such as a finger, a forehead or a fetus' scalp. An attenuated signal is measured after propagation through or reflected from the medium. The pulse oximeter estimates the saturation of oxygenated blood.

[0247] Freshly oxygenated blood is pumped at high pressure from the heart into the arteries for use by the body. The volume of blood in the arteries varies with the heartbeat, giving rise to a variation in absorption of energy at the rate of the heartbeat, or the pulse.

[0248] Oxygen depleted, or deoxygenated, blood is returned to the heart by the veins along with unused oxygenated blood. The volume of blood in the veins varies with the rate of breathing, which is typically much slower than the heartbeat. Thus, when there is no motion induced variation in the thickness of the veins, venous blood causes a low frequency variation in absorption of energy. When there is motion induced variation in the thickness of the veins, the low frequency variation in absorption is coupled with the erratic variation in absorption due to motion artifact.

[0249] In absorption measurements using the transmission of energy through a medium, two light emitting diodes (LED's) are positioned on one side of a portion of the body where blood flows close to the surface, such as a finger, and a photodetector is positioned on the opposite side of the finger. Typically, in pulse oximetry measurements, one LED emits a visible wavelength, preferably red, and the other LED emits an infrared wavelength. However, one skilled in the art will realize that other wavelength combinations could be used. The finger comprises skin, tissue, muscle, both arterial blood and venous blood, fat, etc., each of which absorbs light energy differently due to different absorption coefficients, different concentrations, different thicknesses, and changing optical pathlengths. When the patient is not moving, absorption is substantially constant except for the flow of blood. The constant attenuation can be determined and subtracted from the signal via traditional filtering techniques. When the patient moves, this causes perturbation such as changing optical pathlength due to movement of background fluids (e.g., venous blood having a different saturation than the arterial blood). Therefore, the measured signal becomes erratic. Erratic motion induced noise typically cannot be predetermined and/or subtracted from the measured signal via traditional filtering techniques. Thus, determining the oxygen saturation of arterial blood and venous blood becomes more difficult.

[0250] A schematic of a physiological monitor for pulse oximetry is shown in FIGS.

[0251] The front end analog signal conditioning circuitry

[0252] The signal processing system also provides an emitter current control output

[0253]

[0254] The preferred driver depicted in

[0255] The voltage reference is also chosen as a low noise DC voltage reference for the digital to analog conversion circuit

[0256] In the present embodiment, the output of the voltage to current converters

[0257] In general, the red and infrared light emitters

[0258] It should be understood that in different embodiments of the present invention, one or more of the outputs may be provided. The digital signal processing system

[0259] The light signal is attenuated (amplitude modulated) by the pumping of blood through the finger

[0260] The composite time division signal is provided to the front analog signal conditioning circuitry

[0261] In the present embodiment, the preamplifier

[0262] The output of the preamplifier

[0263] The output of the amplifier

[0264] The programmable gain amplifier is also advantageous in an alternative embodiment in which the emitter drive current is held constant. In the present embodiment, the emitter drive current is adjusted for each patient in order to obtain the proper dynamic range at the input of the analog to digital conversion circuit

[0265] The output of the programmable gain amplifier

[0266] The output of the low-pass filter

[0267] In one advantageous embodiment, the first analog-to-digital converter

[0268] The second analog-to-digital converter

[0269] In addition, by using a single-channel converter, there is no need to tune two or more channels to each other. The delta-sigma converter is also advantageous in that it exhibits noise shaping, for improved noise control. An exemplary analog to digital converter is a Crystal Semiconductor CS5317. In the present embodiment, the second analog to digital converter

[0270] The digital signal processing system

[0271] The microcontroller

[0272] The microcontroller

[0273] FIGS.

[0274] In general, the demodulation operation separates the red and infrared signals from the composite signal and removes the 625 Hz carrier frequency, leaving raw data points. The raw data points are provided at 625 Hz intervals to the decimation operation which reduces the samples by an order of 10 to samples at 62.5 Hz. The decimation operation also provides some filtering on the samples. The resulting data is subjected to statistics and to the saturation transform operations in order to calculate a saturation value which is very tolerant to motion artifacts and other noise in the signal. The saturation value is ascertained in the saturation calculation module

[0275]

[0276] Because the signal processing system

[0277] A sum of the last four samples from each packet is then calculated, as represented in the summing operations

[0278] It should be understood that the 625 Hz carrier frequency has been removed by the demodulation operation

[0279]

[0280]

[0281] As represented in

[0282] Once the DC signal is removed, the signals are subjected to bandpass filtering, as represented in red and infrared Bandpass Filter modules

[0283] After filtering, the last 120 samples from each packet (of now 270 samples in the present embodiment) are selected for further processing as represented in Select Last 120 Samples modules

[0284] Conventional saturation equation calculations are performed on the red and infrared 120-sample packets. In the present embodiment, the conventional saturation calculations are performed in two different ways. For one calculation, the 120-sample packets are processed to obtain their overall RMS value, as represented in the first red and infrared RMS modules _{red}_{IR }

[0285] In addition to the conventional saturation operation

[0286] If the cross correlation is too low, the oximeter

[0287] The red and infrared 120-sample packets are also subjected to a second saturation operation and cross correlation in the same manner as described above, except the 120 samples are divided into 5 equal bins of samples (i.e., 5 bins of 24 samples each). The RMS, ratio, saturation, and cross correlation operations are performed on a bin-by-bin basis. These operations are represented in the Divide Into Five Equal Bins modules

[0288]

[0289] As depicted in

[0290] It should be understood that the scan values could be chosen to provide higher or lower resolution than 117 scan values. The scan values could also be non-uniformly spaced.

[0291] As illustrated in _{n}_{n}

[0292] The ratio “r_{n}_{n}_{n}

[0293] In other words, assuming that the red and infrared sample packets represent the red S_{red}_{IR}

_{ir}_{n}_{red}

[0294] In the present embodiment, the reference signal vectors and the infrared signal are provided as input to the DC removal module

[0295] The bandpass filter

[0296] It should be understood that the red and infrared sample packets may be switched in their use in the reference processor

[0297] The outputs of the reference processor

[0298] The joint process estimator also receives a lambda input ^{−6}

[0299] The joint process estimator

[0300] The joint process estimator

[0301] The Master Power Curve module

[0302] A corresponding transform is completed by the Bin Power Curves module

[0303] In general, in accordance with the signal model of the present invention, there will be two peaks in the power curves, as depicted in _{a}_{a }_{v }_{v }

[0304] In order to obtain arterial oxygen saturation, the peak in the power curves corresponding to the highest saturation value could be selected. However, to improve confidence in the value, further processing is completed.

[0305] The saturation calculation module

[0306] 0.014964670230367

[0307] 0.098294046682706

[0308] 0.204468276324813

[0309] 2.717182664241813

[0310] 5.704485606695227

[0311] 0.000000000000000

[0312] −5.704482606695227

[0313] −2.717182664241813

[0314] −0.204468276324813

[0315] −0.098294046682706

[0316] −0.014964670230367

[0317] This filter performs the differentiation and smoothing. Next, each point in the original power curve in question is evaluated and determined to be a possible peak if the following conditions are met: (1) the point is at least 2% of the maximum value in the power curve; (2) the value of the first derivative changes from greater than zero to less than or equal to zero. For each point that is found to be a possible peak, the neighboring points are examined and the largest of the three points is considered to be the true peak.

[0318] The peak width for these selected peaks is also calculated. The peak width of a power curve in question is computed by summing all the points in the power curve and subtracting the product of the minimum value in the power curve and the number of points in the power curve. In the present embodiment, the peak width calculation is applied to each of the bin power curves. The maximum value is selected as the peak width.

[0319] In addition, the infrared RMS value from the entire snapshot, the red RMS value, the seed saturation value for each bin, and the cross correlation between the red and infrared signals from the statistics module

[0320] If the correlation between the red and infrared signals is too low, the bin is discarded. If the saturation value of the selected peak for a given bin is lower than the seed saturation for the same bin, the peak is replaced with the seed saturation value. If either red or infrared RMS value is below a very small threshold, the bins are all discarded, and no saturation value is provided, because the measured signals are considered to be too small to obtain meaningful data. If no bins contain acceptable data, the exception handling module

[0321] If some bins qualify, those bins that qualify as having acceptable data are selected, and those that do not qualify are replaced with the average of the bins that are accepted. Each bin is given a time stamp in order to maintain the time sequence. A voter operation

[0322] The clip and smooth operation

[0323] In the presently preferred embodiment, the clip and smooth filter

[0324] During high confidence (no motion), the smoothing filter is a simple one-pole or exponential smoothing filter which is computed as follows:

[0325] where x(n) is the clipped new saturation value, and y(n) is the filtered saturation value.

[0326] During motion condition, a three-pole IIR (infinite impulse response) filter is used. Its characteristics are controlled by three time constants t_{a}_{b}_{c }

_{0}

_{1}_{b}_{c}_{a}_{b}

_{2}_{b}_{c}_{a}_{b}_{c}_{a}

_{3}_{b}^{2}_{c}^{2}_{a}

_{0}_{b}_{c}_{a}_{c}_{b}

_{1}_{b}_{c}_{a}

_{2}_{b}_{c}_{b}_{c}_{a}_{b}_{c}_{a}_{a}

[0327]

[0328] As further depicted in

[0329] The average peak width value provides an input to a motion status module

[0330] In the case of motion, motion artifacts are suppressed using the motion artifact suppression module

[0331] In the case of no motion, one of the signals (the infrared signal in the present embodiment) is subjected to DC removal and bandpass filtering as represented in the DC removal and bandpass filter module

[0332] In the present embodiment, the spectral estimation comprises a Chirp Z transform that provides a frequency spectrum of heart rate information. The Chirp Z transform is used rather than a conventional Fourier Transform because a frequency range for the desired output can be designated in a Chirp Z transform. Accordingly, in the present embodiment, a frequency spectrum of the heart rate is provided between 30 and 250 beats/minute. In the present embodiment, the frequency spectrum is provided to a spectrum analysis module

[0333] In the case of motion, a motion artifact suppression is completed on the snapshot with the motion artifact suppression module

[0334] The motion artifact reference processor _{a }

[0335] The motion artifact correlation canceler

[0336] Because only one saturation value is provided to the reference processor, only one output vector of 270 samples results at the output of the motion artifact suppression correlation canceler

[0337] As described above, an alternative joint process estimator uses the QRD least squares lattice approach (

[0338]

[0339] The initialization parameters are referenced in

[0340] Number of Cells=6

[0341] Lambda=0.8

[0342] MinSumErr=10^{−20 }

[0343] GamsInit=10^{−2 }

[0344] SumErrInit=10^{−6}

[0345] The clean waveform output from the motion artifact suppression module

[0346] The output of the spectrum analysis module

[0347] The output filter

[0348] Alternative to Saturation Transform Module—Bank of Filters

[0349] An alternative to the saturation transform of the saturation transform module _{λb}_{λa}

[0350] There are N filter elements in each filter bank. Each of the filter elements in the first filter bank _{1}_{1}_{2}_{2}_{3}_{3}_{4 }_{N−1}_{N }

[0351] It should be understood that the number of filter elements can range from 1 to infinity. However, in the present embodiment, there are approximately 120 separate filter elements with center frequencies spread evenly across a frequency range of 25 beats/minute-250 beats/minute.

[0352] The outputs of the filters contain information about the primary and secondary signals for the first and second measured signals (red and infrared in the present example) at the specified frequencies. The outputs for each pair of matching filters (one in the first filter bank

[0353] The ratio module

[0354] The output of the saturation equation modules

[0355] The results of the histogram provide a power curve similar to the power curve of

[0356] It should be understood that as an alternative to the histogram, the output saturation (not necessarily a peak in the histogram) corresponding to the highest saturation value could be selected as the arterial saturation with the corresponding ratio representing r_{a}_{v}

[0357] Alternative Determination of Coefficients r_{a }_{v }

[0358] As explained above, in accordance with the present invention, primary and secondary signal portions, particularly for pulse oximetry, can be modeled as follows:

_{red}_{1}_{1 }

_{IR}_{2}_{2 }

_{1}_{a}_{2 }_{1}_{v}_{2}

[0359] Substituting Equation (91) into Equation (89) provides the following:

_{red}_{a}_{2}_{v}_{2 (red)}

[0360] Note that S_{red }_{IR }_{red }_{IR }_{1 }_{2 }_{1 }_{2}

[0361] As explained above, determining r_{a }_{v }_{a }_{v }_{a }_{v }_{k }_{k}_{k }_{k }

[0362] where i represents time.

[0363] It should be understood that other correlation functions such as a normalized correlation could also be used.

[0364] Minimizing this quantity often provides a unique pair of r_{a }_{v }_{2 }_{2}_{a }_{v}_{2 }

[0365] inverting the two-by-two matrix provides:

[0366] Thus,

[0367] Preferably, the correlation of equation (93) is enhanced with a user specified window function as follows:

[0368] The Blackman Window is the presently preferred embodiment. It should be understood that there are many additional functions which minimize the correlation between signal and noise. The function above is simply one. Thus,

[0369] In order to implement the minimization on a plurality of discrete data points, the sum of the squares of the red sample points, the sum of the squares of the infrared sample points, and the sum of the product of the red times the infrared sample points are first calculated (including the window function, w_{i}

[0370] These values are used in the correlation equation (93b). Thus, the correlation equation becomes an equation in terms of two variables, r_{a }_{v}_{a }_{v}_{a }_{v }_{a }_{v }_{a }_{v}

[0371] Once r_{a }_{v }_{a }_{v }_{a }_{v}

[0372] In a further implementation to obtain r_{a }_{v}_{a }_{v }_{2 }_{2 }_{2}_{a}_{v }_{a }_{v }_{2 }_{2 }_{2 }_{2 }

[0373] where R_{1 }_{2 }_{1,2 }

[0374] The correlation between s_{2 }_{2 }

[0375] As explained above, the constraint is that s_{k }_{k }

_{1}_{a}_{v}_{12}_{a}_{v}_{2}

[0376] In other words, the goal is to maximize equation (94) under the constraint of equation (98).

[0377] In order to obtain the goal, a cost function is defined (e.g., a Lagrangian optimization in the present embodiment) as follows:

[0378] where μ is the Lagrange multiplier. Finding the value of r_{a}_{v }

[0379] Along the same lines, if we assume that the red and infrared signals S_{red }_{IR }_{1}_{2 }_{12 }_{1 }_{2}

_{1}_{1}_{a}_{v}_{12}_{1}_{a}_{v}_{2}_{1}

_{1}_{2}_{a}_{v}_{12}_{2}_{a}_{v}_{2}_{2}

[0380] Because equations (100) and (101) are non-linear in r_{a }_{v }_{a}_{v}_{a}_{v }

_{12}_{1}_{2}_{1}_{1}_{1}

_{12}_{2}_{2}_{2}_{1}_{2}

[0381] These equation (102) and (103) can be solved for x and y. Then, solving for r_{a }_{v }

[0382] Solving equation (104) results in two values for r_{v}_{v }^{2}_{v}_{v }^{2}_{v}_{v }_{2 }_{2}_{2 }_{v }_{a}_{a}_{v }

[0383] Alternative To Saturation Transform—Complex FFT

[0384] The blood oxygen saturation, pulse rate and a clean plethysmographic waveform of a patient can also be obtained using the signal model of the present invention using a complex FFT, as explained further with reference to FIGS. _{a }_{v }

[0385]

[0386] In this alternative embodiment, the snapshot for red and infrared signals is 562 samples from the decimation module

[0387] The high-pass filter modules

[0388] The window function modules

[0389] The complex FFT modules

[0390] In the first path of processing, the output from the select modules

[0391] The threshold modules

[0392] After thresholding, the data points are forwarded to a point-by-point ratio module

[0393] The phase difference module

[0394] 6. the red sample must pass the red threshold

[0395] 7. the infrared sample must pass the infrared threshold

[0396] 8. the phase between the two points must be less than the predefined threshold as determined in the phase threshold

[0397] For those sample points which qualify, a ratio is taken in the ratio module

[0398] The resulting ratios are provided to a saturation equation module which is the same as the saturation equation modules

[0399] The arterial (and the venous) saturation can then be selected, as represented in the select arterial saturation module

[0400] The fast saturation transform information can also be used to provide the pulse rate and the clean plethysmographic wave form as further illustrated in

[0401] As depicted in

[0402] The window function module performs a windowing function selected to pass those frequencies that significantly correlate to the frequencies which exhibited saturation values very close to the arterial saturation value. In the present embodiment, the following windowing function is selected:

[0403] where SAT_{n }_{art }

[0404] In order to obtain pulse rate, the output points from the window function module

[0405] In order to obtain a clean plethysmographic waveform, the output of the windowing function

[0406] Accordingly, by using a complex FFT and windowing functions, the noise can be suppressed from the plethysmographic waveform in order to obtain the arterial saturation, the pulse rate, and a clean plethysmographic waveform. It should be understood that although the above description relates to operations primarily in the frequency domain, operations that obtain similar results could also be accomplished in the time domain.

[0407] Relation to Generalized Equations

[0408] The measurements described for pulse oximetry above are now related back to the more generalized discussion above. The signals (logarithm converted) transmitted through the finger

_{λa}_{λred1}_{HbO2,λa}^{A}_{HbO2}^{A}_{Hb,λa}^{A}_{Hb}^{A}_{HbO2,λa}^{V}_{HbO2}^{V}_{Hb,λa}^{v}_{Hb}^{v}_{λa}

_{λa}_{HbO2,λa}^{A}_{HbO2}^{A}_{Hb,λa}^{A}_{Hb}^{A}_{λa}

_{λa}_{λa}_{λa}

_{λb}_{λred2}_{HbO2,λb}^{A}_{HbO2}^{A}_{Hb,λb}^{A}_{Hb}^{A}_{HbO2,λb}^{V}_{HbO2}^{V}_{Hb,λb}^{V}_{Hb}^{V}_{λb(}

_{λb}_{HbO2,λb}^{A}_{HbO2}^{A}_{Hb,λb}^{A}_{Hb}^{A}_{λb}

_{λb}_{λb}_{λb(}

[0409] The variables above are best understood as correlated to _{3 }_{4 }_{3 }_{4 }_{5 }_{6 }_{5 }_{6 }^{V}^{V}^{V }_{4}^{A}_{2 }^{A}^{A }_{5 }_{6}

[0410] The wavelengths chosen are typically one in the visible red range, i.e., λa, and one in the infrared range, i.e., λb. Typical wavelength values chosen are λa=660 nm and λb=910 nm. In accordance with the constant saturation method, it is assumed that c^{A}_{HbO2}^{A}_{Hb}_{1 }^{V}_{HbO2}^{V}_{Hb}_{2}_{λa}_{a}_{λb}

_{λa}_{a}_{λb}

_{λa}_{v}_{λb}

_{λa}_{v}_{λb}

[0411] In pulse oximetry, it is typically the case that both equations (108) and (109) can be satisfied simultaneously.

[0412] Multiplying equation (106) by r_{a}

_{λa}_{a}_{λb}

_{HbO2,λa}^{V}_{HbO2}^{V}_{Hb,λa}^{V}_{Hb}^{V}_{λa}_{a}_{HbO2,λb}^{V}_{HbO2}^{V}_{Hb,λb}^{V}_{Hb}^{V}_{λb}

[0413] Multiplying equation (106) by r_{v }

_{λa}_{v}_{λb}

_{λa}_{v}_{λb}

[0414] The constant saturation assumption does not cause the venous contribution to the absorption to be canceled along with the primary signal portions s_{λa}_{λb}

[0415] To illustrate the operation of the oximeter of _{λa}_{λred}_{λb}_{λIR}

[0416] _{λa}_{a}_{λb}_{λa }_{λb}

[0417] It should also be understood that a reference processor could be utilized in order to obtain the primary reference signal s′(t)=s_{λa}_{v}_{λb}

[0418] _{λa}_{λb}_{λa}_{λb}

[0419] It should be understood that approximation n″_{λa}_{λb}_{λa}_{λb}

[0420] Method for Estimating Primary and Secondary Signal Portions of Measured Signals in a Pulse Oximeter

[0421] Implementing the various embodiments of the correlation canceler described above in software is relatively straightforward given the equations set forth above, and the detailed description above. However, a copy of a computer program subroutine, written in the C programming language, which calculates a primary reference s′(t) using the constant saturation method and, using a joint process estimator

[0422] The correspondence of the program variables to the variables defined in equations (54)-(64) in the discussion of the joint process estimator is as follows:

[0423] Δ_{m}

[0424] Γ_{f,m}

[0425] Γ_{b,m}

[0426] f_{m}

[0427] b_{m}

[0428] ℑ_{m}

[0429] β_{m}

[0430] γ_{m}

[0431] ρ_{m,λa}

[0432] ρ_{m,λb}

[0433] e_{m,λa}

[0434] e_{m,λb}

[0435] κ_{m,λa}

[0436] κ_{m,λb}

[0437] A first portion of the program performs the initialization of the registers ^{−1}_{1 }_{2 }_{1}_{2}

[0438] A third portion of the subroutine calculates the primary reference or secondary reference, as in the “CALCULATE PRIMARY OR SECONDARY REFERENCE (s′(t) or n′(t)) FOR TWO MEASURED SIGNAL SAMPLES” action block _{λa}_{λb}_{a}_{v}_{a}_{v}_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0439] A fourth portion of the program performs Z-stage update as in the “ZERO STAGE UPDATE” action block _{o}_{o}_{o }_{0}

[0440] A fifth portion of the program is an iterative loop wherein the loop counter, M, is reset to zero with a maximum of m=NC_CELLS, as in the “m=0” action block _{m}

[0441] A sixth portion of the program calculates the forward and backward reflection coefficient Γ_{m,f}_{m,b}^{th}_{m}_{m}_{m}

[0442] A seventh portion of the program, still within the loop begun in the fifth portion of the program, calculates the regression coefficient register _{m,λa}_{m,λb}^{th }_{m,λa}_{m,λb}_{m,λa}_{m,λb}

[0443] The loop iterates until the test for convergence is passed. The test for convergence of the joint process estimator is performed each time the loop iterates analogously to the “DONE” action block _{m}

[0444] The output of the present subroutine is a good approximation to the primary signals s″_{λa}_{λb}_{λa}_{λb}_{λa}_{λb}

[0445] It should be understood that the subroutine of Appendix B is merely one embodiment which implements the equations (54)-(64). Although implementation of the normalized and QRD-LSL equations is also straightforward, a subroutine for the normalized equations is attached as Appendix C, and a subroutine for the QRD-LSL algorithm is attached as Appendix D.

[0446] While one embodiment of a physiological monitor incorporating a processor of the present invention for determining a reference signal for use in a correlation canceler, such as an adaptive noise canceler, to remove or derive primary and secondary components from a physiological measurement has been described in the form of a pulse oximeter, it will be obvious to one skilled in the art that other types of physiological monitors may also employ the above described techniques.

[0447] Furthermore, the signal processing techniques described in the present invention may be used to compute the arterial and venous blood oxygen saturations of a physiological system on a continuous or nearly continuous time basis. These calculations may be performed, regardless of whether or not the physiological system undergoes voluntary motion.

[0448] Furthermore, it will be understood that transformations of measured signals other than logarithmic conversion and determination of a proportionality factor which allows removal or derivation of the primary or secondary signal portions for determination of a reference signal are possible. Additionally, although the proportionality factor r has been described herein as a ratio of a portion of a first signal to a portion of a second signal, a similar proportionality constant determined as a ratio of a portion of a second signal to a portion of a first signal could equally well be utilized in the processor of the present invention. In the latter case, a secondary reference signal would generally resemble n′(t)=n_{λb}_{m,λa}

[0449] Furthermore, it will be understood that correlation cancellation techniques other than joint process estimation may be used together with the reference signals of the present invention. These may include but are not limited to least mean square algorithms, wavelet transforms, spectral estimation techniques, neural networks, Weiner and Kalman filters among others.

[0450] One skilled in the art will realize that many different types of physiological monitors may employ the teachings of the present invention. Other types of physiological monitors include, but are in not limited to, electro cardiographs, blood pressure monitors, blood constituent monitors (other than oxygen saturation) monitors, capnographs, heart rate monitors, respiration monitors, or depth of anesthesia monitors. Additionally, monitors which measure the pressure and quantity of a substance within the body such as a breathalizer, a drug monitor, a cholesterol monitor, a glucose monitor, a carbon dioxide monitor, a glucose monitor, or a carbon monoxide monitor may also employ the above described techniques.

[0451] Furthermore, one skilled in the art will realize that the above described techniques of primary or secondary signal removal or derivation from a composite signal including both primary and secondary components can also be performed on electrocardiography (ECG) signals which are derived from positions on the body which are close and highly correlated to each other. It should be understood that a tripolar Laplacian electrode sensor such as that depicted in

[0452] Furthermore, one skilled in the art will realize that the above described techniques can also be performed on signals made up of reflected energy, rather than transmitted energy. One skilled in the art will also realize that a primary or secondary portion of a measured signal of any type of energy, including but not limited to sound energy, X-ray energy, gamma ray energy, or light energy can be estimated by the techniques described above. Thus, one skilled in the art will realize that the techniques of the present invention can be applied in such monitors as those using ultrasound where a signal is transmitted through a portion of the body and reflected back from within the body back through this portion of the body. Additionally, monitors such as echo cardiographs may also utilize the techniques of the present invention since they too rely on transmission and reflection.

[0453] While the present invention has been described in terms of a physiological monitor, one skilled in the art will realize that the signal processing techniques of the present invention can be applied in many areas, including but not limited to the processing of a physiological signal. The present invention may be applied in any situation where a signal processor comprising a detector receives a first signal which includes a first primary signal portion and a first secondary signal portion and a second signal which includes a second primary signal portion and a second secondary signal portion. Thus, the signal processor of the present invention is readily applicable to numerous signal processing areas.