MINIMIZING COMMON MODE INTERFERENCE IN A PHYSIOLOGICAL MEASUREMENT DEVICE

Abstract

The present invention relates to a physiological measurement device (11). In order to reduce common mode interference and/or to improve characteristics of the digital input filters (11) of the measurement device (11) without compromising common mode rejection capabilities of the input filters, it is proposed to calculate by means of optimization a set of filter coefficients for at least one digital input filter (31) associated with a specific input channel of the measurement device (11) based on samples (c.sub.i) of multiple input signals (s.sub.i) of the measurement device (11) corresponding to a test signal (TS) applied to multiple input channels of the measurement device (11) and on at least one definition vector (v.sub.j) describing a linear combination of samples (c.sub.i) of at least two input channels.

Claims

1. A method for minimizing common mode interference in a physiological measurement device, the method comprising receiving samples (c.sub.i) of multiple input signals (s.sub.i) of the measurement device corresponding to a test signal (TS) applied to multiple input channels of the measurement device; inputting at least one definition vector (v.sub.i) describing a linear combination of samples (c.sub.i) of at least two input channels; optimizing a metric calculated from at least one vector signal (y.sub.i), the vector signal being based on the samples (c.sub.i) of at least two of the received input signals and the definition vector (v.sub.i); and based on the optimizing, obtaining at least one set filter coefficients (x.sub.i) for at least one digital input filter associated with a specific input channel.

2. The method of claim 1, wherein multiple definition vectors (v.sub.i) are inputted and at least one set of filter coefficients (x.sub.i) is obtained for a specific input vector (v.sub.i).

3. The method of claim 1, wherein optimizing the metric comprises multiple optimizing steps, a set of filter coefficient obtained by one step being kept as constant for a subsequent optimizing step.

4. The method of claim 1, wherein the optimizing is subject to at least one linear equality constraint with respect to a filter coefficient (x.sub.i).

5. The method of claim 4, wherein the method comprises inputting a DC gain value of at least one input channel and determining a corresponding linear equality constraint based on the DC gain value (DC.sub.i).

6. The method of claim 1, wherein the optimizing is subject to at least one nonlinear inequality constraint for limiting a gain of a certain input filter at a given frequency.

7. The method of claim 1, wherein the optimizing is subject to at least one bounds constraint with respect to filter coefficients (x.sub.i).

8. The method of claim 1, wherein the samples (c.sub.i) are received with a calibration sampling rate (f.sub.samp) that is increased compared to an operating sampling rate applied to perform physiological measurements by the measurement device and/or wherein a calibration resolution (res) of the received samples (c.sub.i) differs from a measurement resolution of samples received to perform the physiological measurements.

9. The method according to claim 1, wherein optimizing the metric comprises minimizing a penalty function of output samples (y.sub.i), the output samples corresponding to at least one vector signal (y.sub.i) calculated from the input signal filtered according to the filter coefficients (x.sub.i) and from the definition vector (v.sub.i).

10. The method of claim 9, wherein optimizing comprises minimizing multiple different penalty functions of the output samples (y.sub.i) simultaneously or minimizing a scalar target function of different penalty functions.

11. The method of claim 9, wherein at least one penalty function is a convex or quasi-convex penalty function.

12. A computer program comprising program code means for causing a computer to carry out the steps of the method as claimed in claim 1 when said computer program is carried out on the computer.

13. A calibration device for minimizing common mode interference in a physiological measurement device, the calibration device comprising a processor configured to: receive samples (c.sub.i) of multiple input signals (s.sub.i) of the measurement device corresponding to a test signal (TS) applied to multiple input channels of the measurement device; input at least one definition vector (v.sub.i) describing a linear combination of samples (c.sub.i) of at least two input channels; optimize a metric calculated from at least one vector signal (y.sub.i), the vector signal being based on the samples (c.sub.i) of at least two of the received input signals and the definition vector (v.sub.i); and based on the optimizing, obtain at least one set filter coefficients (x.sub.i) for at least one digital input filter associated with a specific input channel.

