Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
HEART SOUND ANALYZER OF HEART SOUND INFORMATION WITH A PLURALITY OF HEART SOUNDS OF A PLURALITY OF FETUSES
Document Type and Number:
WIPO Patent Application WO/2005/074794
Kind Code:
A1
Abstract:
A heart sound analyzer component (102) of an apparatus (100) in one example extracts a heart sound (124) of a fetus from heart sound information (122) that comprises a plurality of mixtures of a plurality of heart sounds of a plurality of fetuses (304, 306).

Inventors:
FEDER LAURENCE M (US)
PRIEMER ROLAND (US)
NIGAM VIVEK
Application Number:
PCT/US2004/022738
Publication Date:
August 18, 2005
Filing Date:
July 16, 2004
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
FEDER LAURENCE M (US)
PRIEMER ROLAND (US)
NIGAM VIVEK
International Classes:
A61B7/02; A61B7/04; G06K9/00; (IPC1-7): A61B5/024; G06F17/00
Domestic Patent References:
WO2003003905A22003-01-16
Foreign References:
GB2162644A1986-02-05
Other References:
VARADY P; WILDT L; BENYO Z; HEIN A: "An advanced method in fetal phonocardiography", COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE, vol. 71, no. 3, July 2003 (2003-07-01), pages 283 - 296, XP008037413
MOGHAVVEMI M ET AL: "A non-invasive PC-based measurement of fetal phonocardiography", SENSORS AND ACTUATORS A, ELSEVIER SEQUOIA S.A., LAUSANNE, CH, vol. 107, no. 1, 1 October 2003 (2003-10-01), pages 96 - 103, XP004452856, ISSN: 0924-4247
KOMÁROMY B; LAMPÉ L: "Die Bedeutung und die Möglichkeiten der frühzeitigen Erkennung einer Zwillingsschwangerschaft", ZENTRALBLATT FÜR GYNÄKOLOGIE, vol. 91, no. 42, 18 October 1969 (1969-10-18), pages 1378 - 1387, XP008037057
Attorney, Agent or Firm:
Brill, Robert J. (SC 14135 N. Cedarburg Roa, Mequon WI, US)
Download PDF:
Claims:
CLAIMS What is claimed is:
1. An apparatus (100), comprising: a heart sound analyzer component (102) that extracts a heart sound (124) of a fetus from heart sound information (122) that comprises a plurality of mixtures of a plurality of heart sounds of a plurality of fetuses (304,306).
2. The apparatus (100) of claim 1, wherein the plurality of fetuses (304,306) comprises the fetus and one or more other fetuses, wherein the heart sound analyzer component (102) employs a blind source separation algorithm that accounts for propagation delays (318,320, 322, and 324) of the heart sound information (122) to extract the heart sound (124) of the fetus from the heart sound information (122) that comprises one or more other heart sounds of the one or more other fetuses.
3. The apparatus (100) of claim 1, wherein the heart sound analyzer component (102) employs a neural network construct (402) with a plurality of nodes in number at least equal to a number of the plurality of fetuses (304,306) to extract the heart sound (124) of the fetus from the heart sound information (122).
4. The apparatus (100) of claim 1, wherein the heart sound (124) of the fetus comprises a first heart sound from a first fetus, wherein the plurality of heart sounds of the plurality of fetuses (304,306) comprises the first heart sound and a second heart sound from a second fetus; wherein the heart sound information (122) comprises a mixture of the first heart sound, the second heart sound, and background noise; wherein the heart sound analyzer component (102) separates the first heart sound from the second heart sound and the background noise and creates a statistically independent output of the first heart sound.
5. The apparatus (100) of claim 1, further comprising: a plurality of sensor pods (308,310) of a stethoscope sensor head that are locatable over a womb (302) of a woman that carries the plurality of fetuses (304,306) ; wherein the plurality of sensor pods (308,310) comprise a plurality of microphones (104,106, 108,110) in number at least equal to a number of the plurality of fetuses (304, 306); wherein the plurality of microphones (104,106, 108,110) obtain and pass the heart sound information (122) to the heart sound analyzer component (102).
6. The apparatus (100) of claim 1, wherein the heart sound analyzer component (102) sets one or more propagation delay estimates for the heart sound information (122); wherein the heart sound analyzer component (102) iteratively processes the heart sound information (122) to converge the one or more propagation delay estimates closer to one or more true propagation delays of the heart sound information (122) to allow extraction of the heart sound (124) of the fetus from the heart sound information (122).
7. The apparatus (100) of claim 1, wherein the heart sound (124) of the fetus comprises a mixture of a plurality of discrete heart sounds from a plurality of corresponding distinct heart sound sources within a heart of the fetus; wherein the heart sound analyzer component (102) extracts one or more of the plurality of discrete heart sounds from the heart sound (124) of the fetus; wherein the plurality of discrete heart sounds comprise a first discrete heart sound of a first distinct sound source of the plurality of distinct heart sound sources ; wherein the plurality of discrete heart sounds comprise a second discrete heart sound of a second distinct sound source of the plurality of distinct heart sound sources; wherein the first discrete heart sound and the second discrete heart sound occur contemporaneously in the heart sound; wherein the heart sound analyzer component (102) creates statistically independent outputs of one or more of the first discrete heart sound and the second discrete heart sound.
8. The apparatus (100) of claim 1, wherein the heart sound analyzer component (102) collects one or more samples of the heart sound information (122) into a frame; wherein the heart sound analyzer component (102) calculates one or more updates for one or more propagation delay and weight estimates for the frame until a predetermined number of processing iterations have occurred or the one or more updates to the one or more propagation delay and weight estimates drop below a predetermined threshold; wherein the one or more samples comprise one or more first samples, wherein the frame comprises a first frame, wherein the one or more updates comprise one of more first updates; wherein after the predetermined number of processing iterations have occurred or the one or more first updates to the one or more propagation delay and weight estimates drop below the predetermined threshold, the heart sound analyzer component (102) collects one or more second samples of the heart sound information (122) into a second frame; wherein the heart sound analyzer component (102) calculates one or more second updates for the one or more propagation delay and weight estimates for the second frame until a predetermined number of iterations have occurred or the one or more second updates to the one or more propagation delay and weight estimates drop below a predetermined threshold; wherein the heart sound analyzer component (102) processes one or more additional frames until the one or more propagation delay and weight estimates substantially converge to constants that allow separation of the heart sound information (122) into a plurality of independent components.
9. A method, comprising the steps of : obtaining heart sound information (122) that comprises a plurality of mixtures of a plurality of heart sounds of a plurality of fetuses (304,306) ; and extracting a heart sound of a fetus from the heart sound information (122).
10. The method of claim 9, wherein the step of obtaining the heart sound information (122) that comprises the plurality of mixtures of the plurality of heart sounds of the plurality of fetuses (304,306) comprises the step of : moving a plurality of sensor pods (308,310) of a stethoscope sensor head over a womb (302) of a woman that carries the plurality of fetuses (304,306) to obtain the heart sound information (122); wherein the first sensor pod (308) comprises a first plurality of microphones (104, 106), wherein the second sensor pod (310) comprises a second plurality of microphones (108, 110), wherein the first plurality of microphones (104,106) combined with the second plurality of microphones (108,110) comprise a number of microphones in number at least equal to a number of the plurality of fetuses (304,306). * * *.
Description:
HEART SOUND ANALYZER OF HEART SOUND INFORMATION WITH A PLURALITY OF HEART SOUNDS OF A PLURALITY OF FETUSES CROSS-REFERENCE TO RELATED APPLICATIONS This application claims the priority of U. S. provisional Patent Application Serial No.

