動的時間伸縮法を使って2つの心電図信号を比べる
WarpingDistanceを使って,心拍の心電図信号2つを比べる.
In[1]:=
ecg1 = TemporalData[TimeSeries, {CompressedData["
1:eJwl1wk0VO//B3BJZSkktEklZWklW8igDdkSUSlkS6JsWVKypFBkKYpKWrRQ
lCVL3pLs+26Yca9dqaQIlf7P9/d3jvM5M889M3Of87qf5/1Ze+KMiT0nBwfH
HPI/MpeD487t//46ERrw/7VR+P9rk3nC/+qdZ9f/V4/u1Uj4r1YJCN38rw7O
47nxX21IfXfxf9cPG274rz4PWncgntSnydkpN0j9+snr5zlSpZ/lKYiQ6qpX
e744vhM83FYvw2M7kXl0mPvclU4E209yap/sRHLf8f6u1Z3YkRa83j6sA16c
B60V3duxw0rcuP9xK75ZcunYczdjm6eWb8GZelgd4RTV21eBjLTwrWVzSyC8
MVM/x6kQTw/YMaL7C8BZ7WKmbP0epS8enLS/Ww5jMZ7v99bXo/b8/HWGys0o
4VNf4qrcBuWp3X/kr3ZgQmpnp8uKLohHRpjWMXugerbkPZ3bizt6C1+zt9NY
MeY+7c1Vjo4nFUGB7C7G6nCH3puZwwwN+faOaotOhmi9WzdTuREdgt8PntIf
QMKIs6+szSAsxNMtfy7vBVNDZ0D/HQtua2c0tY90YZ737Z3KRkwwbbRKguU7
MYdL5L5MZjsO34gXm85qRWLrySE5pWb0po/vF1NvwKqqJHPXkkpczbH6xPO8
BNf4eZlpjrng0AkNKP7xGJqNgiECuZ6gEku37tkZw+DYY7folF8845KG74Ji
3WCGoG+kmv+/m2gZVQi/IJAJffesnwVrgDdP+niPK5RjzV1Po/KWOkyJ86z4
OdkE/0ranFOoDRzv5cs8H3YgVHe6Ke8JE2rjBi//ru+Gn+dLKW8FFnS/K44s
H2FBfHNEf9ARNkbuXGs9lcTGmhmbj4FFbCyK5aWO17BRIZH/kLOOjacTiv3J
FWyoa7ttOA02mDGZokav2dDWvOF48DEbgu91x3wS2MjKnJBEBBvh9f456wPY
+GQ7FJrmzoan+54vJg5sBCRvrF9Gvi/p2YWlswZsSGyPK+LQZmOFgaidmBIb
l53enNkly4bS9ru62uJs6Bik7D+9mA39VO739XPYYH/ObTv7iQVVh0fB/MUs
cAScvZ5wmoVLAkFzt1Z349aTOH2pZ134OvExd3lTJ8z2h76ZkWmHSdL9itG6
Zlwa9DLW3daAYOeHvM7cVWjjzLFcsr4MTKcnLap1H3GaTnXbLl2JyM1i768d
qMfeO1w7BW2bYW6p9Lvzbxt21p7ytf7ZCbaZw9VEtW4w6L2O9qXkvn87NeS9
o9GQP3H+VnA/mIvr4+b87kWhu6HBIH8dDuldeKV5ooNhmyW+i/41yHh7oZ+p
9rWVIao3lf7ydguMh69s55sehPh2zW+hPMO4zpfTL1rUhzSLHP/EixSi+iqS
eqPZeJSZ9e/xRRZsGmpXtB7uhkWI7kX3BiYWRmskGtZ0wNwgJfjpnjZYO4sv
LN3VjAVBm6dNSuqxe2vKwNM7FeC/x3PU5OZ7UM8WLaXv5eBGi1GY7+Nn4JjU
ihc4lYg1SgHDK/NuYuqznRTr+gNY+io5BMVmQCF7beLg/EJMH1gmuHy2FD4x
ayx042pgORYydmlJE8q62JG7brdiwxGXuKT0DlR8TlriptKF1gnTO5KKLKzy
5Tzdl87Gyp1VrJHsHugfdjkrIk+hbF9ybLIjhW0BFsHulylUyxS8toulYOqe
uVQ7nsKJI+/qeW9SyHUaTrOIohAXHzdceoWCh91hxd1kX671ReR1eVJYvVlF
P8KZgsip8iCTExQKKlz2yR2m0KPisVvKmMJ0we3dCvso6ITQnBYaFHyWBnyN
VqRgVWj4md5EYeexkdQ9khRm/ym4Za+kIB07NSW6hMIa9x9dwrwUJJmurq7/
erDPNPPubG8PDIpkPyS86cGfcW9d2aPktfWHruYsNnwleCc/X2dh4EFB4PzX
XbiXt+7+6nmdkLnQvOow6Yctlgx5Hs4mrLt2eMuzlFqYjHJvlsypxPSzV0Uq
hyoRLTQzKZ1ag5nw/H91to2oPGq1683yVhioi68Wz+iAUubL4MMpXZBIHys2
GWBBNO7al28KFDrM9qw9FtmLNIdnkjImA1iz9GdLSHkfxO3zhxyTGnDxqCHP
gcY2houOo62T3CCjZb6N9U77FobkITqoQrwNiQ/MEp7kDWFKLcNbuWIYW8Ii
mbRfP5Kfafhd3kHD9Y3A39k5FJLtfnQLfWEj7qvLlaRSFvjOqTXeVe5G7Oxt
vqGtTHzf3yyV+7gdxqdKBg8kt2DMptmGodqIese3PL/2VKNxQ2Bp845SCPsa
zXNhFmBTz97tkxyvUbxfNm3xr1RMOVdM+3KRvniqlD2W+AL63bN6edez0W3l
ulrFrxj3vLiD413LEZWi1Hljtg7SavXno/c04+Nb61Cdijb0ef+KvtvZCdX7
ed1eNt1I3nE9V+EoG+lVhn7zG3tgu1nt+RLiZIHYL4GPxRTCWXYF/T8pbLKa
ajRcTJ7fqPNF3ato1IWVuBxeQ8NH3fdHqRi5f7p/aZ4IDc4bEXz1C2lEbNJ1
+z2HRkXxebddvyhkhpxsezJKvHH/27i+l4KZsdjRonYKexzZE261FNrqPnYw
PlB4yyUTLZVHrnv6uVb6FYUWczPRXY8p8K/jFfNMpFDaMPd0VjSFKq/Uuk/E
Oa+2skfxBQovR33y1d0prFQanU8T3y1q6rkL9SgcUPv2LUeMAjefydPeGz0I
f3o3uNqRDb1T0quLw7oh9PuzTjazEw4MjtganzZ8D3FWtm5uQsQf/yoZLXIe
cx37HGJcjZUrRk8FjlYhwaHJ7uyOOnwNEHgxJtsEFdWk9jZ2KzbZbNxg4NGJ
uKgnMjdsu8EdFTsQ/4iNQKvcxffJ8/qB38xAbLQXjxI8Xk2+HMCjuQFSjuv6
0aSSssvCpBGiA9ts5Ve1MX7v2rtRo2OAoX9uG8+rrGaGeOOk1Y3aNuiNxfFv
ERzGoy126fSqERjnP32dXd+PHbtHOAKu07DsvvpZwozCkb2Ln72O7oH/F78F
t8l9Wr/Slfnd0Y0m8B9/Vs+EqvTD8N5dHWCOh579pN4K/2afTzqtjbBfssP/
p08NHB7+k/zg9BGWW/drMxSKwB1SnFVtmgXNDfyFVEwaMrrDRiPinqHUuV9U
XjADWbGf3spV5cLyyXuNRNI3JxRNGqSvVGDvj+qijbwN8BWynLOkvhnufdof
nY3bYWQhKOJgzcSVZe+Fhga7kSHz9ceBHjY0NHPjR2QozFE0cW1+QYHPzXug
h4PG1ErvPZYKNApL+DIljWmI8ExHClnSuPc4u3zsGA0P7eBTOEyjWshDbZkJ
jZojYw5xOjTm/bQPl95Jg57uru2Uo5EhFTKRuoFG/xPe8JiVNBaej3G7LUjD
Ir1XqXAeDVnjN3kzM6TfTZr6mo5R8Ppsu6NygPg0Tzey6CL9Ucgw8l8DhfJa
6bDcMgpiIiPZ9oUU6qPnCWZlUqDGYfOCOP1ga5FYQ/rxp4kz2t7+FNK9Ft58
YEhhiQzL8e2PHiiEzJN9XEnO84P/spI/daMuu134nioT9IGaoPftbahcuE3I
Vq4Z0bfGrh38XY8418aeJu4aBD/Wyud6WA0ug3YX0a912HJkdaFVSBOCdNSm
bUneiwxvr5iu6sSh2W8nLfO78Uq6ec7VGTYGNAo1ZIso3LYP5VIX7cPU49Ht
zzsGcMXVOz1dvx+upmeTlzxqxFN2sPXumlZG+1zNs1syBhjif9ztDts3M5Tt
ji6WFGuHWKKTkYDyMITuTrRb7x7BDU2nlQ0T/RBe/tZG7Q2NBfSOWwF+FDji
l6Skd/XA1MDhQXsGGwuPee4IWMuCobXnS8+lXeD0Mv07Gd6BOFvOVt3gVgzd
fB4wsaIJPxbI645w1uJXfo/q4fGP+KNqI5itDdQM3bsSbp0N68lbT5Y8fgmV
azctLz14Ac/7m2KEhDNhzard6FP9Fh3zQyanhEvgPhNesHteJXR/xWvd0WmA
hdjAXPbCFpSUOZonXm/HZZ00/E4iee9Y8ntKgoVgOWH+jBU9KDb0FDUg598/
g1GD0QoKP65X9b4kfS5n+50FvcRdp5znrM1pGg9XzQvl9qFRtbZk6L4vjenb
oTulvWiUqE1qq7nQ2Hi/XtndlrjrZQ6/N6ehumap81Z9GiccGy/ka9KQ7zv7
9IQijQez/vVSsjSKvAXO86ymcS5UP2KBMI38qUa/tTzEaVyvsPksBcH8goTU
HxRWyPB9Ehyh8DttV2IEm0JAQK/A32Zy/i6ee1CgkkJ8oJVMNfF48Ox3/8w0
Cou0+RrXkr7zqjP4pJA9BZVWre4lghSWV/PoPZwgOWnapGyLKAshMx+3vrNn
Yus33nZ3znY4XxdOTj3aDPbV80lBUg2Y/hiT9upoDbydNi4SmKjGe5XH+kaM
epgcuhUel9YEi4srUuSt2uB6XXfwxHgnEoRKv1vR3VDP8xuSXtuDxQkPzuk2
knxWwBDaKNOH+rBNISc/D8Ay2KlX3aofSv1BUfUVjVgTweH/ObWVYb9rZbD8
vQFGzFXh6EKtZoZ6ygvagdGOnzEOVcE6w1Bx6vd/ajaC0wd9kj7NH0BVv0WQ
bCkNR7lJee0ICjMlxuK9f3oQ0iBrmtvIRrfv2xEebRbeOWS9fqHcBS6HvsfM
5x148vDjTlZKK5hrwnYrqTaB67LeSLlmLa79Ur95c2sZGna3p30uADJmZ+iI
jmxYJHS2Nse8wkn38oiX19KgEJG3JPJqJqYc/GX6NPJQ1+iiuNywBKbMuZfa
dlbi1vywxgS7BhS+XJ8UL90CFbvtdYuetaOO/3i35FsmllwsSyG/AGqfnvwz
Vuwhz5/mPX9TCkFt7otmWimYz3ecuEN8xFp3KqWTPlcXXP9rIfGmcYQ/NSaE
9DntA9K8oTTW/ta47BFIQ91nvZgjWS+WeFxz8Qw5p2PGal8Rh8YnV32ZJQ4l
YxZ+ciYOAwTu7J8kDmNcgifuEofX9GtdT5Dv2VadGL+LOEzflj+sSRyOR0du
O0IclnErP4okDuWr0rd2E4ciT6ds1YhDx6zgu9HE4ea8jC9XiUO/1+9EZcsp
TKYs4pjMoaAk0b9R/y4Fp6/LeiTOEL+vd7/iXUHh1IVrm47z9CD0hoFjrgwL
G3W/LhT2ZSLFT9EnQaQdBhk9PkmuzWgZNNmwT7UBm3Ll9nReqkG4dqUyc3kN
Wr5zRRXa1iP8h0RWUlETrNf8mF7o1oZ9FdfvC3AxYcPZaMYa78b6segmb/ke
XNrqktrdSWHQrqb/iVwfZiXDIup+DODqpsOrIx37sWBxOluvsxE5PPeGY+Jb
GevNThsJ3xpgXEzPjM6Va2YoTs5r0jJqh11MoY+l8TAGo9bccDk2glfVUqa3
+QcgtdHQvrKaxl6/KrophuSPBwGzn+ZR+BXt0x9KkXmw5a1kljELU+/X5W3b
24XZr6sWq+V1wEf5cNOKzFawK36vMNRpQnEIOz3bshZDXCKbNfeXoWXgA5U3
DugIO59q5c+B9EBdivzsK+is/BGw5lca9GUqIze/I/MuZfjut3seTh3cHtTl
WILlDQo9ZuaV2OJ/Z0zFrQEF/UI/RBRb0LXBtdI1ux2RFkX3bEuZ2GOTKjik
y8L4ofzPZlo9eHg+yvnMETJPlKv+YJHzzokVq3V2G42EfzsDPIgjDQMz6sNF
Gh8zHETUwomzsuNpDyNoDLssov8ShytvbNHqI+v6cnKTo57EVfo6CUEn8rrL
PeoA8fvy8aG/GQdo+Me7jG/eS8N9/lKTClUaXy5RKQFbaSg96i89KEk+L1l2
2+7l5Lw38Mgy5qcROWkU5DOXhsyrq2F5UyT/+R45svArBS5NZcNDfSRXPt/q
q9VB4fUW26S2avJ+j8i9WwWkL24T5KAfUPC+zs2d5kHhb8J6vaviJM9uvRSo
LdQDzRQu0cvyLCS9WmxYHchEne/g1R3i7bgueuy6lHcz9ohYehVqN0Dg1Gtm
QnQNknQr1M/I1mD11S5tKc96qAe6rHpX1oTjfIcKl/m24VKLef38hUwMjGYL
FM90Q1k+qm6vag/WKaduymdRyOBpOu+n2IdlP7TOJvwagOyBRccMnfshW5Cw
8y3ViFKplR4qN1oZb/wSg8/GDDA4Y85YB8k2My4+y75md6gdR8wjjPsODuN7
eqjMV+sRTBV9C1wiNIBNm+yPKtbTaDZvvy1/i8LDBKff2jwUYpyv75IcYqMl
QOq2ljkLiqd3PX6p3wX7C2Giz9CB2ScXjH1ySf57E1Rz27AJkjIFhfJOtVA2
94tKMi+DrGT6iaR5xYg7pTM2XyoHml8alGWXZmDb25U7gzekI25MuyCjJRPL
Xj8o0gzPw9kNgZFWHiWoVjmzqeBEJbg1I9hd5xoQ/fKfabpqCz5tCvoznd8O
JdPlXgOVTBQnX3l20IgF4+1nvmft7UFw3G/m0mMUDit5GgeQvrJ1sdwP3u00
Lh1jFi5wIP7Wu3w+QPrcpy960yXXaGzVuvBgfSRxdOlojm8YDU+VhxMWZD3r
SMs5S2/SN2+8Pu3rTOMVf1J2lhWNo4dz+flNyecZhs2GkHzoJbrn5HKSD/9d
3t1URfKhom+S2G2SD1Mm/hQFk3z49kdZWDjJh1sL3F6kkXxoI3JNZJjkQ9va
MXNFkg/PswqqfEk+DMlXCrZikvlb4K/WNzLHcJg1jj57R2ELds379pCCTIlb
Y5YXhauhnf0RayhcOrfMkyFK8tFDIcpbiQX5nhnurMtM7NCTMF4m0Y7E34Y+
f/yaEbo6MunKXtL/glW6V96ugf5fJakeuRoYFXrLh5yvx8uAr9tOVDVBpDbn
WpZ/G1qKi+1uCzDxdcEcX93ZboTlBR4r39mDu1u4fnv1kPmtKa9VTrkPKx1G
zu6fHsDp2dM8f0/3o6hdmWtRfyP8Nq6WHr3WysgLHImvuzHAuDd9uWavVDPD
60MbM/5wO65yeA36mA3jaaJXbsCJEcR05PM8WDKAlzV/JosaaPDrKO9/G0/h
0Mqo6lJeCmVL45TWfGJjw5nGccUjLISKlfskGnVhX9Tzc5ElHTh+zcTkUD7J
f0v7L5490ARWEf200bUWng/SNm0+VoakBpEWcb5ibHPTy/8ml4Op7zM7fq3N
gGBxnMGgWjoapLdM7+3JhMX9OYsK4/IQErunK8qnBIlWbG22YyWm47JYDr4N
EFa50im3swXJl4VW+LxrR5jiMcVjNUwoD4pIfDzAwkXZAbk1uj24NtRnr3ec
zJmMlLqXZL+c9nGdViZzR2OYjIqcI4279WvcfIJoSFxWjB4k85bCWJ3J3iga
a1Q1rsaTfniU0/ybB1n3iYi/70ny4QO/yupIkhelH8mYlFuT93/i2yozGt27
5m2M0aVBySwfWqdB446RxvMmeeLtXKjpPSka33POHwsl83Qr9eFmOJm3PcdH
ZJ/PJ9eX7xbo+03m3b+BTZu/U3hzVqbmzCCF3mI+oYOkX4t6h9tRdRQ0zbc+
SyB5nxnqZtP1iMwhHXea750jffF9T4zHWgqRvU3beZf14PVZHQENFRbq+0Tf
eF1hwj3IK7B1XTt6FWzVsv2bcUr09B81kp89ut59O59UA6snN58oKdaAzoFs
88V6zJZWe86raYLPnUWpBy+24VuRxWGFxUy8WCRkRP3rxkh/iMoRBtnX15wO
OykKbuVOmr+Ivydhe4f4Zwbg/qnmepFLP+6eWrH92kAjNp0SrTwb0co4pihe
J0T8yS/VGWta38zwteTj4D7ajrhHS9dvPTSMdA+TDQq2I3j+Rvz5QeEBcJu+
bD7TSPanKtLeMYHCzamHUu58FPSfPD3v/ZmN5tyi2rCjLNjelxv/bNyF5Xfn
xXd96ICY7h+J5IJWTKUWzs81aUL1pyUqFm61UGhkGaRblSHZTjI+YVExdl8t
TG5TykHS5qJHjRsykByl6Ra0Jx2Sx1QruPszUXjG5b3b7TxY/jqX3e9Xgm2b
d27efKoS3cMtNhF+DZAWsvI8qNECwSKLlltF7TDPzlriW8vEY4nJzWwTFow+
m7kr6fWgLdF7XQvxZ5Buv1mZ7Ne88f1nioi/zzG+fTnE32Tu3Io/xJdk/cZ4
J9L3wh1F+SqIP+ugxZNLyTksR02J/CDrRvq6QePE354LtS58ZB6Z35SVo21D
zu1z08sTiD9Xil+cT49GUAvfk0Tib/19aRVd0mcrwryylkjT6Cp0jpkh/jJZ
9Sl/ib+ighIhsQU0uLQmWg7+IfOtB4fvXeKvg0f2EJv461weEFRM/CUu3Fek
X08hWmv861xQKE4YvbmHzMNtZbHTC7wp9Cj8Tusi/mofppt1EH/jT4QleXew
YDEwAcZVJjRHv0XclWwHv7+izvkLzRAOjRDk0G2A/XtLeY57NUgUO8Gfq1QD
L+bdd0aX6iFU1NzxnvgLEjLYxB/QhsZlAXO+EH+bLD1WxXOwkHJ8Ql1AsweD
Qrxz4sl+7kpjTBip9KGbpVfgQPylhC4KEnHtx63AB9Yqg2T+rV5kUhbeyjDr
rj5dFDXA2LtevXMj8Xe8ot05nfgrWvSaO5D48w8M3XGZ+NudGPflIfHnPrFC
H8Sf/8d7S3KIv+xfWxcVE3+RvD/yFo+yMcOvESxhyUK4e4VE8IEuHCm4tNet
tAM2s9fMFApbITkVCN2DTTjq/fadn3stJv16v9VYl2GhfLplNj/pfx7FAbw7
cqD/s37on3QGim+lspt100EtTf6zcTATSZw+l6MS86B3Un9g7HwJpCulXys5
V8I+fbSR/3wD6ds3b7YSf6XcMc+Xoh2xokseTBF/y+5lnLM8yIJd003fN8Tf
MdmiBUJWFNxneTw9yX5dbKhd/pv4WyTv8mOc+Bvt/aehEkzDpaHk3CPirzD+
nSnnDRpPof7iAPGXartaZRNZn7lCC20kc4jQckPXXcTf/heR0n7EX596/Gwd
8ZfX3nFOi/iT1nXwbCT+ypU9Llwk/oQ0Bdx1iD+hny8qt5C5W79upZmcEI24
ISMtI+Kv/enawFDib6ZyzaN64u/GEa85gkMk90Xte/Wd+Dsm3Oxygfi70+cv
oUD8Gda8MPQi/uTPlKkoE383rBs9eCVIfy/mnVBY3gPG9ZB2e+LvUU5Z2B3i
r0b6nusv4q/qL72nhfhjPz670YL469rQm/qS+DvA8q+wVa6Bb/rPohnib3mY
bYBobRMmduX/tiH+OJvCZxhCTHyuT9YaIf546o7GWhN/x4yP7gbZT675Eo3e
xF/IMmn5SOJPc51VnDrx9+/EjqLLxN+o+MU8X+LvctG1oATib1J0oJSSbGZ4
r1Ap/kX8PZ890fKY+FvXEJP67D9/fXqz5cTfobmxc/qIv9M8Mq6dxN8W2QeL
+oi/DeJh3ObEXyCH2VpH4k9TfrVQOfE3FVxzNYv4U3B7NXqO+BM91pcWQfxR
ei46f4m/nw8YWwxtykCdGj6gKFCMGwalrAWqOQgRulM8I5MBn7CS2NL96fgz
com1digTOr8yfwYl5eHOvCoHJ/8SNOlVcbwk/tRNt9O6xF9s9Al+QUYLBuMf
Rx4g/pwO2n6Xq2MiT32t4k3izzQoJvUz8VdUKNG+mfgr+Xn7TxTZrzVZipai
ZD7dfOfE2KKTNB6F7fhlQnztcZr4VET8XR+/cVmM+BvisLjtTPwJymzs30PW
vYRCL+wm/mbztidYEn9Ta44KxRJ/wxe/PB0i/u6XH+k9TPzF/3vzYIT4s8uV
dblF/M3bWSdpTfx15C/U3Ev8+dJFUTrE3/bDxZoOxN/Hoj977hJ/zxL+bRkk
/sZT+0w2EH8XDFnlC7qJuzEJw1ji772zcLUu8de9T7Q4nPg7JxXroU/8uU9c
9VlN/LUp2sceJ/68nbr5oog/WmK92Qfiz2TD2qbV69vha31R+w/x9zprQ0kA
8Wc3RPP9Jf5UjTXsM4m/P2LPbu4LrEdYV+FVdeKv9P7JR1eIP9W3/dr2xN+2
eWay/HNYuOV8e/Vl4q/XveZoI9nPp/LL88OIv1uT+V/uE3/79pVKGBJ/z0ti
A58Sf08upBmaE398njuEIom//wPA8bVT
"], {{0, 10,
Rational[1, 100]}}, 1, {"Continuous", 1}, {"Discrete", 1}, 1, {
ResamplingMethod -> {"Interpolation", InterpolationOrder -> 1}}},
False, 10.2]["Values"];
ecg2 = TemporalData[TimeSeries, {CompressedData["
1:eJwllwk0Ve3/xTPmrYgiyhAaREVoUOEi4Y0GGqQoRBFNKIrQJKkUIUpSIkUZ
kihsZBYuLtd1h3NM0WxIRdL/eX//u9Zd33XPOuusc59nP/uzt5rrMTt34SlT
pgiR7weRKVPuJPz34YAl8//TpS/uf9Nq+q3/TZMEif9NA4XfYf/NUbFVe/+b
/q75cv9NxY9s/3gyVzcb3Ykh09yhNDqIzF1qq1SmkTnHgSUtf5uD8OpVVwtu
kOeHjGol+ZP7jR49/aHPQZ/CpeDKqA7E80wHz1uy4cuYu+OkTBumDC92dfjS
jLKGK7a7nr9DRkibhMKtKuTsz9SrffkWlcnxn9ZbVuLoS1O9rqI6cO00njXO
akad5Uw5lTQWEvaL0+aH2RARC/UMf87B3x+KVzan82FY87eNkdYFWTFXxQuh
ApSKXly9QP4W/v4K83aT6WKIf5oSuTGHzRgL5x3UXdmCmKqJqZ+u9KJ7xvRy
yZ09uKCNfXr7BFgU903QOtaJ4Ny6v4UCDo4u/vX2YWEHGsRchhtl2Xh0SPx4
bDsLIgGN4tSfZhxo/vl5cmsDlloYP1+2rAoeV49RnNQSrPjwbEx7xwtI6Oot
HNVJxzaBtrD/1keQ3qtrXdiegYWh05g3BPnIv6V/LIBThoFKU5v3trUQnLvZ
mebPRN9hkcHfP1uRdSSo89PHdgzdP+GvtpqD6rJUzTZ+J66vXrihb4wLk6wf
SyOf8qCTsWk0byYf7ZYX7w5s4CPJNPzoiz18+JnE7JV24WNWQeF1wT4+Qm9O
r3xjz0d0zzsh/S18MDNEfcvN+VjzfaXc4fV8aMbuPaGpx0fO4tmbp2jyEZL5
MPHrfD5ETn3zH5nDx4XYnN1SUnyUTs06tVaMj14Fzt+9Ezzox8aUThnhQUHx
xIMD/TwwuFdvNLXyIJx/5yw7nYdVJSeVDpjzYFnaEFF+mIsaj466i00cPImp
75tuxEaNAYcl3N4KJyHL8X83MpG1Njsk07Ae9gPl/klFtTg48TXyePY7uN2J
0HEMaYZEgebLDKM2BFlpzLW70QGvqMboKlEuavVzfFuC+Qgfdx4JWN8FaZkR
SR/zHmTKrJ9ku3UiPKZZOPXuW8beWrmNrUZdDOb6s/0v+osZhmvWF2m3UJgv
4XEjxO89xC+7NcyS70HQHsaCbxYUnto9U5WU4WOoPr7BdRoPx21NDRmFnYh7
IsyJ4XWAFpe4636gHZb2uVlVO1uhWZr7IqSwCS07c2/pRtTAJkR+qv1gGd4Z
eJ7Y9rgQHjGXsz7/eYE89e1K8/kvELpvXrLivkJ8nNembEL0kJsaeKE8swaK
6YVTr85nYrzSQauV3QpGiYt02hI2rr4+eM2X4uCX9g4dhS9cqP0OMk5x4uPR
uXeVH2wEKKd2PWFxBUg8GpcnWEKhbPnI3q02FERnbvvsak9hWL/FU2cnBS+1
TPMN5HqKVrd9ozGF/bO7a71XUCjaI/FNR50CLzXFXV6WQl2TQoq6OIVUvwXR
tr8EiFQ7YJHyUQDNiIlAeb4Av69mTXvWJMC0RTtO7i4XYMu/MSnf8gQ4acF3
zE8T4Hj+w07xWAFclb1cyv0F2LFkx4+d6wQ4MPTD+nQuHztPVlxOSeTBkA5a
0/WsExYHhFm2VWyUTruZFaLPwsDWwJmd9UxIXRTdUuvQgEdqP0Om2jXgoatr
sW81E+8HZ3dvNWZhMF/Sft9nNtRG30WdHO7ENbHzQa5Ez4GxnoeUVWg4mKj9
KvvZjaUnHQoMrHowXaKJ+YzTDIttJpZVcm2MOD/TwEJHPsPPsHGHl0EZ4oaD
RJ4f6kXYg+uvvKPeY8HaPZOqNt1Qjlz9iTpHoUDV0nq2CFnnJVeGFQd5GBqs
p4KPc7FxF31njRMHh00SKn9mtiM2wSZB1KcVXGG1uO7yJpy3/XqhIbUG+eLR
ia5e5ZD+wFI+8O0N3l2u8G6YUoDjjHLprOJCuGW89hn1L8VaXsTZJ2ZVqDRt
Mk843IAbVuIWoedaYPT5Iu2h1I6QH5UK3HUcNKWt8YnL5+JvgVpqYQQfEXY1
qwdqBXh859qLvGgKjUM9QV+FaawNrBEXV6URab1EInYpjcDPV5PnL6fxsjG0
zF6ThntEVuygOo3zn7bj+DwaXit9r4nPojHtVzmV8w+N1Ss27PYRojHn8ZSC
G8MUcnrFNrt1U3hzPK7OhZwjHjWoH15O4eNIaVpLLoVJP0pNL4WCeIYl2+kW
hdBAJeaK8xTKI+7owotCs53RogJzCnpz5z/+d1iAhaKqqf6mAtxckqBXOMDD
d//x39sfdSLm+iOOmA8b54MTptc1tKI1+m0QI5cJp2knIx54NUJSy/BhbVcT
tOLub12d1ALLHYab9Da2g8eYWmxWz4HK/og1dud5yJuceuROoQA5Gs0rN2zv
glzLsIrQ/F7o5d5em5fbhc/yWdVRGgU46NTembSLz0icuNu3d6iFMUdJuExo
oB2m7r4ZHfz3+MW105VM6YPcn9s1a/y7YOZ1McwvmMIWqi+43V6AkGvVhg4a
fJRfPXjHyJ+Lqq1LtnGWc2D77KG67KZ2rLz4ztHkQwt+trApsYRGTLhbRWx9
Wo3h/XH6Op/LIPpyq9+/PkV45Nlj5L77NSQ2Fo1HDRfjE3ak+u17iw/Lw98E
5dWioUrC8OQ1Jp6ucLc1sWVhcXNUf3U1G4LwOg2hnE6EKuPl7WEeWkT2qlTf
FEDWVtjJi+zHqLr4p/2LaKjY86V3bqchJXesiutBY/lVq77Xx2hsS3H+5H2U
RlawsVcPuc6qLhX64kwjwrHUusmeRugB7fqaLTSOLPCw52ykMS96RvQUIxr1
y6Z5Ga+kIZJq5BNN9GV6yz96guhpZZDz8iCip10ugWOSRE+nK2vnpRI9qfZ7
XBUneloZ/SjmwVcKB49m6Nm3UvilZ2dbkkxhirzU5dcrKUg4pSjf3y6A1t+q
vLFsHoLCuFs05TsRFaSpzk9vx8X6e07fpVtRsuzSL8qViWqTpWF5a5ow27/x
YbFUM5jsVR/GtFlQYknN8X3PxibJGxt36HLxwkfHRDaZD4PlHh+Kyfvwg5eu
v3i+G30vvasNw3txU9bdaqyP+P4/lqxdNsWM0pFjtxNFaEby2B/WiHYeozN8
VdZc5S4YZtmFHE3tR4n67T2/zXtByVy86dhIw1fG7Vq3M4UZmsxL9jsE8Elk
c5z/8MAZ5m4fmsdFKjve+p/EDlT3GsqVhbfht2/dItVpLbifp9bI5bzDNDeH
lvGVVTjjz7TxkizD92ZFaUnHIvBc9Rd/WlSEquOnvrFqSmGjni6cMFCJbWbD
pe9y3iF1RT/zi1ILVIIrmEuYbahQPzj0R5GD5X82RjU3cNFrs0PtSjvR58/O
VSvNiJ+Nz9HaMYdGwC2xVPe9NAy2ZH2JDKGRkv/O+1skjUN3b57XiKah4CK6
6y35PeVByUGJcBoPPywMlSL33W7gCj06SSM6OOmmrTeNEkZy25wDNOwM4u6M
ONAwL2Nm9W+j0WjTnTVsScPksupTWQaNtIBNE9arib8sXaMUQ3zo7KkCUfZC
4ktX/zbmK9IoX//zNS1NI74yQHfjFPLcoMvBq0DhA4+ReFqf8O1Yk2zTXAFm
sOTPmWnwILetXnhoPQc9NgGV7qlteJX4Mmf2uhZIyK26GCDPhFVp36iwNBOy
hwJ2d2q1YE9T1ZBOUhvspPKDy604eLv2fJ+LNQ8pMte86MsCvBRW6tKsJuuB
/Lz593sQ/+V33McvPRCbWVtVqdYK11cZk0PMVkZEtqOfhiKPkTfZP9uA5NnS
6hnC1mJ90By7NPXk8n44PNcJu/u+Gz0880SarIe6mFVFvy4F5W2lRV8UBWA9
6BOyT+ChY+kKS6ULnTjednvuYBMbic9muctHs3CgQM0sdZIJT+cr98I766C4
3sQ25E4FOit124zUSsHT2ZP91aoIVjyf1tbxYvSOKMsqR5VDeJWuiM+uGsQ4
2/TJnm7CiZP8gXlhrRhJ3rnUR4WNYhdJXYd1nXgI+4NXXvDwyirqudYlAb6+
WeDeWUaBylISkTen8bl19v0TQTSOni3qsIun0ZIT8CDrIY3hTkpnWgqN43dD
l8xJorE0rn7yUhyNf67syqOu0fhySVjE7gIN1x9DJbzTNG4u2hQUcpzGuL7t
onWHaIx9g8OMfYRHTZe7R3fQWMz9FvPDmobwukzlGRto9OV+ydZeR2Pm2+jC
9bo0puff3nBfg8brUuZhF6KPE1J3ojrFyXXqe6ZFIQX/2QdET5O8sobh7lWR
wceM75ohR925SNSRUlip24GVr/SkpK+yULGk9EfqkWbIvaj1zRM0oXCOiWdM
OhNLy0sDZXa0wrpkbNh4rB0Lmp7Mvna2E2v3lscySZ6YOvW4u4MHhcSwWwZ4
0oVrh9ruPvDrhWnLgu9BU7sRY1Q3ObP6NaTecs2OZfMYc2zeTW+0aWG4d/E3
fo1gI7u7fDRYux8SJl4P+4b7YJfbus6yqguJa0TVgiTJusx/rtlIODgSpHqs
rpSPZCmm75KZPExv89tuWkH6WopT1JigHfvMLbSfBbUiZdyw91NRE/pT6m0t
62tQcHtSYDntLa7NP/0hM70EQelnZayuFMFR1X31pFQpYm44XXU5UQG9PaOK
jyvq0LnvmKbO7GasexrVmEz6TckrG7M26w6kyP29WaLPxRZFDnuaHx/ZPw+e
GJSmYCwTmachQ/Z3yP77CjcauzuPRlyKorFMe33+kTQaOWmzbLMyybldN+BR
+ZTGrFODGx4+opGUMefR7EQauc0apvuJjyQp1MW8Jf4RoCx5zIT4RzAj+kQb
8Y/Qw6/CzhP/MDNxXm9J/EN//Ersgj3Ed2piEubY0mjjhFoqW9EQXC1RXEr8
Q9/enylG/MPHcfaMesKbks+ahfeVaQT9jbRqFKHRoe0yGnWHQlda16BvgQBr
M7M1q2X5+O554/uv052QXeGilzONjTMaoh0NJ1oheHzEsaKciaQJmSMBo014
di8+pcC7GYPjRUKi11i4fuHLZKRZB/QHNNIvXuPivUfF+N5xPgrELWceI+8x
5cvEhbfPu3FPr4g3La8XskuNOQWGfNTYBheXfH/DsNA9/3l5LsVwnrv8s4Eg
lxF5b29iuUUXciQ++T1v7EdgR3AN80gvtL1v1zX+oMH5VJf4M5vkqv0y25NT
BMj0/FossCa8evjWM82Di4lbwovHBzog6FvS/Y3VBtqmzPTazhaclHw7uepY
A+pSlh0R1FRhxxnOk8b7ZejIubP189li2Hz6+l5yYzGcz+s8/LCrDIXHftYZ
n6xC1j2XRTOMGsATffGH79KC4948+rBsO9IMfRdc3c+BP33XKl+aB9USz663
CgKEeglVVFygoDWvbWkt+f+i0ce6LgeSfcxdlLHvPo3qCN2KY89pyE8pP7M7
h4bG8l+mj57ROJly4codohdJJYlmu3uED1+PrL17i4bS0bl+E1do3OC0LDoV
SuNijd48CX8az4RmKOYeofGX01HlS3RXfqRRczPhlNIi7kNDOxrJJ15cMvqX
+IfCRMM6ExqKMt3F7eR9jt1eaVVEdHFohp48j+SPOTJGa7kjFOaVpP9JOERh
x24rym2zAPH0D0m7PTzM9yt4GHyUg01jPTY+7DbcvFt9975XCxZXHs8Js2Ki
p8s457oJEwYZcc++7WlBhlaeSk1DGxr9Y9kCfw6uy87cInaGh9U+e+Ml80kf
s56XKSyg8XNZ5PWSlz1QX3ynfLtoL4qyrpks29IK5+BASZHUVkb0BXp/2jcu
49W9kVce0ytQapFdkzyvD+drZ367ZdaPjMMLVC/97kZUXM3CWJLLCvUfhvbv
ofB3mDXjkrkAMr4OX8LKeLhssrC972knbiYt09rxkw2FiMHN8oUsRDuEjZou
akbwqxrBxOJ6ZF3cJqP7vgLjA9OzdbxLUUHF1+7mFyFZl6mRf70ED6sTrTs5
5RA5WfvueFINqOCfbntLm/C5x9fd40krquM2r/UyZYN7eTz7u2MncOORa0Mb
DzXtxhfkHgvwaGGy1rEeChVSLb6viK9fnWeaPUB48Nl2QWByOsmdGYs+snJJ
fvS2/Cb+guSk4hvhH4kucopjruoRXViV7nFwI34hvzBRNYv4RXJac+Esoovp
q7te/ZdLLE4uFltwikbDCVfjJuIX318f0blF/EJ3anDaUeIXI79OG+wnfnEb
j987Er+4UJYVbE78YnvSVbfzq0jPWSj+TZX0nV9je3L9ZWk4WZ9KX8+iIKPf
l6mpRiF3wf0tulV8iCrzfOJDuHC41+qpSHxxyJ9f/fIRC5d+r7kucZn0xx1l
QkESTGwV4hgElDGhu6X0Wb53K75nOLx5LMPG2RsBtTmxnUgT7M4p2MTHO5Nl
46KhFDpe5PECX3bBv93xi3FoLwTT9qn+mdWNCTnHCmX1N2jXPyoylMBjeD/P
LQlb3cKo6Ps6+3cyG1fGns4+tLYf6TH37uX+6cOBjnDLX8wurJOcck5WgUas
hqRFzVSSo/W6J6+28fFDsDkxbwEP28asM+M7OHBqiHisOtKO3ZZThJJvtuLF
Oj1X794mpBhOHfUcrsH4it0BMrpvUTGbFew+WoLw5tnpCm1FWDdf9/GtlaXI
YI9T9+MqYLYlWSfkQx3Gj+iOexFdxXEP++QOsOB5YXNTvlMHDh9QVfLdyMVO
+cbCqEt8eAcOvdikTkFJ+pRYyVySp6KZg9lk31y+K54fTCA+kWc9u4L4xN97
CU8H80hOdP6YOEL0MVHo6lhEeHJD/J+kuYQjNuv/1m27Q0NPsrI7+SYNsaRD
ef9cprGFHdsVcZb0GGPm9oV+JIesL7zddpiGm9nBl0kuZN73fhK8m8aZG5UK
PltpcKtjq7wtaGid3+quT3rL4JnFdr16NELEZx18TPLnpJz8y4ppNCqaEzUc
0yn4BTYvn1clwG3Tf5Vs1fnw0mNKOF/uRLpZ68wf8mzQkfv0h4Nb0Rjw8XVl
ExN+6Qr70qcy0bEvZGfg6WbkJJ8zTYtnQVyF1TVrawdWla8W64/n4rpSW9BD
cQH6B1WSzhkSjhZXRdbldUPl0vYtk4W9KPDNeeZvxUe4o3CDVtUbRrq5Z/HT
NIpBL5LMjs/OZQSdPDvH26YLnZb+eRtZ/bi3OHfzHt9etJutv28xQfLa9Q+H
m95Q8NjU77PnuQB6s3YdOLuLj/n6Fru1fLlIOFtreGCE8EOjzNyIaoOkrUs1
9rfgb+RmKdcLDfConCjv5FWhMeyeddyLMixZa/9tcWoxMnd/jun3KsZxdkqt
9uEymHsdibMIr4LH5Hvfnq0N6Fkgr6hFfNMgLLojXrkdZ0Z0rJ57cLCy92hx
xVwe7lYUvHijJkBSpV9s2jUKp095vFm+nobRzazoLnK+PWWfdGSSXBkxITee
S3TwcL4+L+ol6Y/tlByLcMT6/JHRUpIvvp1yrzn6gIa/jHJh5m2SS69RcSKk
nyw26BLzvUhjyOFH8/gZkqtGpCvunKCx3lAv0I70Wa9wORH1/UQn4u9fT91F
E52be0/a0OjWrXUfIrnztLbal3Mkd6bgj5uTDskrW2/o+qkQH7o8r3xsnMLg
8pLE1KMU3A7+cnAjvS70Y+JpMxceVuZ6uzqd5KBo9RuKIWjDavuwXEefFqzY
fyZJbysTsyTGopZbMhHadlYskPBVIzi315xwuy13T7tNMAc3TBMy94XyEJpO
pe0qEqB0bmrXwi4aT02WX2sr6MEziwA316m96BnyktYlOVXvIXPfzPutjF1G
G3LxgcsQFUrzlJGvwL4P9TtDlfvQevrh3mMb+8EKHp7vMNkNvw0/6/eTdej5
MWLoS3qpWJ39daFNAuz88OHcjyoeJvXFAi2zO8FZdKIhfYINjxw+va+EhcoK
j7SzWs0Q0bululu3Hhzp/JN+gxUImO3F8zxVCrdJ5+grQ0VQ8KmNffygBPkP
Ho0+6SnH8UiW1Or0Gniumv9VqK4J5bsfvF/0nPT05ZXdUhZsqJzaGXHfpROl
m4dCvTt5GC1zy87IEGCdxD+CX/0UQjlOh0fJ/jBV8jaq3iDn0uCsSQvZd3Pl
H+pTiB4EPqMTy8iUllMpmJ5NY3brHNcthC9SG9jMM6SPdHc9sSiNocFSmbir
cpXGy8krB2LP0Tg3L/vywgAaA+vfPaonPLWvUbgc4U70YcX+7exI+un+pts2
20nv8HztvHETjX27ghSXmNL43fVc12UNjVs3ncfGtUgevRQmtpn0Zk/Tkz95
bNLnXl8sLltIgZ3ydk1pLR8XSt9Yq1zgItDqkGX8lg7IruFXaqazYPy5Eb4R
zbDqeO0rmM6EseoL65YKJhTsG/5MPd4Kp7lIG5RlI/PDL53h+E7MXLPMfWgz
Hwtf3ptcSfKWfJFV38xXXfhoxuY0nuvFbqcY6Suy3ej96SRw0n6DgbT+jIE4
HmPcfDEjTL+F8arLbXIkhY3ZC48/otb1wz4sNkXjbx+e2mcFerZ0wbyyrvEc
8WWxEzpv706jsOJ1yqbdHXw4mwvVnV7Mw+Kzyq93cTnQSpXKY4+2QywkZ8O2
W60YuPcapR+awHFdx3v+swan7Uzc/l39Ft4PwtYX/y1BTYDB9h1UEQL0Vqum
ri3FGR39Ryl3K3BY+XXvya91MOiNY3dpNCN0zhLRBZ9YqL84w1HVuQN1tJE2
y5KL7w3l+79c5iNGNrUuiaxrgHFOQh/JcdGcla1tJA/S/OzuWXdpXEnqzerO
In3DwZY9LZ+GsuOCgH8IR2a+5j1qJLmi0PuO75JU4gtzM746kvvLth7b+5T0
lqNdzHBp0kO2WZxXuBFMI02kb2QR6SFPSmMes7xoGC6dZnjXlYbp46DMANLb
o66Z3ffcRmO+szPbxZLk2L4WscXGNLKND1xh6RO/2vTdIHYRjQIPG63s6YRf
9z66aD+lMD71a0V/tQCMklileQv5+GoWHrTwSifyGQbC+XPZcP2cm/sstBXe
/isag5qZkDv/ZfXef5igz7ZCM7AZa+Zvd952h4V/VucXFGzrgIjD18fhd7go
b5+vYSIhgO9l+c1bCMfK2rMXxb7sxtdczZGq170wV9KdECW5o9R7noFY6RuG
+oFZmxalUgxpI23L6Ke5DPvm8AGtzV1INt54T7itH9jcc0+O9NmowrSyv4Qf
g96/ZI2LKRy4NtejMkuAl4HTy3rt+Xhq9Fk90o+LTRcXe7/93oFblyP+JNBt
6H4XoT6X+Fr0PqcDUmEN+LLzqkcoVYW6WxzVVfllyC5v6bB7WozvHWMXpp4o
hluR9uoC7zI0P/PKr42ogmh+721nuwasPfqL6+bdgt3lCe8bVdphwdtg3+PJ
QdY+g9c/5/EgF7Fdf0BdgKHep/nM6xSkupxRRvjByHP47Et6g9IhWb4RyQu3
D3dlmJJ8ubPh71wVoouJt12tToQnnSMRI5YZND5FKe39TvghrxPzc0U8yRXj
QRfDCD8WPvuZ+p3wQ0i5WyOQ9JniHKnWuT40Kr/csmASfqSZKUYmEX4kzG+T
uED8SUZ7LDRgM8mhjDkbDpvTiJ8fHzdG+DHm0Mpv/Y8fUx8EfSD8+FJ9Rbj7
NwULz8fXrx6jcNFwuZDVTgHk1DbzF7jyYJR48LLhKQ7+jK15p0R4fG6T75E1
vi34oJKvOX0bE+NJWRMipIf0st5bb3VtwZoN22bKt7WhcEGwQCmEg/AHcrnL
zvEQeV/rmXqxAKz0o21XCD9uwur82sIePH76ZXcl4UfwVlmhUcIPafmuvZuT
WhlMy7Uv/w5wGYLy+23XFSqQnSZ/QkD4sSvjQcQ7wo+Ap15BLwk/hKoyFr0h
/EhvU5Vc4kJB1MAsMo/w483zrpLH1Tx4Zs0Y/EH4EZ/78uyuP2z8DUvpmQYW
DPdmTmosbYb6N++ZT/TqUcux9SsbqoCw4lb7F/6lWLHMvFz+exE+j5je0XlU
gslwq9GPveXQ9kmifJ/U4LWGRaphfRPUFi/0fEv4UU/n2mYSfmwd6pdd69qJ
KT+GS4W5PHjH3NYzzhSgXfDsR8QAhSjRe6LG9jQCm9nVhwg/NpzVb1Ug++5c
RkVuJNx4fDfh2HEyPx1kUDaEH2dFmJ6RhB//B6JZ5o8=
"], {{0, 10,
Rational[1, 100]}}, 1, {"Continuous", 1}, {"Discrete", 1}, 1, {
ResamplingMethod -> {"Interpolation", InterpolationOrder -> 1}}},
False, 10.2]["Values"];
ListLinePlot[{ecg1, ecg2}, PlotRange -> Full,
AxesLabel -> {"time (s)", "voltage (mV)"}, DataRange -> {0, 10},
ImageSize -> 400]
Out[1]=
2つの信号間の伸縮距離を計算する.
In[2]:=
WarpingDistance[ecg1, ecg2]
Out[2]=
時間伸縮対応に基づいて,両方の信号を伸縮させ,プロットする.
In[3]:=
{n1, n2} = WarpingCorrespondence[ecg1, ecg2];
In[4]:=
ListLinePlot[{ecg1[[n1]], ecg2[[n2]]}, PlotRange -> Full,
AxesLabel -> {"time (s)", "voltage (mV)"}, DataRange -> {0, 10},
ImageSize -> 400]
Out[4]=