14. A physiological measurement device comprising: multiple input channels for processing multiple input signals (s.sub.i); a memory device configured to store multiple sets of filter coefficients (x.sub.i), at least one of said multiple sets being related to at least one digital input filter of the measurement device, the digital input filter being configured to filter samples (c.sub.i) of one input signal (s.sub.i) to obtain filtered samples of this input signal, and a processor configured to calculate multiple vector signals (y.sub.i), one of said multiple vector signals being calculated from multiple filtered input signals based on a definition vector (v.sub.i); wherein the memory device is configured to store different sets of filter coefficients (x.sub.i) related to one input channel and the processor is configured to select one of the different sets of filter coefficients (x.sub.i) for calculating a specific vector signal (y.sub.i).

15. (canceled)

16. The method of claim 11, wherein the penalty function is a norm.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0024] These and other aspects of the invention will be apparent from and elucidated with reference to the embodiment(s) described hereinafter. In the following drawings:

[0025] FIG. 1 shows a block diagram of a physiological measurement device and a calibration device;

[0026] FIG. 2 shows an exemplary measurement device;

[0027] FIG. 3 shows a signal flow within the measurement device according to an embodiment;

[0028] FIG. 4 shows a flow chart of a calibration method; and

[0029] FIG. 5 shows a further signal flow within the measurement device according to an embodiment.

DETAILED DESCRIPTION OF EMBODIMENTS

[0030] FIG. 1 shows a physiological measurement device 11 connected to a calibration device 13. The physiological measurement device 11 has multiple channels. For the sake of simplicity, only three channels are shown in FIG. 1 but in general N channels may be present. Each channel has an analog input 15 for inputting an analog input signal s.sub.i i=1, . . . , N. The inputs 15 of the individual channels are connected to an analog front end 17. The analog front end 17 includes analog input circuitry 19 associated with the individual channels. The analog input circuitry 19 may comprise one or more amplifiers, an impedance converter, an analog filter, or the like. The input circuitry of the individual channels is connected to a multichannel analog to digital converter 21. The multichannel analog to digital converter 21 is configured to convert the individual input signals s.sub.i preprocessed by the individual analog input circuitries 19 into samples c.sub.i i=1, . . . , N which represent the input signals s.sub.i of the individual channels. In the shown embodiment, the multichannel analog to digital converter 21 comprises multiple analog to digital converters (ADC 23) where each of these ADCs is associated with a specific channel. In another embodiment, a single ADC 23 may be associated with multiple channels and the multichannel analog to digital converter 21 may comprise a multiplexer to selectively connect one channel to the respective ADC. When using a multiplexer, one ADC 33 may be used to convert the input signal of every channel one after the other. Using the multiplexer, several input channels can be sampled sequentially with one ADC. The here described calibration method may compensate for a timing skew that may present in such a multiplexed setup. In a possible implementation, e.g. an ECG device, two ADCs may be present operable to sample e.g. five ECG input channels sequentially per ADC for a total of ten input channels.

[0031] A digital interface of the multichannel analog to digital converter 21 is connected to a processing unit 25 of the measurement device 11 so that the samples c.sub.i of the individual input signals s.sub.i are forwarded to the processing unit 25.

[0032] The processing unit 25 may comprise a first processor 27 and a first memory device 29 that constitute a computing device operable to perform various processing steps needed to perform measurements by the physiological measurement device 11. The processing unit 25 comprises multiple digital input filters 31. The number of input filters may correspond to the number N of channels and one digital input filter 31 may be assigned to one specific channel. The individual input filter 31 may be implemented as Finite Impulse Response (FIR) filters having filter coefficients represented by column vectors x.sub.i i 1, . . . , N. A length K of the individual vectors x.sub.i corresponds to the number of filter coefficients of each digital filter 31. The digital filters 31 are also referred to as Input Channel Specific digital Filters (ICSFs) because they are applied to the samples c.sub.i of the input signals s.sub.i before differential output signals are derived from these samples.

[0033] Note that the present disclosure is not limited to a specific filter topology. Instead of FIR filters other filters like Infinite Impulse Response (IIR) filters may be applied. The filters 31 may be implemented in software, i.e. a computer program stored in the memory device 29 may be programmed so that the first processor 27 executes a filter algorithm based on the filter coefficients x.sub.i. The filter coefficients x.sub.i may be stored in a programmable and non-volatile section of the first memory device 29. The processing unit 25 may be accessible by the calibration device 13 so that the calibration device 13 can read the samples x.sub.i and program the filter coefficients x.sub.i.

