椭球模型下EEG和MEG正问题数值计算

来源 :上海师范大学 | 被引量 : 0次 | 上传用户:xigongdaxigongda
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
EEG(脑电图)/MEG(脑磁图)是脑功能研究及脑疾病临床诊断的两个重要影像技术。与其它影像技术相比,EEG/MEG具有时间分辨率高及无损伤的优势,能够以毫秒级的速度跟踪神经元的生理变化,适于研究人脑功能的动态变化。EEG/MEG是集生物医学工程、超导技术、数学、物理、计算机等领域为一体的交叉学科问题,分为理论、仪器设备及应用等方面。为定量地分析解释EEG/MEG测量数据,理论上我们需要研究生物电磁反问题,即由外部电位势和磁场的测量来确定脑内部神经电流元。因而数学模型的建立,数学方法及数值模拟是必不可少的研究手段。在EEG/MEG问题的理论研究中,电偶极子常用作局部神经电流源的近似源模型,电偶极子定位研究可应用于诊断脑病灶位置和功能区域的划分,已得到广泛关注。本文针对电偶极子源模型,在Maxwell方程的拟稳态的假设下,研究了非球脑模型下EEG/MEG正问题的数值计算,它是研究偶极子定位反问题的一个重要组成部分。EEG和MEG正问题即是在假定的电流源分布条件下分别计算脑皮层上的电位势和脑外部磁场。偶极子的精确定位取决于脑模型及所采用的数学方法。以往的脑模型大多采用的是球模型,近年来研究结果表明球模型用于近似人脑有一些不足之处,医学生物界提出了研究非球模型的需求。目前基于MRI或CT影像所建立真实脑形的数值脑模型开始用于EEG/MEG问题的研究中,但过大的计算量难于适于反问题的研究,而且不利于作理论分析。我们则采用解析的椭球脑模型建立了基于边界元的EEG/MEG正问题的数值计算方法。针对偶极子在三个坐标轴上的不同位置和三个不同的方向,对椭球脑模型及球模型实施了数值模拟,分析比较了两种模型下的EEG/MEG计算结果以及EEG/MEG正问题的解相对偶极子参数变化情况。全文分为四个部分:第一章前言第二章基本理论第三章EEG/MEG正问题的计算第四章数值模拟第一章主要介绍了EEG/MEG影像技术的研究发展历史、测量仪器及EEG/MEG生理基础。第二章给出了在拟稳态Maxwell方程假设下,EEG/MEG各自的场方程。在第三章中,我们利用边界元法分别建立了EEG/MEG正问题计算公式。第四章,分别针对椭球模型和球模型,实施了数值模拟,检验了我们的方法的可靠性和误差,针对偶极子位于三个坐标轴上的不同位置和三个不同方向,比较分析了两种模型下,电位势和磁场与偶极子参数之间的变化情况。
其他文献
曲线的光顺问题是NURBS研究和应用中经常遇到的问题,虽然已经提出了很多方法对NURBS曲线进行处理,但是高度自动化很难实现,且没有统一的光顺准则,计算量大,本文主要完成了以下几个
供应链管理作为全球化市场竞争的一种有效管理模式,强调企业间的合作与协调。然而,当各成员都试图最大化自己的利润时,整个供应链的利润将不可避免地受到损害,这就会出现“双重边际化”。供应链契约是为了尽量减少这种损害而提出的一种解决办法,它通过提供一些激励的措施来协调供应链各成员的关系,使分散决策得到的各方利润之和与一个集中系统决策得到的利润尽量相等。数量折扣就是能够实现这一目标的一种契约。数量折扣是指对
由于广义Pareto分布在金融和保险等领域的广泛应用,对于该分布的统计推断也成为许多学者的研究热点.目前贝叶斯统计推断理论几乎可以作为每一个学科的研究工具之一.为此本文
(λ)Kn表示一个(λ)重的完全图,(λ)重的完全图是指图中任意两点之间恰好有入条边相连.图的圈分解,就是把一个图的边集分解成一些圈的集合,使得这些圈中的边恰好划分图的边集
20世纪90年代中期开始基于偏微分方程(Partial Differential Equations,简写为,PDE)的图像降噪方法在图像处理领域中成为研究热点。显著的去噪效果和广泛的使用范围引起了国
学位
众所周知,对耗散系统而言,吸引子是一个描述系统极限行为的有利工具,因此研究耗散系统的解生成的动力系统的吸引子的存在性是很有意义的.对随机系统而言,亦是如此.但此时,吸引子就
相对同调代数是S.Eilenberg和J.C.Moore于1965年引进的。关于这门学科的理论研究,极大的丰富和发展了同调代数的经典结果,而环与模的相对同调维数理论是相对同调代数这一门学科