A vectorial Matching Pursuit algorithm with an ... - University of Leeds

INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE

A vectorial Matching Pursuit algorithm with an application to the analysis of Vecto-cardiogram signals Catherine Bonnet , Jonathan R. Partington , Michel Sorine

N ˚ 4535 Septembre 2002

ISSN 0249-6399

ISRN INRIA/RR--4535--FR

` THEME 4

apport de recherche

!" ## ! $ %&! ' ' (())* + -,. "0/213 546* 78## $ )# 9;:6@?ACB >EDF6B B >@ v x m 5 5 /§ § 3 §J¯ > v x m 76 5 5 v{hosm})a5¨§z °8[a)p:9; n@s y.¢ hopqa9=< L 3&> ,- ?" D ]*^ adhjsmsa u r y k¦}p~hovÅsr ¢ @h£@as.iz

vm

$

§ {¯ v § {¯ !

§C T ¯

5! 102!

¢ ^ a a;§ ¯ hjv~p ^ akmvwkmnYl hos¦sa u r y km})pris ° ¤asr ¢ {p nVxa5nYv~rik y hj}p{hrismn zp ^ a5hos|Jshp{avqap*rx nxlol` kmsm})pqhor@smv vx ¦ y aJsa y iz § v x ¯-§5!{¯j v x §5!{¯ r ! y

§C D ¯

¢ ^ a ac}Ehovn5vwnYb;ulohos#u¦r@hosip nkmshop£@a}p{r hjs n@s y p ^ av{}nxljn q¥ £nYloka y k¦sm}p{hris vx nYvYho£YaKs#hjs§Vxi¯-° p5 v{p0vqho ^ p hjb;uloabasipqhjs8p ^ atb;nYp{} ^ hjsEum k vwkhop0nxloYr hop ^ b lr|rVivlh4Vxa.n@s¬hosip n@}p{n@loasik¦b#a hj})nYl u r@loab v{hos¦}a5p ^ a y ho})pqhor@smn z«hjvnYl ¢ nWzvhjs|6shopqa0h ¦kp ¢ adsmrYpqadp ^ nxpkmv{hjs§C¯

> vx §S} ¯ v{rp ^ nxp*p{r«bnxhobh s)a M v x r£@a nxlol#² } n@s y hov~a ikmh£nxloasip*p{rp ^ av{hob;uloa 5 §/} ¯ 5 r£Ya nYll`p ^ av{n@b;ulhjs«u6rYhjs@pwv© y } n@s pwn,iV hjs §/} ¯ K 5 §S} ¯ 5 °

§C,Ii¯

"

u r@mlaKb rx b;nhobh s)hjs

]*^ adhopqa nxp{h£@a5vqp{au rx p ^ anYl@r hp ^ b p ^ mk v an y v

, § {¯ § {¯ 3 § vx, /vx § {¯7 §S} ¯ KJ5 §/} ¯ 5 J§ {¯ 3 §-`§/} ¦¯vqhjsm} W§ 3 }6¯ § {¯*n@s y § {¯*n a£@a}p{r v*hjs r aKnY} ^ -°fihjsm})a0p ^ a#nxloYr hop ^ b 5!

! !

5!

!

§C.Si¯

¢ ^ a a0smr ¢ 5! ! ! jh vPv{pqholl nvwu¦aK}hjnxl })nYv{a rx p ^ a«§ a)ljna y ¯b;nxpw} ^ hjs«umk v{khop*nxloYr hop ^ b p ^ a}ris@£@a YaKsm}aÅ rYlolr ¢ v r@b ]*^ a)r a b °oY°

-,. &0/21

7

¤a6a)@hos ¢ hp ^ v{r@ba¦nY}Vi r@kms y ris8alaK}p r})n y horY n@u ^ z.6a r a;nYumuloz|hosm.p ^ ab;nxpw} ^ hjstumk vwkhopnYl@r ¥ hp ^ b p{rp ^ aba y hj})nYl y nYp{nm° hop ^ rikpP@rYhjshjsipqr y a)pwnxholov lapkmvqkmv{p aK})nxlol p ^ nYpPnYsEa)loa})p ho}nxl nY}p{h£|hop z aKv{klop{vÅ r@b ^ an pPam}hop{n ¥ pqhor@s ° ^ hloaPp ^ hjva)m}hop{nxp{hris risipu rY aKv{v{av p ^m r@k ^ p ^ ad})n y hjnY}Åb0kmvw}loa a)loa})p ho}Pu¦r@pqasip{honYlovrx y h ¸a asip br y kljkmvÅn@s yy h aK}pqhor@s¦vn a} anYpqa y ° ]*^ a0a |kho£WnYlaKsip y hju¦r@larY p ^ a ^ an p valaK}p hj})nYlÔnY})pqho£ihop zE§CvqaaB YE ¯ hjvIp ^ aÅv{k¦b£@a})pqr rY`p ^ aKvqaÅu6rYp{asipqhjnxljv° J p£Wn hoavIhos;loas@p ^ n@s y;y h aK}pqhor@s y k hjsp ^ aa)})hpwnxp{hris«}z}loa iy av ¥ } hjhosmn})k £Ya~hjsU§\p ^ aF !5P ¯ ¢ ^ ho} ^ bho ^ p6ah y asipqh¦a y iz#uljnY})hosv{a)£@a nxlalaK}p r y avIr@s«umnYpqhoasip{vnxp nYu¦u riu honYpqaPu6rYhjsip{v*hjstp ^ a6r y zY° J s. nY})p uljnY})hosm0 rik a)loa})p r y av ¢ hlolÔnYllor ¢ r@smap{r;r@pwnxhjs.p ^m aa5vqho@smnYlov ¢ ^ hj} ^ })n@s6aP}rismv{h y a a y nYv¦ahos#p ^ ar p ^ rYYrismnxlmu rYa}p{hrismvIrY p ^ a0F !P h p ^ aPa)loa}p r y aKvn aPv{khop{n@loz uljnY})a y ° J s¤kmv{k¦nxl }lohosmho}nxlIv{p{k y hoav rislzp ¢ ra)loa})p r y avdn aumlon@}a y pqrr@mp{nxhjs¤ntv{hosmYloa)R !#P ° ]*^ a«nYsmnYlzv{hovrx y })n hon@}« a |kas¦}z¬hov;n@}p{k¦nxlolz y risa. r@b p ^ hjv)R !" P iz³h y asipqh z|hos¤p ^ a br@v{phob;u6r p{nYsip;u¦aKn,|V v#rxÅp ^ a v{hismnxlI§Cn hov{hjs.nxpÅanY} ^^ an p6anxp-¯*n@s y }rismvqp k¦}pqhjsp ^ a#v{r ¥ })nYlloa y hjsipqa £nxlC° ]*^ a0}n y hjnY}5u6a hr y hjv })nYlo}kljnxpqa y n@v~p ^ abaKnYstrx p ^ a5hjsipqa £nxljv*¦ap ¢ aastp ¢ r«u6aniV v° K i£|hrikmvqloz nYs³n@smnxlozvqhjvrYp ^ a "vqho@s¦nxlv ^ rikl y @h£@a.br ahjs| r b;nxpqhor@s ° J s©umn p{ho}klon })r@b;umn a y pqr nE¼n@smnxlozvqhjv hopv ^ rikl y nxlolr ¢ n@san@vqhoa y a)p{a}p{hris8rYv{b;nYll u¦aKn,|V vnYs y nt¦apqp{a av{pqhjb;nxpqhor@sErxv{r@ba hjs@p{a q¥ u6anV y k nxp{hrismv° ©© jÒ ª«¬®¬

