مرگ ناگهانی قلبی (SCD) نتیجه تخریب شدید عملکرد قلبی است که سبب از بین رفتن سیستم قلبی در افراد می شود. وقتی این اتفاق رخ می دهد ، خون دیگر نمی تواند برای مدتی به قسمت‌های مختلف بدن پمپ شود [1] این رخداد، هر ساله، عامل مرگ تقریبا 250000 تا300000  انسان در آمریکاست [2]. همچنین در آمریکای شمالی و اروپا این تعداد از میان هر 100000 نفر به 50 الی 100 نفر می رسد [3-8]. با استفاده از تجهیزات پزشکی، از قبیل دیفیبریلاتور تعداد این نوع مرگ‌ها می تواند کاهش یابد. با وجود این، راه‌های مناسبی برای پیش بینی مرگ ناگهانی قلبی که پزشکان بتوانند از طریق آن تصمیمات مناسبی را برای بیماران در معرض خطر بگیرند، وجود ندارد. تنها 2-1% از بیمارانی که دچار SCD شده اند می‌توانند در بیرون بیمارستان جان سالم به در ببرند [9]. این اتفاق، غالبا با یک تاکی آریتمی شدید، از قبیل: تاکیکاردی بطنی (VT)، لرزش بطنی (VFL) ویا فیبریلاسیون بطنی، شروع می شود. در موارد نادری نیز بر اثر برادی آریتمی است [1،10]. مطالعات وتحقیقاتی که در زمینه تشخیص SCD با الکترو کاردیوگرام صورت گرفته، به تعداد اندکی محدود می شود که غالبا از طریق استخراج ویژگی‌های کلاسیک که تمایز درخور توجهی نیز ایجاد نمی کنند، به بررسی این موضوع پرداخته اند [1،11]. آنها شواهدی دارند که نشان می دهد HRV در افرادی که تجربه SCD دارند، کم و در افراد جوان سالم زیاد است [1،11]. در تحقیق دیگری نشان داده شده است که پیک های تنفسی در تغییرات ضربان قلبی (HRV) در بیماران SCD، در یک هفته قبل از مرگشان ناپدید شده اند[12]. اگر مرگ ناگهانی قلبی، قبل از رخداد تشخیص داده شود، می تواند به نجات جان بیمار قلبی بینجامد.

منظور از HRV تغییرات در فواصل زمانی بین ضربان‌های متوالی قلب است. به عبارت دیگر، زمان بین دو موج R در دو سیکل متوالی قلبی در سیگنال الکترو کاردیوگرام را بازه زمانی RR گوئیم. تغییرات در این سری زمانی (HRV) می‌تواند به عنوان یک علامت خطر برای مرگ و میر پس از سکته قلبی در نظر گرفته شود [1]. در این مقاله از سیگنال تغیرات ضربان قلبی (HRV) به منظور تشخیص SCD استفاده شده است. در این تحقیق به منظور پیش بینی مرگ ناگهانی قلبی به استخراج ویژگی‌های معمول خطی و همچنین، حوزه زمان فرکانس پرداخته شده است. از طرفی، با توجه به اینکه تاکنون روش‌های خطی مختلفی برای تحلیل سیگنال HRV به کار گرفته شده است، اخیراً اثبات شده که روش‌های پردازش غیر خطی می توانند اطلاعات بیشتری نسبت به روش‌های خطی معمول (تحلیل‌های حوزه زمان و حوزه فرکانس ) فراهم کنند و مکمل خوبی برای آنها باشند. تحلیل غیر خطی سیگنال HRV به دو دلیل عمده بسیار مورد توجه قرار گرفته است: دلیل اول ذات غیر خطی سیگنال مشاهده شده از قلب است که به صورت یک نوسانگر دینامیک غیر خطی عمل می کند و دلیل دوم، لزوم به دست آوردن دانش کافی درباره این پدیده واقعی است. از این رو، به تحلیل‌های غیر خطی نیز پرداخته شده است .

