Wavelet analysis of the blood pressure and pulse frequency measuements with Maple
\302\251 -2013, Prof. Gennady P. Chuiko and Irina A. Danishewska,
Petro Mohyla BlackSea Satate University, Mykolayiv, Ukraine
gp47@mail.ru
<Text-field style="Heading 2" layout="Heading 2">1. Introduction and measurement details</Text-field>
A significant part of medical signals, or observations, is non-stationary, discrete time sequences. Thus, the computer methods analysis, as well as refinement and compression, are very helpful as for the problems of recognition and detection of their key diagnostic features. We are going to illustrate here this statement with examples of very common, and even routine medical measurements of blood pressure as well as pulse rate and with possibilities of Maple.
The package of Discrete Wavelet transforms (DWT) within Maple 16 [1] was recently added as new research software just for such tasks. The practical testing of this package was additional goal of present study.
Results of measurments of the blood pressure and the pulse frequency of the patient were obtained by automatic blood pressure monitor "Micrilife BP 3AG1" (Microlife AG, Switzerland, http://www.microlife.com.ua). The technical parameters of this device are presented in table 1 in accordance with the accompanying documentation.
Table 1. Technical parameters of device
The results of observations are presented as data tables and vectors below
restart;
Data_s:=convert(Data_s,Vector[row], datatype=float[8]): whattype(Data_s);# systolic blood pressures (mmHg)
JkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiNJJHJvd0dGJw==
Data_d:=convert(Data_d,Vector[row], datatype=float[8]):# systolic blood pressures (mmHg)
Data_p:=convert(Data_p,Vector[row], datatype=float[8]):# pulse frequences (beats per minute - bpm)
<Text-field style="Heading 2" layout="Heading 2">2. Primary results and their statistical characteristics</Text-field>
Let us show the prymary results grafically (table 2)
gs:=plots[listplot](Data_s,view=[0..80,125..185],thickness=2, gridlines=true,tickmarks=[15,16],color=black, title=`fig.1.Systola`,labeldirections=[horizontal,vertical],labels=[`counts`,`Blood pressure, mmHg.`],font=[helvetica,14]):
gd:=plots[listplot](Data_d,view=[0..80,70..110],thickness=2, gridlines=true,tickmarks=[15,16],color=black, title=`fig.2.Diastola`,labeldirections=[horizontal,vertical],labels=[`counts`,`Blood pressure, mmHg.`],font=[helvetica,14]):
gp:=plots[listplot](Data_p,view=[0..80,65..105],thickness=2, gridlines=true,tickmarks=[15,16],color=black, title=`fig.3.Pilse`,labeldirections=[horizontal,vertical],labels=[`counts`,`Pulse, bpm.`],font=[helvetica,14]):
Table 2
The next table shows the statistical characteristics of results
Table 3
The results of table 3 as for standard deviations and accuracies are in certain disagreements with values declared by table 1.
Other visual characteristics of primary measurements are statistical frequency plots (table 4).
Table 4
The readers can decide themself how convincing look the graphs of table 4 for the solving of question of closeness of these distributions of probabilities to the normal law.
A significant positive correlation coefficient has been observed only between systolic pressure and pulse frequency.
Statistics:-Correlation(Data_d,Data_s);
Statistics:-Correlation(Data_d,Data_p);
Statistics:-Correlation(Data_p,Data_s);
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyMkITAlPUhcO2otSCEjOzcjXUJGOURCOTEzMjAzREUxNEY=
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyMkIjA3WmszKTQnMyIhIzo3I10zRkJCQ0REQTM4Mjc2MkZD
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyMkIjBySkA+MkZmJSEjOjcjXTNGREQ2NEIxMDNDMUZCNUQ=
However, this correlation also looks not too compelling (see the scatter plot below).
Statistics:-ScatterPlot(Data_p,Data_s,title = "fig.7 Systola-Pulse correlation", legend = "Point data",view=[76..88,140..170],symbol=solidcircle,symbolsize=18,color=black,gridlines=true,font=[helvetica,16],tickmarks=[12,10],labeldirections=[horizontal,vertical],labels=[`Pulse, bpm`,` Systolic pressure, mmHg`], titlefont=[helvetica,16]);
LSUlUExPVEc2LC0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJD8oISIiLSUpQk9VTkRTX1lHNiMkIiQhUiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJXFfISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJVNGISIiLSUpQ0hJTERSRU5HNiItJSdDVVJWRVNHNiY3XHA3JCQiJCF5ISIiJCIlKzkhIiI3JCQiJFMpISIiJCIlXTshIiI3JCQiJCspISIiJCIlXTkhIiI3JCQiJCspISIiJCIlKzshIiI3JCQiJCspISIiJCIlKzkhIiI3JCQiJF0oISIiJCIlNTohIiI3JCQiJCF5ISIiJCIlPzkhIiI3JCQiJDUpISIiJCIlPzohIiI3JCQiJGcpISIiJCIlUzshIiI3JCQiJEkpISIiJCIlPzohIiI3JCQiJGcoISIiJCIlXTkhIiI3JCQiJDUoISIiJCIlZzkhIiI3JCQiJCF6ISIiJCIlPzohIiI3JCQiJGcoISIiJCIlNTohIiI3JCQiJCF6ISIiJCIlPzkhIiI3JCQiJFMpISIiJCIlIWYiISIiNyQkIiQ/KSEiIiQiJT86ISIiNyQkIiQrKSEiIiQiJVM6ISIiNyQkIiRTKSEiIiQiJSFbIiEiIjckJCIkXSkhIiIkIiU/OiEiIjckJCIkKykhIiIkIiVdOSEiIjckJCIkPykhIiIkIiVdOSEiIjckJCIkKykhIiIkIiVdOSEiIjckJCIkUykhIiIkIiUrOiEiIjckJCIkIXohIiIkIiU/OiEiIjckJCIkcSkhIiIkIiUrOiEiIjckJCIkcSchIiIkIiVdOyEiIjckJCIkIXohIiIkIiVTOiEiIjckJCIkXSkhIiIkIiUrOyEiIjckJCIkIXohIiIkIiU1OSEiIjckJCIkZykhIiIkIiU/OiEiIjckJCIkZykhIiIkIiVxOiEiIjckJCIkZykhIiIkIiU1OiEiIjckJCIkPykhIiIkIiUrOiEiIjckJCIkSSkhIiIkIiU/OiEiIjckJCIkISopISIiJCIlNTshIiI3JCQiJHEpISIiJCIlZzshIiI3JCQiJDUpISIiJCIlNTohIiI3JCQiJFMpISIiJCIlUzohIiI3JCQiJCspISIiJCIlcTkhIiI3JCQiJGcpISIiJCIlIVsiISIiNyQkIiQheiEiIiQiJSFcIiEiIjckJCIkNSkhIiIkIiUrOCEiIjckJCIkIXohIiIkIiUhWyIhIiI3JCQiJD8oISIiJCIlSTghIiI3JCQiJCF6ISIiJCIlKzkhIiI3JCQiJFMpISIiJCIlPzohIiI3JCQiJCF6ISIiJCIlIVwiISIiNyQkIiQ1KCEiIiQiJSs6ISIiNyQkIiRxKCEiIiQiJWc5ISIiNyQkIiQheiEiIiQiJSs8ISIiNyQkIiRdKSEiIiQiJVM7ISIiNyQkIiRxKSEiIiQiJSFwIiEiIjckJCIkUykhIiIkIiVTOiEiIjckJCIkSSkhIiIkIiVxOiEiIjckJCIkUykhIiIkIiUhXCIhIiI3JCQiJD8qISIiJCIlSTwhIiI3JCQiJCEqKSEiIiQiJUk8ISIiNyQkIiRJKiEiIiQiJSFvIiEiIjckJCIkZykhIiIkIiVJPCEiIjckJCIkXSkhIiIkIiU/OiEiIjckJCIkISopISIiJCIlcTshIiI3JCQiJF0pISIiJCIlNTohIiI3JCQiJCF5ISIiJCIlXTkhIiI3JCQiJCspISIiJCIlPzkhIiI3JCQiJF0pISIiJCIlSTohIiI3JCQiJF0pISIiJCIlKzohIiI3JCQiJGcpISIiJCIlcTohIiI3JCQiJHEoISIiJCIlPzohIiI3JCQiJFMpISIiJCIlZzshIiI3JCQiJEkpISIiJCIlPzwhIiI3JCQiJEkpISIiJCIlUzwhIiI3JCQiJEkpISIiJCIlKzwhIiI3JCQiJEkpISIiJCIlXTshIiI3JCQiJCEpKSEiIiQiJXE5ISIiNyQkIiRJKSEiIiQiJSF5IiEiIjckJCIkIXohIiIkIiUrOyEiIjckJCIlPzUhIiIkIiVTOyEiIjckJCIkKyohIiIkIiU/PSEiIjckJCIkPykhIiIkIiUrPSEiIi0lJ0xFR0VOREc2Iy0lKV9UWVBFU0VURzYjLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1LUkjbXNHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2I1ErUG9pbnR+ZGF0YTYiLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuNiIvJSVzaXplR1EjMTI2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EmZmFsc2U2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdub3JtYWw2Ii0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUmU1RZTEVHNiMlJlBPSU5URy0lJVZJRVdHNiQ7JCIkZyghIiIkIiQhKSkhIiI7JCIlKzkhIiIkIiUrPCEiIi0mJSZfQVhJU0c2IyIiIjYnLSUrX0dSSURMSU5FU0c2IyUoREVGQVVMVEctJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIi0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Jy0lK19HUklETElORVNHNiMlKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiItJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUqQVhFU1RJQ0tTRzYkIiM3IiM1LSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRK1B1bHNlLH5icG02Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxNjYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1UTl+U3lzdG9saWN+cHJlc3N1cmUsfm1tSGc2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxNjYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzslK0hPUklaT05UQUxHJSlWRVJUSUNBTEctJSZUSVRMRUc2JS0lKV9UWVBFU0VURzYjLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY2LUkjbXNHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2I1FAZmlnLjd+U3lzdG9sYS1QdWxzZX5jb3JyZWxhdGlvbjYiLyUnZmFtaWx5R1EqaGVsdmV0aWNhNiIvJSVzaXplR1EjMTY2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EmZmFsc2U2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUwZm9udF9zdHlsZV9uYW1lR1ElVGV4dDYiLyUsbWF0aHZhcmlhbnRHUSdub3JtYWw2Ii0lJUZPTlRHNiQlKmhlbHZldGljYUciIzstJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUnU1lNQk9MRzYkJS1fU09MSURDSVJDTEVHIiM9LSUlRk9OVEc2JCUqaGVsdmV0aWNhRyIjOw==
The need to clean up the primary data from input noises is the main conclusion of this section. The package of Discrete Wavelet transforms (DWT) looks as suitable tool for this aim [1,2].
JSFH
<Text-field style="Heading 2" layout="Heading 2">3. VectorDWT and InverseVectorDWT procedures</Text-field>
Write down the two procedures for direct Discrete Wavelet transformation (VectorDWT) and for inverse transforms with recovering of the original signal (InverseVectorDWT), using samples of [1]
with(DiscreteTransforms);
NyhJOURpc2NyZXRlV2F2ZWxldFRyYW5zZm9ybUc2IkkxRm91cmllclRyYW5zZm9ybUdGJElASW52ZXJzZURpc2NyZXRlV2F2ZWxldFRyYW5zZm9ybUdGJEk4SW52ZXJzZUZvdXJpZXJUcmFuc2Zvcm1HRiRJNFdhdmVsZXRDb2VmZmljaWVudHNHRiRJLFdhdmVsZXRQbG90R0Yk
VectorDWT := proc (V, iters, w1, w2) local i, s, temp1, temp2; description `Direct wavelet transformation`;s := LinearAlgebra:-Dimensions(V); if modp(s, 2^iters) <> 0 then error "the length of V is not divisible by enough powers of 2" end if; temp1 := Vector(V); temp2 := Vector(V); for i to iters do DiscreteWaveletTransform(s/2^(i-1), temp1, w2, w1, storagetype = singlearray, temp2); ArrayTools:-Copy(temp2, temp1) end do; return temp1 end proc; # direct Discrete Wavelet transformation (VectorDWT)
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSpWZWN0b3JEV1RHRihmKjYmSSJWR0YoSSZpdGVyc0dGKEkjdzFHRihJI3cyR0YoNiZJImlHRihJInNHRihJJnRlbXAxR0YoSSZ0ZW1wMkdGKEYoRihDKD5GNi1fSS5MaW5lYXJBbGdlYnJhR0YlSStEaW1lbnNpb25zR0YoNiNGMEAkMC1JJW1vZHBHRiY2JEY2KSIiI0YxIiIhWVFXdGhlfmxlbmd0aH5vZn5Wfmlzfm5vdH5kaXZpc2libGV+Ynl+ZW5vdWdofnBvd2Vyc35vZn4yRig+RjctSSdWZWN0b3JHRiVGPz5GOEZLPyhGNSIiIkZPRjFJJXRydWVHRiZDJC1fSTNEaXNjcmV0ZVRyYW5zZm9ybXNHRidJOURpc2NyZXRlV2F2ZWxldFRyYW5zZm9ybUdGKDYoKiZGNkZPKUZGLCZGNUZPRk8hIiJGWkY3RjNGMi9JLHN0b3JhZ2V0eXBlR0YoSSxzaW5nbGVhcnJheUdGKEY4LV9JK0FycmF5VG9vbHNHRiVJJUNvcHlHRig2JEY4RjdPRjc2I0k+RGlyZWN0fndhdmVsZXR+dHJhbnNmb3JtYXRpb25HRihGKEYoNyNGLg==
InverseVectorDWT := proc (V, iters, w1, w2) local s, i, temp1, temp2;description `Inverse Wavelet Transformation`; s := LinearAlgebra:-Dimensions(V); if modp(s, 2^iters) <> 0 then error "the length of V is not divisible by enough powers of 2" end if; temp1 := Vector(V); temp2 := Vector(V); for i from iters by -1 to 1 do InverseDiscreteWaveletTransform(s/2^(i-1), temp1, w2, w1, storagetype = singlearray, temp2); ArrayTools:-Copy(temp2, temp1) end do; return temp1 end proc;# inverse Discrete Wavelet transformation (InverseVectorDWT)
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+STFJbnZlcnNlVmVjdG9yRFdUR0YoZio2JkkiVkdGKEkmaXRlcnNHRihJI3cxR0YoSSN3MkdGKDYmSSJzR0YoSSJpR0YoSSZ0ZW1wMUdGKEkmdGVtcDJHRihGKEYoQyg+RjUtX0kuTGluZWFyQWxnZWJyYUdGJUkrRGltZW5zaW9uc0dGKDYjRjBAJDAtSSVtb2RwR0YmNiRGNSkiIiNGMSIiIVlRV3RoZX5sZW5ndGh+b2Z+Vn5pc35ub3R+ZGl2aXNpYmxlfmJ5fmVub3VnaH5wb3dlcnN+b2Z+MkYoPkY3LUknVmVjdG9yR0YlRj8+RjhGSz8oRjZGMSEiIiIiIkkldHJ1ZUdGJkMkLV9JM0Rpc2NyZXRlVHJhbnNmb3Jtc0dGJ0lASW52ZXJzZURpc2NyZXRlV2F2ZWxldFRyYW5zZm9ybUdGKDYoKiZGNUZQKUZGLCZGNkZQRlBGT0ZPRjdGM0YyL0ksc3RvcmFnZXR5cGVHRihJLHNpbmdsZWFycmF5R0YoRjgtX0krQXJyYXlUb29sc0dGJUklQ29weUdGKDYkRjhGN09GNzYjST9JbnZlcnNlfldhdmVsZXR+VHJhbnNmb3JtYXRpb25HRihGKEYoNyNGLg==
The input parameters of both procedures (such as V, iters, w1, w2 etc.) will be explained below.
<Text-field style="Heading 2" layout="Heading 2">4. Wavelet analysis, de-noising , compressin and reconstruction of sygnals</Text-field>
Let us determine the number of scale levels for beginning. This number (SDlevels) does not exceed a maximal power of two, which is the divisor of counts number (80). Because LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUklbXN1cEdGJDYlLUkjbW5HRiQ2JFEiMkYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy1GIzYlLUYvNiRRIjNGJ0YyLyUnaXRhbGljR1EldHJ1ZUYnL0YzUSdpdGFsaWNGJy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictSSNtb0dGJDYtUSJ+RidGMi8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGSC8lKXN0cmV0Y2h5R0ZILyUqc3ltbWV0cmljR0ZILyUobGFyZ2VvcEdGSC8lLm1vdmFibGVsaW1pdHNHRkgvJSdhY2NlbnRHRkgvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZXRjI= and LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUklbXN1cEdGJDYlLUkjbW5HRiQ2JFEiMkYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy1GIzYlLUYvNiRRIjRGJ0YyLyUnaXRhbGljR1EldHJ1ZUYnL0YzUSdpdGFsaWNGJy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictSSNtb0dGJDYtUSJ+RidGMi8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGSC8lKXN0cmV0Y2h5R0ZILyUqc3ltbWV0cmljR0ZILyUobGFyZ2VvcEdGSC8lLm1vdmFibGVsaW1pdHNHRkgvJSdhY2NlbnRHRkgvJSdsc3BhY2VHUSYwLjBlbUYnLyUncnNwYWNlR0ZXRjI= both are the divisors of counts number (80) we can choose SDlevels=3 or SDlevels=4. If we select SDlevels=4, the highest detail scale, as well as the scale of the approximation would be relatively narrow (only 5 counts, because LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkmbWZyYWNHRiQ2KC1JI21uR0YkNiRRIzgwRicvJSxtYXRodmFyaWFudEdRJ25vcm1hbEYnLUYjNiUtSSVtc3VwR0YkNiUtRi82JFEiMkYnRjItRiM2JS1GLzYkUSI0RidGMi8lJ2l0YWxpY0dRJXRydWVGJy9GM1EnaXRhbGljRicvJTFzdXBlcnNjcmlwdHNoaWZ0R1EiMEYnRkJGRS8lLmxpbmV0aGlja25lc3NHUSIxRicvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGTy8lKWJldmVsbGVkR1EmZmFsc2VGJ0Yy =5). So it looks better to choose one level less
SDlevels:=3;
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSlTRGxldmVsc0dGKCIiJDcjRi4=
Now we can choose the wavelet family and the lenght (10) of both digital filters ( SDw1,SDw2) [1]:
SDw1,SDw2:=WaveletCoefficients(Symlet,10);
LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYqLUkjbWlHRiQ2JVElU0R3MUYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RKCZjb21tYTtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRjEvJSlzdHJldGNoeUdGPS8lKnN5bW1ldHJpY0dGPS8lKGxhcmdlb3BHRj0vJS5tb3ZhYmxlbGltaXRzR0Y9LyUnYWNjZW50R0Y9LyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdRLDAuMzMzMzMzM2VtRictRiw2JVElU0R3MkYnRi9GMi1GNjYtUSkmQXNzaWduO0YnRjlGOy9GP0Y9RkBGQkZERkZGSC9GS1EsMC4yNzc3Nzc4ZW1GJy9GTkZYLUkobWFjdGlvbkdGJDYlLUkobWZlbmNlZEdGJDYoLUYjNiYtSSdtdGFibGVHRiQ2Pi1JJG10ckdGJDYmLUkkbXRkR0YkNigtSSNtbkdGJDYkUTUwLjAxOTUzODg4MjczNTI0OTgyN0YnRjkvJSlyb3dhbGlnbkdRIUYnLyUsY29sdW1uYWxpZ25HRltwLyUrZ3JvdXBhbGlnbkdGW3AvJShyb3dzcGFuR1EiMUYnLyUrY29sdW1uc3BhbkdGYnBGaW9GXHBGXnAtRmBvNiYtRmNvNigtRmZvNiRRNTAuMDIxMTAxODM0MDI0Njg5MDQyRidGOUZpb0ZccEZecEZgcEZjcEZpb0ZccEZecC1GYG82Ji1GY282KC1GIzYlLUY2Ni1RKiZ1bWludXMwO0YnRjlGO0ZWRkBGQkZERkZGSC9GS1EsMC4yMjIyMjIyZW1GJy9GTkZmcS1GZm82JFE0MC4xNzUzMjgwODk5MDgwNTYyM0YnRjlGOUZpb0ZccEZecEZgcEZjcEZpb0ZccEZecC1GYG82Ji1GY282KC1GIzYlRmJxLUZmbzYkUTQwLjAxNjYwMjEwNTc2NDUxMDg1RidGOUY5RmlvRlxwRl5wRmBwRmNwRmlvRlxwRl5wLUZgbzYmLUZjbzYoLUZmbzYkUTIwLjYzMzk3ODk2MzQ1Njc5MkYnRjlGaW9GXHBGXnBGYHBGY3BGaW9GXHBGXnAtRmBvNiYtRmNvNigtRiM2JUZicS1GZm82JFEzMC43MjM0MDc2OTA0MDQwNDA3RidGOUY5RmlvRlxwRl5wRmBwRmNwRmlvRlxwRl5wLUZgbzYmLUZjbzYoLUZmbzYkUTQwLjE5OTM5NzUzMzk3Njg1NTU4RidGOUZpb0ZccEZecEZgcEZjcEZpb0ZccEZecC1GYG82Ji1GY282KC1GZm82JFE1MC4wMzkxMzQyNDkzMDIzMTM4NDRGJ0Y5RmlvRlxwRl5wRmBwRmNwRmlvRlxwRl5wLUZgbzYmLUZjbzYoLUZmbzYkUTQwLjAyOTUxOTQ5MDkyNTcwNjI2RidGOUZpb0ZccEZecEZgcEZjcEZpb0ZccEZecC1GYG82Ji1GY282KC1GIzYlRmJxLUZmbzYkUTUwLjAyNzMzMzA2ODM0NDk5ODc2OEYnRjlGOUZpb0ZccEZecEZgcEZjcEZpb0ZccEZecC8lJmFsaWduR1ElYXhpc0YnL0Zqb1EpYmFzZWxpbmVGJy9GXXBRJ2NlbnRlckYnL0ZfcFEnfGZybGVmdHxockYnLyUvYWxpZ25tZW50c2NvcGVHRjEvJSxjb2x1bW53aWR0aEdRJWF1dG9GJy8lJndpZHRoR0Zfdi8lK3Jvd3NwYWNpbmdHUSYxLjBleEYnLyUuY29sdW1uc3BhY2luZ0dRJjAuOGVtRicvJSlyb3dsaW5lc0dRJW5vbmVGJy8lLGNvbHVtbmxpbmVzR0Zqdi8lJmZyYW1lR0Zqdi8lLWZyYW1lc3BhY2luZ0dRLDAuNGVtfjAuNWV4RicvJSplcXVhbHJvd3NHRj0vJS1lcXVhbGNvbHVtbnNHRj0vJS1kaXNwbGF5c3R5bGVHRj0vJSVzaWRlR1EmcmlnaHRGJy8lMG1pbmxhYmVsc3BhY2luZ0dGZ3YvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXUYnLyUpcmVhZG9ubHlHRj1GOUY5L0krbXNlbWFudGljc0dGJFEqQ29sVmVjdG9yRicvJSVvcGVuR1EiW0YnLyUmY2xvc2VHUSJdRidGYngvJSthY3Rpb250eXBlR1EucnRhYmxlYWRkcmVzc0YnLyUpcnRhYmxlaWRHUSs0NDI5MjA3Mzc4RidGNS1GZW42JS1GaG42KC1GIzYmLUZdbzY+LUZgbzYmLUZjbzYoRl91RmlvRlxwRl5wRmBwRmNwRmlvRlxwRl5wRmJ0LUZgbzYmLUZjbzYoLUYjNiVGYnFGX3RGOUZpb0ZccEZecEZgcEZjcEZpb0ZccEZecEZkcy1GYG82Ji1GY282KEZhc0Zpb0ZccEZecEZgcEZjcEZpb0ZccEZecEZkci1GYG82Ji1GY282KEZhckZpb0ZccEZecEZgcEZjcEZpb0ZccEZecEZccS1GYG82Ji1GY282KC1GIzYlRmJxRmlwRjlGaW9GXHBGXnBGYHBGY3BGaW9GXHBGXnBGX29GYnVGZXVGZ3VGaXVGW3ZGXXZGYHZGYnZGZXZGaHZGW3dGXXdGX3dGYndGZHdGZndGaHdGW3hGXXhGYHhGOUY5RmJ4RmV4Rmh4RmJ4Rlt5L0ZfeVErNDQyOTIyMDYwMkYnRjk=
WaveletPlot(SDw1,SDw2,output=both);# wawelet digital filters
LSUlUExPVEc2Jy0lJ0NVUlZFU0c2JTdkXGw3JCQiIiEhIiIkIjBteU5UL2NfIiEjQTckJCIlREohIiYkIi5laEtTd2siISM/NyQkIiREJyEiJSQhMGw5UXc0bE8mISNCNyQkIiV2JCohIiYkIjBcJCo9J3kpM0giISNANyQkIiREIiEiJCQiL0BFMXc6W1AhIz83JCQiJkRjIiEiJiQiMCRIYGp2ZSNvJSEjQDckJCIldj0hIiUkIjAnKkghUng0QGUhI0A3JCQiJnY9IyEiJiQiMC9ESlQpb1hRISNANyQkIiNEISIjJCIvM24hPiUpKj1ZISNANyQkIiZEIkchIiYkIjAoR1YicGRzNyIhIz83JCQiJURKISIlJCIwLSEqZjFcQEwkISM/NyQkIiZ2ViQhIiYkIjBCLCgqUSdSIUcnISM/NyQkIiR2JCEiJCQiLzpqMC9EOSoqISM+NyQkIiZEMSUhIiYkIi9NV28mWzE+IiEjPTckJCIldlYhIiUkIjBVZ21LRHM/IiEjPjckJCImdm8lISImJCIwIltZXG00QDghIz43JCQiIiYhIiIkIjA0MGAoenAqWyIhIz43JCQiJkRKJiEiJiQiMFlzbzJWQE8iISM+NyQkIiVEYyEiJSQiMEY2YS5NSjciISM+NyQkIiZ2JGYhIiYkIjBvL0ttJzN2cSEjPzckJCIkRCchIiQkIjB0dT1OYVtzIyEjPzckJCImRGMnISImJCIwNmQrcyFHJTMiISM+NyQkIiV2byEiJSQiMHdKNSV5XjJKISM+NyQkIiZ2PSghIiYkIi8jSHFJLjtzJiEjPTckJCIjdiEiIyQiL3Y9Iz5FaCgqKSEjPTckJCImRCJ5ISImJCIwJXo9Q1g4bTchIz03JCQiJUQiKSEiJSQiMEtCYjtbM2wiISM9NyQkIiZ2VikhIiYkIjBDOjJac0c1IyEjPTckJCIkdikhIiQkIjBadm9GPG1lIyEjPTckJCImRDEqISImJCIwcHhdXSwmPkghIz03JCQiJXYkKiEiJSQiMCg+dThZQEFKISM9NyQkIiZ2byohIiYkIjBBbkFpa1I+JCEjPTckJCIjNSEiIiQiMDRtZmtEMDokISM9NyQkIidESjUhIiYkIjBeXSxDKnBKSyEjPTckJCImRDEiISIlJCIwZEJ5IWYmM1MkISM9NyQkIid2JDQiISImJCIwJ1s9SjpWMk8hIz03JCQiJUQ2ISIkJCIwI1EkZVd2SCVRISM9NyQkIidEYzYhIiYkIjA4TDkmKjNUIVEhIz03JCQiJnY9IiEiJSQiMEZTekZfc1wkISM9NyQkIid2PTchIiYkIjBuZlhbKSopKj0kISM9NyQkIiREIiEiIyQiMCMpek9QNWQnRyEjPTckJCInRCJHIiEiJiQiMCJSaUxOPnJCISM9NyQkIiZESiIhIiUkIjBtXm9DMSFwPSEjPTckJCIndlY4ISImJCIwWSNvIT1HZ0UiISM9NyQkIiV2OCEiJCQiMF8nPSlSQChlcyEjPjckJCInRDE5ISImJCIwInBFIVsoW2M3ISM9NyQkIiZ2ViIhIiUkIjAwQ1h0KltHRyEjPTckJCIndm85ISImJCIwPClvYC8mSDgmISM9NyQkIiM6ISIiJCIwKWZiWFUkbz0pISM9NyQkIidESjohIiYkIi8iW1pYKnpfNiEjOzckJCImRGMiISIlJCIwL2UvbXQyXiIhIzw3JCQiJ3YkZiIhIiYkIjAqUXgtKVElPj4hIzw3JCQiJUQ7ISIkJCIwJXokb1oiR25CISM8NyQkIidEYzshIiYkIjBCXlU5OXUjRyEjPDckJCImdm8iISIlJCIwO0pYSDo8SSQhIzw3JCQiJ3Y9PCEiJiQiMGZ1TkwpXHJQISM8NyQkIiR2IiEiIyQiMCQ+VFZwV1BVISM8NyQkIidEInkiISImJCIwKFErVEx2ZFohIzw3JCQiJkQiPSEiJSQiMDEkM1YhKVE2YCEjPDckJCIndlY9ISImJCIwd10qUVdEJSllISM8NyQkIiV2PSEiJCQiMDxZW0AqKTRZJyEjPDckJCInRDE+ISImJCIwRjFIJCpwLiNwISM8NyQkIiZ2JD4hIiUkIi8lZmlDMW9EKCEjOzckJCIndm8+ISImJCIvIlE+aiRbU3YhIzs3JCQiIz8hIiIkIjBiaFZGTiZleCEjPDckJCInREo/ISImJCIwWVdfJyk9YCh5ISM8NyQkIiZEMSMhIiUkIjAiXEYlcDtCI3ohIzw3JCQiJ3YkNCMhIiYkIjAlW0hdYjtzeSEjPDckJCIlREAhIiQkIjBSKlInKm85YXghIzw3JCQiJ0RjQCEiJiQiMGojPiJ6Mm16KCEjPDckJCImdj0jISIlJCIvIno5VnpbKnohIzs3JCQiJ3Y9QSEiJiQiMFQjKWVSMClIIykhIzw3JCQiJEQjISIjJCIwI3pkNm5GNyYpISM8NyQkIidEIkcjISImJCIwdi1nWy5FJCkpISM8NyQkIiZESiMhIiUkIjBCVTZ4X3g4KiEjPDckJCIndlZCISImJCIwYiJ6c0lULiYqISM8NyQkIiV2QiEiJCQiMDpxKmU2U28pKiEjPDckJCInRDFDISImJCIwbCYzYEkjcCgpKiEjPDckJCImdlYjISIlJCIwZlRPMC1KYiohIzw3JCQiJ3ZvQyEiJiQiMGRNJWZbaHcqKSEjPDckJCIjRCEiIiQiMEFycCtfdDopISM8NyQkIidESkQhIiYkIjA9YTx4U3ROKCEjPDckJCImRGMjISIlJCIwRVFYYEJEZCchIzw3JCQiJ3YkZiMhIiYkIi9KNHonZTdyJiEjOzckJCIlREUhIiQkIjBham1MPkghWyEjPDckJCInRGNFISImJCIwTTA9UktWdSQhIzw3JCQiJnZvIyEiJSQiMCMzJlEjPjtURCEjPDckJCIndj1GISImJCIweFFML0AnKVIiISM8NyQkIiR2IyEiIyQiL1phNz0tR0ohIzw3JCQiJ0QieSMhIiYkITBVeD9hKVxoJCohIz03JCQiJkQiRyEiJSQhL0V0MkEqcD4jISM7NyQkIid2VkchIiYkITAkKilwZiozOlskISM8NyQkIiV2RyEiJCQhMEg4X2wnel1ZISM8NyQkIidEMUghIiYkITBwWHMib3FHWyEjPDckJCImdiRIISIlJCEwck1rX0MsLCUhIzw3JCQiJ3ZvSCEiJiQhMDQkZSopb0xIRCEjPDckJCIjSSEiIiQhMCVIKFFpZDlRJCEjPTckJCInREpJISImJCIwSzVcXShbI1IjISM8NyQkIiZEMSQhIiUkIjBaSUFhSzlkJiEjPDckJCIndiQ0JCEiJiQiMHQiUmt3QDclKiEjPDckJCIlREohIiQkIjAlcDJeQVZ6OCEjOzckJCInRGNKISImJCIwJXBzL2I1MT0hIzs3JCQiJnY9JCEiJSQiMDk3KT4lKSlIQSMhIzs3JCQiJ3Y9SyEiJiQiME4iR0coPTlsIyEjOzckJCIkRCQhIiMkIjByalo+eSkpMyQhIzs3JCQiJ0QiRyQhIiYkIjBuXnA0aiRlTiEjOzckJCImREokISIlJCIuLVsoMzRrUyEjOTckJCIndlZMISImJCIwKj4hb013W2UlISM7NyQkIiV2TCEiJCQiMCQ9eDNxLUdeISM7NyQkIidEMU0hIiYkIjA+WSYqZiZSSGQhIzs3JCQiJnZWJCEiJSQiMD0qZjFpOCRRJyEjOzckJCIndm9NISImJCIwVFo+RUNyMyghIzs3JCQiI04hIiIkIjBqemQwWGwkeSEjOzckJCInREpOISImJCIwQyoqbylIRyRlKSEjOzckJCImRGMkISIlJCIvOC9XKD45TCohIzo3JCQiJ3YkZiQhIiYkIjBlIillJikqKiozNSEjOjckJCIlRE8hIiQkIjAkZWlkYnQmMyIhIzo3JCQiJ0RjTyEiJiQiLzMtJXlNI282ISM5NyQkIiZ2byQhIiUkIi9lZCZmJT1jNyEjOTckJCIndj1QISImJCIwcCpwIT5JV00iISM6NyQkIiR2JCEiIyQiLz03LkMtTDkhIzk3JCQiJ0QieSQhIiYkIjBtLD5dIzREOiEjOjckJCImRCJRISIlJCIwOEYyZDN0aCIhIzo3JCQiJ3ZWUSEiJiQiMCVwSnk9czU8ISM6NyQkIiV2USEiJCQiMHBBaCNlMi09ISM6NyQkIidEMVIhIiYkIjB4UjUhb0ZzPSEjOjckJCImdiRSISIlJCIvXkVyX0tAPiEjOTckJCIndm9SISImJCIwKnkmKm8qeWImPiEjOjckJCIjUyEiIiQiMDArNFNTVCg+ISM6NyQkIidESlMhIiYkIjB4QDlHeUUpPiEjOjckJCImRDElISIlJCIvVyZSbmpDKT4hIzk3JCQiJ3YkNCUhIiYkIjBobCFlIT4obz4hIzo3JCQiJURUISIkJCIvLmRLdWJWPiEjOTckJCInRGNUISImJCIwKDQ6eXdeOz4hIzo3JCQiJnY9JSEiJSQiMC0oemNqQigpPSEjOjckJCIndj1VISImJCIwOyRwTV0lUiY9ISM6NyQkIiREJSEiIyQiMCVcY1QtKW8iPSEjOjckJCInRCJHJSEiJiQiMF81X3dQL3giISM6NyQkIiZESiUhIiUkIi8xJ3p1Ul5yIiEjOTckJCIndlZWISImJCIwUSo9YFoiUmwiISM6NyQkIiV2ViEiJCQiMEE4JiopZXcnZSIhIzo3JCQiJ0QxVyEiJiQiMC5LY3pgY14iISM6NyQkIiZ2ViUhIiUkIjBbSDNiNDdXIiEjOjckJCIndm9XISImJCIwYCZvJEhuOk8iISM6NyQkIiNYISIiJCIwayg+R3JfeDchIzo3JCQiJ0RKWCEiJiQiMFAxPCN6ayM+IiEjOjckJCImRGMlISIlJCIwKTNDJlElMzE2ISM6NyQkIid2JGYlISImJCIwbXQzLC4kPTUhIzo3JCQiJURZISIkJCIwQnFDLHkkKUcqISM7NyQkIidEY1khIiYkIjApeUlRKlFZSCkhIzs3JCQiJnZvJSEiJSQiMGQudyR6bzFzISM7NyQkIid2PVohIiYkIjBFPCJbJUheMychIzs3JCQiJHYlISIjJCIwPkFPRnYjSFwhIzs3JCQiJ0QieSUhIiYkIjB1bF0qcEFEUCEjOzckJCImRCJbISIlJCIwPnIhSFEtMkQhIzs3JCQiJ3ZWWyEiJiQiMEMwUD1TX0QiISM7NyQkIiV2WyEiJCQiMCg9bS4tQGZcISM+NyQkIidEMVwhIiYkITAqeUU2YExcNSEjOzckJCImdiRcISIlJCEwaSEqXDspSDM+ISM7NyQkIid2b1whIiYkITA4QjwiZV1KRSEjOzckJCIjXSEiIiQhMFxJLyQqKmU1SyEjOzckJCInREpdISImJCEwbkAmZSxHOFAhIzs3JCQiJkQxJiEiJSQhMCQqXFcpRyUpW1QhIzs3JCQiJ3YkNCYhIiYkIS90I29vWXhZJSEjOjckJCIlRF4hIiQkITBYV0lzLykpbyUhIzs3JCQiJ0RjXiEiJiQhMC1icU93YClbISM7NyQkIiZ2PSYhIiUkITAlPUkjNElQMCYhIzs3JCQiJ3Y9XyEiJiQhMG1EVC49JSo9JiEjOzckJCIkRCYhIiMkITBGN15vKT4kSCYhIzs3JCQiJ0QiRyYhIiYkITBybjJpaGZJJiEjOzckJCImREomISIlJCEwOldrWU8lUl8hIzs3JCQiJ3ZWYCEiJiQhMFxbIVIoR102JiEjOzckJCIldmAhIiQkITBfQU8pRygqUlwhIzs3JCQiJ0QxYSEiJiQhMD5oZUBXKXlaISM7NyQkIiZ2ViYhIiUkITBCKlIiPTlbaiUhIzs3JCQiJ3ZvYSEiJiQhMHNOIj4kbyR5VyEjOzckJCIjYiEiIiQhMHYiXCJHZmhKJSEjOzckJCInREpiISImJCEwJjQiZSopKlxeVCEjOzckJCImRGMmISIlJCEwLCxbIzQvd1IhIzs3JCQiJ3YkZiYhIiYkITB5Nm9CPGAhUSEjOzckJCIlRGMhIiQkITAzKnkoW2syaiQhIzs3JCQiJ0RjYyEiJiQhMC1Ycmp1PVEkISM7NyQkIiZ2byYhIiUkITBIbCVSaUJpSSEjOzckJCIndj1kISImJCEwQ0pLSjwmM0YhIzs3JCQiJHYmISIjJCEwO1BTcWopPkIhIzs3JCQiJ0QieSYhIiYkITBwYnUtJ1swPiEjOzckJCImRCJlISIlJCEwel5hJTRrIlsiISM7NyQkIid2VmUhIiYkITAsJnAkeUw0LiIhIzs3JCQiJXZlISIkJCEwZVkuTEVnciYhIzw3JCQiJ0QxZiEiJiQhMCkpPk5UW1goPiEjPDckJCImdiRmISIlJCIwcyoqZXFtckIqISM9NyQkIid2b2YhIiYkIjA0RzVfdkpDJCEjPDckJCIjZyEiIiQiMDhwb1RVKFtcISM8NyQkIidESmchIiYkIjBAVTciSGM3ayEjPDckJCImRDEnISIlJCIuUDBuTWltKCEjOjckJCIndiQ0JyEiJiQiMG04dyZbbG0lKSEjPDckJCIlRGghIiQkIjBeX0pxMmEqKSkhIzw3JCQiJ0RjaCEiJiQiMCQpb3NJTmhAKiEjPDckJCImdj0nISIlJCIwPi9qaXo+VCohIzw3JCQiJ3Y9aSEiJiQiMDg4LjYkcDYmKiEjPDckJCIkRCchIiMkIjA9WnNWLlReKiEjPDckJCInRCJHJyEiJiQiMDdoJSo0J3pLIiohIzw3JCQiJkRKJyEiJSQiMCo9YlYoPm9XKSEjPDckJCIndlZqISImJCIweDJ0PylmSXYhIzw3JCQiJXZqISIkJCIwcm4iZidHVlcnISM8NyQkIidEMWshIiYkIjBabCl6PCoqW2MhIzw3JCQiJnZWJyEiJSQiMHgrZkgzMjomISM8NyQkIid2b2shIiYkIjAtazxyTDd4JSEjPDckJCIjbCEiIiQiMClmKWZgIVxSWCEjPDckJCInREpsISImJCIwbyYpbzRdJXBWISM8NyQkIiZEYychIiUkIjBaY2VaWTVBJSEjPDckJCIndiRmJyEiJiQiMG4pPSIqKilbMFUhIzw3JCQiJURtISIkJCIwMzRgZiUqKm9VISM8NyQkIidEY20hIiYkIjBtSXBsSEs1JSEjPDckJCImdm8nISIlJCIwWzZPQSR5QFAhIzw3JCQiJ3Y9biEiJiQiMDw8NDw+N0IkISM8NyQkIiR2JyEiIyQiMEJxRj1VbmkjISM8NyQkIidEInknISImJCIwNyM+I29EPS0jISM8NyQkIiZEIm8hIiUkIjBWJypmQW4+VyIhIzw3JCQiJ3ZWbyEiJiQiMCxyI1xUXSM0KSEjPTckJCIldm8hIiQkIjBSWVBxdzpoIiEjPTckJCInRDFwISImJCEwaVopZUBMKVskISM9NyQkIiZ2JHAhIiUkITA8Y3hXRS1EKCEjPTckJCIndm9wISImJCEwa25xelg2KyIhIzw3JCQiI3EhIiIkITAsUD1MUFM8IiEjPDckJCInREpxISImJCEwMy8iPjZiQTghIzw3JCQiJkQxKCEiJSQhMF9IKEhfJXlXIiEjPDckJCIndiQ0KCEiJiQhMGJwWXhdIjM6ISM8NyQkIiVEciEiJCQhL3RsPEJEOTohIzs3JCQiJ0RjciEiJiQhMHRqKmVQTiFbIiEjPDckJCImdj0oISIlJCEwVTBGNTxYUyIhIzw3JCQiJ3Y9cyEiJiQhMCR6XHQlW29JIiEjPDckJCIkRCghIiMkITA0UHI3Qmk9IiEjPDckJCInRCJHKCEiJiQhMEEhPkdpXSIpKiohIz03JCQiJkRKKCEiJSQhMCk0SWphbT13ISM9NyQkIid2VnQhIiYkIS9OcDhTQUBbISM8NyQkIiV2dCEiJCQhMG09WikpXCdlPCEjPTckJCInRDF1ISImJCIwKjRHJFFcR0olISM+NyQkIiZ2VighIiUkIi8sQUIrS1w8ISM8NyQkIid2b3UhIiYkIjBIVilvMm8qZSMhIz03JCQiI3YhIiIkIjBOeTMmKT0mKSpHISM9NyQkIidESnYhIiYkIjAkcFdLYm8xSSEjPTckJCImRGMoISIlJCIwOSRmKEcpZigpSCEjPTckJCIndiRmKCEiJiQiMDxUZWchR2lEISM9NyQkIiVEdyEiJCQiMCZIT0dHeV09ISM9NyQkIidEY3chIiYkIjBKXnc2TnlSIiEjPTckJCImdm8oISIlJCIwQ1V6N2kpejYhIz03JCQiJ3Y9eCEiJiQiMEsjeSJcdSoqMyIhIz03JCQiJHYoISIjJCIwZCFbIz5uZTgiISM9NyQkIidEInkoISImJCIwLWNqKjNHKioqKiEjPjckJCImRCJ5ISIlJCIwLXZJKTQ6Y3EhIz43JCQiJ3ZWeSEiJiQiMEFVIWZCLjtSISM+NyQkIiV2eSEiJCQiMC49SilcJlxgJiEjPzckJCInRDF6ISImJCEwI3lZVUpMJ3kiISM+NyQkIiZ2JHohIiUkITBoWCEpcFFaLyQhIz43JCQiJ3ZveiEiJiQhMEpkaS9UOyhRISM+NyQkIiMhKSEiIiQhMDNTNW4hemdUISM+NyQkIidESiEpISImJCEwLyZbLkJlTFIhIz43JCQiJkQxKSEiJSQhMGV3KillOmlNJCEjPjckJCIndiQ0KSEiJiQhMDpKJz0iPS46IyEjPjckJCIlRCIpISIkJCEwOjBzJkgpeUgmISM/NyQkIidEYyIpISImJCIwNio9QjloblYhIz83JCQiJnY9KSEiJSQiMDgqbyY9XHEjeSEjPzckJCIndj0jKSEiJiQiMFctKG9kVmciKSEjPzckJCIkRCkhIiMkIjBubTN0JXosXiEjPzckJCInRCJHKSEiJiQiMEA0bS9TJTNNISM/NyQkIiZESikhIiUkIjA3ZUcpW25QSyEjPzckJCIndlYkKSEiJiQiMCR6KGUocCE9KD4hIz83JCQiJXYkKSEiJCQiMF9SKHp3KnopPiEjQDckJCInRDElKSEiJiQhME9IbjZSbSd6ISNANyQkIiZ2VikhIiUkITBKIz0lKm9oVjYhIz83JCQiJ3ZvJSkhIiYkITBkSykpeSpmaCIqISNANyQkIiMmKSEiIiQhMD07K2ckKTRdIiEjQDckJCInREomKSEiJiQiMHM4Mjt4IT4+ISNANyQkIiZEYykhIiUkIjBQZCEqeix1TCIhI0A3JCQiJ3YkZikhIiYkIjBsNzVwIypvOCohI0E3JCQiJUQnKSEiJCQiMCo9J2YmUUNUZCEjQjckJCInRGMnKSEiJiQhMFcoUXdDIUgnRyEjQTckJCImdm8pISIlJCEwSDdmIlJfdkkhI0I3JCQiJ3Y9KCkhIiYkIjBKZzYkKUd4JUchI0ItJSdMRUdFTkRHNiMtJSlfVFlQRVNFVEc2Iy1JJW1yb3dHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NS1JI21zRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNiNRL0ZhdGhlcn5XYXZlbGV0NiIvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW42Ii8lJXNpemVHUSMxMjYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSZmYWxzZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ25vcm1hbDYiLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIzUhIiItJSdDVVJWRVNHNiU3ZFxsNyQkIiIhISIiJCIwV0IjPiRwMDQiISNBNyQkIiVESiEiJiQiMEZRPChmIXk8IiEjQTckJCIkRCchIiUkITBWKHBneUBPUSEjQjckJCIldiQqISImJCIwUUJGd055QSohI0E3JCQiJEQiISIkJCIwQChlaiRbJHpFISNANyQkIiZEYyIhIiYkIi9uaDMsS1pMISM/NyQkIiV2PSEiJSQiMEN0VVJ4NjslISNANyQkIiZ2PSMhIiYkIjBmd2RvbiFcRiEjQDckJCIjRCEiIyQiMGFuRGhgPUkkISNBNyQkIiZEIkchIiYkIjBeVXg4TiJlISkhI0A3JCQiJURKISIlJCIwWyVmNnInPlEjISM/NyQkIiZ2ViQhIiYkIjAlPlVFVl0qWyUhIz83JCQiJHYkISIkJCIwIyllQ1dWcjMoISM/NyQkIiZEMSUhIiYkIjAqKm9UMyJHNiYpISM/NyQkIiV2ViEiJSQiMEtDRzF6KEgnKSEjPzckJCImdm8lISImJCIwaSdSKlwieVYlKiEjPzckJCIiJiEiIiQiMDxbTiI9IVwxIiEjPjckJCImREomISImJCIwZSJSbCUqPlAoKiEjPzckJCIlRGMhIiUkIi9oI2VVZCdHISkhIz43JCQiJnYkZiEiJiQiMEMlXCMpUmVkXSEjPzckJCIkRCchIiQkIi5oaDFZeSU+ISM9NyQkIiZEYychIiYkIjBMRj96OzR2KCEjPzckJCIldm8hIiUkIjBNJilcITNSQEEhIz43JCQiJnY9KCEiJiQiMHRLaVFhKzQlISM+NyQkIiN2ISIjJCIwTyg9ISkzYDtrISM+NyQkIiZEInkhIiYkIjBjSngpeikzMCohIz43JCQiJUQiKSEiJSQiME4kenchKjQhPSIhIz03JCQiJnZWKSEiJiQiMHlwZilmQS46ISM9NyQkIiR2KSEiJCQiMFpAIzMiRyFcPSEjPTckJCImRDEqISImJCIwJSl5ZVYpKXAzIyEjPTckJCIldiQqISIlJCIweSdIMW4qPUIjISM9NyQkIiZ2byohIiYkIjBlN294Jz0kRyMhIz03JCQiIzUhIiIkIjBeIylRQTc4RCMhIz03JCQiJ0RKNSEiJiQiMHYoPioqR0Y0QiEjPTckJCImRDEiISIlJCIvYTctNlBKQyEjPDckJCIndiQ0IiEiJiQiMD5cQiY+enJEISM9NyQkIiVENiEiJCQiL3chSCQ+I3BzIyEjPDckJCInRGM2ISImJCIwJiopMzo0NSVwIyEjPTckJCImdj0iISIlJCIwIykpeidcMydvQyEjPTckJCIndj03ISImJCIvKT1zI3phZkEhIzw3JCQiJEQiISIjJCIwdVJ4MFlnLyMhIz03JCQiJ0QiRyIhIiYkIjB1Jj4kKkhFTTshIz03JCQiJkRKIiEiJSQiMFBGJHA2VGM2ISM9NyQkIid2VjghIiYkIi8vci07TmtjISM9NyQkIiV2OCEiJCQhME52NCY0VGY6ISM/NyQkIidEMTkhIiYkIjBVdCNSRjZqRCEjPjckJCImdlYiISIlJCIwWFV0QzY2UCIhIz03JCQiJ3ZvOSEiJiQiMD0jKjQ2YHEmSCEjPTckJCIjOiEiIiQiLyZRIm8wQFxdISM8NyQkIidESjohIiYkIjBPbilbV1ExdiEjPTckJCImRGMiISIlJCIweEp3QT4lPjUhIzw3JCQiJ3YkZiIhIiYkIjBQNFhiZFJMIiEjPDckJCIlRDshIiQkIjAiZjtfdWF4OyEjPDckJCInRGM7ISImJCIwKSk9XV0xRic+ISM8NyQkIiZ2byIhIiUkIjByL19xJW8jPiMhIzw3JCQiJ3Y9PCEiJiQiMCg0XipRIWUoUSMhIzw3JCQiJHYiISIjJCIwOHBiIyozX2EjISM8NyQkIidEInkiISImJCIwUC9nKGVaPUYhIzw3JCQiJkQiPSEiJSQiMDc+Zk9RbyFIISM8NyQkIid2Vj0hIiYkIjBQT1c9aEUyJCEjPDckJCIldj0hIiQkIjBiJEgncFtUQSQhIzw3JCQiJ0QxPiEiJiQiL3p6UGMydEwhIzs3JCQiJnYkPiEiJSQiMDZady8jSC9OISM8NyQkIid2bz4hIiYkIjBWOillbFJvTyEjPDckJCIjPyEiIiQiL2xjNF1qWlEhIzs3JCQiJ0RKPyEiJiQiMHpsbSo+TigpUSEjPDckJCImRDEjISIlJCIwUms2ZUkpSFEhIzw3JCQiJ3YkNCMhIiYkIjBUTFQmPjojbyQhIzw3JCQiJURAISIkJCIwQHlMaGUqcE0hIzw3JCQiJ0RjQCEiJiQiMDwiKikqeVU0XyQhIzw3JCQiJnY9IyEiJSQiMC9FMUA2eCNRISM8NyQkIid2PUEhIiYkIjB1bDUrQT87JSEjPDckJCIkRCMhIiMkIjBMdF9vISkpUlghIzw3JCQiJ0QiRyMhIiYkIjBjblAtUDwuJiEjPDckJCImREojISIlJCIwYnRwQDRJXiYhIzw3JCQiJ3ZWQiEiJiQiME5pOEI5JCozJyEjPDckJCIldkIhIiQkIjA6P25TNSpHbSEjPDckJCInRDFDISImJCIwXCVHTmkzVWohIzw3JCQiJnZWIyEiJSQiMFgnKWYmUWNpXyEjPDckJCIndm9DISImJCIwZCp5Y1U/L08hIzw3JCQiI0QhIiIkIjBbJ3pMd1JtOCEjPDckJCInREpEISImJCEwaUs8cT9CKyIhIzw3JCQiJkRjIyEiJSQhMHJLQSgpZlxbJCEjPDckJCIndiRmIyEiJiQhMFVtTWZUJSpHJyEjPDckJCIlREUhIiQkITBSVSVvSzZRJCohIzw3JCQiJ0RjRSEiJiQhMFFmNnFQLkUiISM7NyQkIiZ2byMhIiUkITBqcDdGNCE0OyEjOzckJCIndj1GISImJCEwInkleWNYSCY+ISM7NyQkIiR2IyEiIyQhMD09LURsSEgjISM7NyQkIidEInkjISImJCEwPENbKFFWdkUhIzs3JCQiJkQiRyEiJSQhMFBhQVF4czIkISM7NyQkIid2VkchIiYkITBbVCplKSlcJFwkISM7NyQkIiV2RyEiJCQhMG9ydSc0bS9SISM7NyQkIidEMUghIiYkITBIIXl3V1x3VCEjOzckJCImdiRIISIlJCEwOiwnKlJLakklISM7NyQkIid2b0ghIiYkITA0JlJNNCllTiUhIzs3JCQiI0khIiIkITApKmUiejEzOlYhIzs3JCQiJ0RKSSEiJiQhMHNbKVsoPWM9JSEjOzckJCImRDEkISIlJCEwX2lEKDQmKSopUiEjOzckJCIndiQ0JCEiJiQhMEU5ciRwVSRwJCEjOzckJCIlREohIiQkITAzMzUjeTVATCEjOzckJCInRGNKISImJCEva3IyVElPSSEjOjckJCImdj0kISIlJCEwJEhGQlEoSCRHISM7NyQkIid2PUshIiYkITBSJjNmYGpYRSEjOzckJCIkRCQhIiMkITBGU2ZGSip5QyEjOzckJCInRCJHJCEiJiQhL2ooUixdTEcjISM6NyQkIiZESiQhIiUkITBeXV1maXIuIyEjOzckJCIndlZMISImJCEvJ2U8WzcnKXoiISM6NyQkIiV2TCEiJCQhMDRoLikzJ2ZgIiEjOzckJCInRDFNISImJCEwKnp6Xzg4KzYhIzs3JCQiJnZWJCEiJSQhL18iXF8hcCo0JiEjOzckJCIndm9NISImJCIwVFIqeiFRPFQjISM8NyQkIiNOISIiJCIwMkB3JlJ5UDYhIzs3JCQiJ0RKTiEiJiQiMCd5d3B6UHA+ISM7NyQkIiZEYyQhIiUkIi8jW3kjeldwRiEjOjckJCIndiRmJCEiJiQiMDknZmQoKj52TiEjOzckJCIlRE8hIiQkIi9kaiozNGpSJSEjOjckJCInRGNPISImJCIwRiVSWCxwOmIhIzs3JCQiJnZvJCEiJSQiLi9tbCpHQXAhIzk3JCQiJ3Y9UCEiJiQiMHdfMVE8YFEpISM7NyQkIiR2JCEiIyQiMFhTSi16VCIqKiEjOzckJCInRCJ5JCEiJiQiMHc9QVFUKmY2ISM6NyQkIiZEIlEhIiUkIjAubSc0VFNJOCEjOjckJCIndlZRISImJCIvUSMqZTEyNzohIzk3JCQiJXZRISIkJCIvI2ViLDIwcCIhIzk3JCQiJ0QxUiEiJiQiMEVNJFs7ZyN5IiEjOjckJCImdiRSISIlJCIwTmk6ZWo3eiIhIzo3JCQiJ3ZvUiEiJiQiMC5yJzN5X088ISM6NyQkIiNTISIiJCIwMF5eIVsjemgiISM6NyQkIidESlMhIiYkIjAuJ3lHK2ElWyIhIzo3JCQiJkQxJSEiJSQiMCNINXdreU44ISM6NyQkIid2JDQlISImJCIwSi5jRzY0OiIhIzo3JCQiJURUISIkJCIwWEUpPVJiYiQqISM7NyQkIidEY1QhIiYkIjBoLGM7aVN6JyEjOzckJCImdj0lISIlJCIvZjAxI1soR1EhIzo3JCQiJ3Y9VSEiJiQiLyYqR11QUShcKCEjOzckJCIkRCUhIiMkITA0blQwMipcQyEjOzckJCInRCJHJSEiJiQhMDdUUHVtYjUnISM7NyQkIiZESiUhIiUkITAqPSdSQj88KioqISM7NyQkIid2VlYhIiYkITAycV8rW1hUIiEjOjckJCIldlYhIiQkITBtUnIhKUhdJD0hIzo3JCQiJ0QxVyEiJiQhMFpfIlxdZUZAISM6NyQkIiZ2ViUhIiUkITBkNyM9OEQlSCMhIzo3JCQiJ3ZvVyEiJiQhMDxHO2ddVFAjISM6NyQkIiNYISIiJCEwWS1qIXBcakIhIzo3JCQiJ0RKWCEiJiQhMEUpKXlJLyo+QiEjOjckJCImRGMlISIlJCEwKFIubDFPWkEhIzo3JCQiJ3YkZiUhIiYkITBOKkczWSNHNiMhIzo3JCQiJURZISIkJCEwPCg9RWMtRz4hIzo3JCQiJ0RjWSEiJiQhMD5OW28xJz08ISM6NyQkIiZ2byUhIiUkIS8iUmtbJVIkWyIhIzk3JCQiJ3Y9WiEiJiQhMHhPY3JDKFE3ISM6NyQkIiR2JSEiIyQhMDsxL3ZzZCUpKiEjOzckJCInRCJ5JSEiJiQhMD5HZV5bTXUnISM7NyQkIiZEIlshIiUkITBOKkdhOVN1SyEjOzckJCIndlZbISImJCIwaC4pWz91b1whIzw3JCQiJXZbISIkJCIwenlHJSp6M1MlISM7NyQkIidEMVwhIiYkIjBKajEnNFI0dCEjOzckJCImdiRcISIlJCIwNDh0W2ZuQSohIzs3JCQiJ3ZvXCEiJiQiL14hcHNoPzAiISM5NyQkIiNdISIiJCIwOjc1R1pXNiIhIzo3JCQiJ0RKXSEiJiQiME1VKltQKipbNiEjOjckJCImRDEmISIlJCIwJj0sUUc/aDYhIzo3JCQiJ3YkNCYhIiYkIi92RFpzTUI2ISM5NyQkIiVEXiEiJCQiMFYzLzFlay8iISM6NyQkIidEY14hIiYkIjBSIXo5TypScSohIzs3JCQiJnY9JiEiJSQiMCVHNzwneU4kKikhIzs3JCQiJ3Y9XyEiJiQiMClmZiUqKVs4OykhIzs3JCQiJEQmISIjJCIwPnAlRzYuIlIoISM7NyQkIidEIkcmISImJCIvJz4pKj5yKmVpISM6NyQkIiZESiYhIiUkIjAnW1ImWyx3JlshIzs3JCQiJ3ZWYCEiJiQiL2opM0ZcKilHJCEjOjckJCIldmAhIiQkIjBNPiRvU0VAOyEjOzckJCInRDFhISImJCIwOXhUTFJuJ1EhIzw3JCQiJnZWJiEiJSQhMCIzcFwqejI3JSEjPDckJCIndm9hISImJCEwS2RTXjpObCohIzw3JCQiI2IhIiIkITBMRD9mI1FZNyEjOzckJCInREpiISImJCEvJD5zNy1dUyIhIzo3JCQiJkRjJiEiJSQhMDxiJmYxKVxaIiEjOzckJCIndiRmJiEiJiQhMHdnXSFvY0c4ISM7NyQkIiVEYyEiJCQhMG1XLipIOEE1ISM7NyQkIidEY2MhIiYkITBFY1FUKSkzLikhIzw3JCQiJnZvJiEiJSQhMGNzJkdwbjJtISM8NyQkIid2PWQhIiYkITBWTC5jNXNnJiEjPDckJCIkdiYhIiMkIS9PKFEnRyE9MSYhIzs3JCQiJ0QieSYhIiYkITBePDtwek5JJCEjPDckJCImRCJlISIlJCEwIj0hUktVUHYmISM9NyQkIid2VmUhIiYkIjBLZkRtNFpdIyEjPDckJCIldmUhIiQkIjApZVknejokPWUhIzw3JCQiJ0QxZiEiJiQiMDQ1I3knKSlSKXohIzw3JCQiJnYkZiEiJSQiMExaMylmTyUpKikhIzw3JCQiJ3ZvZiEiJiQiMHk5TiZ5Tm0kKiEjPDckJCIjZyEiIiQiMHlPWl1kUjAqISM8NyQkIidESmchIiYkIjBKbCVwRiRwUCkhIzw3JCQiJkQxJyEiJSQiLzhBXzM1UXUhIzs3JCQiJ3YkNCchIiYkIjB0IkhRVWBBZiEjPDckJCIlRGghIiQkIjA3NXNbJ2UpKVIhIzw3JCQiJ0RjaCEiJiQiLzc6ci06W0IhIzs3JCQiJnY9JyEiJSQiMEBOLkYrY3YqISM9NyQkIid2PWkhIiYkITBUN091ZzRCIyEjPTckJCIkRCchIiMkITA4ZmwoKiozTjchIzw3JCQiJ0QiRychIiYkITB3KT0oUiZvUkQhIzw3JCQiJkRKJyEiJSQhMCJbSEdlNWVTISM8NyQkIid2VmohIiYkITBKVjwhZlU9YyEjPDckJCIldmohIiQkITBOMTFgNW88KCEjPDckJCInRDFrISImJCEvWHA2JWYnKUcpISM7NyQkIiZ2VichIiUkIS9hI1JXekAmKikhIzs3JCQiJ3ZvayEiJiQhMCU+JlFJS3ZLKiEjPDckJCIjbCEiIiQhL1pIJVF5X1IqISM7NyQkIidESmwhIiYkITAxMywrUTVJKiEjPDckJCImRGMnISIlJCEwJWVgWSV5I28hKiEjPDckJCIndiRmJyEiJiQhMDYiKkhBX1RmKSEjPDckJCIlRG0hIiQkITA8KW9JdihcI3ohIzw3JCQiJ0RjbSEiJiQhL2w+OWlFN3MhIzs3JCQiJnZvJyEiJSQhMDlRMEIvNFgnISM8NyQkIid2PW4hIiYkITA4Tig+O0xxYyEjPDckJCIkdichIiMkITBLKVwjWyYqRyhbISM8NyQkIidEInknISImJCEvO2ZnJWZLKVEhIzs3JCQiJkQibyEiJSQhL1hWalZmaEYhIzs3JCQiJ3ZWbyEiJiQhMFlfWEQuX2EiISM8NyQkIiV2byEiJCQhMGA3J0h1UFFHISM9NyQkIidEMXAhIiYkIjBHN25jOVd3JyEjPTckJCImdiRwISIlJCIwJj54U2VYTjghIzw3JCQiJ3ZvcCEiJiQiMGVbVHBqKjM9ISM8NyQkIiNxISIiJCIwJCkqKSp5eEsiMyMhIzw3JCQiJ0RKcSEiJiQiMEc4TFM5XkUjISM8NyQkIiZEMSghIiUkIjBISUZeQSV5QiEjPDckJCIndiQ0KCEiJiQiMD52JSpHPW1MIyEjPDckJCIlRHIhIiQkIjAyclRWKj11QCEjPDckJCInRGNyISImJCIwY11RWyl6Qz8hIzw3JCQiJnY9KCEiJSQiMHJfQk84QSk9ISM8NyQkIid2PXMhIiYkIjBIJD4lKm8xVTwhIzw3JCQiJEQoISIjJCIvTlgkcCJmMDshIzs3JCQiJ0QiRyghIiYkIjBtMjNmZy5PIiEjPDckJCImREooISIlJCIwIj1BUEdpSjUhIzw3JCQiJ3ZWdCEiJiQiMEJWOVc8a2AnISM9NyQkIiV2dCEiJCQiMFxTYlQ4I1JDISM9NyQkIidEMXUhIiYkITAjWyNvJykqeiM+JiEjPjckJCImdlYoISIlJCEwJDNMNlhbRUIhIz03JCQiJ3ZvdSEiJiQhMCQqW00ubmdfJCEjPTckJCIjdiEiIiQhMDAkPnQ4IypRUyEjPTckJCInREp2ISImJCEwLDQzdjxqQSUhIz03JCQiJkRjKCEiJSQhMChcR2kpek0+JSEjPTckJCIndiRmKCEiJiQhMF4oXC4sLiVmJCEjPTckJCIlRHchIiQkITBMYz1teCcqZSMhIz03JCQiJ0RjdyEiJiQhLzNkdSxVXz4hIzw3JCQiJnZvKCEiJSQhMDtKOUIjPl07ISM9NyQkIid2PXghIiYkITAyPiRRRjJEOiEjPTckJCIkdighIiMkITB5TiopKj0oKillIiEjPTckJCInRCJ5KCEiJiQhMDwsaiJlISkpUiIhIz03JCQiJkQieSEiJSQhMGFxYSZbKjMoKSohIz43JCQiJ3ZWeSEiJiQhL01MeEY7eWEhIz03JCQiJXZ5ISIkJCEvPSdwLCUpR3UoISM+NyQkIidEMXohIiYkIjAnPUYoKUciKilcIyEjPjckJCImdiR6ISIlJCIwRiIzKVsvJGZVISM+NyQkIid2b3ohIiYkIjAqUUImW2pnVCYhIz43JCQiIyEpISIiJCIwTz4lXG5iP2UhIz43JCQiJ0RKISkhIiYkIi94OTJOci1iISM9NyQkIiZEMSkhIiUkIjAqZiNSV1U1byUhIz43JCQiJ3YkNCkhIiYkIjBtZ1YhUTQzSSEjPjckJCIlRCIpISIkJCIwKmV3UkNDNnUhIz83JCQiJ0RjIikhIiYkITBvcWRpeik0aCEjPzckJCImdj0pISIlJCEwJlFVSDQkXDQiISM+NyQkIid2PSMpISImJCEweCVII2ZvOjkiISM+NyQkIiREKSEiIyQhMHUqZWlMJHA4KCEjPzckJCInRCJHKSEiJiQhMDktbmopM29aISM/NyQkIiZESikhIiUkIS9gOmFUP0hYISM+NyQkIid2ViQpISImJCEvIz5AVXQkZUYhIz43JCQiJXYkKSEiJCQhMHVvOnFBNXkjISNANyQkIidEMSUpISImJCIwXyVRYSNlVzYiISM/NyQkIiZ2VikhIiUkIjB1UWorOCkqZiIhIz83JCQiJ3ZvJSkhIiYkIjBFTSFwP2kiRyIhIz83JCQiIyYpISIiJCIwPkVFZU4oKjQjISNANyQkIidESiYpISImJCEwUnRiTTRZbyMhI0A3JCQiJkRjKSEiJSQhMFVdQiYqKiozKD0hI0A3JCQiJ3YkZikhIiYkITBpeHN4bCJ5NyEjQDckJCIlRCcpISIkJCEwaFVSJUdZSiEpISNCNyQkIidEYycpISImJCIvKFE7aktcKyUhI0A3JCQiJnZvKSEiJSQiMC9ibjdxQkklISNCNyQkIid2PSgpISImJCEwWU0/NzFQKVIhI0ItJSdMRUdFTkRHNiMtJSlfVFlQRVNFVEc2Iy1JJW1yb3dHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NS1JI21zRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNiNRL01vdGhlcn5XYXZlbGV0NiIvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW42Ii8lJXNpemVHUSMxMjYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSZmYWxzZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ25vcm1hbDYiLSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJiUmX0FYSVNHNiMiIiI2Jy0lK19HUklETElORVNHNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Jy0lK19HUklETElORVNHNiYtJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIS0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiRxJCEiIi0lKUJPVU5EU19ZRzYjJCIjXSEiIi0lLUJPVU5EU19XSURUSEc2IyQiJSFlJCEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUrTyEiIi0lKUNISUxEUkVORzYi
Let us transform all primary data by choosed wavelets:
SData_d := VectorDWT(Data_d, SDlevels, SDw1, SDw2):
SData_s := VectorDWT(Data_s, SDlevels, SDw1, SDw2):
SData_p := VectorDWT(Data_p, SDlevels, SDw1, SDw2):
Table 5.
plots[listplot](SData_s,thickness=3,color=red,gridlines=true,title=`fig.9 Wavelet transformation of Systola (Symlet,10)`,font=[helvetica,14],titlefont=[Arial,14],view=[0..80,-50..500],tickmarks=[8,20],labeldirections=[horizontal,vertical],labels=[`counts`,`w-coeffs`]);
LSUlUExPVEc2LS0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJCFvISIiLSUpQk9VTkRTX1lHNiMkIiRnJCEiIi0lLUJPVU5EU19XSURUSEc2IyQiJV1QISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJTVNISIiLSUpQ0hJTERSRU5HNiItJSdDVVJWRVNHNiM3XHA3JCQiIzUhIiIkIjFyRHh0VWFpViEjODckJCIjPyEiIiQiMid6JUcjPkVFIkglISM5NyQkIiNJISIiJCIxTzgoKVssJEgxJSEjODckJCIjUyEiIiQiMS5OXyZ6ISkqSFkhIzg3JCQiI10hIiIkIjFbQ3hSPnNnViEjODckJCIjZyEiIiQiMTN3VSxjY3JYISM4NyQkIiNxISIiJCIyLlBIJGU8TjRaISM5NyQkIiMhKSEiIiQiMSJbeiQpSC8kKkclISM4NyQkIiMhKiEiIiQiMFVpKWYkcFlGJSEjNzckJCIkKyIhIiIkIjFKPjRHK2xDVSEjODckJCIkNSIhIiIkIjIqKjMkKVsoSFVrOyEjOzckJCIkPyIhIiIkIjJHJGYqKVxRXCdHIiEjOjckJCIkSSIhIiIkITJQJypRMm9JeXMiISM6NyQkIiRTIiEiIiQhMnRUR1IzSEo/IiEjOjckJCIkXSIhIiIkIjJkdT46WyR5NUAhIzo3JCQiJGciISIiJCIyYiEqKSlRJkh0IVIlISM7NyQkIiRxIiEiIiQiMWtMKnpJT09vKSEjOjckJCIkIT0hIiIkITJbaz1kb1R5QCIhIzo3JCQiJCE+ISIiJCEvQlsjbz1vQCchIzg3JCQiJCsjISIiJCEyUm1sWCt4aC8iISM6NyQkIiQ1IyEiIiQhMSI0UVxHJ293TyEjOzckJCIkPyMhIiIkIjFUJCkpR11cKHkmKSEjOzckJCIkSSMhIiIkITI5JHlMJ2Y1WzQiISM7NyQkIiRTIyEiIiQhMU4nZnFFTT47KiEjOjckJCIkXSMhIiIkITIpKm8iUVQqeUpvIyEjOzckJCIkZyMhIiIkITJPbzlBRyUqXHojISM7NyQkIiRxIyEiIiQiMk4qSGMwLFEkKm8hIzs3JCQiJCFHISIiJCIxKSlcPSI9aXg5IyEjOjckJCIkIUghIiIkITIydjhsYyE0LjUhIzo3JCQiJCskISIiJCEyMGVpT3dCJDQ4ISM6NyQkIiQ1JCEiIiQiMWNOcjslZVd5IiEjOTckJCIkPyQhIiIkITEoeWohPUE9L0EhIzo3JCQiJEkkISIiJCIxKls+TG1ZbTMlISM6NyQkIiRTJCEiIiQhMXZrUVo+OyQqUiEjOjckJCIkXSQhIiIkITIsKVFKXXFqMDchIzo3JCQiJGckISIiJCIxWHBcQlEuPGchIzo3JCQiJHEkISIiJCExXiVmY0VtaTcjISM5NyQkIiQhUSEiIiQhMT5kWyFvcF5cJiEjOjckJCIkIVIhIiIkIjIvUHhKbkNcKD4hIzs3JCQiJCslISIiJCIyelk0OGRaJXo4ISM6NyQkIiQ1JSEiIiQhMXIjKipwR1coKlEoISM6NyQkIiQ/JSEiIiQhMGxWYmR2cSQ9ISM5NyQkIiRJJSEiIiQiMi9dMiVvWz8tdiEjOzckJCIkUyUhIiIkITAlbypcOSRSQUghIzk3JCQiJF0lISIiJCIxQzlHZSI9MzUmISM7NyQkIiRnJSEiIiQhMm5aLTVTQiFwNiEjOjckJCIkcSUhIiIkITFLO3crOllQayEjOzckJCIkIVshIiIkITEjKip6cyQpKVtRVyEjOjckJCIkIVwhIiIkITEqUkZTTSl6JylcISM8NyQkIiQrJiEiIiQhMkMjSCl5eGlbWiIhIzs3JCQiJDUmISIiJCIxOSR5UUsoM0ljISM6NyQkIiQ/JiEiIiQiMTp2QSdweHJuKCEjOjckJCIkSSYhIiIkIjJGOipmOWhhKjMiISM6NyQkIiRTJiEiIiQhMChRJj5eJVxrXSEjOTckJCIkXSYhIiIkIjIneVs7aDQkb1ciISM7NyQkIiRnJiEiIiQhMiNwZDcjPiNcW0shIzs3JCQiJHEmISIiJCIwQjMqUW4qXGIqISM5NyQkIiQhZSEiIiQiMjoieVNAQiQ9PSIhIzs3JCQiJCFmISIiJCExLSQpW2hAMGReISM6NyQkIiQrJyEiIiQhMVs4VjQlUVhYIiEjOTckJCIkNSchIiIkIjEzSSx2LTldUyEjOzckJCIkPychIiIkIjEpPT9edXNodiUhIzo3JCQiJEknISIiJCIxXmp6a2J3LnEhIzo3JCQiJFMnISIiJCIxImZ4Ymc6YCJvISM6NyQkIiRdJyEiIiQiMiYpeiIqcDgyYnQnISM7NyQkIiRnJyEiIiQiMWMoKls5SHVReSEjOjckJCIkcSchIiIkIjI8J0hASTVycjchIzs3JCQiJCFvISIiJCExL1JuYHZAJkgpISM6NyQkIiQhcCEiIiQhMkExUT8jPnUlUiIhIzo3JCQiJCsoISIiJCIyR21Wcl1lZ3kjISM7NyQkIiQ1KCEiIiQhMmw2OUV3ODQvKCEjOzckJCIkPyghIiIkITFpI3A4XUYneUYhIzo3JCQiJEkoISIiJCExdV1pJipRVTRdISM6NyQkIiRTKCEiIiQhMnd1IykpW1BXPDshIzw3JCQiJF0oISIiJCEyOFdCQzsqUlM7ISM7NyQkIiRnKCEiIiQhMDJGIW9mLEpBISM4NyQkIiRxKCEiIiQiMS04P3UrMklRISM6NyQkIiQheSEiIiQhMSdSOV5IcGlwJSEjOjckJCIkIXohIiIkITFyJ1EicERvdEIhIzk3JCQiJCspISIiJCEwcC9dSCczKCkpKiEjOS0lJVZJRVdHNiQ7JCIiISEiIiQiJCspISIiOyQhJCsmISIiJCIlK10hIiItJiUmX0FYSVNHNiMiIiI2Jy0lK19HUklETElORVNHNiMlKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiItJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNictJStfR1JJRExJTkVTRzYjJShERUZBVUxURy0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIiLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lKkFYRVNUSUNLU0c2JCIiKSIjPy0lK0FYRVNMQUJFTFNHNictSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USdjb3VudHM2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxNDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USl3LWNvZWZmczYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzE0NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii0lJUZPTlRHNiUlKERFRkFVTFRHJShERUZBVUxURyIjOSUrSE9SSVpPTlRBTEclKVZFUlRJQ0FMRy0lJlRJVExFRzYlLSUpX1RZUEVTRVRHNiMtSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY2UVRmaWcuOX5XYXZlbGV0fnRyYW5zZm9ybWF0aW9ufm9mflN5c3RvbGF+KFN5bWxldCwxMCk2Ii8lJ2ZhbWlseUdRJkFyaWFsNiIvJSVzaXplR1EjMTQ2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJTBmb250X3N0eWxlX25hbWVHUSVUZXh0NiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLSUlRk9OVEc2JCUmQXJpYWxHIiM5LSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUqVEhJQ0tORVNTRzYjIiIkLSUlRk9OVEc2JCUqaGVsdmV0aWNhRyIjOQ==
plots[listplot](SData_d,thickness=3,color=red,gridlines=true,title=`fig.10 Wavelet transformation of Diastola (Symlet,10)`,font=[helvetica,14],titlefont=[Arial,14],view=[0..80,-50..500],tickmarks=[8,20],labeldirections=[horizontal,vertical],labels=[`counts`,`w-coeffs`]);
LSUlUExPVEc2LS0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJCFvISIiLSUpQk9VTkRTX1lHNiMkIiRnJCEiIi0lLUJPVU5EU19XSURUSEc2IyQiJV1QISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJTVNISIiLSUpQ0hJTERSRU5HNiItJSdDVVJWRVNHNiM3XHA3JCQiIzUhIiIkIjFiSSFcb20mcEIhIzg3JCQiIz8hIiIkIjJkZyMpKkdSSStEISM5NyQkIiNJISIiJCIya3ZCNidHWmNCISM5NyQkIiNTISIiJCIyeCcqKXlVMiRRRSMhIzk3JCQiI10hIiIkIjItbGosYClvSEIhIzk3JCQiI2chIiIkIjIxTEolUT9KXEIhIzk3JCQiI3EhIiIkIjJDMFImM3UwS0MhIzk3JCQiIyEpISIiJCIyN0V4V1EjSFlBISM5NyQkIiMhKiEiIiQiMmB5eXNPYltQIyEjOTckJCIkKyIhIiIkIjIoW1BHI0dOKD5DISM5NyQkIiQ1IiEiIiQhMiJmSmhHWERMOSEjOjckJCIkPyIhIiIkITJEX11kZCoqRzYiISM7NyQkIiRJIiEiIiQiMjFha0FWJ1IoMyQhIzs3JCQiJFMiISIiJCEyYVRGLV9aOXcjISM7NyQkIiRdIiEiIiQhMSNHenckSCkpM08hIzo3JCQiJGciISIiJCIyJkhGMlckem1ZIyEjOzckJCIkcSIhIiIkIjJzVW4yZ2E5LSMhIzs3JCQiJCE9ISIiJCExPj0tJD5pWXkjISM6NyQkIiQhPiEiIiQhMlkqKTQneT8uKz8hIzo3JCQiJCsjISIiJCIxN05SWk0jPj8mISM6NyQkIiQ1IyEiIiQiMicpKUc+WnNvInkkISM7NyQkIiQ/IyEiIiQhMEh6Iz1QKSk0KiohIzk3JCQiJEkjISIiJCEybFA8R3YoPW1vISM7NyQkIiRTIyEiIiQiMk88KnBQdiI+IT0hIzs3JCQiJF0jISIiJCEyUHBkMCU+MnlGISM7NyQkIiRnIyEiIiQhMUJYKVEmKSkpKVE8ISM5NyQkIiRxIyEiIiQhMm0rKkdBIkgoKSk9ISM7NyQkIiQhRyEiIiQhMS93N28rVyh6IyEjOjckJCIkIUghIiIkIjFDTjAqZmgoeiMqISM6NyQkIiQrJCEiIiQiMmFwdWwmeScpNDkhIzs3JCQiJDUkISIiJCExPmFxI1tGMUklISM6NyQkIiQ/JCEiIiQiMWBfT0xtTVdgISM6NyQkIiRJJCEiIiQiMWE8TSs6O10/ISM7NyQkIiRTJCEiIiQhMXByNVMoeSVcXCEjOjckJCIkXSQhIiIkITJEcG5gJj1CaF0hIzs3JCQiJGckISIiJCIxaHk6QiFlaClwISM6NyQkIiRxJCEiIiQhMlF3ZCtEcTUsIiEjOjckJCIkIVEhIiIkITIyaFZqKFtgSzshIzs3JCQiJCFSISIiJCIyJ3BcWjMrUi05ISM8NyQkIiQrJSEiIiQhMSc9azBnRGZBJCEjOjckJCIkNSUhIiIkIjEkXDY3ITNyKnokISM6NyQkIiQ/JSEiIiQiMihbJmZAZGdwcSMhIzs3JCQiJEklISIiJCEyVWFWRVknSHg9ISM7NyQkIiRTJSEiIiQiMlo4RnghKTNtIUchIzs3JCQiJF0lISIiJCIySk5oOiVmOHg2ISM7NyQkIiRnJSEiIiQiMVApSGtUNWxVJSEjOjckJCIkcSUhIiIkITJiWShwcGlWWDshIzo3JCQiJCFbISIiJCIxbmE7aHVZRTwhIzk3JCQiJCFcISIiJCExLXBDVCVScVclISM6NyQkIiQrJiEiIiQiMSQ9I2UyOk4pUiQhIzo3JCQiJDUmISIiJCIyM2pnaGV1NDQiISM7NyQkIiQ/JiEiIiQiMTE8LFFkKzVVISM6NyQkIiRJJiEiIiQiMUYzK3d0dCFRJSEjOjckJCIkUyYhIiIkIjEmNEdNaSUpZVwjISM6NyQkIiRdJiEiIiQiMSZSVGdqNEJgJiEjOjckJCIkZyYhIiIkITJQUVR2alxaUiMhIzs3JCQiJHEmISIiJCExLHhudHpPJSlIISM6NyQkIiQhZSEiIiQiMShRZSopUiNIYlEhIzo3JCQiJCFmISIiJCEyLjVlWik9L3U+ISM7NyQkIiQrJyEiIiQhMXA/KjN1RipwYSEjOzckJCIkNSchIiIkIjJhKyk0cENYJ1wjISM7NyQkIiQ/JyEiIiQiMU84cTIiendBJyEjOzckJCIkSSchIiIkIjImZnMnPV1CSzIiISM6NyQkIiRTJyEiIiQhMTtFIyllcXYxVCEjOjckJCIkXSchIiIkITJzZTFiNEJzNiIhIzs3JCQiJGcnISIiJCEwZiVHN3pdOWkhIzo3JCQiJHEnISIiJCIwJHonUSNmUHU5ISM5NyQkIiQhbyEiIiQiMSgqNDYtOXJiQyEjOjckJCIkIXAhIiIkIjIlUiY9ViIzSUZHISM7NyQkIiQrKCEiIiQhMFhOXkh1dXMkISM5NyQkIiQ1KCEiIiQiMnVFVllhWSkpKT0hIzs3JCQiJD8oISIiJCEyO3Z0OylvYCVwJCEjOzckJCIkSSghIiIkITIoKSlHMD5DcWE2ISM8NyQkIiRTKCEiIiQiMXR5bytKIzQvKiEjOjckJCIkXSghIiIkITIkR08zPCU0RjwiISM8NyQkIiRnKCEiIiQhMVowXyhcZjwmWyEjOjckJCIkcSghIiIkITFKbXheU0lGZCEjOjckJCIkIXkhIiIkITJsKkhiMSJbPy8kISM7NyQkIiQheiEiIiQiMnVBRkYpKTMiZk0hIzs3JCQiJCspISIiJCIyJSlbRD06RVI4JCEjOy0lJVZJRVdHNiQ7JCIiISEiIiQiJCspISIiOyQhJCsmISIiJCIlK10hIiItJiUmX0FYSVNHNiMiIiI2Jy0lK19HUklETElORVNHNiMlKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiItJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNictJStfR1JJRExJTkVTRzYjJShERUZBVUxURy0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIiLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lKkFYRVNUSUNLU0c2JCIiKSIjPy0lK0FYRVNMQUJFTFNHNictSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USdjb3VudHM2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxNDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USl3LWNvZWZmczYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzE0NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii0lJUZPTlRHNiUlKERFRkFVTFRHJShERUZBVUxURyIjOSUrSE9SSVpPTlRBTEclKVZFUlRJQ0FMRy0lJlRJVExFRzYlLSUpX1RZUEVTRVRHNiMtSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY2UVZmaWcuMTB+V2F2ZWxldH50cmFuc2Zvcm1hdGlvbn5vZn5EaWFzdG9sYX4oU3ltbGV0LDEwKTYiLyUnZmFtaWx5R1EmQXJpYWw2Ii8lJXNpemVHUSMxNDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lMGZvbnRfc3R5bGVfbmFtZUdRJVRleHQ2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItJSVGT05URzYkJSZBcmlhbEciIzktJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSpUSElDS05FU1NHNiMiIiQtJSVGT05URzYkJSpoZWx2ZXRpY2FHIiM5
plots[listplot](SData_p,thickness=3,color=red,gridlines=true,title=`fig.11 Wavelet transformation of Pulse (Symlet,10)`,font=[helvetica,14],titlefont=[Arial,14],view=[0..80,-50..500],tickmarks=[8,20],labeldirections=[horizontal,vertical],labels=[`counts`,`w-coeffs`]);
LSUlUExPVEc2LS0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJCFvISIiLSUpQk9VTkRTX1lHNiMkIiRnJCEiIi0lLUJPVU5EU19XSURUSEc2IyQiJV1QISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJTVNISIiLSUpQ0hJTERSRU5HNiItJSdDVVJWRVNHNiM3XHA3JCQiIzUhIiIkIjFfbjopNC5STCMhIzg3JCQiIz8hIiIkIjE8LDUsKnpIUSMhIzg3JCQiI0khIiIkIjI5UTdeb3NuPCMhIzk3JCQiI1MhIiIkIjJWP1QkKW9QM1QjISM5NyQkIiNdISIiJCIyMyJ6YndqNTRDISM5NyQkIiNnISIiJCIyJW9CJHAocCdITSMhIzk3JCQiI3EhIiIkIjJ1WjQqZlNCU0MhIzk3JCQiIyEpISIiJCIxYCMqeSJIMEFEIyEjODckJCIjISohIiIkIjIoW201Y2IyZkEhIzk3JCQiJCsiISIiJCIyMypvWXFZJXBGIyEjOTckJCIkNSIhIiIkITJtOCpvUSopKltrJCEjOzckJCIkPyIhIiIkIjJgOVc5M25OJj4hIzs3JCQiJEkiISIiJCIxWE11SEgpXCJIISM6NyQkIiRTIiEiIiQhMVx3LThMJ2U3JyEjOjckJCIkXSIhIiIkIjJESD48MkVoUiMhIzs3JCQiJGciISIiJCIyInlQPXpIWlA5ISM7NyQkIiRxIiEiIiQiMSNbXG0xRCQ9cSEjOjckJCIkIT0hIiIkITJiaDh0KXA9ZV4hIzs3JCQiJCE+ISIiJCExYSZwSz0wTiJlISM6NyQkIiQrIyEiIiQiMUc1Vl0wUS11ISM6NyQkIiQ1IyEiIiQhMURZc0FTYEZbISM6NyQkIiQ/IyEiIiQhMiYqeSZSdCY0bWoiISM7NyQkIiRJIyEiIiQhMUZrMzIsO0NjISM7NyQkIiRTIyEiIiQiMEhKSSY+Z3EjKiEjOTckJCIkXSMhIiIkITI8bnlXLDkxeiMhIzs3JCQiJGcjISIiJCEyM0l5MlcmUls7ISM7NyQkIiRxIyEiIiQiMiUpXFQqcGJQUD8hIzs3JCQiJCFHISIiJCIyMGY/JXlNdElFISM8NyQkIiQhSCEiIiQhMUYtMiFHIzR2IykhIzo3JCQiJCskISIiJCExKVJnWDJLaHMlISM6NyQkIiQ1JCEiIiQiMXYiNF8hR3JldCEjOjckJCIkPyQhIiIkIjFYckt2QSRmYSUhIzs3JCQiJEkkISIiJCIxWyopcDZoIXBfIyEjOjckJCIkUyQhIiIkITE6PGRLI3lnRiYhIzo3JCQiJF0kISIiJCEyTypbVT0mNEIiPSEjOzckJCIkZyQhIiIkIjBVanJ2NyR5SCEjOTckJCIkcSQhIiIkITEiejwlSC5kayEpISM7NyQkIiQhUSEiIiQhMiRIQ1ZFX3ZvSyEjOzckJCIkIVIhIiIkIjFrQyJcKzV0UCkhIzs3JCQiJCslISIiJCIxN2dmQ09NNyoqISM6NyQkIiQ1JSEiIiQiMEo7WCVcI3lSJCEjOTckJCIkPyUhIiIkITE8JFtLJltabVkhIzs3JCQiJEklISIiJCIwRkdHOytAI3ohIzo3JCQiJFMlISIiJCIvKEhoRjstTyMhIzg3JCQiJF0lISIiJCIybWRUS04rckkkISM7NyQkIiRnJSEiIiQhMTsmUSRmdnR0SiEjOjckJCIkcSUhIiIkIjJXZmUxaUU6QCMhIzs3JCQiJCFbISIiJCExNzVVKjR1elIiISM6NyQkIiQhXCEiIiQhMnV6QmZadWpEIyEjOzckJCIkKyYhIiIkITBlSVlQW10tJCEjOTckJCIkNSYhIiIkITFXO0s+IzQleSIpISM6NyQkIiQ/JiEiIiQhMTNHIkc+RVFTJyEjOjckJCIkSSYhIiIkIjF3U0JZXyk0KHkhIzo3JCQiJFMmISIiJCEwL1YhXFguWGIhIzo3JCQiJF0mISIiJCIyZTtOJD5TeDZDISM7NyQkIiRnJiEiIiQhMShlU0cjZTEoKlEhIzo3JCQiJHEmISIiJCIwRmNrJSo9MD8lISM5NyQkIiQhZSEiIiQiMVpQJ3B5YFl2IyEjOjckJCIkIWYhIiIkIjEqKnlIXU5HTlIhIzo3JCQiJCsnISIiJCEyd0teenRGKCkzIiEjOzckJCIkNSchIiIkITF1IylwbClRaSRRISM6NyQkIiQ/JyEiIiQiMllSQDsuQSkpZiMhIzs3JCQiJEknISIiJCExWChRUyhHOC9WISM6NyQkIiRTJyEiIiQhMi5HXXdsZVkuIiEjOzckJCIkXSchIiIkIjApKT1AKCpcLkUjISM5NyQkIiRnJyEiIiQiMW1ETnIpXElsKCEjOzckJCIkcSchIiIkIjJPKFxNWiVRV1gkISM7NyQkIiQhbyEiIiQiMThcJzQleVNVTCEjOjckJCIkIXAhIiIkITFhVjlyX0A4WCEjOjckJCIkKyghIiIkIjFid1dudmBdVCEjOjckJCIkNSghIiIkITJuJz5mXTgnKlFKISM7NyQkIiQ/KCEiIiQhMl86OmBXbkpqIiEjOzckJCIkSSghIiIkITE0M1RbcDpXXCEjOjckJCIkUyghIiIkIjJLI3AmZTpaYEciISM7NyQkIiRdKCEiIiQiMk9zJ1JfJSpcTUUhIzw3JCQiJGcoISIiJCIxTmElSCNcIylcTCEjOjckJCIkcSghIiIkITJ5P1tiV2V2ZiIhIzo3JCQiJCF5ISIiJCIxby0xNVhUIW8lISM6NyQkIiQheiEiIiQhMT5rZyIqcHFtYCEjOjckJCIkKykhIiIkIjJXbj5YYnB6UiMhIzwtJSVWSUVXRzYkOyQiIiEhIiIkIiQrKSEiIjskISQrJiEiIiQiJStdISIiLSYlJl9BWElTRzYjIiIiNictJStfR1JJRExJTkVTRzYjJShERUZBVUxURy0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqTElORVNUWUxFRzYjIiIiLSUqVEhJQ0tORVNTRzYjIiIhLSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0mJSZfQVhJU0c2IyIiIzYnLSUrX0dSSURMSU5FU0c2IyUoREVGQVVMVEctJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIi0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJSpBWEVTVElDS1NHNiQiIikiIz8tJStBWEVTTEFCRUxTRzYnLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEnY291bnRzNiIvJSdmYW1pbHlHUShERUZBVUxUNiIvJSVzaXplR1EjMTQ2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEpdy1jb2VmZnM2Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxNDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzklK0hPUklaT05UQUxHJSlWRVJUSUNBTEctJSZUSVRMRUc2JS0lKV9UWVBFU0VURzYjLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NlFUZmlnLjExfldhdmVsZXR+dHJhbnNmb3JtYXRpb25+b2Z+flB1bHNlfihTeW1sZXQsMTApNiIvJSdmYW1pbHlHUSZBcmlhbDYiLyUlc2l6ZUdRIzE0NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUwZm9udF9zdHlsZV9uYW1lR1ElVGV4dDYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii0lJUZPTlRHNiQlJkFyaWFsRyIjOS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJUZPTlRHNiQlKmhlbHZldGljYUciIzk=
It is noticeable from the graphs of table 5, that thewavelet coefficients (w-coeffs) differ significantly from zero only on a scale of approximation (first 10 counts). At the same time for all three scales of detail (11-20, 21-40, and 41-80 counts) the w-coeffs fluctuate around zero. We can then clear the images of signals from high-frequency noises, using the method of so-called "hard" threshold [3,4]. Let us find t first the noise levels at different scale counts, according to the recommendations of [3,4].
sigma1s:=(Statistics[Median](map(abs,(LinearAlgebra:-SubVector(SData_s,41..80))))-Statistics[Mean](LinearAlgebra:-SubVector(SData_s,41..80)))/0.6745;
sigma2s:=(Statistics[Median](map(abs,(LinearAlgebra:-SubVector(SData_s,21..40))))-Statistics[Mean](LinearAlgebra:-SubVector(SData_s,21..40)))/0.6745;
sigma3s:=(Statistics[Median](map(abs,(LinearAlgebra:-SubVector(SData_s,11..20))))-Statistics[Mean](LinearAlgebra:-SubVector(SData_s,11..20)))/0.6745;
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SShzaWdtYTFzR0YoJCIwLiFvKmUnNDk1ISM4NyNdNDAyNDQ4MkNBRTk5Mzc5Qg==
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SShzaWdtYTJzR0YoJCIwNmkjbzMmR1AqISM5NyNdNDAyMkJFRTY0RjFBRjEyMg==
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SShzaWdtYTNzR0YoJCIwaj4pR2ljMj0hIzg3I100MDMyMTM1RTlBOEQ1N0VD
Now these values would be corrected by numbers of counts [3,4]
SDthresh1s:=evalf(sigma1s*sqrt(2*ln(40))):
SDthresh2s:=evalf(sigma2s*sqrt(2*ln(20))):
SDthresh3s:=evalf(sigma3s*sqrt(2*ln(10))):
T:=[SDthresh1s,SDthresh2s,SDthresh3s];
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSJUR0YoNyUkIjBtOjFCI1xhRiEjOCQiMFVBNWdPVUgjRjEkIjBIJSlccnYqeVFGMTcjNyVdNDAzQjhCODAwNzNCQ0RCNV00MDM2RjEzRUU2MUFDNjVCXTQwNDM2NTE2QzMyNTJEQzM=
One can see that the thresholds (4,.4) are slightly differ for different scales. We are guaranteed to get zero w-coeffs on all three scales, choosing the highest of these thresholds
SDthresh_s:=max(T);
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SStTRHRocmVzaF9zR0YoJCIwSCUpXHJ2KnlRISM4NyNdNDA0MzY1MTZDMzI1MkRDMw==
Thresholds for diastolic pressure an pulse frequencies were obtained in the same way
SDthresh_d := 22.2790061630243;
SDthresh_p := %:
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SStTRHRocmVzaF9kR0YoJCIwVi1qaCt6QSMhIzg3I0Yu
Previous cleaning from fluctuations with \302\253hard threshold\302\273 method uses the obtained thresholds SDRaf_1s:= Vector(map( z -> `if`(abs(z) <= SDthresh_s, 0, z),SData_s),datatype = float[8]):
SDRaf_1d:= Vector(map( z -> `if`(abs(z) <= SDthresh_d, 0, z),SData_d),datatype = float[8]):
SDRaf_1p:= Vector(map( z -> `if`(abs(z) <= SDthresh_p, 0, z),SData_p),datatype = float[8]):
About the same way we can de-noise w-coeffs on the approximation scale, using the statistical standard deviation as threshold and the mathematical expectation as a level for the w-coefficients on these scales (1-10 counts)
avs:=Statistics[Mean](LinearAlgebra:-SubVector(SDRaf_1s,1..10));
var_s:=sqrt(Statistics[Variance](LinearAlgebra:-SubVector(SDRaf_1s,1..10)));
avd:=Statistics[Mean](LinearAlgebra:-SubVector(SDRaf_1d,1..10));
var_d:=sqrt(Statistics[Variance](LinearAlgebra:-SubVector(SDRaf_1d,1..10)));
avp:=Statistics[Mean](LinearAlgebra:-SubVector(SDRaf_1p,1..10));
var_p:=sqrt(Statistics[Variance](LinearAlgebra:-SubVector(SDRaf_1p,1..10)));
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSRhdnNHRigkIjAiZksjMylweFYhIzc3I100MDdCNUM1MTIyNkZGMkMw
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSZ2YXJfc0dGKCQiMGRQXkcuLysjISM4NyNdNDAzNDAxMDg0QzA0QzRBMw==
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSRhdmRHRigkIjBBKCpHXzZVTyMhIzc3I100MDZEOEQ3QTE0NURERkJF
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSZ2YXJfZEdGKCQiMHREIzM3JG9lKCEjOTcjXTQwMUU1OEVBNDdCRTY4MjA=
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSRhdnBHRigkIjBJWi9qLSZHQiEjNzcjXTQwNkQxQjM1NUFEODQzOUQ=
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSZ2YXJfcEdGKCQiMDBNN0xfcl4pISM5NyNdNDAyMTA4QzgzMEJCNjY5Qg==
SDRaf_s:= Vector(map( z -> `if`(abs(z-avs) <= var_s, avs, z),SDRaf_1s),datatype = float[8]):
SDRaf_d:= Vector(map( z -> `if`(abs(z-avd) <= var_d, avd, z),SDRaf_1d),datatype = float[8]):
SDRaf_p:= Vector(map( z -> `if`(abs(z-avp) <= var_p, avp, z),SDRaf_1p),datatype = float[8]):
Table 6 shows the de-noised wavelet-images of sygnals for a comparicon with table 5
Table 6
plots[listplot](SDRaf_s,thickness=3,color=red,gridlines=true,title=`fig.12 Wavelet transformation of Systola (Symlet,10)`,font=[helvetica,14],titlefont=[Arial,14],view=[0..80,-50..500],tickmarks=[8,20],labeldirections=[horizontal,vertical],labels=[`counts`,`w-coeffs`]);
LSUlUExPVEc2LS0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJCFvISIiLSUpQk9VTkRTX1lHNiMkIiRxJiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJSE9JCEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUhPSQhIiItJSlDSElMRFJFTkc2Ii0lJ0NVUlZFU0c2IzdccDckJCIjNSEiIiQiMlkiZksjMylweFYhIzk3JCQiIz8hIiIkIjJZImZLIzMpcHhWISM5NyQkIiNJISIiJCIxTzgoKVssJEgxJSEjODckJCIjUyEiIiQiMS5OXyZ6ISkqSFkhIzg3JCQiI10hIiIkIjJZImZLIzMpcHhWISM5NyQkIiNnISIiJCIyWSJmSyMzKXB4ViEjOTckJCIjcSEiIiQiMi5QSCRlPE40WiEjOTckJCIjISkhIiIkIjJZImZLIzMpcHhWISM5NyQkIiMhKiEiIiQiMlkiZksjMylweFYhIzk3JCQiJCsiISIiJCIyWSJmSyMzKXB4ViEjOTckJCIkNSIhIiIkIiIhISIiNyQkIiQ/IiEiIiQiIiEhIiI3JCQiJEkiISIiJCIiISEiIjckJCIkUyIhIiIkIiIhISIiNyQkIiRdIiEiIiQiIiEhIiI3JCQiJGciISIiJCIiISEiIjckJCIkcSIhIiIkIiIhISIiNyQkIiQhPSEiIiQiIiEhIiI3JCQiJCE+ISIiJCIiISEiIjckJCIkKyMhIiIkIiIhISIiNyQkIiQ1IyEiIiQiIiEhIiI3JCQiJD8jISIiJCIiISEiIjckJCIkSSMhIiIkIiIhISIiNyQkIiRTIyEiIiQiIiEhIiI3JCQiJF0jISIiJCIiISEiIjckJCIkZyMhIiIkIiIhISIiNyQkIiRxIyEiIiQiIiEhIiI3JCQiJCFHISIiJCIiISEiIjckJCIkIUghIiIkIiIhISIiNyQkIiQrJCEiIiQiIiEhIiI3JCQiJDUkISIiJCIiISEiIjckJCIkPyQhIiIkIiIhISIiNyQkIiRJJCEiIiQiIiEhIiI3JCQiJFMkISIiJCIiISEiIjckJCIkXSQhIiIkIiIhISIiNyQkIiRnJCEiIiQiIiEhIiI3JCQiJHEkISIiJCIiISEiIjckJCIkIVEhIiIkIiIhISIiNyQkIiQhUiEiIiQiIiEhIiI3JCQiJCslISIiJCIiISEiIjckJCIkNSUhIiIkIiIhISIiNyQkIiQ/JSEiIiQiIiEhIiI3JCQiJEklISIiJCIiISEiIjckJCIkUyUhIiIkIiIhISIiNyQkIiRdJSEiIiQiIiEhIiI3JCQiJGclISIiJCIiISEiIjckJCIkcSUhIiIkIiIhISIiNyQkIiQhWyEiIiQiIiEhIiI3JCQiJCFcISIiJCIiISEiIjckJCIkKyYhIiIkIiIhISIiNyQkIiQ1JiEiIiQiIiEhIiI3JCQiJD8mISIiJCIiISEiIjckJCIkSSYhIiIkIiIhISIiNyQkIiRTJiEiIiQiIiEhIiI3JCQiJF0mISIiJCIiISEiIjckJCIkZyYhIiIkIiIhISIiNyQkIiRxJiEiIiQiIiEhIiI3JCQiJCFlISIiJCIiISEiIjckJCIkIWYhIiIkIiIhISIiNyQkIiQrJyEiIiQiIiEhIiI3JCQiJDUnISIiJCIiISEiIjckJCIkPychIiIkIiIhISIiNyQkIiRJJyEiIiQiIiEhIiI3JCQiJFMnISIiJCIiISEiIjckJCIkXSchIiIkIiIhISIiNyQkIiRnJyEiIiQiIiEhIiI3JCQiJHEnISIiJCIiISEiIjckJCIkIW8hIiIkIiIhISIiNyQkIiQhcCEiIiQiIiEhIiI3JCQiJCsoISIiJCIiISEiIjckJCIkNSghIiIkIiIhISIiNyQkIiQ/KCEiIiQiIiEhIiI3JCQiJEkoISIiJCIiISEiIjckJCIkUyghIiIkIiIhISIiNyQkIiRdKCEiIiQiIiEhIiI3JCQiJGcoISIiJCIiISEiIjckJCIkcSghIiIkIiIhISIiNyQkIiQheSEiIiQiIiEhIiI3JCQiJCF6ISIiJCIiISEiIjckJCIkKykhIiIkIiIhISIiLSUlVklFV0c2JDskIiIhISIiJCIkKykhIiI7JCEkKyYhIiIkIiUrXSEiIi0mJSZfQVhJU0c2IyIiIjYnLSUrX0dSSURMSU5FU0c2IyUoREVGQVVMVEctJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIi0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Jy0lK19HUklETElORVNHNiMlKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiItJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUqQVhFU1RJQ0tTRzYkIiIpIiM/LSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRJ2NvdW50czYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzE0NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRKXctY29lZmZzNiIvJSdmYW1pbHlHUShERUZBVUxUNiIvJSVzaXplR1EjMTQ2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM5JStIT1JJWk9OVEFMRyUpVkVSVElDQUxHLSUmVElUTEVHNiUtJSlfVFlQRVNFVEc2Iy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjZRVWZpZy4xMn5XYXZlbGV0fnRyYW5zZm9ybWF0aW9ufm9mflN5c3RvbGF+KFN5bWxldCwxMCk2Ii8lJ2ZhbWlseUdRJkFyaWFsNiIvJSVzaXplR1EjMTQ2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJTBmb250X3N0eWxlX25hbWVHUSVUZXh0NiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLSUlRk9OVEc2JCUmQXJpYWxHIiM5LSUtVFJBTlNQQVJFTkNZRzYjJCIiISEiIi0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUqVEhJQ0tORVNTRzYjIiIkLSUlRk9OVEc2JCUqaGVsdmV0aWNhRyIjOQ==
plots[listplot](SDRaf_d,thickness=3,color=red,gridlines=true,title=`fig.13 Wavelet transformation of Diastola (Symlet,10)`,font=[helvetica,14],titlefont=[Arial,14],view=[0..80,-50..500],tickmarks=[8,20],labeldirections=[horizontal,vertical],labels=[`counts`,`w-coeffs`]);
LSUlUExPVEc2LS0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJCFvISIiLSUpQk9VTkRTX1lHNiMkIiRxJiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJSE9JCEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUhPSQhIiItJSlDSElMRFJFTkc2Ii0lJ0NVUlZFU0c2IzdccDckJCIjNSEiIiQiMUBzKkdfNlVPIyEjODckJCIjPyEiIiQiMmRnIykqR1JJK0QhIzk3JCQiI0khIiIkIjFAcypHXzZVTyMhIzg3JCQiI1MhIiIkIjJ4JyopeVUyJFFFIyEjOTckJCIjXSEiIiQiMUBzKkdfNlVPIyEjODckJCIjZyEiIiQiMUBzKkdfNlVPIyEjODckJCIjcSEiIiQiMUBzKkdfNlVPIyEjODckJCIjISkhIiIkIjI3RXhXUSNIWUEhIzk3JCQiIyEqISIiJCIxQHMqR182VU8jISM4NyQkIiQrIiEiIiQiMUBzKkdfNlVPIyEjODckJCIkNSIhIiIkIiIhISIiNyQkIiQ/IiEiIiQiIiEhIiI3JCQiJEkiISIiJCIiISEiIjckJCIkUyIhIiIkIiIhISIiNyQkIiRdIiEiIiQiIiEhIiI3JCQiJGciISIiJCIiISEiIjckJCIkcSIhIiIkIiIhISIiNyQkIiQhPSEiIiQiIiEhIiI3JCQiJCE+ISIiJCIiISEiIjckJCIkKyMhIiIkIiIhISIiNyQkIiQ1IyEiIiQiIiEhIiI3JCQiJD8jISIiJCIiISEiIjckJCIkSSMhIiIkIiIhISIiNyQkIiRTIyEiIiQiIiEhIiI3JCQiJF0jISIiJCIiISEiIjckJCIkZyMhIiIkIiIhISIiNyQkIiRxIyEiIiQiIiEhIiI3JCQiJCFHISIiJCIiISEiIjckJCIkIUghIiIkIiIhISIiNyQkIiQrJCEiIiQiIiEhIiI3JCQiJDUkISIiJCIiISEiIjckJCIkPyQhIiIkIiIhISIiNyQkIiRJJCEiIiQiIiEhIiI3JCQiJFMkISIiJCIiISEiIjckJCIkXSQhIiIkIiIhISIiNyQkIiRnJCEiIiQiIiEhIiI3JCQiJHEkISIiJCIiISEiIjckJCIkIVEhIiIkIiIhISIiNyQkIiQhUiEiIiQiIiEhIiI3JCQiJCslISIiJCIiISEiIjckJCIkNSUhIiIkIiIhISIiNyQkIiQ/JSEiIiQiIiEhIiI3JCQiJEklISIiJCIiISEiIjckJCIkUyUhIiIkIiIhISIiNyQkIiRdJSEiIiQiIiEhIiI3JCQiJGclISIiJCIiISEiIjckJCIkcSUhIiIkIiIhISIiNyQkIiQhWyEiIiQiIiEhIiI3JCQiJCFcISIiJCIiISEiIjckJCIkKyYhIiIkIiIhISIiNyQkIiQ1JiEiIiQiIiEhIiI3JCQiJD8mISIiJCIiISEiIjckJCIkSSYhIiIkIiIhISIiNyQkIiRTJiEiIiQiIiEhIiI3JCQiJF0mISIiJCIiISEiIjckJCIkZyYhIiIkIiIhISIiNyQkIiRxJiEiIiQiIiEhIiI3JCQiJCFlISIiJCIiISEiIjckJCIkIWYhIiIkIiIhISIiNyQkIiQrJyEiIiQiIiEhIiI3JCQiJDUnISIiJCIiISEiIjckJCIkPychIiIkIiIhISIiNyQkIiRJJyEiIiQiIiEhIiI3JCQiJFMnISIiJCIiISEiIjckJCIkXSchIiIkIiIhISIiNyQkIiRnJyEiIiQiIiEhIiI3JCQiJHEnISIiJCIiISEiIjckJCIkIW8hIiIkIiIhISIiNyQkIiQhcCEiIiQiIiEhIiI3JCQiJCsoISIiJCIiISEiIjckJCIkNSghIiIkIiIhISIiNyQkIiQ/KCEiIiQiIiEhIiI3JCQiJEkoISIiJCIiISEiIjckJCIkUyghIiIkIiIhISIiNyQkIiRdKCEiIiQiIiEhIiI3JCQiJGcoISIiJCIiISEiIjckJCIkcSghIiIkIiIhISIiNyQkIiQheSEiIiQiIiEhIiI3JCQiJCF6ISIiJCIiISEiIjckJCIkKykhIiIkIiIhISIiLSUlVklFV0c2JDskIiIhISIiJCIkKykhIiI7JCEkKyYhIiIkIiUrXSEiIi0mJSZfQVhJU0c2IyIiIjYnLSUrX0dSSURMSU5FU0c2IyUoREVGQVVMVEctJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIi0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Jy0lK19HUklETElORVNHNiMlKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiItJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUqQVhFU1RJQ0tTRzYkIiIpIiM/LSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRJ2NvdW50czYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzE0NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRKXctY29lZmZzNiIvJSdmYW1pbHlHUShERUZBVUxUNiIvJSVzaXplR1EjMTQ2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM5JStIT1JJWk9OVEFMRyUpVkVSVElDQUxHLSUmVElUTEVHNiUtJSlfVFlQRVNFVEc2Iy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjZRVmZpZy4xM35XYXZlbGV0fnRyYW5zZm9ybWF0aW9ufm9mfkRpYXN0b2xhfihTeW1sZXQsMTApNiIvJSdmYW1pbHlHUSZBcmlhbDYiLyUlc2l6ZUdRIzE0NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUwZm9udF9zdHlsZV9uYW1lR1ElVGV4dDYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii0lJUZPTlRHNiQlJkFyaWFsRyIjOS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJSZDT0xPUkc2JiUkUkdCRyQiIzUhIiIkIiIhISIiJCIiISEiIi0lKlRISUNLTkVTU0c2IyIiJC0lJUZPTlRHNiQlKmhlbHZldGljYUciIzk=
plots[listplot](SDRaf_p,thickness=3,color=red,gridlines=true,title=`fig.14 Wavelet transformation of Pulse (Symlet,10)`,font=[helvetica,14],titlefont=[Arial,14],view=[0..80,-50..500],tickmarks=[8,20],labeldirections=[horizontal,vertical],labels=[`counts`,`w-coeffs`]);
LSUlUExPVEc2LS0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJCFvISIiLSUpQk9VTkRTX1lHNiMkIiRxJiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJSE9JCEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUhPSQhIiItJSlDSElMRFJFTkc2Ii0lJ0NVUlZFU0c2IzdccDckJCIjNSEiIiQiMi9JWi9qLSZHQiEjOTckJCIjPyEiIiQiMi9JWi9qLSZHQiEjOTckJCIjSSEiIiQiMjlRN15vc248IyEjOTckJCIjUyEiIiQiMi9JWi9qLSZHQiEjOTckJCIjXSEiIiQiMi9JWi9qLSZHQiEjOTckJCIjZyEiIiQiMi9JWi9qLSZHQiEjOTckJCIjcSEiIiQiMnVaNCpmU0JTQyEjOTckJCIjISkhIiIkIjIvSVovai0mR0IhIzk3JCQiIyEqISIiJCIyL0laL2otJkdCISM5NyQkIiQrIiEiIiQiMi9JWi9qLSZHQiEjOTckJCIkNSIhIiIkIiIhISIiNyQkIiQ/IiEiIiQiIiEhIiI3JCQiJEkiISIiJCIiISEiIjckJCIkUyIhIiIkIiIhISIiNyQkIiRdIiEiIiQiIiEhIiI3JCQiJGciISIiJCIiISEiIjckJCIkcSIhIiIkIiIhISIiNyQkIiQhPSEiIiQiIiEhIiI3JCQiJCE+ISIiJCIiISEiIjckJCIkKyMhIiIkIiIhISIiNyQkIiQ1IyEiIiQiIiEhIiI3JCQiJD8jISIiJCIiISEiIjckJCIkSSMhIiIkIiIhISIiNyQkIiRTIyEiIiQiIiEhIiI3JCQiJF0jISIiJCIiISEiIjckJCIkZyMhIiIkIiIhISIiNyQkIiRxIyEiIiQiIiEhIiI3JCQiJCFHISIiJCIiISEiIjckJCIkIUghIiIkIiIhISIiNyQkIiQrJCEiIiQiIiEhIiI3JCQiJDUkISIiJCIiISEiIjckJCIkPyQhIiIkIiIhISIiNyQkIiRJJCEiIiQiIiEhIiI3JCQiJFMkISIiJCIiISEiIjckJCIkXSQhIiIkIiIhISIiNyQkIiRnJCEiIiQiIiEhIiI3JCQiJHEkISIiJCIiISEiIjckJCIkIVEhIiIkIiIhISIiNyQkIiQhUiEiIiQiIiEhIiI3JCQiJCslISIiJCIiISEiIjckJCIkNSUhIiIkIiIhISIiNyQkIiQ/JSEiIiQiIiEhIiI3JCQiJEklISIiJCIiISEiIjckJCIkUyUhIiIkIiIhISIiNyQkIiRdJSEiIiQiIiEhIiI3JCQiJGclISIiJCIiISEiIjckJCIkcSUhIiIkIiIhISIiNyQkIiQhWyEiIiQiIiEhIiI3JCQiJCFcISIiJCIiISEiIjckJCIkKyYhIiIkIiIhISIiNyQkIiQ1JiEiIiQiIiEhIiI3JCQiJD8mISIiJCIiISEiIjckJCIkSSYhIiIkIiIhISIiNyQkIiRTJiEiIiQiIiEhIiI3JCQiJF0mISIiJCIiISEiIjckJCIkZyYhIiIkIiIhISIiNyQkIiRxJiEiIiQiIiEhIiI3JCQiJCFlISIiJCIiISEiIjckJCIkIWYhIiIkIiIhISIiNyQkIiQrJyEiIiQiIiEhIiI3JCQiJDUnISIiJCIiISEiIjckJCIkPychIiIkIiIhISIiNyQkIiRJJyEiIiQiIiEhIiI3JCQiJFMnISIiJCIiISEiIjckJCIkXSchIiIkIiIhISIiNyQkIiRnJyEiIiQiIiEhIiI3JCQiJHEnISIiJCIiISEiIjckJCIkIW8hIiIkIiIhISIiNyQkIiQhcCEiIiQiIiEhIiI3JCQiJCsoISIiJCIiISEiIjckJCIkNSghIiIkIiIhISIiNyQkIiQ/KCEiIiQiIiEhIiI3JCQiJEkoISIiJCIiISEiIjckJCIkUyghIiIkIiIhISIiNyQkIiRdKCEiIiQiIiEhIiI3JCQiJGcoISIiJCIiISEiIjckJCIkcSghIiIkIiIhISIiNyQkIiQheSEiIiQiIiEhIiI3JCQiJCF6ISIiJCIiISEiIjckJCIkKykhIiIkIiIhISIiLSUlVklFV0c2JDskIiIhISIiJCIkKykhIiI7JCEkKyYhIiIkIiUrXSEiIi0mJSZfQVhJU0c2IyIiIjYnLSUrX0dSSURMSU5FU0c2IyUoREVGQVVMVEctJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIi0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Jy0lK19HUklETElORVNHNiMlKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiItJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUqQVhFU1RJQ0tTRzYkIiIpIiM/LSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRJ2NvdW50czYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzE0NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRKXctY29lZmZzNiIvJSdmYW1pbHlHUShERUZBVUxUNiIvJSVzaXplR1EjMTQ2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM5JStIT1JJWk9OVEFMRyUpVkVSVElDQUxHLSUmVElUTEVHNiUtJSlfVFlQRVNFVEc2Iy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjZRVGZpZy4xNH5XYXZlbGV0fnRyYW5zZm9ybWF0aW9ufm9mfn5QdWxzZX4oU3ltbGV0LDEwKTYiLyUnZmFtaWx5R1EmQXJpYWw2Ii8lJXNpemVHUSMxNDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lMGZvbnRfc3R5bGVfbmFtZUdRJVRleHQ2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItJSVGT05URzYkJSZBcmlhbEciIzktJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSpUSElDS05FU1NHNiMiIiQtJSVGT05URzYkJSpoZWx2ZXRpY2FHIiM5
Wavelet-images now are describing literally a few different from zero coefficients, i.e. they are effectively compressed and de-niised. Now we can reconstruct all sygnals by inverse transformations
SDataNew_s := InverseVectorDWT(SDRaf_s, SDlevels, SDw1, SDw2):
SDataNew_d := InverseVectorDWT(SDRaf_d, SDlevels, SDw1, SDw2):
SDataNew_p := InverseVectorDWT(SDRaf_p, SDlevels, SDw1, SDw2):
Let us present the results graphically
g2s:=plots[listplot](SDataNew_s,thickness=5,color=[red],gridlines=true,title=`\320\240\320\270\321\201.15 Systolic pressure`,font=[helvetica,14],titlef ont=[Arial,14],view=[0..80,125..185],tickmarks=[15,16],labeldirections=[horizontal,vertical],labels=[`counts`,`pressure, mm Hg`]):
g2d:=plots[listplot](SDataNew_d,thickness=5,color=[red],gridlines=true,title=`\320\240\320\270\321\201.15 Diastolic pressure`,font=[helvetica,14],titlef ont=[Arial,14],view=[0..80,65..105],tickmarks=[15,16],labeldirections=[horizontal,vertical],labels=[`counts`,`pressure, mm Hg`]):
g2p:=plots[listplot](SDataNew_p,thickness=5,color=[red],gridlines=true,title=`\320\240\320\270\321\201.15 Pulse`,font=[helvetica,14],titlef ont=[Arial,14],view=[0..80,65..105],tickmarks=[15,16],labeldirections=[horizontal,vertical],labels=[`counts`,`pulse, bpm`]):
Table 7
All curves are characterized by minimum about count\342\211\210(46-52). This feature is associated with the course of treatment of the patient in a sanatorium (Truskvavec\342\200\262, Ukraine). The effect looks quite short-living.
Comparison with the data of table 2 shows the efficiency of cleaning of small-scale fluctuations with high frequencies. De-noised and compressed sygnals of Table 7 include only changes with relative long periods. Such long periodical cyclical nature usually presents itself from autocorrelation functions (Table 8).
Statistics:-AutoCorrelation(SDataNew_s):
acfs:=plots[listplot](%,gridlines=true,thickness=3,title=`fig17. Systolic pressure autocorrelation `, font=[helvetica,14],color=red):
Statistics:-AutoCorrelation(SDataNew_d):
acfd:=plots[listplot](%,gridlines=true,thickness=3,title=`fig.18 Diastolic pressure autocorrelation`, font=[helvetica,14],color=red):
Statistics:-AutoCorrelation(SDataNew_p):
acfp:=plots[listplot](%,gridlines=true,thickness=3,title=`fig.19 Pulse autocorrelation`, font=[helvetica,14],color=red):
Table 8
These autocorrelations are quite expressive, very similar and demonstrating the main period about 12 counts.
<Text-field style="Heading 2" layout="Heading 2">5. Some comparisons and conclusions</Text-field>
Let us recalculate some statistical characteristics of de-noised sygnals
Table 9
It is easy to see that de-noising of primary signals has no influence on the mathematical expectations, but markedly increases the accuracy of the results, comparing the characteristics of table 3 and table 4.
The correlation coefficient between systolic pressure and the pulse frequency, which was essentially masked by noises within primary signals, also increased noticeably (almost twice)
Statistics:-Correlation(SDataNew_p,SDataNew_s);
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyMkIjA2Jm9XPiZwJykpISM6NyNdM0ZFQzVGQ0U5OUI4Njc1Ng==
Correlation between reconstructed data also looks much more convincing visually, as is demonstrated in table 10
Table 10
Statistics:-ScatterPlot(Data_p,Data_s,title = "fig. 20 Systola-Pulse correlation (primary data)", legend = "Point data",view=[76..88,140..170],symbol=solidcircle,symbolsize=18,color=black,gridlines=true,font=[helvetica,16],tickmarks=[12,10],labeldirections=[horizontal,vertical],labels=[`pulse, bpm`,`pressure, mmHg`]);
LSUlUExPVEc2LC0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJD8oISIiLSUpQk9VTkRTX1lHNiMkIiQhUiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJWdpISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJTVKISIiLSUpQ0hJTERSRU5HNiItJSdDVVJWRVNHNiY3XHA3JCQiJCF5ISIiJCIlKzkhIiI3JCQiJFMpISIiJCIlXTshIiI3JCQiJCspISIiJCIlXTkhIiI3JCQiJCspISIiJCIlKzshIiI3JCQiJCspISIiJCIlKzkhIiI3JCQiJF0oISIiJCIlNTohIiI3JCQiJCF5ISIiJCIlPzkhIiI3JCQiJDUpISIiJCIlPzohIiI3JCQiJGcpISIiJCIlUzshIiI3JCQiJEkpISIiJCIlPzohIiI3JCQiJGcoISIiJCIlXTkhIiI3JCQiJDUoISIiJCIlZzkhIiI3JCQiJCF6ISIiJCIlPzohIiI3JCQiJGcoISIiJCIlNTohIiI3JCQiJCF6ISIiJCIlPzkhIiI3JCQiJFMpISIiJCIlIWYiISIiNyQkIiQ/KSEiIiQiJT86ISIiNyQkIiQrKSEiIiQiJVM6ISIiNyQkIiRTKSEiIiQiJSFbIiEiIjckJCIkXSkhIiIkIiU/OiEiIjckJCIkKykhIiIkIiVdOSEiIjckJCIkPykhIiIkIiVdOSEiIjckJCIkKykhIiIkIiVdOSEiIjckJCIkUykhIiIkIiUrOiEiIjckJCIkIXohIiIkIiU/OiEiIjckJCIkcSkhIiIkIiUrOiEiIjckJCIkcSchIiIkIiVdOyEiIjckJCIkIXohIiIkIiVTOiEiIjckJCIkXSkhIiIkIiUrOyEiIjckJCIkIXohIiIkIiU1OSEiIjckJCIkZykhIiIkIiU/OiEiIjckJCIkZykhIiIkIiVxOiEiIjckJCIkZykhIiIkIiU1OiEiIjckJCIkPykhIiIkIiUrOiEiIjckJCIkSSkhIiIkIiU/OiEiIjckJCIkISopISIiJCIlNTshIiI3JCQiJHEpISIiJCIlZzshIiI3JCQiJDUpISIiJCIlNTohIiI3JCQiJFMpISIiJCIlUzohIiI3JCQiJCspISIiJCIlcTkhIiI3JCQiJGcpISIiJCIlIVsiISIiNyQkIiQheiEiIiQiJSFcIiEiIjckJCIkNSkhIiIkIiUrOCEiIjckJCIkIXohIiIkIiUhWyIhIiI3JCQiJD8oISIiJCIlSTghIiI3JCQiJCF6ISIiJCIlKzkhIiI3JCQiJFMpISIiJCIlPzohIiI3JCQiJCF6ISIiJCIlIVwiISIiNyQkIiQ1KCEiIiQiJSs6ISIiNyQkIiRxKCEiIiQiJWc5ISIiNyQkIiQheiEiIiQiJSs8ISIiNyQkIiRdKSEiIiQiJVM7ISIiNyQkIiRxKSEiIiQiJSFwIiEiIjckJCIkUykhIiIkIiVTOiEiIjckJCIkSSkhIiIkIiVxOiEiIjckJCIkUykhIiIkIiUhXCIhIiI3JCQiJD8qISIiJCIlSTwhIiI3JCQiJCEqKSEiIiQiJUk8ISIiNyQkIiRJKiEiIiQiJSFvIiEiIjckJCIkZykhIiIkIiVJPCEiIjckJCIkXSkhIiIkIiU/OiEiIjckJCIkISopISIiJCIlcTshIiI3JCQiJF0pISIiJCIlNTohIiI3JCQiJCF5ISIiJCIlXTkhIiI3JCQiJCspISIiJCIlPzkhIiI3JCQiJF0pISIiJCIlSTohIiI3JCQiJF0pISIiJCIlKzohIiI3JCQiJGcpISIiJCIlcTohIiI3JCQiJHEoISIiJCIlPzohIiI3JCQiJFMpISIiJCIlZzshIiI3JCQiJEkpISIiJCIlPzwhIiI3JCQiJEkpISIiJCIlUzwhIiI3JCQiJEkpISIiJCIlKzwhIiI3JCQiJEkpISIiJCIlXTshIiI3JCQiJCEpKSEiIiQiJXE5ISIiNyQkIiRJKSEiIiQiJSF5IiEiIjckJCIkIXohIiIkIiUrOyEiIjckJCIlPzUhIiIkIiVTOyEiIjckJCIkKyohIiIkIiU/PSEiIjckJCIkPykhIiIkIiUrPSEiIi0lJ0xFR0VOREc2Iy0lKV9UWVBFU0VURzYjLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1LUkjbXNHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2I1ErUG9pbnR+ZGF0YTYiLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuNiIvJSVzaXplR1EjMTI2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EmZmFsc2U2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdub3JtYWw2Ii0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUmU1RZTEVHNiMlJlBPSU5URy0lJVZJRVdHNiQ7JCIkZyghIiIkIiQhKSkhIiI7JCIlKzkhIiIkIiUrPCEiIi0mJSZfQVhJU0c2IyIiIjYnLSUrX0dSSURMSU5FU0c2IyUoREVGQVVMVEctJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIi0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Jy0lK19HUklETElORVNHNiMlKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiItJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUqQVhFU1RJQ0tTRzYkIiM3IiM1LSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRK3B1bHNlLH5icG02Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxNjYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1US9wcmVzc3VyZSx+bW1IZzYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzE2NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii0lJUZPTlRHNiUlKERFRkFVTFRHJShERUZBVUxURyIjOyUrSE9SSVpPTlRBTEclKVZFUlRJQ0FMRy0lJlRJVExFRzYkLSUpX1RZUEVTRVRHNiMtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjYtSSNtc0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjYjUVFmaWcufjIwflN5c3RvbGEtUHVsc2V+Y29ycmVsYXRpb25+KHByaW1hcnl+ZGF0YSk2Ii8lJ2ZhbWlseUdRKmhlbHZldGljYTYiLyUlc2l6ZUdRIzE2NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJmZhbHNlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lMGZvbnRfc3R5bGVfbmFtZUdRJVRleHQ2Ii8lLG1hdGh2YXJpYW50R1Enbm9ybWFsNiItJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUnU1lNQk9MRzYkJS1fU09MSURDSVJDTEVHIiM9LSUlRk9OVEc2JCUqaGVsdmV0aWNhRyIjOw==
Statistics:-ScatterPlot(SDataNew_p,SDataNew_s,title = "fig. 21 Systola-Pulse correlation (reconstructed data)", legend = "Point data",view=[76..88,140..170],symbol=solidcircle,symbolsize=18,color=red,gridlines=true,font=[helvetica,16],tickmarks=[12,10],labeldirections=[horizontal,vertical],labels=[`pulse, bpm`,`pressure, mmHg`]);
LSUlUExPVEc2LC0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJD8oISIiLSUpQk9VTkRTX1lHNiMkIiQhUiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJXFlISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJTVLISIiLSUpQ0hJTERSRU5HNiItJSdDVVJWRVNHNiY3XHA3JCQiMXNxOS9pdCJmKSEjOSQiMiMqZWZISSRHYTshIzk3JCQiMSQpNEM8dF1DJikhIzkkIjFZPXojb2BYaiIhIzg3JCQiMXkkcFxSaHhXKSEjOSQiMlkmKXlXbU08aCIhIzk3JCQiMScpKVEqPlZcXCQpISM5JCIyTkpyayR6YCNlIiEjOTckJCIxKj4qM3lOJClRIykhIzkkIjJDI288KHpMJ1w6ISM5NyQkIjE9MU1CIjRoOykhIzkkIjFLLWEsVS1HOiEjODckJCIxdVpccSkpUUkiKSEjOSQiMkFAOVk/UXVeIiEjOTckJCIxXGFHN1MyOiIpISM5JCIyazBRPz4nKkdeIiEjOTckJCIwNlw+I1EjPjcpISM4JCIyRkQkKnBhRlxeIiEjOTckJCIxVS4kPlNeVzgpISM5JCIxI29VRzdZJz06ISM4NyQkIjE4KlxbKyM+XCIpISM5JCIybzc5R2ZASV8iISM5NyQkIjFoWGchbyp6eSIpISM5JCIxTGknSDI1PWAiISM4NyQkIjEtP2FiOyJ5QCkhIzkkIjIoKUduPi8hUlY6ISM5NyQkIjEjRyc9YVcjPUMpISM5JCIyKSpmbiQpeTwwYiIhIzk3JCQiMGslM3Y6cF4jKSEjOCQiMXlTNmRtV2A6ISM4NyQkIjFtPmdOMmVgIykhIzkkIjEvV3k9dSthOiEjODckJCIxKVteI0hmI3BDKSEjOSQiMnReQFYkWy5fOiEjOTckJCIxLC07WiZbSUMpISM5JCIyVVImNEdUKTNiIiEjOTckJCIxMVI2eCM0Q0MpISM5JCIyZUkqKXlTIXBdOiEjOTckJCIxNEFIJlsoSFEjKSEjOSQiMlleMWxJJVxcOiEjOTckJCIxLHBFWihRQkIpISM5JCIyQ2NLQDxxeGEiISM5NyQkIjFCJSp6RHA3SCMpISM5JCIyLXN5YHdMb2EiISM5NyQkIjEiUWorMGh6QSkhIzkkIjJRVTpOTzNsYSIhIzk3JCQiMSQ+JUdlbS9IIykhIzkkIjJ1NScqenUlelk6ISM5NyQkIjFGc2dqLitLIykhIzkkIjJ4Uyt4Qjd3YSIhIzk3JCQiMTs5QjglMzNCKSEjOSQiMEQhPVJbWFo6ISM3NyQkIjFkKlsnKkc3aEEpISM5JCIxJipRJlEtZWthIiEjODckJCIxKj5ZbV4mSD8jKSEjOSQiMTFOPlFjR1g6ISM4NyQkIjBcKycqeS5LQCkhIzgkIjFSaVE+XShRYSIhIzg3JCQiMSI0InpvP0U0IykhIzkkIjJCJ3k6aSl6SWEiISM5NyQkIi9DKkgiKlwjNCMpISM3JCIyJDRyVnFDN1Y6ISM5NyQkIjFfYFg/OF8xIykhIzkkIjJ5MCwlZjdcVTohIzk3JCQiMFFuMmUqSC0jKSEjOCQiMnZkXVIoW1xUOiEjOTckJCIxVFwpXGNkcz8pISM5JCIyaXhHZS4pKUdhIiEjOTckJCIxdCk0Mj46cUApISM5JCIyKD0iPV0hW21YOiEjOTckJCIxbW1xRiQzOkIpISM5JCIyMyIqKT4oKXpoXDohIzk3JCQiMWpEXlZoYlojKSEjOSQiMmBCVypIWTNhOiEjOTckJCIxPUZlKXpdXkIpISM5JCIxJjReXnhcQGIiISM4NyQkIjFvQ0JyRzIlPikhIzkkIjF1IzREXVVPYSIhIzg3JCQiMGtFPFEqM1ciKSEjOCQiMUBjUGUpSFBgIiEjODckJCIxIkg6IipbZUkzKSEjOSQiMmMneiopKVJwPF8iISM5NyQkIjFYSSNRTjRAKykhIzkkIjF6N041PjovOiEjODckJCIwcERhQXI5InohIzgkIjJWRikqSFBFUFsiISM5NyQkIjEkZncoUlJLL3khIzkkIjAwcSEpZSkzZjkhIzc3JCQiMWNMKWZpajBwKCEjOSQiMkFvLVhUP0dWIiEjOTckJCIxODJJRiY0bWooISM5JCIybioqUWo4IXBCOSEjOTckJCIwNncmKm9IPGsoISM4JCIxSyc9NlIjZUo5ISM4NyQkIjEqPSopKlwqeW5uKCEjOSQiMmVhV11Ea3JXIiEjOTckJCIxOE4qZU0lb1d4ISM5JCIyYkAkMzo2U3I5ISM5NyQkIjBMOGM+IylmJHkhIzgkIjJaai05ditRXSIhIzk3JCQiMXRKW3QmcCwleiEjOSQiMihSMENxQVtTOiEjOTckJCIxQiEzKzw7UTIpISM5JCIxXjQ3cDUvJ2UiISM4NyQkIjFeKVskSCdmV0EpISM5JCIyJFtGU1QuQ087ISM5NyQkIjAiKTN0SmVMSykhIzgkIjI9PUg4alNkbSIhIzk3JCQiMWRJY2lpLnMkKSEjOSQiMnUoKXohKilHK3Y7ISM5NyQkIjBVbWhBUURSKSEjOCQiMnU6OzcrLU1uIiEjOTckJCIxUEZVbGJxI1EpISM5JCIxVWhPZnktZzshIzg3JCQiMTxzdXMhM3RPKSEjOSQiMkNXJEdtISl6VDshIzk3JCQiMW0kM3hIQjFOKSEjOSQiMlcsOVopZUlAOyEjOTckJCIxOGQib1szWUopISM5JCIyczU4LTRAP2YiISM5NyQkIjFjNWUhUW9tRSkhIzkkIjI9JXprRiZSdmIiISM5NyQkIjEteCI9JHkpb0IpISM5JCIxLj09WU4+TjohIzg3JCQiMWUkeUdVX05BKSEjOSQiMk86aHRRblZfIiEjOTckJCIxQCQpZkAhNElBKSEjOSQiMnQ3IVElUXg0XyIhIzk3JCQiMSU0PmF4MV5CKSEjOSQiMiNbQk5hNUpEOiEjOTckJCIxLnNTcXlyTyMpISM5JCIxMm9lN3c5RzohIzg3JCQiMS84MHozWUkjKSEjOSQiMiRmQXovZGBIOiEjOTckJCIxOVprNjgrRCMpISM5JCIxZShIRFdNVGAiISM4NyQkIjFmRGxdIilmPyMpISM5JCIxPnh5Vk0oNGEiISM4NyQkIjFNPXBWciRRQikhIzkkIjIoSG9sSiRlKlw6ISM5NyQkIjFHOSopKVxiYEUpISM5JCIxKlJ2SVhHOWMiISM4NyQkIjEzbFMnM0w3SSkhIzkkIjIkbyIqM3F3ZXM6ISM5NyQkIjApPlA6RzlWJCkhIzgkIjJsKlthJz4uUWUiISM5NyQkIjFDaSFROiRSLCUpISM5JCIyNyFvITNGUy5nIiEjOTckJCIxUVYpenZvJG8lKSEjOSQiMUIwUCRwWCtpIiEjODckJCIwTDZEbFF1YSkhIzgkIjIkbyYqKikpR01FayIhIzk3JCQiMTF6TyE+XTpqKSEjOSQiMXI7QUshUWptIiEjODckJCIxJ3BELHE8OW4pISM5JCIyOmZuP1F4dW4iISM5NyQkIjF1I29ydFh2bSkhIzkkIjJFdlc4NUVobiIhIzk3JCQiMWMzVilcNzxrKSEjOSQiMmxtV0dqOSdvOyEjOS0lJ0xFR0VOREc2Iy0lKV9UWVBFU0VURzYjLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1LUkjbXNHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2I1ErUG9pbnR+ZGF0YTYiLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuNiIvJSVzaXplR1EjMTI2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EmZmFsc2U2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdub3JtYWw2Ii0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUmU1RZTEVHNiMlJlBPSU5URy0lJVZJRVdHNiQ7JCIkZyghIiIkIiQhKSkhIiI7JCIlKzkhIiIkIiUrPCEiIi0mJSZfQVhJU0c2IyIiIjYnLSUrX0dSSURMSU5FU0c2IyUoREVGQVVMVEctJSZDT0xPUkc2JiUkUkdCRyQiIiEhIiIkIiIhISIiJCIiISEiIi0lKkxJTkVTVFlMRUc2IyIiIi0lKlRISUNLTkVTU0c2IyIiIS0lLVRSQU5TUEFSRU5DWUc2IyQiIiEhIiItJiUmX0FYSVNHNiMiIiM2Jy0lK19HUklETElORVNHNiMlKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiItJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUqQVhFU1RJQ0tTRzYkIiM3IiM1LSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRK3B1bHNlLH5icG02Ii8lJ2ZhbWlseUdRKERFRkFVTFQ2Ii8lJXNpemVHUSMxNjYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiItSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1US9wcmVzc3VyZSx+bW1IZzYiLyUnZmFtaWx5R1EoREVGQVVMVDYiLyUlc2l6ZUdRIzE2NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2Ii0lJUZPTlRHNiUlKERFRkFVTFRHJShERUZBVUxURyIjOyUrSE9SSVpPTlRBTEclKVZFUlRJQ0FMRy0lJlRJVExFRzYkLSUpX1RZUEVTRVRHNiMtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjYtSSNtc0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjYjUVdmaWcufjIxflN5c3RvbGEtUHVsc2V+Y29ycmVsYXRpb25+KHJlY29uc3RydWN0ZWR+ZGF0YSk2Ii8lJ2ZhbWlseUdRKmhlbHZldGljYTYiLyUlc2l6ZUdRIzE2NiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJmZhbHNlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lMGZvbnRfc3R5bGVfbmFtZUdRJVRleHQ2Ii8lLG1hdGh2YXJpYW50R1Enbm9ybWFsNiItJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUnU1lNQk9MRzYkJS1fU09MSURDSVJDTEVHIiM9LSUlRk9OVEc2JCUqaGVsdmV0aWNhRyIjOw==
The readers may experiment with different lengths of digital filters in the expression (4.2) and make sure that the selected length=10 is optimal in certain sense. Someone could even change the family of wavelets, since Maple offers a choice of four such families: Daubechies, Symlet, Coiflet, and Battle-Lemarie [1].
The general conclusion that we do here, this is the statement that the Maple package Discrete Wavelet Transforms (DWT) is a perfect tool for analyzing, cleaning and compression of medical signals, if you know how and where to use it.
Legal Notice: The copyright for this application is owned by the author(s). Neither Maplesoft nor the author are responsible for any errors contained within and are not liable for any damages resulting from the use of this material. This application is intended for non-commercial, non-profit use only. Contact the author for permission if you wish to use this application in for-profit activities.
<Text-field style="Heading 2" layout="Heading 2">6. References</Text-field>
[1] Wavelets Examples and Applications. \302\251 - 2012, Maplesoft
[2] Y. Jin, E. Angelini , and A. Laine. "Wavelets in Medical Image Processing: Denoising, Segmentation, and Registration," in Handbook of Biomedical Image Analysis Vol 1- Segmentation Models - Part a, Ed.: D. L. W. Jasjit Suri, Swamy Laximinarayan Kluwer Academic Plenum Publishers, pp. 305-358, (2005).
[3] D.Donoho, and I. Johnstone. Adapting to Unknown Smoothness via Wavelet Shrinkage., Journal of American Statistics Association, Vol. 90, No. 432, pp. 200-1224, (1995).
[4] R. Coifman and D. Donoho.Translation-invariant De-noising. in Wavelets and Statistics, Antoniadis, A. and Oppenheim, G., Eds. New York NY: Springer-Verlag, 1995.
JSFH TTdSMApJN1JUQUJMRV9TQVZFLzQ0MjkxNTg0NDJYKiUpYW55dGhpbmdHNiI2IltnbCckJSEhISJbcCJbcDQwNjE4MDAwMDAwMDAwMDA0MDY0QTAwMDAwMDAwMDAwNDA2MjIwMDAwMDAwMDAwMDQwNjQwMDAwMDAwMDAwMDA0MDYxODAwMDAwMDAwMDAwNDA2MkUwMDAwMDAwMDAwMDQwNjFDMDAwMDAwMDAwMDA0MDYzMDAwMDAwMDAwMDAwNDA2NDgwMDAwMDAwMDAwMDQwNjMwMDAwMDAwMDAwMDA0MDYyMjAwMDAwMDAwMDAwNDA2MjQwMDAwMDAwMDAwMDQwNjMwMDAwMDAwMDAwMDA0MDYyRTAwMDAwMDAwMDAwNDA2MUMwMDAwMDAwMDAwMDQwNjNFMDAwMDAwMDAwMDA0MDYzMDAwMDAwMDAwMDAwNDA2MzQwMDAwMDAwMDAwMDQwNjI4MDAwMDAwMDAwMDA0MDYzMDAwMDAwMDAwMDAwNDA2MjIwMDAwMDAwMDAwMDQwNjIyMDAwMDAwMDAwMDA0MDYyMjAwMDAwMDAwMDAwNDA2MkMwMDAwMDAwMDAwMDQwNjMwMDAwMDAwMDAwMDA0MDYyQzAwMDAwMDAwMDAwNDA2NEEwMDAwMDAwMDAwMDQwNjM0MDAwMDAwMDAwMDA0MDY0MDAwMDAwMDAwMDAwNDA2MUEwMDAwMDAwMDAwMDQwNjMwMDAwMDAwMDAwMDA0MDYzQTAwMDAwMDAwMDAwNDA2MkUwMDAwMDAwMDAwMDQwNjJDMDAwMDAwMDAwMDA0MDYzMDAwMDAwMDAwMDAwNDA2NDIwMDAwMDAwMDAwMDQwNjRDMDAwMDAwMDAwMDA0MDYyRTAwMDAwMDAwMDAwNDA2MzQwMDAwMDAwMDAwMDQwNjI2MDAwMDAwMDAwMDA0MDYyODAwMDAwMDAwMDAwNDA2MkEwMDAwMDAwMDAwMDQwNjA0MDAwMDAwMDAwMDA0MDYyODAwMDAwMDAwMDAwNDA2MEEwMDAwMDAwMDAwMDQwNjE4MDAwMDAwMDAwMDA0MDYzMDAwMDAwMDAwMDAwNDA2MkEwMDAwMDAwMDAwMDQwNjJDMDAwMDAwMDAwMDA0MDYyNDAwMDAwMDAwMDAwNDA2NTQwMDAwMDAwMDAwMDQwNjQ4MDAwMDAwMDAwMDA0MDY1MjAwMDAwMDAwMDAwNDA2MzQwMDAwMDAwMDAwMDQwNjNBMDAwMDAwMDAwMDA0MDYyQTAwMDAwMDAwMDAwNDA2NUEwMDAwMDAwMDAwMDQwNjVBMDAwMDAwMDAwMDA0MDY1MDAwMDAwMDAwMDAwNDA2NUEwMDAwMDAwMDAwMDQwNjMwMDAwMDAwMDAwMDA0MDY0RTAwMDAwMDAwMDAwNDA2MkUwMDAwMDAwMDAwMDQwNjIyMDAwMDAwMDAwMDA0MDYxQzAwMDAwMDAwMDAwNDA2MzIwMDAwMDAwMDAwMDQwNjJDMDAwMDAwMDAwMDA0MDYzQTAwMDAwMDAwMDAwNDA2MzAwMDAwMDAwMDAwMDQwNjRDMDAwMDAwMDAwMDA0MDY1ODAwMDAwMDAwMDAwNDA2NUMwMDAwMDAwMDAwMDQwNjU0MDAwMDAwMDAwMDA0MDY0QTAwMDAwMDAwMDAwNDA2MjYwMDAwMDAwMDAwMDQwNjY0MDAwMDAwMDAwMDA0MDY0MDAwMDAwMDAwMDAwNDA2NDgwMDAwMDAwMDAwMDQwNjZDMDAwMDAwMDAwMDA0MDY2ODAwMDAwMDAwMDAwRiU= TTdSMApJN1JUQUJMRV9TQVZFLzQ0MjkyMDcyMzRYKiUpYW55dGhpbmdHNiI2IltnbCckJSEhISJbcCJbcDQwNTRDMDAwMDAwMDAwMDA0MDUzNDAwMDAwMDAwMDAwNDA1NTAwMDAwMDAwMDAwMDQwNTQwMDAwMDAwMDAwMDA0MDU0QzAwMDAwMDAwMDAwNDA1MzQwMDAwMDAwMDAwMDQwNTRDMDAwMDAwMDAwMDA0MDUzODAwMDAwMDAwMDAwNDA1MzgwMDAwMDAwMDAwMDQwNTQ0MDAwMDAwMDAwMDA0MDU0QzAwMDAwMDAwMDAwNDA1M0MwMDAwMDAwMDAwMDQwNTNDMDAwMDAwMDAwMDA0MDUzNDAwMDAwMDAwMDAwNDA1NTQwMDAwMDAwMDAwMDQwNTQwMDAwMDAwMDAwMDA0MDUzMDAwMDAwMDAwMDAwNDA1QTgwMDAwMDAwMDAwMDQwNUI0MDAwMDAwMDAwMDA0MDU0QzAwMDAwMDAwMDAwNDA1NDgwMDAwMDAwMDAwMDQwNTUwMDAwMDAwMDAwMDA0MDU2MDAwMDAwMDAwMDAwNDA1NTAwMDAwMDAwMDAwMDQwNTUwMDAwMDAwMDAwMDA0MDU0NDAwMDAwMDAwMDAwNDA1NEMwMDAwMDAwMDAwMDQwNTM0MDAwMDAwMDAwMDA0MDU0QzAwMDAwMDAwMDAwNDA1MzgwMDAwMDAwMDAwMDQwNTU4MDAwMDAwMDAwMDA0MDU0QzAwMDAwMDAwMDAwNDA1NTgwMDAwMDAwMDAwMDQwNTRDMDAwMDAwMDAwMDA0MDU5NDAwMDAwMDAwMDAwNDA1OTQwMDAwMDAwMDAwMDQwNTQ0MDAwMDAwMDAwMDA0MDU0ODAwMDAwMDAwMDAwNDA1NjAwMDAwMDAwMDAwMDQwNTUwMDAwMDAwMDAwMDA0MDU1MDAwMDAwMDAwMDAwNDA1NUMwMDAwMDAwMDAwMDQwNTVDMDAwMDAwMDAwMDA0MDU2NDAwMDAwMDAwMDAwNDA1NjgwMDAwMDAwMDAwMDQwNTRDMDAwMDAwMDAwMDA0MDUzMDAwMDAwMDAwMDAwNDA1MzAwMDAwMDAwMDAwMDQwNTY4MDAwMDAwMDAwMDA0MDUyQzAwMDAwMDAwMDAwNDA1M0MwMDAwMDAwMDAwMDQwNTQ0MDAwMDAwMDAwMDA0MDUzMDAwMDAwMDAwMDAwNDA1M0MwMDAwMDAwMDAwMDQwNTQ0MDAwMDAwMDAwMDA0MDU1MDAwMDAwMDAwMDAwNDA1NTAwMDAwMDAwMDAwMDQwNTQ0MDAwMDAwMDAwMDA0MDU0MDAwMDAwMDAwMDAwNDA1MzQwMDAwMDAwMDAwMDQwNTU0MDAwMDAwMDAwMDA0MDU0NDAwMDAwMDAwMDAwNDA1NDAwMDAwMDAwMDAwMDQwNTU0MDAwMDAwMDAwMDA0MDU1MDAwMDAwMDAwMDAwNDA1NEMwMDAwMDAwMDAwMDQwNTU4MDAwMDAwMDAwMDA0MDU2NDAwMDAwMDAwMDAwNDA1MzQwMDAwMDAwMDAwMDQwNTNDMDAwMDAwMDAwMDA0MDU2ODAwMDAwMDAwMDAwNDA1NDAwMDAwMDAwMDAwMDQwNTY0MDAwMDAwMDAwMDA0MDU1NDAwMDAwMDAwMDAwNDA1M0MwMDAwMDAwMDAwMDQwNTUwMDAwMDAwMDAwMDA0MDUzQzAwMDAwMDAwMDAwNDA1NkMwMDAwMDAwMDAwMDQwNTY0MDAwMDAwMDAwMDA0MDU3NDAwMDAwMDAwMDAwRiU= TTdSMApJN1JUQUJMRV9TQVZFLzQ0MjkyMDczMDZYKiUpYW55dGhpbmdHNiI2IltnbCckJSEhISJbcCJbcDQwNTM4MDAwMDAwMDAwMDA0MDU1MDAwMDAwMDAwMDAwNDA1NDAwMDAwMDAwMDAwMDQwNTQwMDAwMDAwMDAwMDA0MDU0MDAwMDAwMDAwMDAwNDA1MkMwMDAwMDAwMDAwMDQwNTM4MDAwMDAwMDAwMDA0MDU0NDAwMDAwMDAwMDAwNDA1NTgwMDAwMDAwMDAwMDQwNTRDMDAwMDAwMDAwMDA0MDUzMDAwMDAwMDAwMDAwNDA1MUMwMDAwMDAwMDAwMDQwNTNDMDAwMDAwMDAwMDA0MDUzMDAwMDAwMDAwMDAwNDA1M0MwMDAwMDAwMDAwMDQwNTUwMDAwMDAwMDAwMDA0MDU0ODAwMDAwMDAwMDAwNDA1NDAwMDAwMDAwMDAwMDQwNTUwMDAwMDAwMDAwMDA0MDU1NDAwMDAwMDAwMDAwNDA1NDAwMDAwMDAwMDAwMDQwNTQ4MDAwMDAwMDAwMDA0MDU0MDAwMDAwMDAwMDAwNDA1NTAwMDAwMDAwMDAwMDQwNTNDMDAwMDAwMDAwMDA0MDU1QzAwMDAwMDAwMDAwNDA1MEMwMDAwMDAwMDAwMDQwNTNDMDAwMDAwMDAwMDA0MDU1NDAwMDAwMDAwMDAwNDA1M0MwMDAwMDAwMDAwMDQwNTU4MDAwMDAwMDAwMDA0MDU1ODAwMDAwMDAwMDAwNDA1NTgwMDAwMDAwMDAwMDQwNTQ4MDAwMDAwMDAwMDA0MDU0QzAwMDAwMDAwMDAwNDA1NjQwMDAwMDAwMDAwMDQwNTVDMDAwMDAwMDAwMDA0MDU0NDAwMDAwMDAwMDAwNDA1NTAwMDAwMDAwMDAwMDQwNTQwMDAwMDAwMDAwMDA0MDU1ODAwMDAwMDAwMDAwNDA1M0MwMDAwMDAwMDAwMDQwNTQ0MDAwMDAwMDAwMDA0MDUzQzAwMDAwMDAwMDAwNDA1MjAwMDAwMDAwMDAwMDQwNTNDMDAwMDAwMDAwMDA0MDU1MDAwMDAwMDAwMDAwNDA1M0MwMDAwMDAwMDAwMDQwNTFDMDAwMDAwMDAwMDA0MDUzNDAwMDAwMDAwMDAwNDA1M0MwMDAwMDAwMDAwMDQwNTU0MDAwMDAwMDAwMDA0MDU1QzAwMDAwMDAwMDAwNDA1NTAwMDAwMDAwMDAwMDQwNTRDMDAwMDAwMDAwMDA0MDU1MDAwMDAwMDAwMDAwNDA1NzAwMDAwMDAwMDAwMDQwNTY0MDAwMDAwMDAwMDA0MDU3NDAwMDAwMDAwMDAwNDA1NTgwMDAwMDAwMDAwMDQwNTU0MDAwMDAwMDAwMDA0MDU2NDAwMDAwMDAwMDAwNDA1NTQwMDAwMDAwMDAwMDQwNTM4MDAwMDAwMDAwMDA0MDU0MDAwMDAwMDAwMDAwNDA1NTQwMDAwMDAwMDAwMDQwNTU0MDAwMDAwMDAwMDA0MDU1ODAwMDAwMDAwMDAwNDA1MzQwMDAwMDAwMDAwMDQwNTUwMDAwMDAwMDAwMDA0MDU0QzAwMDAwMDAwMDAwNDA1NEMwMDAwMDAwMDAwMDQwNTRDMDAwMDAwMDAwMDA0MDU0QzAwMDAwMDAwMDAwNDA1NjAwMDAwMDAwMDAwMDQwNTRDMDAwMDAwMDAwMDA0MDUzQzAwMDAwMDAwMDAwNDA1OTgwMDAwMDAwMDAwMDQwNTY4MDAwMDAwMDAwMDA0MDU0ODAwMDAwMDAwMDAwRiU= TTdSMApJN1JUQUJMRV9TQVZFLzQ0MjkyMDczNzhYKiUpYW55dGhpbmdHNiI2IltnbCcjJSEhISIrIiszRjk0MDIwMDM5NjRBMkJDM0Y5NTlCQjgxQzE2MkRFOEJGQzY3MTI2OUUxQ0Q0RjdCRjkxMDAyNDc1MzY5NjI4M0ZFNDQ5OEU0MDRDQzNDQ0JGRTcyNjI3RTI3RUJBNUEzRkM5ODVEQkJGQUFCNEQ1M0ZBNDA5Njc4MUQwNUI5NzNGOUUzQTVCODA3QjVFNURCRjlCRkQzMzJBOTA3RDMwRiU= TTdSMApJN1JUQUJMRV9TQVZFLzQ0MjkyMjA2MDJYKiUpYW55dGhpbmdHNiI2IltnbCcjJSEhISIrIiszRjlCRkQzMzJBOTA3RDMwM0Y5RTNBNUI4MDdCNUU1REJGQTQwOTY3ODFEMDVCOTczRkM5ODVEQkJGQUFCNEQ1M0ZFNzI2MjdFMjdFQkE1QTNGRTQ0OThFNDA0Q0MzQ0MzRjkxMDAyNDc1MzY5NjI4QkZDNjcxMjY5RTFDRDRGN0JGOTU5QkI4MUMxNjJERTgzRjk0MDIwMDM5NjRBMkJDRiU=