60/537, 941 (by Roland Priemer, filed January 21,2004, and entitled"EXTRACTION OF ONE OR MORE DISCRETE HEART SOUNDS FROM HEART SOUND INFORMATION").

This application contains subject matter that is related to the subject matter of the following applications, which are assigned to the same assignee as this application. The below-listed applications are hereby incorporated herein by reference in their entireties.

"EXTRACTION OF ONE OR MORE DISCRETE HEART SOUNDS FROM HEART SOUND INFORMATION, "by Roland Priemer, co-filed herewith.

"IDENTIFICATION OF ONE OR MORE DISTINCT HEART SOUND SOURCES THAT PRODUCE ONE OR MORE DISCRETE HEART SOUNDS,"by Roland Priemer, co-filed herewith.

"EXTRACTION OF HEART SOUND OF FETUS FROM HEART SOUND INFORMATION THAT COMPRISES PLURALITY OF MIXTURES OF PLURALITY OF HEART SOUNDS OF PLURALITY OF FETUSES,"by Roland Priemer, co-filed herewith.

TECHNICAL FIELD This invention relates generally to the medical arts and more particularly to sensing of heart sound information.

BACKGROUND A multiple gestation is a pregnancy in which a woman carries more than one fetus.

Multiple gestations are concerning because women who are expecting more than one baby are at increased risk of certain pregnancy complications, including preterm delivery. Some of the complications associated with multiple gestation can be minimized if they are detected and diagnosed early.