هدف ما در این مقاله، انجام تحلیل‌های خطی کلاسیک، زمان- فرکانس و غیر خطی روی سیگنال HRV  به دست آمده از افراد سالم و افراد در معرض مرگ ناگهانی قلبی است تا بتوانیم با استفاده از یک بردار ویژگی ترکیبی ، دو گروه افراد سالم وافراد در آستانه مرگ قلبی را با استفاده از طبقه بندی کننده شبکه عصبی با دقت خوبی از یکدیگر تفکیک کرده، از طریق آشکار کردنِ افزایش احتمال خطر، به پیش بینی مرگ ناگهانی قلبی بپردازیم .

تحقیقات انجام شده در این زمینه نشان داده اند که روش‌های خطی زمان و فرکانس در جهت پیش بینی مرگ ناگهانی قلبی، فاقد توانمندی لازم هستند[1،13].

در تحقیقات اولیه در این زمینه، ما نشان دادیم که روش‌های زمان- فرکانس نسبت به روش‌های خطی کلاسیک توانمندی بمراتب بیشتری در تفکیک بین افراد سالم و ریسک پذیر دارند [13].

در ارزیابی صورت گرفته بین دو روش زمان-فرکانس و روش کلاسیک به ترتیب صحت تفکیک پذیری 99.16% و 74.36% برای 1 دقیقه سیگنال ECG قبل از واقعه از طریق طبقه بندی کننده MLP به دست آمده است [13].

در ادامه، پس از استخراج سیگنال HRV  از سیگنال ECG  ویژگی های خطی استخراج شده و سپس تبدیل ویگنرویل اعمال کرده و به استخراج ویژگی در حوزه زمان – فرکانس پرداخته شده و در نهایت نیز ویژگی‌های غیر خطی استخراج شده است. در مرحله بعد بهترین ترکیب ویژگی ها بر اساس بیشترین ایجاد تمایز بین دو کلاس انتخاب شده، سپس با اعمال PCA  به بردار ویژگی ترکیبی، ابعاد ویژگی کاهش یافته و در نهایت، از طریق شبکه عصبی MLP افراد سالم و افرادی که دچار  SCD شده اند، دسته بندی می شوند. این دسته بندی در چهار مرحله ، برای دقایق اول و دوم و سوم و چهارم پیش از SCD  صورت گرفته است. در شکل (1) خلاصه مراحل الگوریتم نشان داده شده است .

شکل(1) :فلوچارتکلیالگوریتمارائهشده.

1- روش‌ها و ابزارها

الگوریتم ارائه شده بر روی مجموعه دادهایی متشکل از 35 نفر از افرادی که دچار مرگ ناگهانی قلبی شده اند (شامل 19 مرد و 16 زن و ECG با فرکانس نمونه برداری 256 نمونه در ثانیه) و 35 فرد سالم بدون سابقه قلبی (شامل 17 زن و 18 مرد و ECG با فرکانس نمونه برداری 128نمونه در ثانیه )، ارزیابی

شده است . این داده ها از پایگاه داده  MIT-BIH Sudden Cardiac Death Holter database  , Normal Sinus Rhythm database   به دست آمده است [14،15]. در مواردی که از هر بیمار دو کانال در دسترس بوده است، هر کانال به عنوان یک مشاهده (بیمار) در نظر گرفته شده است.  

1-1- پیش پردازش سیگنال ECG

از افرادی که خطر بالای مرگ ناگهانی قلبی داشته‌اند، به مدت طولانی سیگنالECG  ثبت شده است .این افراد غالباً  کسانی بوده اند که یا تجربه سکته قلبی را داشته اند و یا  به علت یک تاکی آریتمی شدید در سیستم قلبی - عروقی مستعد برایSCD  بوده اند، که در نهایت نیز به این حادثه دچار شده اند. از سیگنال ECG این افراد دقیقا پیش از رخ دادن واقعه بازه های زمانی 1 دقیقه ای جدا کرده و به ترتیب به نام دقیقه اول، دقیقه دوم و دقیقه سوم و دقیقه چهارم نام گذاری شده‌اند. در شکل (2) سیگنال الکتروکاردیوگرام یک مرد 34 ساله که در معرض مرگ ناگهانی قلبی است، از 3 دقیقه قبل تا لحظاتی پس از وقوع مرگ قلبی، ملاحظه می کنید .