[0034] The processing unit 25 is operable for calculating one or more output signals y.sub.j j=1, . . . , M from the samples c.sub.i of the individual input signals filtered by the individual digital input filters 31. The output signals y.sub.j, also referred to as vector signals or just to as vectors, may be calculated by forming a linear combination of filtered samples {tilde over (c)}.sub.l of multiple channels, thereby performing a differential measurement with respect to at least two input signals s.sub.i. This linear combination may be specified by a definition vector v.sub.j of a specific output signal y.sub.j, where 1.sup.T*v.sub.j=0 and 1 stands for a column vector of ones, i.e. 1.sup.T=[1, . . . , 1] and * denotes the scalar product. The output signal is y.sub.j=y.sub.j={tilde over (c)}.sub.l*v.sub.j. The linear combination described by the definition vector specifies a differential measurement. The differential measurement may be based on two or more channels. In some applications like ECG or EEG, the vector signals are often also referred to as leads.

[0035] In physiology, voltage measurements are almost always differential. Other physical quantities like pressure or temperature are usually measured as absolute values (e.g. a temperature of 37? C. or an intra-arterial pressure of 100 mmHg), but in some cases, differential measurements can be of interest (cardiac output can be measured by thermodilution, which requires measuring the temperature difference at two points in a blood vessel). The here described technology can thus be applied to differential measurements that are not voltage measurements, too.

[0036] For the sake of simplicity, FIG. 1 shows only one digital filter 31 per input channel. However, it is possible to provide multiple digital filters 31 for one input channel i, each of them processing the samples c.sub.i. The different digital filters 31 may be associated with different output signals. Accordingly, a specific variant of the filtered signal {tilde over (c)}.sub.l.sup.j may be used for calculating a certain output signal y.sub.j.

[0037] The measurement device 11 may comprise an output device 33 coupled with the processing unit 25 so that the output device 33 can output the output signals y.sub.j. For example, the output device 33 may include elements for visualizing the output signals y.sub.j, for example a display or a printer. One or more post-processing steps may be present for processing at least one output signal y.sub.i before outputting it by the output device 33.

[0038] The physiological measurement device 11 may be an electrophysiological measurement device such as an ECG or an EEG. Various configurations, in particular regarding the number N of channels and the number M of output signals are possible.

[0039] A comparatively simple example is an ECG or EEG device having a single output signal (single vector) and two electrodes connected to two inputs 15 (two channels). The only definition vector needed for such a measurement device 11 may be v.sub.1=[1,?1].

[0040] According to another example shown in FIG. 2, a 12-lead ECG is provided having nine measurement electrodes connected to N=9 inputs 15 of the measurement device 11. Such measurement device thus has nine channels. The electrodes may include three limb electrodes (right arm RA, left arm LA, left leg LL), six chest electrodes (V1-V6), which are connected via a measurement cable 35 to nine different inputs 15 of the physiological measurement device 11.

[0041] The output vector signals y.sub.1, . . . , y.sub.2 may be defined according to twelve leads visualized by means of the output device 33 to medical personnel. The leads and respective vector signals include lead I (LA-RA), lead 11 (LL-RA), lead III (LL-LA), aVR (RA-0,5*(LA+LL)), aVL (LA-0,5*(RA+LL)), aVF (LL-(RA+LA)), and (Vn-WCT); n=1, . . . , 6. WCT corresponds to the Wilson's central terminal voltage which is approximated as (RA+LA+LL)/3.

[0042] The present disclosure is not limited to the above two exemplary measurement devices. For example a 15-lead or an 18-lead ECG can be provided based on the present disclosure as well. It is also possible to provide a high-density EEG based on the present disclosure having 32 or more channels.