The embryonic heart starts beating twenty-two days after conception. The heart at this stage may be too small to hear. By the ninth or tenth week after the last menstrual period of the expecting mother, one is more likely to be able to hear the heartbeat of the fetus. By the twelve week, the heartbeat can usually be heard consistently by using a Doppler instrument for amplification.

Although the human ear is unequalled in detecting sounds over a particular frequency range, the spectrum occupied by fetal heart sounds is on the threshold of audibility. This aural processing has always caused problems relating to degree and accuracy of fetal auscultation. Electronically processing the fetal phonocardiogram ("FPCG") presents a way of avoiding these problems. The greater low frequency sensitivity of the phonocardiographic transducer captures all the acoustic information generated by the fetal heart. Then, additional visual examination of the FPCG provides an effective way of determining temporal measures of fetal cardiac function.

The FPCG allows the measurement of the instantaneous fetal heart rate, beat-to-beat differences and duration of systolic and diastolic phases. These measures are sensitive indicators of cardiac function, reflecting fetal well-being. With the record of various critical cardiac activities readily available, some abnormalities of the fetal heart can be detected in the early stages of the pregnancy. A problem arises with the detection of the abnormalities of the fetal heart in the event of a multiple gestation. When more than one fetus is present in a

mother's womb, the acoustic transducer (acoustic or electronic stethoscope) captures multiple FPCGs. As one shortcoming, it is difficult to clearly discern the heart sounds from one fetus in the mother's womb from the heart sounds of another fetus in the mother's womb.

Thus, a need exists for a capability to discern a heart sound of one fetus from the heart sound of another fetus.

SUMMARY The invention in one implementation encompasses an apparatus. The apparatus comprises a heart sound analyzer component that extracts a heart sound of a fetus from heart sound information that comprises a plurality of mixtures of a plurality of heart sounds of a plurality of fetuses.

Another implementation of the invention encompasses a method. Heart sound information is obtained that comprises a plurality of mixtures of a plurality of heart sounds of a plurality of fetuses. A heart sound of a fetus extracted from the heart sound information.

DESCRIPTION OF THE DRAWINGS Features of exemplary implementations of the invention will become apparent from the description, the claims, and the accompanying drawings in which: FIG. 1 is a representation of one exemplary implementation of an apparatus that comprises a heart sound analyzer, a sensor array of a plurality of sensors, a plurality of signal paths, and one or more speaker components.

FIG. 2 is a representation of the heart sound analyzer of the apparatus of FIG. 1, and illustrates one or more analog to digital converters, one or more data buffer components, one or more software memory components, one or more processor components, one or more digital to analog converters, one or more filters, one or more amplifiers, one or more storage devices, one or more input devices, and one or more output devices.

FIG. 3 is a representation of a womb of a woman that carries a plurality of fetuses that each contributes a heart sound to heart sound information captured by the plurality of sensors of the apparatus of FIG. 1.

FIG. 4 is a representation of an exemplary neural network that extracts two statistically independent outputs from heart sound information captured by the plurality of sensors of the apparatus of FIG. 1.

FIG. 5 is a representation of an exemplary interaction between a neural network and a recursive parameter adjustment algorithm of the heart sound analyzer of the apparatus of FIG.

1.

FIG. 6 is a representation of one exemplary plot of two fetal phonocardiograms that comprise heart sounds of two fetal hearts.

FIG. 7 is a representation of another exemplary plot of heart sound information representing the output of two sensors of the plurality of sensors of the apparatus of FIG. 1.

Two fetal heart sounds are delayed by propagation delays between the fetal hearts and the plurality of sensors and then mixed by each of the plurality of sensors.

FIG. 8 is a representation of another exemplary plot of two fetal phonocardiograms extracted from heart sound information by the heart sound analyzer of the apparatus of FIG.

1.