شکل(2) : سیگنال ECG فرد در آستانه مرگ ناگهانی قلبی از 3 دقیقه قبل از واقعه تا لحظاتی بعد از واقعه .

شکل (3) : سیگنال ECG بیمار در لحظه مرگ ناگهانی قلبی [15]

شکل ظاهری سیگنال برای فردی که دچار SCD شده (قبل از حادثه) تفاوت چندانی با سیگنال فرد سالم ندارد. در شکل (3) یک سیگنال ECG فردی را که دچار SCD شده است، از چند ثانیه قبل تا هنگام وقوع حادثه ملاحظه می کنید [15].

در این مرحله از  ECG افرادی که دچار مرگ ناگهانی قلبی شده اند، دقیقاً دقیقه قبل از رخ دادن حادثه، وهمچنین افراد سالم، به مدت 1 دقیقه سیگنال ECG  اخذ می شود. سپس به آشکارسازی موج R و درنهایت استخراج HRV پرداخته شده است. در این مرحله با استفاده از الگوریتم [16] Pan-Tampkins موج R را به دست آورده و سپس از کنار هم قرار دادن فواصل RR متوالی، سیگنال HRV را تشکیل می دهیم. به این ترتیب، سیگنال های HRV پیش پردازش شده برای استخراج ویژگی آماده می‌شوند. در شکل‌های(4)و (5) سیگنال ECG  و HRV فرد سالم وSCD  نشان داده شده است. چنانچه قبلاً اشاره شد و در شکل (4)و (5) مشخص است سیگنال ECG فرد سالم و ریسک پذیر تفاوت چندانی با یکدیگر ندارند .

(الف)

(ب)

شکل (4) : الف) یک دقیقه سیگنالECG  فرد سالم ب) سیگنال HRV استخراج شده از شکل الف .

(الف)

(ب)

شکل (5) : الف) یک دقیقه سیگنالECG  فرد در آستانه مرگ قلبی دقیقاً قبل از رخداداً واقعه؛ ب) سیگنال HRV استخراج شده از شکل الف .

1-2- استخراج ویژگی های کلاسیک

در این مرحله ویژگی های معمول خطی در حوزه زمان و فرکانس استخراج شده اند. این ویژگی ها شامل پنج ویژگی در حوزه زمان و چهار ویژگی در حوزه فرکانس است.

1-2-1- ویژگی های حوزه زمان

ویژگی های حوزه زمان شامل مجموعه ای از ویژگی‌های آماری است [17،18،19،20]. این ویژگی‌ها عبارتند از : MNN ( میانگین فواصل RR موجود در هر قطعه از سیگنال HRV) - SDNN (انحراف معیار فواصل RR موجود در هر قطعه از سیگنال HRV) - RMSSD (جذر میانگین مربعات اختلاف فواصل RR موجود در هر قطعه از سیگنال HRV) - SDSD (انحراف معیار اختلاف بین فواصل RR موجود در هر قطعه از سیگنال HRV) - PNN50 (درصد تعداد فواصل RR موجود در هر قطعه از سیگنال HRV که اختلافشان بیشتر از 50 میلی ثانیه است). توزیع فضای برخی از این ویژگی‌ها در شکل (6) نشان داده شده است. به طوری که در این شکل دیده می‌شود، این ویژگی ها می توانند تفکیک پذیری مناسبی بین دو گروه سالم وSCD  ایجاد کنند.

شکل (6) : توزیع فضای ویژگی میانگین (MNN)در محور افقی و انحراف معیار(Std NN)در محور عمودی نمایش داده شده است .

1-2-2- ویژگی های حوزه فرکانس