[0043] FIG. 3 shows an example where an output signal y.sub.i is generated by subtracting a filtered input signal s.sub.2 from another filtered input signal s.sub.i (definition vector [1, ?1]). Accordingly, the output signal y.sub.1 corresponds to a signal obtained by differential measurement based on the two input signals s.sub.1 and s.sub.2. When performing physiological measurements, the input signals s.sub.1, s.sub.2 are subject to common mode interference. In an ideal situation, subtracting the two signals s.sub.1, s.sub.2 from each other would eliminate the common mode part entirely. However, in practical implementation of the measurement device 11, the signal y.sub.1 includes a common mode part to some extent due to differences of electrical characteristics e.g. of components 19, 21, 23 of a signal path related to the individual channels.

[0044] In electrophysiology voltages are measured that occur between points on the surface of a patient due to activity of muscle cells or nerve cells. The signals of interest are usually in the microvolt (EEG) to millivolt (ECG) range, while the common-mode component of the signal can be tens or hundreds of millivolts. In a medical context, there are many possible sources of interference (including common mode interference) ranging from other pieces of medical equipment, electronic communication or networking devices and power line noise. In order for the output signal y.sub.j, to be useful for diagnosis and therapy, electrophysiological measurement devices 11 are required to have a high common mode rejection ratio (CMRR). Although there are standardized test procedures for determining the CMRR in the range of power line frequency of 50 Hz or 60 Hz, common-mode interference should also be mitigated at other frequencies. The frequencies where common-mode interference should be mitigated may include frequencies outside the ECG/EEG bands, especially if the recording device performs some kind of out-of-band processing (for example for electrode impedance measurement or ECG pacemaker pulse detection).

[0045] To mitigate common-mode interference, the digital filters 31 may be designed such that a common-mode part of the respective output signals y.sub.j is minimized. To this end the calibration device 13 shown in FIG. 1 can be applied to determine the filter coefficients x.sub.i so that the common-mode component of the individual output signals y.sub.j are minimized. The calibration device 13 comprises a further processing unit 37 including a second processor 39 and a second memory device 41. The further processing unit 37 includes Input/Output (IO) circuitry configured for coupling the further processing unit 37 with a test signal generator 33 of the calibration device 13. The IO circuitry is configured to access the first memory device 29, in particular the non-volatile section thereof, so that the calibration device 13 can read the samples x.sub.i and program the filter coefficients x.sub.i. The IO circuitry is configured for performing IO operations to access data present in the measurement device 11 and/or to control the measurement device to carry out a herein described calibration method 45. The further processing unit 37 is coupled with the test signal generator 33 of the calibration device 13. The test signal generator 43 can be controlled by the processing unit 37, via the IO circuitry, to generate a test signal TS fed into the individual inputs 15 of the calibration device when performing a calibration method. A computer comprising the further processing unit 37 may be programmed to perform a herein described method for calibrating the physiological measurement device 13. A respective program can be stored e.g. in the first memory device 13.

[0046] FIG. 4 shows a method for calibrating the physiological measurement device 11. The method 45 may be carried out by the calibration device 13. Accordingly, the further processing unit 37 can be configured to execute the method 45. In particular, the second memory device 41 may include a computer program that is programmed to execute the method 45 when run on a computer such as the further processing unit 37. Alternatively, at least some operations of the method 45 may carried out by the measurement device 11, in particular the processing unit 25.

[0047] The method 45 may be executed as a part of a production final test related to the measurement device 11. During the production final test, or in general when performing a calibration, the calibration device 13 is connected to the measurement device 11. There may be an electrical connection between the test signal generator 43 and the inputs 15 configured for applying the test signal TS to the inputs 15 and a communication connection between the processing unit 37 of the calibration device and the processing unit 25 of the measurement device configured for reading the samples c.sub.i and for writing back the filter coefficients x.sub.i calculated by the calibration device 13 as described below.

[0048] After a start 47 of the method 45, the test signal TS is applied to the inputs 15 in step 49. To this end, the further processing unit 37 may control the test signal generator 43 to output the test signal TS. The test signal should cover the frequency range of interest of the measurement device and the amplitude should be large enough to cover a significant part of the measurement devices ADCs 23 input range in order to minimize error due to quantization. Different types of signals can be used subsequently, for example chirp or sweep signals, pulse train signals, periodic signals with sharp transitions like square or saw tooth waves, or even wide or band-limited noise.