DETAILED DESCRIPTION Turning to FIG. 1, an apparatus 100 in one example comprises a heart sound analyzer 102, a sensor array 103 of a plurality of sensors 104,106, 108, and 110, a plurality of signal paths 112,114, 116, and 118, and one or more speaker components 120. The sensors 104, 106, 108, and 110 in one example obtain heart sound information 122 from a plurality of fetal hearts. The heart sound information 122 comprises a plurality of mixtures of a plurality of heart sounds of a plurality of fetuses. For example, the sensors 104,106, 108, and 110 each capture fetal phonocardiograms ("FPCGs") from the hearts of the plurality of fetuses within a womb of a woman. The heart sound information 122 comprises mixtures of heart sounds from the plurality of fetuses. The sensors 104,106, 108, and 110 pass the heart sound information 122 through the signal paths 112, 114, 116, and 118 to the heart sound analyzer 102. The heart sound analyzer 102 extracts a heart sound 124 of a fetus from the heart sound information 122 for aural examination by a doctor. For example, the heart sound 124 represents a heartbeat of a single fetal heart substantially without interference from heartbeats of the other fetal hearts within the womb. The heart sound analyzer 102 separates the heart sound 124 from one or more other sounds of the heart sound information 122. The heart sound analyzer 102 in one example comprises a fetal phonocardiogram analyzer.

The plurality of fetuses may comprise twins, triplets, quadruplets, quintuplets, or more. For example, the heart sound information 122 may comprise a mixture of two, three, four, five, or more fetal hearts. To extract the heart sound 124 associated with a single fetal heart, the apparatus 100 comprises a number of the plurality of sensors 104,106, 108, and 110 at least equal to a number of the plurality of fetuses. In one example, the apparatus 100 comprises a number of the plurality of sensors 104,106, 108, and 110 greater than a number of the plurality of fetuses. For example, the plurality of fetuses may comprise two fetuses and the apparatus 100 may comprise three sensors to allow separation of background noise,

or the like, from the output of the heart sound 124 of one of the plurality of fetuses. The sensors in one example comprise microphones. The microphones capture the heart sound information 122.

Referring to FIGS. 1-2, the heart sound analyzer 102 comprises one or more analog to digital converters 204, one or more data buffer components 206 and 212, one or more software memory components 208, one or more processor components 210, one or more digital to analog converters 214, one or more filters 216, one or more amplifiers 218, one or more storage devices 220, one or more input devices 222, and one or more output devices 224. Upon receipt of the heart sound information 122 from the sensor array 103, the heart sound analyzer 102 processes the heart sound information 122 and in one example outputs the heart sound 124.

The processor component 210 of the heart sound analyzer 102 extracts heartbeats of one or more of the plurality of fetuses from the heart sound information 122. For example, the heart sound 124 comprises a heartbeat of a first fetus of the plurality of fetuses within the womb. In one example, the heart sound information 122 comprises a plurality of mixtures of a first heartbeat of a first fetus and a second heartbeat of a second fetus. The heart sound analyzer 102 employs a blind source separation algorithm that accounts for propagation delays of the heart sound information 122 to create statistically independent outputs of one or more of the first heartbeat and the second heartbeat. The heart sound analyzer 102 in one example employs an independent component analysis procedure to extract the heart sound 124 from the heart sound information 122.

The analog to digital converter 204 in one example obtains the heart sound information 122 from the sensor array 103. The analog to digital converter 204 in one example digitizes the heart sound information 122. For example, the processor component 210 employs the analog to digital converter 204 to digitize the heart sound information 122.

The analog to digital converter 204 in one example outputs an MxK matrix of data of the heart sound information 122 to the data buffer component 206, where, for example, M is the number of sensors (e. g. , the sensors 104,106, 108, and 110) in the sensor array 103 and K is the number of times the analog to digital converter 204 samples the heart sound information 122. The data buffer component 206 in one example comprises MxK point double buffered data memory.

The data buffer component 206 in one example provides the MxK matrix of data to the processor component 210. The data buffer component 206 in one example comprises a demultiplexer, first and second data buffers, and a multiplexer. The processor component 210 in one example controls the demultiplexer to steer the output of the analog to digital converter 204 to the first data buffer with an MxK matrix of data. While the first data buffer is filling with data, the processor component 210 in one example controls the multiplexer to steer data from the second data buffer to the processor component 210 until the processor component 210 has received an MxK matrix of data. The processor component 210 processes the MxK matrix of data with software stored in the software memory component 208.