در این مرحله با استفاده از روش تخمین طیف ((PSD، انرژی سیگنال در باند فرکانسی خیلی پائین (0.04-0.003هرتز) و باند فرکانسی پایین (0.15-0.04هرتز) و باند فرکانسی بالا (0.4-0.15هرتز) استخراج شده است [21,22,18]. مطالعات قبلی نشان داده است که فرکانس‌های بالا در طیف توان سیگنال HRV نشان دهنده فعالیت بخش پاراسمپاتیک سیستم عصبی و همچنین، فرکانس‌های پایین نشان دهنده فعالیت بخش سمپاتیک سیستم عصبی خودکار کنترل کننده نرخ ضربان قلب است [17،20]. از این رو، نسبت انرژی سیگنال در باند فرکانسی پایین (LF) به انرژی سیگنال در باند فرکانسی بالا (HF) می تواند به عنوان یک ویژگی که تعادل سمپاتیک و پاراسمپاتیک را ارزیابی می کند، استفاده شود [17،20]. طیف توان سیگنال(PSD) با استفاده از روش پارامتریک Burg محاسبه شده است. در شکل (7) طیف توان سیگنال و نیز باندهای فرکانسی مختلف آن نشان داده شده است.

شکل (7) : الف) سیگنال HRV استخراج شده؛ ب) سیگنال HRV در حوزه فرکانس و تفکیک انرژی در هر باند فرکانسی.

1-3- آنالیز زمان – فرکانس

یکی از روش های تحلیل HRV استفاده از تبدیل های زمان- فرکانس است . که به سه دسته مهم طبقه بندی می‌شوند:

1-روش های زمان –فرکانس غیر پارامتری، شامل تبدیل فوریه زمان کوتاه [21،22] و ویولت [23،24].

2- نمایندگان زمان – فرکانس درجه دوم غیر پارامتری، شامل توزیع ویگنرویل[25،28].

3 – روش‌های تغییر زمان زمان [29،31]. پارامتری براساس مدل های AR با ضرایب تغییر.

در این تحقیق، روش توزیع ویگنرویل (SPWVD) برای تحلیل HRV انتخاب شده است .این روش تفکیک پذیری زمان- فرکانس بهتری نسبت به روش های غیر پارامتریک ارائه می دهد و دارای کنترل های مستقلی برای زمان و فرکانس است و همچنین، هنگام وقوع تغییرات سریع، تخمین توان با واریانس کمتری نسبت به روش‌های پارامتریک ارائه می شود [29].

توزیع ویگنرویل از سیگنال گسسته x[n] این گونه تعریف می شود :

که  m,n به ترتیب شاخص های زمان و فرکانس گسسته هستند.h(k)  پنجره متحرک فرکانسی به طول 2N-1 و g(p) پنجره متحرک زمانی به طول 2M-1 و rx(n,k) تابع همبستگی لحظه ای است که به فرم زیر تعریف می شود:

در شکل (8) سیگنال HRV فرد پس از اعمال تبدیل ویگنرویل نمایش داده شده است .

شکل(8) : نمایش HRV استخراج شده یک فرد SCD، در حوزه زمان - فرکانس .

1-3-1- استخراج ویژگی های   TF

در این مرحله سیگنال HRV را که بر حسب نمونه است، به پنج بازه زمانی مساوی تقسیم می کنیم که هر پنجره بیانگر 15 ثانیه در حوزه زمان است و میانگین انرژی را در هر پنجره زمانی به دست می آوریم.

Max_w : مقدار ماکزیمم انرژی بین پنجره های زمانی

Min_w : مقدار مینیمم انرژی بین پنجره های زمانی

: Dif_w  تفاوت میان بیشترین مقدار انرژی و کمترین آنها در بین پنجره ها

Std_w: انحراف معیار بین انرژی پنجره های زمانی

به همین ترتیب سیگنال به دست آمده در حوزه TF رااز سمت محور فرکانس به 3 بازه تقسیم می‌کنیم .

 : E_vlf مجموع انرژی سیگنال در باند فرکانسی خیلی پایین (0.04-0.003) هرتز تقسیم بر (0.003-0.04).

