# How to remove heartbeat signal from blood pressure signal?

I am working on a medical engineering project that involves processing signals received from a human body. We use a sensor to record the blood pressure, $$B(t)$$.

However, if you put your finger on many parts of your skin, you can notice your heartbeat. Accordingly, the recorded signal also has some weak heartbeat signal added, $$A'(t)$$. To obtain accurate blood pressure signal, we have a sensor that records the heartbeat signal, $$A(t)$$. However, $$A'(t)$$ is a delayed and weak version of $$A(t)$$, where the phase lag is unknown.

Knowing $$A(t)$$, how can we find accurate $$B(t)$$?

• @rodrigo-de-azevedo , as mentioned we have recorded h(t) for a period of time and have some recorded signal like B'(t) in the same period, which has B(t) and h'(t). h'(t) has lower amplitudes than h(t) (as the sensor captures it weakly). However, it is not something like multiplying h(t) by some constant. Usually, h(t) has higher frequencies than B(t) so a low pass filter may be used. However, as h'(t) may change from person to person and B(t) may have high frequency components, I am looking for a better solution. for instance some sort of decomposing B'(t) and finding components related to.. Sep 12 at 14:28
• h'(t) is unknown and we wish to find it, If you mean B'(t), I dont know if cross correlation would help for finding delay. besides, while h' is much similar to h at each point t we have something like h'(t)=ah(t-v) for some noisy a. so finding only the delay would not solve the problem Sep 12 at 14:48
• Are you trying to record the blood pressure averaged over time, or are you trying to record systolic and diastolic blood pressure? If the latter, I'm not sure how you're going to separate these out, because the instantaneous blood pressure varies in sync with the heartbeat. Sep 12 at 14:55
• Just as a note -- you may want to review your variable names. In signal processing, $h(t)$ is almost universally used to denote a system impulse response, so unless $h(t)$ is already used to denote variations due to heartbeat in biomechanical circles, you'll probably confuse more than you clarify by that choice of signal name. I'd see if I can find what people tend to use in papers; if that doesn't work I'm not sure what to suggest -- $x_h(t)$, or $p_h(t)$, maybe. Sep 12 at 14:57
• @TimWescott , I am not only recording the blood pressure, It was just one sample. There are other types of signals which are mixed with heart beat. The choice of h(t) was poor and I am changing it. Sep 12 at 15:53

If $$A(t)$$ is known, it can be zeroed in the synchrosqueezed representation - the remainder is then $$B(t)$$, recovered by inversion. $$A(t)$$ need not be known perfectly - just enough to indentify its time-frequency ridges.