After the processor 210 has completed processing the MxK matrix of data and the first data buffer has filled with data, the processor component 210 in one example controls the demultiplexer to steer data to the second data buffer with an MxK matrix of data. While the second data buffer is filling with data, the processor component 210 in one example controls the multiplexer to steer data from the first data buffer to the processor component 210 until the processor component 210 has received an MxK matrix of data. The processor component 210 again employs the software stored in the software memory component 208 to process the MxK matrix of data. By continuously reversing the roles of the first and second data buffers, the heart sound analyzer 102 can work in real time.

The software memory component 208 in one example stores software for use by the processor component 210. The processor component 210 employs the software to extract the heart sound 124 from the heart sound information 122. The software memory component 208 in one example comprises an instance of the recordable data storage medium 146. The software memory component 208 in one example stores a blind source separation algorithm that accounts for propagation delays within the heart sound information 122. The processor component 210 employs the blind source separation algorithm to extract the heart sound 124 from the heart sound information 122.

After extraction of the heart sound 124 from the heart sound information 122, the processor component 210 stores the heart sound 124. For example, the processor component 210 may store the heart sound 124 in the storage device 220. The processor component 210 may then access the heart sound 124 upon request from a user of the heart sound analyzer 102. For example, the user may employ the input device 222 to cause the processor 210 to output the heart sound 124 to the output device 224 and/or the speaker component 120. The input device 222 in one example comprise one or more of a button, a dial, a mouse, a keyboard, and/or a touch-screen. The output device 224 in one example comprise a liquid crystal display ("LCD").

The user in one example chooses with the input device 222 to listen to the heart sound 124 with the speaker component 120. In another example, the user chooses with the input device 222 to output the heart sound 124 to the output device 224. The user may also choose to filter, amplify, and/or shift a spectral content of the heart sound 124 from a first frequency range to a second frequency range with the processor component 210. In one example, the processor component 210 employs digital filtering software to filter out one or more frequency ranges of the heart sound 124. In another example, the processor component 210 amplifies one or more regions of the heart sound 124. In still another example, the processor

component 210 may shift the spectral content of the heart sound 124 from a first frequency range to a second frequency range. For example, the user may hear better in the second frequency range than in the first frequency range.

The data buffer component 212 obtains data from the processor component 210. In one example, the data buffer component 212 obtains the heart sound 124 from the processor component 210. In another example, the data buffer component 212 obtains the heart sound 124 which has been spectrally modified, spectrally shifted, and/or amplified. The data buffer component 212 in one example comprises K point double buffered data memory.

The digital to analog converter 214 in one example converts a digital representation of the heart sound 124 into an analog representation of the heart sound 124. The digital to analog converter 214 in one example outputs the heart sound 124 to one or more of the speaker component 120, the filter 216, and/or the amplifier 218. The filter 216 in one example filters one or more frequency ranges of the analog representation of the heart sound 124. The filter 216 in one example comprises one or more low pass filters. The amplifier 218 in one example amplifies an output of the filter 216 to drive the speaker component 120.

The user may employ one or more of the speaker component 120 and/or the output device 224 to listen to and/or view the heart sound 124 and/or a representation of the heart sound 124. The user may employ the speaker component 120 and/or the output device 224 for examination of the heart sound 124. The user may employ the heart sound analyzer 102 to automatically diagnose one or more fetal heart dysfunctions. For example, the processor component 210 compares one or more features (e. g. , signatures) of the heart sound 124 to a normal range of the features (e. g. , signatures) of the heart sound 124. If the heart sound 124 is outside of the normal range, then the heart sound 124 may indicate one or more dysfunctions.

The heart sound 124 in one example is a mixture of a plurality of discrete heart sounds from a plurality of corresponding distinct heart sound sources (e. g. , a heart valve) of a single fetal heart. The heart sound analyzer 102 may extract from the heart sound 124 one or more discrete heart sounds of the corresponding distinct heart sound sources, as described in "EXTRACTION OF ONE OR MORE DISCRETE HEART SOUNDS FROM HEART SOUND INFORMATION, "by Priemer, co-filed herewith. For example, the heart sound 124 may comprise a first discrete heart sound from a first distinct sound source of the single fetal heart. The heart sound 124 may also comprise a second discrete heart sound from a second distinct sound source of the single fetal heart. The first and second discrete heart sounds in one example occur contemporaneously in the heart sound 124. The heart sound analyzer 102 in one example creates statistically independent outputs of one or more of the first discrete heart sound and the second discrete heart sound to allow aural examination of the individual first and second discrete heart sounds by a doctor.