E_lf : مجموع انرژی سیگنال در باند فرکانسی پایین (0.15-0.04)هرتز تقسیم بر (0.15-0.04).

E_hf : مجموع انرژی سیگنال در باند فرکانسی بالا (0.4-0.15) هرتز تقسیم بر(0.15-0.4).

: F_vlf میانگین انرژی سیگنال در باند فرکانسی خیلی پایین (0.04-0.003هرتز) .

F_lf : میانگین انرژی سیگنال در باند فرکانسی پایین (0.15-0.04 هرتز) .

F_hf : میانگین انرژی سیگنال در باند فرکانسی بالا (0.4-0.15 هرتز) .

همچنین، ویژگی دیگری به نام مشتق مرتبه اول تعریف می‌کنیم که حاصل تفاضل میانگین انرژی هرپنجره با پنجره قبلی است، که در مورد پنجره اول (15 ثانیه اول) اختلافش با آخرین پنجره در بازه دقیقه دوم سنجیده می شود   W_dif = |w_n - w_n-1|. نتایج  بررسی ویژگی ها در بازه های زمانی 15 ثانیه ای نشان می دهد که برای موارد SCD تغییرات ویژگی ها از یک پنجره به پنجره بعدی بسیار بارزتر است. به همین خاطر، ویژگی مشتق مرتبه اول را تعریف کرده ایم.

1-4- استخراج ویژگی های غیر خطی

از آنجا که سیستم قلبی - عروقی خیلی پیچیده تر از یک سیستم خطی عمل کرده، رفتار غیر ایستا از خود نشان می‌دهد، علاوه بر ویژگی های زمان-فرکانس، دو تحلیل غیر خطی که نشان دهنده خصوصیات دینامیکی آشوبناک در سیگنال HRV هستند، برای طبقه بندی بین افراد سالم و افراد در معرض مرگ قلبی به منظور پیش بینی مرگ ناگهانی قلبی استفاده شده‌اند. این دو تحلیل عبارتند: از معیار پوانکاره که سه ویژگی SD1,SD2,SD1/SD2 از آن استخراج شده است و تحلیل نوسان دترند شده (DFA)که ویژگی α از آن استخراج شده است.

1-4-1- معیار پوآنکاره  

نمودار پوانکاره تکنیک نسبتا جدیدی برای تحلیل دینامیک‌های غیر خطی همچون HRV است. هر نقطه به روی این نمودار به صورت i = 1,2,3,…,N  (RR(i),RR(i+1)) , کهN  طول سیگنال است. مشخص می‌شود[32]. از نظر آماری این نمودار همبستگی بین فواصل متوالی R-R  را به صورت گرافیکی نمایش می‌دهد، اما مفهوم اصلی و مهم آن، این است که این نمودار یک فضای حالت دو بعدی ساخته شده از فواصل متوالی است، که بیانگر دینامیک غیر خطی آن است. این نمودار اطلاعات مفیدی هم از نوسان‌های کوتاه مدت و هم از نوسان‌های بلند مدت در اختیار ما قرار می دهد. نمودار پوانکاره به صورت کمی با محاسبه انحراف معیار فواصل RR(i) با خطوط y=x و y = -x+2RR_m تحلیل می شود که RR_m متوسط RR(i) ها است. انحراف معیارها SD1و SD2نامیده می شوند که SD1 مربوط به تغییرات سریع ضربان به ضربان در دادگان است که عمدتا به RSA مربوطمی شود، در حالی که SD2 تغییرات بلند مدت RR(i) ها را توصیف می کند. نسبت SD2/SD1 را نیز می توان برای توصیف ارتباط بین این مؤلفه ها محاسبه کرد[23]. در شکل(9) نمودار پوانکاره را برای یک فرد ریسک پذیر و یک فرد سالم مشاهده می کنید.

(الف)

(ب)

شکل (9): الف)نمودار پوانکاره مربوط به فرد SCD؛ ب) نموداره پوانکاره مربوط به فرد سالم.