[0049] In a step 51, the method 45 inputs or otherwise determines one or more definition vectors v.sub.1, . . . , v.sub.M, which may correspond to the leads supported by the measurement device 11. As shown in FIG. 4, the definition vectors can be arranged in a matrix V=[v.sub.1, . . . , v.sub.M].

[0050] Step 53 of the method 45 inputs at least one DC gain value DC.sub.i of a specific channel i. The DC gain values DC.sub.i may be determined in a separate measurement procedure, which may be carried out during the production final test, and describe the gain at frequency 0 within a specific channel, in particular within the analog circuitry 19 of that channel.

[0051] In a step 55 of the method 45 a sampling frequency f.sub.SAMP and a resolution res of the ADCs 23 are configured by the calibration device 13. The sampling rate f.sub.SAMP may be increased compared to a sampling rate used during normal operation of the measurement device 11. The resolution res used during the calibration method 45 may be reduced compared to a resolution used during normal operation. For example, the sampling rate during calibration may be at least 2000 Hz, wherein a typical ECG recorder usually operates at sampling rates between 500 Hz and 2000 Hz when performing ECG measurements. To configure the special sampling rate and resolution during the calibration, the measurement device, in particular the multichannel analog digital converter 21 and/or the processing unit 25 may be switched to a special operating mode.

[0052] While the test signal TS is applied, the measurement device 11 records the unfiltered ADC samples c.sub.i of all input channels. In step 57, the method 45 receives the samples c.sub.1, . . . , c.sub.N related to the individual channels and representing the test signal TS distorted by the analog front end 17. As shown in FIG. 4, the samples may be represented by an L?N matrix C=[c.sub.1, . . . , c.sub.N]. L is the number of recorded samples.

[0053] The recorded samples C and the definition vectors V may be used to establish an optimization problem to minimize a target function TF in step 59. The target function TF may depend on an output vector X including at least some, preferably all, filter coefficients x.sub.i. The filter coefficients may be represented by a column vector X=[x.sub.1.sup.T, . . . , x.sub.NT.sup.T].sup.T which has a length of N*K elements, K corresponding to the number of filter coefficients per channel. Accordingly, all filters 31 of the shown embodiment have the same number of coefficients. However, the individual filters 31 may also differ in terms of number of coefficients and the calibration method 45 may be adapted accordingly.

[0054] In general, when considering cases of more than one digital filter 31 per channel, the length of the coefficient vector X is 1 . . . imax, N?imax?N*M. The minimum value for imax represents the case of only one matching filter per input channel, and the maximum value represents the case of one matching filter per input channel and possible output vector. Without loss of generality, the following will assume that imax=N.

[0055] The target function TF may include a convex or quasi-convex penalty function, in particular a norm, preferably the L-infinity norm (maximum norm). Other norms like the L2 norm (Euclidean norm) can be used as well. The optimization problem to be solved in step 59 is comparatively large due to the number of samples but can be handled with existing solvers for optimization problems.

[0056] An argument of the target function may include the unfiltered samples c.sub.i, the definition vectors v.sub.j and the filter coefficients x.sub.i. The filter coefficients x.sub.i are target variables to be determined during optimization. Based on the optimization optimal or suboptimal values of the filter coefficients are obtained as a result of the optimization. The optimization may comprise applying a numerical solver for optimization problems.

[0057] The argument of the target function may comprise a column vector of all samples of the vector signal and can be formally modelled as follows.

[00001]$Y=\left[\begin{array}{c}\begin{array}{c}{C}_{v1}\\ .Math.\end{array}\\ C_{\mathrm{vM}}\end{array}\right].Math.x$

where

[0058] To(c.sub.i) is the Toeplitz matrix of ci elements, size (L?K+1, K). The first column contains the last (L?K+1) elements of ci; the first row contains the first K elements of ci in reverse order. Example for c=[1 2 3 4 5 6 7 8 9]T, L=9, K=4:

[00002]$\mathrm{To}\left(c\right)=\left[\begin{array}{cccc}4& 3& 2& 1\\ 5& 4& 3& 2\\ 6& 5& 4& 3\\ 7& 6& 5& 4\\ 8& 7& 6& 5\\ 9& 8& 7& 6\end{array}\right]$