Referring to FIGS. 1 and 3, a womb 302 of a woman may carry a plurality of fetuses 304 and 306. Each of the fetuses 304 and 306 contribute heart sounds to the heart sound information 122 captured by the sensors 104,106, 108, and 110. The sensors 104,106, 108, and 110 in one example are organized into a plurality of sensor pods 308 and 310. In one example, each of the sensor pods 308 and 310 comprise a single sensor or microphone. In another example, the sensor pod 308 comprises a plurality of sensors (e. g. , the sensors 104 and 106) and the sensor pod 310 comprises a plurality of sensors (e. g. , the sensors 108 and 110).

The sensor pods 308 and 310 are independently movable over an abdomen 312 of the woman to capture the heart sound information 122 of the fetuses 304 and 306 within the womb 302. For example, a doctor may move the sensor pod 308 over the abdomen 312 to a

location that is near the fetus 304 and the doctor may move the sensor pod 310 over the abdomen 312 to a location that is near the fetus 306.

The heart sound information 122 comprises a mixture of first heart sounds from a heart 314 of the fetus 304 and second heart sounds from a heart 316 of the fetus 306. The first heart sounds from the fetus 304 experience a propagation delay (d22) 318 traveling to the sensor pod 308 and a propagation delay (dl2) 320 traveling to the sensor pod 310. The second heart sounds from the fetus 306 experience a propagation delay (d21) 322 traveling to the sensor pod 308 and a propagation delay (dll) 324 traveling to the sensor pod 310. To extract the heart sound 124 from the heart sound information 122, the blind source separation algorithm accounts for the propagation delays (d22, dl2, d2i, and dll) 318,320, 322, and 324 of the first and second heart sounds from the hearts 314 and 316.

Referring to FIGS. 1-4, a neural network 402 illustrates an exemplary extraction of two statistically independent outputs from the heart sound information 122 captured by two sensors (e. g. , the sensors 104 and 106). The processor component 210 passes the heart sound information 122 through a plurality of external nodes of the neural network 402 in number at least equal to the number of the plurality of fetuses 304 and 306. At each external node of the plurality of external nodes, the processor component 210 passes the heart sound information 122 through a plurality of internal nodes of the neural network 402 equal to the number of external nodes to determine one or more weights and one or more delays of each of the heart sounds from the plurality of fetuses 304 and 306. For example, the neural network 406 outputs an MxK matrix of data that represents the heart sounds of the plurality of fetuses 304 and 306. The heart sounds of the plurality of fetuses 304 and 306 in one example are statistically independent.

The heart sound information 122 comprises a plurality of composite heart sounds.

For example, the plurality of composite heart sounds comprise a composite heart sound signal

(xi) 404 captured by the sensor 104 and a composite heart sound signal (x2) 406 captured by the sensor 104. The outputs of the sensors 104 and 106 are sampled at discrete time points t.

The discrete time can be written as t = nT, where n is an integer and T is the time increment between samples. The sampling rate is 1/T samples/sec, and n is called the discrete time index. The composite heart sound signals (xi and x2) 404 and 406 in one example comprise delayed mixtures of fetal phonocardiogram (si) produced by the fetus 306 and fetal phonocardiogram (s2) produced by the fetus 304. The composite heart sound signals (xl and x2) 404 and 406 in one example are modeled by: x1(t) = a11s1(t-d11)+a12s2(t-d12) x2(t) = a21s1(t-d21)+a22s2(t-d22) where d22, dl2, d2l, and du represent the propagation delays 318, 320, 322, and 324 of the first and second heart sounds. Also, where all, al2, a2l, and a22 represent attenuations of the first and second heart sounds due to signal strength losses traveling through the abdomen 312 of the woman to the sensors (e. g. , the sensors 104 and 106).

The neural network 402 serves to separate the composite heart sound signals (xi and x2) 404 and 406 to obtain estimates (ul and u2) 408 and 410 of the fetal phonocardiograms si and s2. The neural network computes: <BR> <BR> <BR> <BR> <BR> <BR> u1(t) = w11x1(t)+w12u2(t-(b12))+w01<BR> u2(t) = w22x2(t)+w21u1(t-(b21))+w02 where bl2 and b2l represent delay differences 414 and 412 between the propagation delays 318,320, 322, and 324. Also, where wol, w11, w12, w02, w22 and w21 represent weights 416, 420,424, 418,422, and 426 of the neural network 402.