1-4-2- تحلیل نوسان دترند شده (DFA)

از تحلیل DFA برای کمی کردن خصوصیات فرکتال سیگنال فواصل RR کوتاه مدت استفاده می شود [34,35,18]. در ابتدا ، از سری‌های زمانی RR با استفاده از رابطه (3) انتگرال (سیگما) گرفته می شود :

که در آن y(k)، k  امین مقدار سری‌های سیگما گرفته شده است، RR(i) ، i امین فاصله  RR است و RR_avg میانگین فواصل RR روی کل سری است. سپس سری زمانی سیگما گرفته شده به پنجره های با طول مساوی n تقسیم می شود. در هر پنجره با طول n ، یک خط با معیار حداقل مربع خطا روی دیتای RR منطبق می شود که ترند را در آن پنجره نشان می دهد. اگر این تکه خط‌های راست را با y_n(k) نشان دهیم، طبق رابطه (4) با کم کردن آن از سری زمانی سیگما گرفته شده در هر پنجره یک سری زمانی ترند شده به دست می آید :

این محاسبه روی طول های مختلف پنجره تکرار می‌شود تا ارتباط بین F(n)و طول پنجره یعنی n به دست آید. شیب خط مربوط به  log(F(n))نسبت به log(n) با عنوان فاکتور خود همانندی به عنوان یک ویژگی در طبقه بندی استفاده می‌شود.

1-5- انتخاب ویژگی

در این مرحله سعی شده است بهترین ترکیب ویژگی برای ایجاد بیشترین تمایز بین کلاس‌ها برگزیده شود. به این ترتیب که در ابتدا با هر یک از ویژگی ها به طور جداگانه طبقه بندی انجام و براساس بیشترین میزان صحت تفکیک پذیری بهترین ویژگی انتخاب می شود. سپس این ویژگی را با تک تک ویژگی ها ترکیب کرده، بهترین ترکیب دوتایی حاصل می شود.به همین ترتیب، این فرایند ادامه پیدا می کند تا بهترین ترکیب حاصله که باعث بیشترین تمایز بین کلاس‌ها می‌شود، به دست آید. ترکیب مذکور شامل چهارده ویژگی است (سه ویژگی زمان،دو ویژگی فرکانسی، شش ویژگی زمان- فرکانس، سه ویژگی غیر خطی).

1-6- کاهش ابعاد ویژگی

در ادامه به منظور بهتر شدن کارایی دسته بندی کننده و کاهش زمان یادگیری از آنالیز مؤلفه اساسی (PCA) استفاده شده است. آنالیز مؤلفه اساسی (PCA) یک روش کاهش بعد با سرپرستی است. هدف از اعمال PCA پیدا کردن یک ماتریس تبدیل است که بردارهای داده، را که متعلق به کلاس‌های مختلف هستند، به یک فضای با بعد کمتر چنان تصویر کند که نسبت پراکندگی بین کلاسی به پراکندگی داخل کلاسی ماکزیمم باشد [36]. با اعمال PCA، بردار ویژگی ترکیبی منتجه از مرحله استخراج ویژگی با چهارده ویژگی به یازده ویژگی ترکیبی کاهش یافته با حفظ 99% انرژی و برای دسته بندی استفاده می‌شوند.

1-7- شبکه عصبی