(m - 12® 2(m 2/ , o m u 2 I

] r hololokmv{p nxp{a;p ^ aba)p ^ r y` ¢ a ¢ r V ¢ hp ^ Wx@ y nYp{nu6rYhjs@pwv an@} ^ })r@b;u oh v{hosm})rir y hosmnYpqaKv ¢ ^ ho} ^ ¢ a Qa Vihos y lz5u r£|h y a y iz ° 8 nYk asm})aIEn@sYhjs]2p ^ av{a au av{asipÔPvqaK}r@s y v wF !P¨ban@v{k abaKs@pwvÔvwnYb;uloa y x@0p{hobaKvu¦a v{a}ris y °

r++U< wN $ S ]6)" > 2G-;* "" +m/$

Jp

¢ Yn v; rikms y })r@si£YaKshaKsip«pqr²vwkmp n@}p.n³})r@smv{p{n@s@p rib±p ^ a y nYp{n p{r©Yho£Yahoptban@s£WnYlokma s)a r hjs r y a pqr©nW£Yr@h y p ^ aE¦a ^ nW£|hrik rx S¦ahos y r@bhjsmnxp{a y iz²hop{v £WnYlokma©§ p ^ a8£nxljka8nxp|¯ mkp.p ^ hjvhjv sr@pÅavwv{asipqhjnxlC° ljvqr.p ^ a¦P n@b V¦f|} ^ b' h y pÅu r}avwv ¢ nYvÅkmvqa y pqrtu a ¥ u r}aKv{vp ^ a y nxpwnnYs y bnVxap ^ a0p ^m a)a } ^ n@smsa)ljv~r p ^ r@sr b;nxlC° hik a# au av{asip{v~nr@smad6anYp!F !P n@s y n;vw})nYlon u rYa}p{hris { ¢ a5})nYs v{a)adp ^ nYp*p ^ aF !P ^ nYv~p ^m a)a b;nxhjslr|r@u¦v ¢ hp ^ y hon@ba)pqa v })r w aKv{u6r@s y hosdn p{a u rxaK}p{hris#p{rdp ^ aÅk¦v{kmnYl¦u6an,V|v rYJp ^ aR !P ¯!° 6r amnYb;uloa hopb;nWz6a#aKnYv{ha pqr y apqa})p hosEp ^ a F !Pp ^ n@s8hjs8p ^ a §

¢ R !P nxlolr hjsan@vqhoa aKvqp{hob;nxp{hris.rY p ^ a y k nxp{hris.rY[p ^ a ] hjsipqa £nxl §\£Yasip hj})kljn * au6rYljn h snxp{hris6¯!° hik a~dv ^ r ¢ v ndulr@p rxJp ^ a~ klol y nYp{nP r p ^ ad}r|r y hjsmnxp{av nYs y ho@k a*v ^ r ¢ v ndp ^¦ aa ¥ y hobaKsmvqhor@s¦nxl ulorYp*rY[p ^ a y nxpwnm° hik aÅc#v ^ r ¢ vn0ulorYprxÔp ^ aP kllJ a |kaKsm}z«vwu¦aK}p kmb¡rY p ^ advqho@s¦nxl §\v{p ho})pqlozv{u6an V@hjs p ^ au6rYhjsip ¢ hov{a sr b rx p ^ 9a ¦r@k hoa p n@smvq r b ¦¢ ^ ho} ^ hjvÅnYlov{r«£YaK}p{r q¥ £nxljka y ¯!° ]*^ hov ¢ n@v y apqa bhosma y iz.pw,n VihosmnmnYs y ¥ ¢ h y p ^ rY !Sx@# "tn@s y hjs@p{a u¦r@lonYpqhjs ¢ hop ^ p ^ aPv{p{nYs y n y }n y hjsmnxl¦v{a hoav° ho@k aÅdv ^ r ¢ vInYs;aKslon @abasip rx p ^ a5lr ¢ ¥ a |kasm})zu¦r p{hris.rY[p ^ av{u6a})p kmb ° ] ,n VihosmnEmnYs ym¢ h y p ^ rx ¦ YYnYs y § § aKv{u6a})pqho£Ya)lo z § I|¯ p ^ a.n@umu rWhob;nxp{hris u r y km})a y iz KYhopqa nxp{hrismv*rY p ^ a0£@a})pqr hjnxl § avwu6a}p{h£@a)loz.£@a}p{r honYl alonxa y ¯ nxp{} ^ hos k v{khopnxloYr hop ^ b Yho£YaKvn a)ljnxpqho£Ya«£@a})pqr hjnxl 3 ¥ a w r sr b a |kmnxlpqrE¦° @. S¬§ avwu6a}p{h£@a)loz¤¦° @Y|¯nYv5p ^ atn@umu rW|hjb;nxp{hris¤u r y km})a y kmv{hosmp ^ av{}nxljn vwulohosmaPnxloYr hop ^ b §\nW£nYhljnYloahjs«v{})hljnYr b;nxp{lon@6¯ ris;an@} ^ }r|r y hjsmnxp{a*r@sKYda |kh y hjvqpwnYsip u6rYhjs@pwvnYbrisp ^ a« kmll y nxpwnYho£Yav5n a)ljnxp{h£@a;a w r a |kmnxlIp{r8m° T Y° r@pqa«p ^ nxph ¢ an@umu rWhob;nxp{a#aKnY} ^ }r|r y hjsmnYpqaEizM@¬hp{a nxpqhor@s¦vrxp ^ a¬v{}nxljn a)ljna y b;nYp{} ^ hjs®u¦k vwkhoptnYl@r hp ^ b ¢ hop ^ § I p ^ a £YaK}p{r honYl 3 a { r smr b+hjv. a |kmnxlp{r®¦° @ D v ^ r ¢ hos²p ^ aKs¨p ^ nYp.p ^ aE£YaK}p{r honYlÅnxloYr hop ^ b })n@s¨hob;u r£Ya nYu¦u rWhjb;nxpqhor@s ° hik a T V S.n a;mn@vqa y ris8p ^ anYumu rWhob;nYpqhor@s8u r y km})a y @zp ^ a«£YaK}p{r honYl a)ljna y Enxpw} ^ hjs k v{khop nxloYr hop ^ b ¢ hp ^ § 2 I { hik a T v ^ r ¢ v ndulr@p rx¸p ^ a*£@a})pqr hjnxlnYumu rWhjbnYpqhor@s ¢ hp ^ §2 InYs y ho@k a D v ^ r ¢ v~np ^¦ aa ¥ y hobaKsmvqhor@s¦nxl¸ulr@p)'° ho@k a I;v ^ r ¢ vnumlr@prYÔp ^ anYumu rWhob;nYpqa y v{u6a})p kmb nYs y hik a S v ^ r ¢ v~n«ulr@p~rx p ^ an@umu rW|hjb;nxp{a y lor ¢ ¥ a ikmasm})zv{u6a})p kmb ° hik aK0@h£@av~p ^ anYumu rWhjbnYpqhor@s«rip{nYhosma y £|honp ^ av{}nxljn v{umlhjsa5nxloYr hop ^ b ° Å%"7$ Å; ]*^ an@kp ^ r v nYpqa kmlloz.nY}V|sr ¢ loa y Ya ^ a)lju| kl y hovw})kmvwv{hrismv ¢ hp ^ kloha~aKvqp{a)l nYs y a n@s ljnxh nYb0¦nYklop)°

/xP Å*

B Ed²° VWnWz.n@s y RÅ° ° p0°Ekml y a J s@£@av{pqho@nYpqhjs#p ^ a5a¸a})p*rx bnYpqa smnYl nxlj}r ^ rYl`hjsip{n,VYaPris ^ k¦bn@s ap{nYl aKnxp ^ hos nYpqa~kmv{hosn y n@upqho£Yapqhjba ¥ a ikmasm})znYsmnYlzv{hov ¥ bap ^ r y v .30$ ;,M + c T §q SIS T ¯ Wx=VJ T c¦° B E ° °~a Ya fJ° |V v{a)l r y` °,PP r y r@sn@s y ° °r ^ aKs x s#a® }hoasip nxloYr hop ^ bS r vwu¦aK}p nYl|n@smnxlozvqhjv rY ^ an p nYpqad£n hon@holhop z ! (m 3 2.3" 4" -;*" > 2'- .3" 4" ÅK | r y sho£Ya v{hop z avwv §{ SIS T ¯-° B E f¸°?#P °6EnxlolonYpnYs y ° ~ ° ^ nYs Enxpw} ^ hjsumk v{kmhpwv ¢ hp ^ pqhjba ¥ a ikmasm})z y hj}p{hrismn haKv 1! mE" - > 225< ®"E" c6§{ SISYi¯ Y.S D Vxc¦W° B YE ° #F ° nYs¦lr£n@s y ° 0F °mPp rYl y as¤§\a y v°Ú¯ $0+U/ - # M ,3 ]6(mm r ^ s holoa)z8§{ S.S T ¯-° B KE ° ° n pqhjs@pqr@s ! 1® +m / [ ;,1/ [/ -;*"$0+U. O ÅK | r y smh£@a vqhop z avwv§q S.S D ¯!° =

=

©© jÒ ª«¬®¬

(m - 12® 2(m 2/ , o m u 2

0.21

Rx

0.17

0.13

0.09

0.05

Tx Px

0.01

Qx

−0.03

Ux

Sx

−0.07

−0.11

−0.15 0

0.4

Z 0.10

]

−0.01

0.8

1.2

1.6

2.0

2.4

2.8

3.2

f

−0.12 −0.078

0.028

Y

−0.064 0.050 X

0.134 0.163

hik a;{ r@sa56anxp!F !P3nYs y n«vw})nYlon u rYa}p{hris8§SR !P5¯

dfe © dhg

0.19

0.15

0.11

0.07

0.03

−0.01

−0.05

−0.09

−0.13 0

0.4

0.8

1.2

1.6

2.0

2.4

2.8

3.2

[ >