The blocks 428 and 430 that convert the estimates (ul and u2) 408 and 410 to signals (yl and y2) 432 and 434 represent nonlinearities described by the equations below. A logistic

function relates the estimate (ul) 408 to the signal (y1) 432 and the estimate (u2) 410 to the signal (y2) 434. For example a monotonically increasing nonlinear function converts the estimate (ul) 408 to the signal (yi) 432 and the estimate (u2) 410 to the signal (y2) 434.

Below are representative logistic functions that relate the estimate (ul) 408 to the signal (y1) 432 and the estimate (u2) 410 to the signal (y2) 434: Yl (t) =g(u1(t))=1/(1+exp(-u1(t))) Y2 (t) =g(u2(t))=1/(1+exp(-u2(t))) Assuming that the fetal phonocardiograms (sl and s2) are statistically independent, while composite heart sound signals (xi and x2) 404 and 406 cannot be statistically independent, the parameters (e. g. , the weights (wolf wn, w12, wo2, W2z and w21) 416,420, 424, 418,422, and 426, and the delay differences (bl2 and b21) 414 and 412) of the neural network 402 are calculated to make the signals (y1 and y2) 432 and 434 statistically independent.

Then, the estimates (ul and u2) 408 and 410 must also be statistically independent. If the conditions below are achieved: <BR> <BR> <BR> <BR> <BR> b12 # d12 - d22, w12w22a22 # -w11a12<BR> <BR> <BR> <BR> <BR> <BR> <BR> <BR> b21 d21 - d11, w11w21a11 # -w22a21 then the estimate (ul) 408 becomes proportional to the fetal phonocardiogram (s1) and the estimate (u2) 410 becomes proportional to the fetal phonocardiogram (s2), delayed by the propagation delays (dll and d22) 324 and 318, respectively, which means that the estimates (ul and u2) 408 and 410 are statistically independent, resulting in: <BR> <BR> <BR> u1(t) # a11w11s1(t-d11)#s1(t-d11)<BR> <BR> <BR> <BR> <BR> <BR> <BR> <BR> <BR> u2(t) # a22w22s2(t-d22)#s2(t-d22)

Determining the weights (w01, w11, w12, won, w2 and w21) 416,420, 424,418, 422, and 426, and the delay differences (bl2 and b21) 412 and 414 that will make the signals (yl and y2) 432 and 434 statistically independent is the goal of the phonocardiogram separator depicted in FIG. 5. The recursive parameter adjustment algorithm 502 to estimate weights and delays of the neural network 402 is based on the gradients given by: Awoi oc (1-2yi) #wii#(1-2yi)xi+1/wii for i = 1, 2 <BR> <BR> <BR> <BR> <BR> <BR> <BR> #wij#(1-2yi)uj(t-bij)<BR> <BR> <BR> <BR> <BR> <BR> <BR> <BR> Abij #-(1-2yi)wij#j(t-bij) uj (t) = uj (t)-uj (t-1) for i, j = 1, 2, iNj The formulas above give the gradients of entropy of the signals H (y1, y2) with respect to each of the parameters of the neural network 402. For example, consider the determination of a parameter of the neural network 402, such as the weight (w11) 420. The gradient of H with respect to the weight (wl l) 420 is denoted by A wl l.

For each input sample xi (t) and x2 (t), y1 (t) and y2 (t) and ul (t) and u2 (t) can be calculated and then #w11 (t) can also be calculated. Awl, can be calculated for many (e. g., K) input samples, and then the resulting average of #w11 can be calculated. Then, the average of #w11, denoted by #w11 (K), is used to update the value of the weight (wll) 420 with: w11 (new) = w, 1 (old) + h Awol (K) where h, a small number, is called the learning rate, that is selected by trial and error. Each of the other weights (w01, wn, w12, wo2, w22 and w21) 416,420, 424,418, 422, and 426, and the delay differences (bl2 and b21) 414 and 412 are updated analogously.

Before moving on to the next set of K input samples, all the updated weights and delay differences are used to calculate #w11 (K) and the other gradients again, and then again the above equation may be used to update wll and all the other weights (won, wn, wl2, w02, W22 and w2l) 416,420, 424,418, 422, and 426, and delay differences (bl2 and b21) 414 and <BR> <BR> 412. This loop is continued until the absolute value of lw (new) -w (old)# for all of the weights (poli wn, w12, w02, w22 and w21) 416,420, 424,418, 422, and 426 and b (new) -b (old) l for all of the delay differences (bl2 and b21) 414 and 412 becomes smaller than some predetermined threshold or the loop has been executed more than a predetermined number of times.