در این مرحله برای نشان دادن اختلاف بین سیگنال ECG فرد سالم و فردی که دچار مرگ ناگهانی قلبی شده، از طبقه بندی کننده MLP استفاده شده است. به این ترتیب که سیگنالECG مربوط به 1 دقیقه قبل از وقوع حادثه، با یک دقیقه سیگنال ECG فرد سالم بر اساس ویژگی های استخراج شده مقایسه شده است. شبکه MLP استفاده شده دارای دو لایه پیشخور است که با استفاده از الگوریتم پس انتشار خطا با نرخ یادگیری متغییر آموزش داده می شود. تعداد نورون های لایه اول برابر با تعداد ویژگی ها (11نورون) در نظر گرفته شده و لایه خروجی شامل 1 نورون است. که مقادیری بین 0و1 اختیار می کنند. همچنین، تابع فعالیت نورون ها از نوع تانژانت سیگموئید انتخاب و با تغییر تعداد نورون های لایه میانی سعی در بهینه سازی معماری شبکه عصبی می گردد. تعدادی از انواع ساختار های معمول شبکه عصبی مصنوعی بررسی و بهترین شبکه منتجه ، یک شبکه عصبی دو لایه دارای 5 نورون در لایه پنهان با تابع فعالیت سیگموئید استاندارد، است. آموزش شبکه تا زمانی که مربع خطا کمتر از 0.01 شود یا تعداد تکرار های آموزش به 1000 برسد، ادامه می یابد. از آنجا که تعداد مشاهدات به 70 عدد محدود می شود، این شبکه با توجه به روش Leave One Out در هر مرحله یکی از مشاهدات به عنوان دادهtest  و 69 داده دیگر به عنوان دادگان train ودر مرحله بعد داده دوم به عنوان test  و 69 داده دیگر به عنوان دادگان train انتخاب می شوند. به این ترتیب 70 مرتبه شبکه آموزش می بیند[37]. در هر مرحله خطای شبکه استخراج شده و در نهایت میانگین آن محاسبه می شود.یکی از مزایای این روش این است که  همه دادگان به طور غیر همزمان هم در train وهم در test حضور داشته وشبکه تمام قدرت خود را نشان داده است .

دراین مرحله ابتدا ویژگی‌های استخراج شده به منظور ارزیابی میزان تفکیک پذیری، هم به صورت جداگانه(ویژگی های خطی و غیر خطی و زمان-فرکانس) و هم ترکیبی با یکدیگر مقایسه شده‌اند. جدول (1) میزان صحت تفکیک پذیری با توجه به ویژگی های خطی ، زمان-فرکانس ، غیر خطی و همچنین ترکیبی را در 3 دقیقه (180 ثانیه) قبل از وقوع SCD نمایش داده است .

جدول (1) : میزان صحت تفکیک پذیری افراد سالم و ریسک پذیر با توجه به روش‌های تحلیلی مختلف به صورت مجزا و ترکیبی

همان طور که مشاهده می شود، ویژگی های ترکیبی توانمندی بیشتری در تفکیک افراد دارند. از این رو، با توجه به بردار ویژگی مذکور از طریق بررسی صحت تفکیک پذیری در بازه‌های یک دقیقه‌ای(دقیقه اول، دوم ،سوم، چهارم) قبل از وقوع حادثه، سعی در پیش بینی مرگ ناگهانی قلبی شده است. نتایج نشان می دهند در بررسی یک دقیقه قبل از وقوع حادثه ، روش پیشنهادی بردار ترکیبی، توانسته است افراد سالم و افراد در معرض خطر را با صحت 99.37% دسته بندی کند. این نتیجه برای دقیقه دوم برابر 97.86% است که در جدول (2) قابل مشاهده است. این صحت حاکی از بروز علایم شدید مرگ ناگهانی قلبی در افراد ریسک پذیر در 2 دقیقه قبل از وقوع حادثه است.

جدول(2) : تفکیک پذیری دو دقیقه اول قبل از SCD

در بررسی دقیقه سوم ، سیگنال HRV استخراج شده به سه بازه زمانی مساوی تقسیم شده است که هر کدام معرف 20 ثانیه هستند. سپس مطابق قبل از هر پنجره 20 ثانیه ای استخراج ویژگی و انتخاب بهترین ترکیب ویژگی صورت گرفته و از طریق PCA، کاهش بعد یافته اند. در نهایت، از طریق شبکه MLP به دسته بندی بین افراد ریسک پذیر و سالم در هر بازه زمانی 20 ثانیه ای پرداخته شده است. نتایج این دسته بندی در جدول (3) مشاهده می‌شود.