[0059] CA=[To(c1) To(c2) . . . To(cN)] is augmented input sample matrix for convenient calculation of all output samples. CA can also be constructed from Hankel matrices; this merely changes the order of elements and hence the order of the matching filter coefficients in the output. [0060] ? is the Hadamard (element-wise) matrix product [0061] .Math.: Kronecker product [0062] C.sub.vi=C.sub.A?(v.sub.i.sup.T.Math.1.sub.L?K+1,K): Sample matrix that only retains input channel samples used in vector i, multiplied by the appropriate sign and coefficient. Note the use of Hadamard/element-wise product and Kronecker product. [0063] C.sub.vi.Math.X is a column vector of output samples of vector signal i

[0064] Based on the target function argument, the optimization can be formulated as a multi-objective optimization: Minimize 1-infinity norm of the vector of output samples and also minimize another norm (or norm-like penalty function like the Huber loss function that has the robust behavior to outliers of the 11-norm while being a smooth function which the 11-norm is not), usually the 12-norm of the output sample vector (with lower priority), i.e.,

[00003]$\mathrm{minimize}{.Math.Y.Math.}_{?},{.Math.Y.Math.}_2\mathrm{with}\mathrm{regard}\mathrm{to}X.$

[0065] In general, the optimization problem should minimize the 1-infinity norm (i.e. the maximum absolute deviation from zero when a pure common-mode signal is applied to the inputs) if the subsequent processing of the differential signal is sensitive to outliers (e.g. pace pulse detection in ECG) instead of signal energy. However, some 1-infinity norm optimality can be traded for finding a solution that has a small 12-norm (or norm other than 1-infinity). Depending on the used solvers, the solution of the 1-infinity norm optimization might already be biased towards also minimizing other norms. If this is not the case, the multi-objective optimization problem can either be scalarized to

minimize ?y??+??y?.sub.2 with regard to X

or it can be split into a first 1-infinity norm optimization, followed by the minimization of a different norm while keeping the 1-infinity norm constant (possibly allow some deviation from the value of the first optimization to simplify the job of the optimization algorithm).

[0066] The optimization problem may include constraints. For example, a linear equality constraint of the form 1.sup.T*x.sub.i=b.sub.i can be applied in order to compensate differences between the individual channel as regards the DC gain. The value b; may be derived from the DC gain values DC.sub.i obtained in step 53. Similarly, the gain of at least one digital filter 31 at the Nyquist frequency can be constraint using a linear equality constraint.

[0067] At least one non-linear equality constraint may be applied to the optimization problem, e.g. to constrain the gain of a digital filter 31 to be equal to a value b.sub.i. This constrain has the form (sinf.sup.T*x.sub.i).sup.2+(cosf.sup.T*x.sub.i).sup.2=(b.sub.i).sup.2. This type of constraint can also be formulated using complex exponentials instead of real valued vectors of sine and cosine values. In this case, scalar products may be used to produce real-valued output.

[0068] Moreover, a non-linear inequality constraint may be applied to the optimization problem, e.g. in order to constrain the gain of at least one of the digital filters 31 at a given frequency to be equal or less than a value b.sub.i. This constraint has the form (sinf.sup.T*x.sub.i).sup.2+(cosf.sup.T*x.sub.i).sup.2?(b.sub.i).sup.2. This type of constraint can also be formulated using complex exponentials instead of real valued vectors of sine and cosine values. In this case, scalar products may be used to produce real-valued output.

[0069] Bound constraints, which are a special case of linear inequality constraints, may be applied to constrain values of coefficients of the individual digital filters 31. These constraints have the form b.sub.min?x?b.sub.max. These constraints can be used to avoid numeric overflows when storing the coefficients and to avoid numeric overflows when calculating the filter output. Moreover high-pass behavior of the individual filters 31 may be limited.

[0070] A further constraint ??x.sub.i?b.sub.i related to a norm calculated from the coefficients x.sub.i may be applied. The norm may be the 12-norm of the filter coefficients. Such constraints limit the average gain of the filters and the amount of high pass behavior of the filter. This constraint corresponds to a non-linear convex inequality constraint.