The blind source separation algorithm that accounts for the propagation delays of the heart sound information 122 comprises a number of steps. The first step is for the heart sound analyzer 102 to collect a fixed number (e. g. , K) of incoming mixture samples in a single frame. The second step is for the heart sound analyzer 102 to calculate the updates to be made to the weights (w01, wll, wl2, w02, W22 and w21) 416,420, 424,418, 422, and 426, and propagation delay differences (bl2 and b21) 414 and 412 by iterating over the frame once.

This is equivalent to averaging K gradients and updating the weights (w01, w11, w12, w02, w22 and w21) 416,420, 424,418, 422, and 426, and propagation delay differences (bl2 and b2i) 414 and 412. The third step is for the heart sound analyzer 102 to repeat the first and second steps until a predetermined number of iterations over a frame have been completed, or the updates to weights (w01, w11, w12, w02, w22 and w21) 416,420, 424,418, 422, and 426, and propagation delay differences (bl2 and b21) 414 and 412 have become smaller than a predetermined threshold. The fourth step is for the heart sound analyzer 102 to input the next frame, which comprises new incoming data, and discard the last frame. By overlapping successive frames, convergence of the weights (w01, wll, wl2, w02, W22 and w21) 416,420, 424,418, 422, and 426, and propagation delay differences (bl2 and b21) 414 and 412 may be smoother. The fifth step is for the heart sound analyzer 102 to repeat the first, second, third,

and fourth steps until the weights (won, wl, w12, W02, W22 and w21) 416,420, 424,418, 422, and 426, and propagation delay differences (bl2 and b21) 414 and 412 have converged to constants.

Referring to FIG. 5, the neural network 402 receives the composite heart sound signals (xl and x2) 404 and 406 as inputs. A recursive parameter adjustment algorithm 502 also receives the heart sound signals (xi and x2) 404 and 406 as inputs. The neural network 402 produces the estimates (ul and u2) 408 and 410 and the signals (yl and y2) 432 and 434.

The recursive parameter adjustment algorithm 502 receives the estimates (ul and u2) 408 and 410 and the signals (yl and y2) 432 and 434 as inputs. The recursive parameter adjustment algorithm 502 serves to determine the parameters (e. g. , the weights and delay differences) of the neural network 402.

Referring to FIGS. 1 and 6-8, fetal phonocardiograms 602 and 604 illustrate heart sounds of two fetal hearts. FIG. 7 depicts heart sound information 702 and 704 comprised of mixtures of heart sounds of the two fetal hearts delayed by propagation delays between the fetal hearts and the sensors. Fetal phonocardiograms 802 and 804 illustrate a separation of the heart sounds of the two fetal hearts from the mixed phonocardiograms 702 and 704. For example, the heart sound information 122 comprises the two phonocardiogram mixtures 702 and 704. The heart sound analyzer 102 extracts the heart sound of each of the two fetuses into statistically independent outputs.

The apparatus 100 in one example comprises a plurality of components such as one or more of electronic components, hardware components, and computer software components.

A number of such components can be combined or divided in the apparatus 100. An exemplary component of the apparatus 100 employs and/or comprises a set and/or series of computer instructions written in or implemented with any of a number of programming languages, as will be appreciated by those skilled in the art.

The apparatus 100 in one example employs one or more computer-readable signal-bearing media. Examples of a computer-readable signal-bearing medium for the apparatus 100 comprise the recordable data storage medium 146 of the heart sound analyzer 102. For example, the computer-readable signal-bearing medium for the apparatus 100 comprises one or more of a magnetic, electrical, optical, biological, and atomic data storage medium. In one example, the computer-readable signal-bearing medium comprises a modulated carrier signal transmitted over a network comprising or coupled with the apparatus 100, for instance, one or more of a telephone network, a local area network ("LAN"), the internet, and a wireless network.

The steps or operations described herein are just exemplary. There may be many variations to these steps or operations without departing from the spirit of the invention. For instance, the steps may be performed in a differing order, or steps may be added, deleted, or modified.

Although exemplary implementations of the invention have been depicted and described in detail herein, it will be apparent to those skilled in the relevant art that various modifications, additions, substitutions, and the like can be made without departing from the spirit of the invention and these are therefore considered to be within the scope of the invention as defined in the following claims.