جدول(3) : تفکیک پذیری دقیقه سوم قبل از SCD با بازه‌های 20 ثانیه‌ای.

برای دقیقه چهارم نیز به همین ترتیب  HRV استخراج شده به سه بازه مساوی تقسیم شده که هر بازه معرف 20 ثانیه است. سپس مطابق قبل، پس از انتخاب بهترین ترکیب ویژگی و اعمال PCA به دسته بندی پرداخته شده است . نتایج این دسته بندی در جدول شماره (4) مشاهده می‌شود.

جدول (4) : تفکیک پذیری دقیقه چهارم قبل از SCD با بازه‌های 20 ثانیه ای

همچنانکه مشاهده می شود؛ شخص در معرض مرگ ناگهانی قلبی تقریبا از 4 دقیقه قبل از وقوع حادثه، علایم شدید بروز این اتفاق را در سیگنال HRV خود به همراه دارد. روش پیشنهادی بردار ترکیبی توانایی لازم را برای آشکار کردن این علایم دارد. این در حالی است که سیگنال ECG  فرد ریسک پذیر تفاوت چندانی با فرد سالم نمی کند و پزشک قادر به تشخیص علائم بروز مرگ با توجه به مشاهده سیگنال ECG نیست. از دیگر توانایی‌های روش پیشنهادی ،آشکار کردن روند رو به رشد علایم خطر است، به طوری که هرچه به لحظه SCD نزدیکتر می شود، این علایم افزایش می‌یابد.

2- نتیجه‌گیری

نتایج نشان می دهند درسیگنال HRV افرادی که در آستانه مرگ ناگهانی (SCD) قلبی قرار دارند، ویژگی‌هایی وجود دارد که با افراد سالم تفاوت های بارزی دارد و روش ترکیبی توانایی درخور توجهی در آشکار کردن این اختلاف دارد. نتایج نشان می دهد از 4 دقیقه قبل از وقوع واقعه آشکار سازی احتمال خطر وجود دارد؛ به طوری که هرچه به محل وقوع حادثه نزدیک می‌شویم، احتمال آشکار سازی افزایش می یابد، به این ترتیب که 2 دقیقه سیگنال ECG قبل از SCD  فرد در معرض مرگ قلبی با فرد سالم کاملا متمایز است . در بررسی دقیقه سوم قبل از رخ دادن واقعه، مشاهده می شود که خطر بالای مرگ ناگهانی قلبی  همچنان وجود داشته و قابل آشکار کردن است. از طرفی مشاهده می‌شود که این روند احتمال خطر هرچه به واقعه نزدیکتر می‌شویم، بیشتر می شود. در بررسی دقیقه چهارم، همچنان خطر احساس می شود؛ به طوری که تقریباً بعد از20 ثانیه اول می توان این اختلاف را مشاهده کرد. همچنین، روند روبه رشد احتمال خطر نیز، خود را نمایان کرده است. نکته دیگری که در بررسی‌های انجام شده در این تحقیق قابل ذکر است، این است که فرد سالم و فرد ریسک پذیر به علت وجود تاکی آریتمی های موجود و یا تجربه سکته قلبی، قابل تفکیک هستند، ولی نکته مهم این است که درصد تفکیک پذیری از تقریباً چهار دقیقه قبل به طرز فزآینده ای رشد می کند و روش پیشنهادی توانایی لازم برای آشکار کردن این تغییرات را داراست. پیش بینی مرگ ناگهانی قلبی تا 4 دقیقه زودتر، این امکان را به پزشک و مراکز درمانی می دهد که برای جلوگیری از این واقعه راهکارهایی اتخاذ نمایند که به نجات جان بیمار منجر می‌شود. در جدول (5) میزان تفکیک پذیری را از 4 دقیقه قبل از مرگ قلبی مشاهده می‌کنید.

جدول (5) : میانگین درصد تفکیک پذیری بین افراد سالم و در آستانه مرگ ناگهانی از 4 دقیقه قبل از وقوع حادثه با استفاده از روش پیشنهادی بردار ترکیبی