[0071] The present disclosure is not limited to the above described constraints applied to the optimization problem. Other types of constraints may also be included, for example constraining the value of at least some coefficient vector x.sub.i to have decreasing magnitudes, which forces the filter to be closer to a minimum-phase (and minimum-delay) system.

[0072] The optimization problem can be varied, e.g. to reduce its computational complexity. One approach is splitting the optimization problem into sub-problems than are, in sum, smaller than the original problem. For example, in an ECG application, the calibration could first calculate the matching filter coefficients for the vectors formed by limb leads (I, II, III, aVR, aVL, aVF), and then, keeping the first set of coefficients constant, calculate the coefficients for each chest electrode individually (as no vectors with two chest electrodes are used in practice, the chest electrode filter coefficients can be calculated individually for each chest electrode). This results in seven smaller optimization problems.

[0073] Alternatively, the calculation of the Wilson's Central Terminal (WCT) vector signal (RA+LA+LL)/3 could use its own coefficients, in which case the filter coefficients of channels related to chest electrodes V1, . . . , V6 must be matched simultaneously (i.e. during one optimization run) as they all require the WCT. Calculating filter coefficients specially for the WCT vector signal splits the original problem into two smaller sub-problemsone optimization run for the limb lead vector signals RA, LA, LL and one optimization run for the chest lead vectors V1, . . . , V6.

[0074] When using multiple matching filters per electrode, the optimization problem can be split into entirely independent sub-problems that can each be solved individually without adversely affecting the overall optimality of the total solution.

[0075] After completion of the optimization run(s), step 61 of the method 45 obtains the filter coefficients X of the individual digital filters.

[0076] The obtained filter coefficients X may be stored, in step 63, into an appropriate storage region of the first memory device 29 of the measurement device.

[0077] In step 65, the method 45 is terminated.

[0078] In one embodiment, each channel has exactly one digital filter 31 the filter samples {tilde over (c)}.sub.l generated by these filters 31 is used to calculate every output signal y.sub.j corresponding to a specific vector or lead. In a different embodiment illustrated in FIG. 5 at least one channel has at least two sets of filter coefficients corresponding to different digital filters 31a, 31b. These filters may be used to calculate one or more specific output values. In the example shown in FIG. 5 the digital filter 31a is used to calculate the output signal y.sub.i and the digital filter 31b is used to calculate the output signal y.sub.2. Both digital filters 31a, 31b are configured to filter the input signal s.sub.2 of one channel. In an embodiment, each output signal y.sub.j has a dedicated set of digital filters 31. In other words, the set of filter coefficients corresponding to the individual digital filters 31 are obtained for a specific vector corresponding to the output signal y.sub.j.

[0079] To sum up, the present calibration method obtains filter coefficients of the digital filters 31 by means of optimization of a target function based on both samples of a test signal applied to inputs 15 of different channels as well as the definition vectors. This allows for determining the filter coefficients directly without selecting a reference signal. Moreover, constraints related to characteristics of the filters 31 that have nothing to do with common-mode interference mitigation can easily be incorporated into the calibration method. The calibration method 45 may be used in combination with measurement devices 11 that filter the sampled input signal s.sub.i in the digital domain and then calculate the vector signals from the filtered samples corresponding to a filtered input signal. Moreover, a measurement device 11 is described that comprises more than one digital input filter 31 associated with a certain channel. The digital input filters 31 associated with the same channel may be associated with different vector signals.

[0080] While the invention has been illustrated and described in detail in the drawings and foregoing description, such illustration and description are to be considered illustrative or exemplary and not restrictive; the invention is not limited to the disclosed embodiments. Other variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed invention, from a study of the drawings, the disclosure, and the appended claims.

[0081] In the claims, the word comprising does not exclude other elements or steps, and the indefinite article a or an does not exclude a plurality. A single element or other unit may fulfill the functions of several items recited in the claims. The mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measures cannot be used to advantage.

[0082] A computer program may be stored/distributed on a suitable non-transitory medium, such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems.

[0083] Any reference signs in the claims should not be construed as limiting the scope.