چکیده

این نرم‌افزار شبیه‌سازی با هدف فراهم‌کردن ابزاری کارآمد برای دانشجویان جهت مطالعه پتانسیل‌های عمل منتشرشونده در رشته‌های عصبی بدون میلین توسعه یافته است. این نرم‌افزار به کاربر اجازه می‌دهد تمامی پارامترهای آکسون، غشاء سلولی، محیط اطراف، الکترودهای تحریک و الکترودهای ثبت را به دلخواه تنظیم کند.

ابزارهای گرافیکی و ترسیمی متنوعی در اختیار کاربر قرار گرفته است تا متغیرهای مختلف غشایی مانند پتانسیل غشاء، مؤلفه‌های جریان غشایی و سایر پارامترهای فیزیولوژیکی را به صورت پویا یا ایستا، در قالب نمودارهای سه‌بعدی و نمودارهای پارامتری دوبعدی مشاهده کند. همچنین نرم‌افزار امکان انتخاب روش‌های عددی مختلف برای حل معادلات هاجکین–هاکسلی را که فرآیند فیزیولوژیکی تولید و انتشار پتانسیل عمل را مدل می‌کنند، فراهم می‌سازد.

نتایج ارزیابی روش‌های عددی نشان می‌دهد که هیچ روش واحدی برای تمامی شرایط بهترین عملکرد را ندارد. کاربر می‌تواند اندازه گام‌های گسسته‌سازی زمانی و مکانی (dt و dz) را تغییر داده و متناسب با معیار مورد نظر خود مانند سرعت اجرا یا میزان مصرف حافظه، روش عددی مناسب را انتخاب کند.

آزمایش‌های انجام‌شده نشان دادند که نرم‌افزار قادر است نتایجی منطبق با مطالعات معتبر پیشین تولید کند. برای نمونه، با استفاده از پارامترهای پیش‌فرض مدل هاجکین–هاکسلی، سرعت هدایت پتانسیل عمل برابر با 18.75 متر بر ثانیه به دست آمد که تنها حدود 12 درصد با نتایج گزارش‌شده توسط Cooley و Dodge در سال 1966 اختلاف دارد.

فهرست مطالب

فصل 1: مقدمه .......................................................... 11
    1-1 پیشینه ........................................................ 11
    1-2 نمای کلی نرم‌افزار ............................................ 12

فصل 2: مبانی نظری ..................................................... 15
    2-1 معادلات هادی مرکزی (Core Conductor) .......................... 15
    2-2 چگالی جریان غشایی ............................................ 17
    2-3 رسانایی‌های غشایی ............................................ 18
    2-4 مدل‌سازی الکترودهای تحریک .................................... 19
    2-5 جریان‌های طولی ............................................... 21
    2-6 اختلاف پتانسیل‌های داخل و خارج سلول .......................... 21
    2-7 مقادیر پیش‌فرض پارامترها ..................................... 22

فصل 3: روش‌های عددی .................................................. 23
    3-1 روش اویلر پیشرو .............................................. 25
    3-2 روش اویلر پسرو ............................................... 27
    3-3 روش کرانک–نیکلسون ............................................ 29
    3-4 روش کرانک–نیکلسون با گام‌بندی پلکانی ........................ 32

فصل 4: نرم‌افزار MATLAB ............................................... 35
    4-1 بخش PAP Control .............................................. 36
    4-2 بخش PAP Workspace ............................................ 37
    4-3 بخش PAP Parameters ........................................... 39
    4-4 بخش PAP Stimulus ............................................. 43
    4-5 بخش PAP Numerics ............................................. 45
    4-6 بخش PAP Voltage Recorder ..................................... 46
    4-7 بخش PAP Variable Summary ..................................... 47
    4-8 بخش PAP Space-Time Evolution ................................. 48
    4-9 بخش PAP 3D Plots ............................................. 49
    4-10 بخش PAP Comparison Plots .................................... 50
    4-11 راهنمای گام‌به‌گام ایجاد شبیه‌سازی سفارشی .................. 56

فصل 5: ارزیابی نرم‌افزار ............................................. 59
    5-1 آزمون‌های نرم‌افزار .......................................... 59
    5-2 ارزیابی روش‌های عددی ......................................... 61
        5-2-1 پایداری ................................................ 62
        5-2-2 کارایی محاسباتی ........................................ 63
        5-2-3 مصرف حافظه ............................................. 66
        5-2-4 دقت .................................................... 69

فصل 6: نتیجه‌گیری .................................................... 73

پیوست A: نمونه اسکریپت ارزیابی ....................................... 75
پیوست B: کد MATLAB روش کرانک–نیکلسون پلکانی .......................... 79
منابع ................................................................ 83

خلاصه پایان‌نامه

این پایان‌نامه با هدف توسعه یک محیط آموزشی و پژوهشی برای مطالعه بیوفیزیک سلولی و به‌طور خاص بررسی پدیده انتشار پتانسیل عمل در فیبرهای عصبی بدون میلین انجام شده است. یکی از چالش‌های مهم در آموزش نوروفیزیولوژی و بیوفیزیک سلولی، درک رفتار دینامیکی غشاء سلول‌های عصبی و نحوه انتشار سیگنال‌های الکتریکی در طول آکسون است. اگرچه مدل هاجکین–هاکسلی از دهه 1950 به عنوان یکی از بنیادی‌ترین مدل‌های ریاضی در توصیف فعالیت الکتریکی نورون‌ها شناخته می‌شود، اما پیچیدگی معادلات آن باعث می‌شود دانشجویان و پژوهشگران درک عمیقی از تأثیر پارامترهای مختلف بر رفتار سیستم نداشته باشند.

در این تحقیق تلاش شده است با طراحی و پیاده‌سازی یک نرم‌افزار آموزشی در محیط MATLAB، امکان مشاهده و تحلیل مستقیم رفتار مدل هاجکین–هاکسلی فراهم شود. این نرم‌افزار به کاربران اجازه می‌دهد ساختار آکسون، ویژگی‌های غشاء، شرایط محیطی، الکترودهای تحریک و الکترودهای ثبت را به صورت کامل تنظیم کرده و اثر هر تغییر را بر انتشار پتانسیل عمل مشاهده کنند.

بخش نظری پایان‌نامه به بررسی مبانی فیزیولوژیکی و ریاضی انتشار پتانسیل عمل اختصاص دارد. ابتدا معادلات Core Conductor معرفی می‌شوند که جریان‌های الکتریکی داخل و خارج سلول را توصیف می‌کنند. سپس نحوه مدل‌سازی جریان‌های یونی سدیم، پتاسیم و جریان نشتی بر اساس مدل هاجکین–هاکسلی تشریح می‌شود. همچنین وابستگی رسانایی‌های یونی به متغیرهای فعال‌سازی و غیرفعال‌سازی کانال‌های یونی به تفصیل مورد بررسی قرار گرفته است.

یکی از مهم‌ترین اهداف این پایان‌نامه مقایسه روش‌های مختلف حل عددی معادلات دیفرانسیل جزئی حاکم بر انتشار پتانسیل عمل است. به همین دلیل چهار روش عددی شامل اویلر پیشرو، اویلر پسرو، کرانک–نیکلسون و کرانک–نیکلسون پلکانی در نرم‌افزار MATLAB پیاده‌سازی شده‌اند. تمامی شبیه‌سازی‌ها، محاسبات و تحلیل‌های ارائه‌شده در این پژوهش با استفاده از نرم‌افزار MATLAB انجام شده است و این محیط نقش اصلی را در توسعه مدل و اجرای آزمایش‌ها ایفا می‌کند.

نویسنده نشان می‌دهد که هر روش عددی دارای نقاط قوت و ضعف خاص خود است. روش اویلر پیشرو از نظر سادگی و سرعت پیاده‌سازی مناسب است اما محدودیت‌های پایداری شدیدی دارد. در مقابل، روش اویلر پسرو پایداری بالاتری ارائه می‌دهد اما نیازمند محاسبات بیشتری است. روش کرانک–نیکلسون تعادل مناسبی بین دقت و پایداری ایجاد می‌کند و روش کرانک–نیکلسون پلکانی نیز با هدف افزایش بهره‌وری محاسباتی توسعه یافته است.

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

تمامی این قابلیت‌ها بر پایه شبیه‌سازی‌های انجام‌شده در MATLAB طراحی شده‌اند. در واقع MATLAB نه تنها برای حل معادلات هاجکین–هاکسلی مورد استفاده قرار گرفته، بلکه تمامی ابزارهای نمایش داده‌ها، رسم نمودارها، تحلیل نتایج و ارزیابی روش‌های عددی نیز در همین محیط توسعه یافته‌اند.

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

در نهایت این پایان‌نامه به این نتیجه می‌رسد که ترکیب مدل هاجکین–هاکسلی با یک محیط تعاملی مبتنی بر MATLAB می‌تواند ابزار قدرتمندی برای آموزش و پژوهش در حوزه بیوفیزیک سلولی باشد. نرم‌افزار توسعه‌یافته نه تنها امکان مطالعه عمیق‌تر رفتار نورون‌ها را فراهم می‌کند، بلکه بستری مناسب برای آزمایش روش‌های عددی مختلف، بررسی اثر پارامترهای فیزیولوژیکی و درک بهتر فرآیند انتشار پتانسیل عمل در سیستم عصبی ارائه می‌دهد. بنابراین مسئله اصلی این پژوهش، یعنی دشواری آموزش و تحلیل رفتار دینامیکی نورون‌ها، از طریق توسعه یک سامانه شبیه‌سازی جامع در محیط MATLAB تا حد زیادی برطرف شده است.

“`