《多参数电刺激下神经元放电节律的研究.doc》由会员分享，可在线阅读，更多相关《多参数电刺激下神经元放电节律的研究.doc（41页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 多参数电刺激下神经元放电节律的研究多参数电刺激下神经元放电节律的研究摘 要非线性动力学系统有复杂和简单之分，其中复杂的非线性动力学系统之一就是神经元，神经元在不同参数的电刺激的时候，神经元就会产生比较复杂的动作电位的不同序列模式，有周期峰放电、混沌放电、交替放电、周期簇放电等多种丰富的放电节律，不同的刺激信号被神经元不同的放电节律所代表。在分析和讨论外电场下细胞膜电生理特性的通常基础上，本文建立了一个有用模型， 该模型是基于神经元建立的，各种振幅和不一样的频率都会对神经元产生明显的影响。多方研究之后可以发现神经元被多种多样的不同交流外电场刺激时神经元的动态行为是不一样的：放电规律会随着频率的
2、变化而变化，频率不一样（比如增大），放电节律也会不一样（簇放电到峰放电）；峰与峰的间期序列根据两种周期放电模式的变换呈现的特征是整数倍关系。我们对其放电节律的深入研究和探讨都对于神经信息进行不同的编码有着十分重要的意义，我们对神经元的研究结果，进一步揭示了交流外电场对神经元的不同放电节律的作用的不同影响，该题目的研究在外电场和神经系统兴奋程度的联系和开发相关疾病的治疗等方面意义重大。关键词：神经元；放电节律；交流外电场；HH模型AbstractNeurons are complex nonlinear dynamic systems ,A complex action potential se
3、quence pattern can be generated under different parameters of electrical sti-mulation ,(Alternating discharge、chaotic discharge、periodic peak discharge、period-ic cluster discharge, etc.)Different firing rhythms of neurons represent differe-nt stimuli.On the basis of analyzing the electrophysiologica
4、l characteristics ofcell membrane under external electric field, the neure model under external electric field was established .The dynamic behavior of HHneurons underthe action of different amplitude and frequency external electric field was studied.With the increase of frequency,the reverse additi
5、on period bifurcation sequence from cluster discharge to peak discharge and chaotic discharge is shown in the discharge rhythm.There is a new alternating rhythm in the process of periodic rhythm transition. The discharge sequence is the alternation of two periodic discharge modes, and the interpeak
6、sequence has the characteristic ofinteger multiple.The study of its discharge rhythmis of great significance forthe study of neural information coding, The study of this topic is of great significance for the excitement of external electric fields and nervous systems and the development of related d
7、isease treatments.Keywords: neure ; Discharge rhythm ; Ac external electric field ;HH model2目录摘 要1Abstract2第1章 绪 论31.1 课题背景31.1.1国内发展现状31.1.2国外发展现状41.2研究的目的51.3研究的意义5第2章 放电节律研究的理论基础72.1 生物神经元72.2动作电位的基本概念82.2.1 动作电位的产生82.2.2 动作电位的特点102.3神经元的静息电位112.4 神经元模型建立112.4.1 HH经典模型112.4.2外电场作用下的细胞模型142.4.3外电场
8、下HH模型152.5神经元的放电节律162.6本章总结17第三章 MATLAB基本知识183.1 MATLAB组成及功能特点183.1.1 MATLAB平台的组成183.1.2 MATLAB功能193.1.3 MATLAB特点203.2 Simulink优点及建立223.2.1 Simulink优点223.2.2 Simulink建立233.2.3 Simulink建立的基本步骤253.3 本章小结25第四章 多参数电刺激下神经元放电节律的研究264.1刺激频率对神经元放电节律的影响264.2刺激振幅对神经元放电节律的影响28设计总结与存在的问题33致谢34心得体会35参考文献35第1章 绪
9、论1.1 课题背景目前，社会竞争压力越来越大，人口和家庭结构变化越来越明显、抑郁和痴呆来源于老年人的比例变多、非法药品的出售和滥用、多数人口的自杀现象增多以及重多受到灾害的人群的心理危机感凸显等各个方面的问题日益增多和突出，有大部分大学生的各种心理问题也会不断发生，且趋势也在不断上升。因此，我们必须通过有效的干预以及合理的进行控制这种神经精神类疾病的发生，这种情况已经刻不容缓，它不仅可以改善我国人民的健康状况和生活质量，而且可以降低我国的医疗支出的费用。虽然利用药物来预防和控制神经精神疾病是一种比较传统的方法，但是实践证明：药物治疗能给神经精神疾病带来一定程度上的疗效，但也存在着很多的不足之处
10、。例如，对大脑用药的时候其余的正常脑区域也会受到药物的负面影响，这些药物会干扰到部分正常的脑区域进行的工作；然而药物也会对于多数具有抗药性的有关各种具有运动障碍的疾病，例如原发性震颤、帕金森症以及肌张力的异常等药物的治疗效果都表现为不佳的状态。研究学者们想要通过克服药物的治疗过程中所出现的各种各样的问题，人们则需要通过研究使用非药物性的神经刺激的不同方法来对具有神经功能性的疾病来进行相关的控制和具有针对性的治疗。1.1.1国内发展现状这些年来技术不断进步，我们国家的很多研究人员在建立关于神经细胞的数学模型方面颇有建树，不仅在国内，在国际上也得到了很多同行的认可。我国学者惠萍在20年前首次发现了
11、具有非线性特征的HH模型，科学家们利用HH模型具有很多种状态的特性解释清楚了为什么生物的细胞膜会表现出激发现象。2004年，学者刘青云等人针对白噪声设计实施了很多实验，最终表明如果噪音的分贝数达到一定的阀值，自共振现象就会相应而生。学者韩晓鹏等人在这一年基于HH模型也做了很多研究，他们得出的结论是单个的神经细胞到达临界值后会产生随机共振现象，但这只是单个神经细胞自己的作用，其对神经元网络（双层结构）判别临界值的能力影响甚微。八年后，研究人员刘少宝研究了噪声对K离子和Na离子的影响，结果表明：当K离子和Na离子的通道收到噪声影响后会对HH模型产生影响，不仅影响其空间模式，还会影响模型的时间模式。
12、学者邓斌在2014年仿真模拟了主神经元和丛神经元的协同控制问题并最终得到了多种图像。主神经元在适当的算法的控制下可以表现出和丛神经元一样的放电情况。HH模型中参数众多，细胞膜的特性可以从模型中找到相应的参数和它对应，默认状态下这些参数对应的是能够使细胞膜活跃的动态性质。1.1.2国外发展现状过去的一段时间，国际上在研究神经组织与神经细胞的过程中，取得了越来越多的研究成果。最早在1791年，意大利科学家加尔瓦尼首次发现生物体中的带电现象。60年后，科学家们发现不仅生物体带电，神经组织中的神经元也会出现电现象。之后科学家对神经组织对研究在上个世纪中期有来长足对进步。英国生物学家赫胥黎和霍奇金首次对
13、生物放电现象进行了详细的解释，并提出来一个能准确描述其规律对数学模型，即HH模型。他们以大西洋枪乌贼为研究对象，重点研究其神经元的放电状况，然后对他们获得的数据进行拟合，并将这一结果发在了国际期刊上。他们的发现极大的鼓舞了后来的生物学家，这也成为生物学发展的一个里程碑，从此以后，很多生物学家进行研究的同时，不仅要做出相应的实验，还要找到可以反映实验情况的数学模型，这极大地改变了科学研究的方式。到了20世纪末，科学家们基于HH模型做的大量的研究，他们在研究中发现，只要有轻微的耦合作，神经细胞中的神经元就会产生放电现象，而且这种放电现象呈现一定的规律。而且研究人员发现，如果这种耦合是兴奋性的，那么
14、在这种情况下神经元的放电现象就会减弱。科学家在1998年又发现基于HH模型的神经元在噪声刺激下会产生相干共振，这又是一个新发现。转眼来到21世纪，外国科学家在2005年发现当外界有电信号刺激到神经细胞以后，神经组织中的信号传输速率就会受到一定程度上的影响。2013年又有两个科学家发现，。当给定一个刺激信号，如果这个刺激信号没有达到它的阀值，那我在这种情况下HH模型的性质是没有任何变化的。因此科学界得到如下结论：如果外界的刺激频率没有达到神经网络中的阀值，那么神经网络中就会产生一个和这个刺激频率以相匹配的一个信号，二者之间有明显的相关性。到了第二年，也就是2014年，科学家们研究了鲨鱼的神经系统
15、，科学家发现鲨鱼能够对外界的噪音造成产生反应，即使这个刺激很微弱。这也说明鲨鱼的神经系统相当发达，灵敏。这也帮助鲨鱼驰骋大海，轻易的捕获猎物。然后科学家们又根据HH模型研究分析来鲨鱼的神经系统，对鲨鱼大脑的工作机制有了深入的了解。这种研究方法对以后人类的研究也大有裨益。1.2研究的目的神经细胞和神经胶质共同构成了人体的神经系统，其中，在人体中主要起主导作用的是神经系统。神经系统在人体内发挥着巨大的作用，对人体各种信号之间传递有着不可或缺的作用。神经系统的最基本单位是神经元，神经元是神经系统中信号传递的最基本单位，当神经元传递信号时神经元自身会产生放电行为，传递的信号不同，放电行为也会有所差异。
16、因此就可以通过实验分析神经元的放电规律进而解析出神经元传递的信号。这对于在不同刺激状态下神经信息在神经系统中的编码规律的研究意义重大。而且最重要的是，这项研究不仅可以应用于生物学，还可以用于临床医学，对现代神经学疾病的研究也有很大的作用。现如今，人们越来越多的疾病和疼痛来自于电场，神经元的生理变化和外电厂之间的联系如果能研究明白，那对于这类疾病的治愈将很有帮助。大量的实验已经证明，英国科学家赫胥黎提出的HH模型可以将神经元的膜电位变化形象的的描述出来，HH模型给神经系统的研究指出来一条可行之路。因此本论文基于HH模型分析探讨了细胞在外电场作用下的工作原理和神经元的电位规律，本文使用的外电场为大
17、脑的电刺激。这种治疗方式会对病人产生一定的损伤，因为在治疗的时候会在病人的脑部植入一个电极，当然为了安全起见这个电极很微小。这个电极会产生不同频率和不同幅度的电信号来刺激细胞。虽然现在已经有方案（DBS）用于治疗神经疾病，但这种治疗方法的具体工作原理等尚未明朗，还需要进一步探索。利用外电场的刺激之所以能够治疗神经类疾病是因为外电场的控制器能够发出和脑刺激一样的电信号来代替脑刺激，对大脑深部的神经元产生影响，这样就能够控制和治愈疾病。所以这就为抑郁症等精神类疾病的治疗提供类解决方法。1.3研究的意义通过技术的进步和时代的不断更新，世界上当前已经发展了很多种包含有纯数学的或者具有不同生理意义的相关
18、的神经元的具体模型，尽管HH模型具有较高的维数和非线性，以至于HH模型不能够适用于大规模的神经网络的模拟，但是，我们要详细研究单个的神经元具有的不同的特性，因为，HH模型的参数规定的生理意义比较具体和清晰，具有十分典型的代表性。所以我将以HH模型为我的研究基础，于是便建立了在外电场的作用下通过改进后的神经元的具体模型，通过修改不同参数的电刺激来研究神经元的放电节律，并给出了仿真结果和理论解释。如果将精神疾病的病因放到了细胞水平来看的话，那就是因为神经细胞的兴奋阀值变化不规律导致的。神经细胞的兴奋阀值，不论是减小还是变大都会引起疾病，因此治病要治本，要从根本上抑制这些精神类的疾病，必须要改变神经
19、细胞兴奋阀值的兴奋程度。世界卫生组织对中国精神疾病的负担做了一下预算，根据他的推测，中国精神疾病的负担会逐年上升，到2020年会增长到疾病总负担的四分之一，这也给中国医疗卫生行业敲响了警钟，我们必须要对精神类疾病加以控制和治疗，此任务迫在眉睫，刻不容缓。精神类的疾病如果能治好，不仅能为国家节省大笔的财政，还可以极大的提升人们的生活质量和幸福指数，让人们活得越好。市面上已经有了很多精神类疾病的药物，但是尽管这样的药物给疾病带来一定程度上的缓解，却无法根治，而且使用不当还会有很大的副作用，因为大脑毕竟是十分敏感的区域，使用不当会影响大脑的正常功能，因此，药物治疗不是长久之计。科学家一直以来也在想办
20、法克服这个问题，他们进行了很多研究和实验，最后也开发出了一些可以控制和治疗精神疾病的神经刺激药物。神经元网络的研究也被列为21世纪以脑科学主要研究目标之一。因为在神经系统中神经元是神经系统最基本的单位，想要研究神经系统，必须先从神经元入手，通过对神经元的研究，我们最终也会搞明白信息是如何在神经系统中进行传递的，神经元在此过程中具有不可或缺的地位和作用。38第2章 放电节律研究的理论基础2.1 生物神经元生物神经元作为最基本的人类大脑的神经系统。每个神经元与一万多个其他神经元相连接，而一百多亿个神经元构成了生物神经元。人类大脑的神经元细胞可以产生、传递和处理相关的电信号，这就是人类大脑神经元与身
21、体神经元所不同之处。细胞核、细胞膜和细胞质这三个重要的部分共同构成了神经元。神经元轴突投影长度长短不一，从最短几微米到最长1米。轴突的表面十分光滑，厚度也很均匀。神经元的信息通过胞体向四周扩散，通过树突传输到其他神经元。神经元之间的连接方式是通过突触实现的，有轴突和细胞体，轴突和轴突，轴突和树突，树突和树突还有许多其他的连接形式。同样的脉冲串的宽度和振幅是神经元单元，输出脉冲激发它的轴突是神经元的高频率；另一方面，突触的输出脉冲减少了神经元的频率，即没有脉冲。科学家将突触分为两类，活跃性突触和克制性突触，这是突触对下一个神经元将要进行活动的不同影响。活跃性突触能够让下一个神经元变得很活跃，相应
22、的，克制性突触就会抑制神经元的活动。一般情况下在以下两种情况发生时神经元的静息电位就会被触发：第一种情况是由于细胞膜具有特殊的成分和结构，导致它有特殊的功能，比如细胞膜对不同对离子对通过能力是不同的，有的离子通过的多，有的离子就不能通过，这样也起到一定的保护作用。细胞膜所具有的特性是对离子的选择通透性；二是细胞膜两侧的离子的浓度有所不同。在未受到外界的刺激时，神经元的细胞膜上会大批量的掀开钾离子通过所使用的相应通道，神经元细胞内部的钾离子大量流出，那钠离子相应的通道就只有很少被开启，细胞内部很多带有负电的大分子因为不能像钾离子一样经过细胞膜，就留在来细胞内部。这样细胞膜两面的带电性就会不同，细
23、胞膜外面带正电，细胞膜内带负点。当内外达到平衡的时候就是一种稳定状态，也叫静息态。通过以上分析可以知道，影响静息电位的最重要的因素是钾离子的浓度，如果细胞膜两侧的浓度突然发生变化，稳定态，也就是静息态就会发生改变，科学家经过科学的实验，也证明了这一点。如果系统处在一种平衡的状态，也就是静息态，因为细胞膜的选择通过性这一特性就会导致细胞膜两侧产生浓度差，也就会产生外正内负的电位差，也就是我们经常说的静息电位。人类最早发现静息电位是在上个世纪30年代末期，英国科学家赫胥黎等人以枪乌贼为研究对象，对其神经轴进行研究的时候偶然发现了一个电位差。于是他们进行了更进一步的研究，如果把检测单位的装置都放在细
24、胞膜外，装置就检测不到电位，如果把其中的一个缓慢的穿过细胞膜就会出现电位差，这就是我们后来所说的外正内负的静息电位，它的大小介于-10mV到100mV之间。细胞膜把神经元的胞体分成两个部分，即外部和内部，由于外正内负，外部电位高于内部电路，这种情况下就会形成电位差，最终产生膜电位。2.2动作电位的基本概念当没有外部刺激时，神经元的细胞膜的状态是静息电位，当外部信号刺激到神经元时，稳定状态不再持续，静息电位的状态就会发生改变，产生动作电位。形成动作电位主要有以下这些基本的条件：（1) 细胞膜对带电离子的的选择通过性（2)外部刺激。（3)细胞膜两侧存在浓度差。2.2.1 动作电位的产生神经元细胞的
25、活动规律的表现形式为放电行为，也就是说神经元相互之间的放电行为，这种放电行为也是神经元之间相互传递信息的一种方式，大脑活动和神经网络之间的联系也靠放电行为来完成。当没有外界环境的刺激时，神经元处在稳定状态，该电位称为静息电位。如果有外界环境刺激到神经元的细胞，神经元的稳定状态就会被破坏，静息电位发生改变，会产生一个脉冲，这就是动作电位。神经元细胞膜两侧电位不一样，表现为外正内负。如果这种电位差情况加剧，就会成为超极化。当有外界环境的刺激时，随之产生的动作电位会削弱这种极化，也叫去极化。神经元的电位经常在这两种模式之间切换。后来霍奇金等科学家仔细研究了动作电位的产生机制，并结合神经元的离子模型阐
26、述了这种机制：(l) 因为细胞膜对离子具有选择透过性，这就会导致细胞膜内外出现离子浓度差，这样就会形成静息电位。(2) 如果有神经元细胞受到外界的信号刺激，去极化就会产生静息电位转变成动作电位，这时候细胞膜的离子通透性也会发生改变。首先会打通钠离子的通道，这样钠离子就可以大量通过，从细胞膜外流入细胞内，随着钠离子的大量进入，细胞膜内的电位也会发生变化，细胞膜的原来是负电位现在变成了正电位，去极化也会加速进行，如此循环，细胞膜内的电位越来越高，而且会因为此时会产生相当大的内向电流，最终导致出现超射现象，动作电位进一步上升。(3) 细胞膜内部的电位不会盲目的增长，当电位增加到一定程度就会趋于饱和，
27、细胞膜就会关闭一部分钠离子通道，相应的内向电流也会逐渐降低。(4) 细胞膜上的钠离子通道逐渐关闭以后，细胞膜内部的动作电位达到最高，此时钾离子开始向外流动形成外向电流，细胞膜内的动作电位开始降低。(5) 钾离子外流后细胞膜内的离子浓度重新趋于平衡，重新形成新一轮的静息电位。神经元膜上的钠泵能吸收钾离子并排泄钠离子，从而维持神经元膜内外的离子浓度差。通过前面的描述我们已经知道，静息电位和动作电位的产生原理是不一样的，但是在研究二者的时候可以采用类比的方法去研究。正常情况下因为细胞膜对离子的选择通过性会导致细胞膜内的哪里的浓度很低，细胞膜外的钠离子浓度比较高，但即使存在这样的浓度差，钠离子也不能随
28、意进入到细胞膜内部，因为细胞膜特殊的结构，正常状态下细胞膜上的钠离子通道是不会开启的，这样钠离子就进不来。如果此时有外界环境的刺激信号，细胞膜上的钠离子通道就会打开一部分，少数钠离子就可以通过细胞膜上的钠离子通道进入细胞内部，这时候细胞膜内外的离子浓度差就会发生变化，有离子浓度差导致的电位差在这种情况下也会变小。如果细胞膜两侧的电位差继续减小的话，这样会导致细胞膜上钠离子通道更多的被打开，钠离子从细胞外进入细胞内的速度加快，细胞内的钠离子浓度迅速增加，细胞内的正电荷因为钠离子的涌入越来越多，细胞内的电位也越来越高，当电位高到一定程度之后就会趋于饱和，阻碍钠离子的流动，钠离子的数量也会趋于饱和，
29、这个时候进入到钠离子的平衡电位，细胞膜上的钠离子通道这个时候就会关闭，钠离子不再从细胞外流向细胞内部。这时候虽然钠离子通道已经关闭，但是细胞内部过高的钠离子会刺激到钾离子的通道开启，由于细胞内部的钾离子本身就比外部的钾离子浓度高，因此在这个时候钾离子就会通过细胞膜上的钾离子通道从细胞内流向细胞外，如果钾离子外流的速度和钠离子内流的速度达到平衡的时候就会出现峰值电位，随着钾离子继续外流，这个峰值电位的状态就会发生改变，随着细胞内部大量的钾离子流到细胞外部，细胞内部的正电荷迅速减少，细胞内部的膜电位也迅速降低，即动作电位降低。动作电位降低以后，就会接近于静息电位，但这并不是真实的静息电位，因为本来
30、细胞膜内部应该有大量的钾离子，而这时候由于钾离子大量流出和钠离子的大量涌入导致细胞内部主要是钠离子。这时候就需要细胞膜上的ATP酶发挥作用了。细胞膜上有两种特殊的ATP酶，一种是钠ATP酶，一种是钾ATP酶。钠ATP酶会把细胞内部的钠离子输送到细胞外，钾ATP酶则把钾离子从细胞外重新带回到细胞内部，这就会回到钾离子和钠离子在细胞膜两侧具有浓度差的初始状态。这种状态也为以后的兴奋状态提供了有利的条件。下面我们来做一下总结，细胞膜上的钠离子通道打开，钠离子从细胞的外部迅速流向细胞内部，这就会导致动作电位上升。然后钾离子通过细胞膜上的钾离子通道流向细胞外部，导致动作电位降低。动作电位的变化主要是因为
31、钠离子在细胞外和细胞内的浓度差所导致的，如果钠离子向细胞内流动，细胞内部的钠离子浓度就会升高，相应的，动作电位也会在一定程度上升高。因此，如果想要限制动作电位的变化，就需要改变细胞膜对钠离子的通透性，具体方法可以向细胞膜上注射药物，达到阻断钠离子开启的效果。2.2.2 动作电位的特点（1） 动作电位的第一大特性就是在传递的时候不会衰减，这就是说，从某一个神经元产生的动作电位在不同的神经元之间进行传递的时候，动作电位不会发生衰减，也就是说，这个动作电位的幅值和它的形态都不会发生任何变化。（2） 动作电位的第二大特性就是有或者无：神经元的工作原理并不是我们通常想象的那样，一有刺激就会兴奋，只有外界
32、的刺激，超过神经元自身的阀值时，神经元才能够达到兴奋状态，然后动作电位才会出现，但是这个动作电位变化与否，这个动作电位的状态性质，与外界刺激的类型没有关系，也就是说无论外界的刺激类型如何发生改变，动作电位的特性也不会变化，比如持续时间或者是幅值。（3） 动作电位的第三个特点是不能够叠加，因为动作电位不会因为外界刺激的改变而发生变化，也就是说，只要动作电位产生了，它就不会再发生改变，而且不能叠加，也就是说增加刺激也不会增加动作电位，在任何维度都不可以。2.3神经元的静息电位神经元能够产生静息电位主要有两个原因，首先一个就是细胞膜自身的原因，细胞膜对不同离子的离子通透性是不一样的，也就是说，它对离
33、子具有选择透过性。第二个就是由于细胞模对离子具有选择透过性，会导致细胞膜两侧的离子浓度有差异。正常情况下，细胞膜上的钾离子通道开启，大量的钾离子从细胞内部经过细胞膜上的钾离子通道流到细胞外部，但细胞膜上的钠离子通道开启的很少，钠离子就进不去，这就会导致细胞膜内外出现浓度差。除此之外，由于细胞膜上没有可供大分子通过的通道，细胞膜内部的带负电的大分子只能留在细胞内部，而由于带正电的钾离子已经流出，这样就会导致细胞膜两侧的电位出现外正内负的情况，因为细胞膜外呈现出正电状态，正电荷趋于饱和，钾离子外流的情况就会受到抑制，直到最终钾离子不再向细胞膜外流动，这个时候就会达到一种平衡稳定的状态，这种状态就叫
34、做静息态，这种电位叫做静息电位。因此我们可以得到这样的结论，钾离子是影响静息电位的最重要的原因，当钾离子外流到细胞膜外以后，细胞膜内钾离子浓度就会发生变化，这样就会改变静息电位，科学家柯蒂斯等人也通过实验研究证明来这一点。由于神经元的细胞膜对不同的离子具有选择透过性，所以如果没有外部刺激的话神经系统细胞膜的两侧就会存在浓度差，钾离子在细胞膜内部含量低而在细胞膜外部含量高，由于钾离子带正电，此时就会产生电位差，呈现一种外正内负的状态，这种状态也叫做静息电位，上个世纪30年代末期，英国生物学家赫胥黎和霍奇金首次对静息电位进行了详细的解释。他们以枪乌贼为研究对象，对其神经轴进行研究的时候偶然发现了一
35、个电位差。于是他们进行了更进一步的研究，如果把检测单位的装置都放在细胞膜外，装置就检测不到电位，如果把其中的一个缓慢的穿过细胞膜就会出现电位差，这就是我们后来所说的外正内负的静息电位。2.4 神经元模型建立2.4.1 HH经典模型在HH神经元模型建模过程中，用等效电路来模拟神经元具有的电学特性。取一小段乌贼的巨大轴突，其等效的电路如图一所示。电容则表示细胞膜内存储电荷的能力大小，存在差异的离子通道由电阻进行表示，细胞膜的膜内和膜外因为离子的浓度差异产生的电位差由电动势进行表示。图一 HH模型的神经元等效电路图在图一中跨膜电流主要包括两个部分：（1） 通过细胞膜的膜电容电流。（2） 通过各种不同
36、的离子通道的电流（钠电流INa、钾电流IK和漏电流IL）。根据基尔霍夫电流定律表示等效电路的微分方程为： (2-1)式中，Cm表示膜电容；V表示细胞膜内外的电位差；Iion = INa+IK+IL表示Na离子和钾离子穿过细胞时产生的电流，其中Iext表示外界环境的刺激信号的电流大小。以下微分方程共同组成了经典的HH模型： (2-2)式中，Cm表示膜电容；Iext为外部施加的刺激电流；、分别为钠离子通道、钾离子通道和漏离子通道的最大电导；m，h和n分是跟离子有关的变量，n是跟钾离子有关的激活变量，m和h则和钠离子有关，一个控制激活，一个控制抑制；VNa、VK和VL是一种平衡电位，对着着各自离子的
37、电流；=3(T-6.3)/10是一个常数，这个常数和电流电压无关，只跟温度有关，在本文中这个温度T设定为6.3。和分别表示离子通道开通和关闭的转换率。跨膜电压V的非线性函数： (2-3)HH神经元的模型中各参数值如下表一所示，为了便于研究和观察，我们仿真时将真实的膜电位数值向上平移了约65 mV，使该模型的静息电位变为零。设定V、m、h、n的初始数值分别为0、0.053、0.596、0.317。表一 HH神经元模型的典型参数HH模型参数参数数值单位符号Cm1.0F/cm2120mS/cm236mS/cm2gL0.3mS/cm2VNa115mVVK-12mVVL10.59mVT6.3将（22）改
38、写成标准形式： 其中 （2-4）2.4.2外电场作用下的细胞模型 图二 外电场作用下的细胞模型图图二中E(t)为外加电场场强,为外电场与其在细胞上某点处感应产生的V的夹角，如下表所示为各参数的典型值：表二细胞模型中各参数的典型值参数表示符号典型值细胞半径R10m膜厚度d0.5nm细胞质电导率i0.3S/m细胞膜电导率m3107S/m膜外介质电导率o1.2S/m细胞介质电常数i6.410-10As/Vm细胞膜介电常数m4.410-11As/Vm膜外介质介电常数o6.410-10As/Vm细胞质的电容因为外加来低频电场它的容抗值不小因此它的电容效应a影响不大，可以被忽略掉，这是允许的。也就是说我们
39、不再考虑这两个介电常数 i 和 o ;细胞质的电容是可以忽略掉的，因为它影响不大，但是细胞膜电容万万不可被忽略，因为细胞膜不同于细胞质，它的电阻远远超过来细胞质。前人针对VE的计算已经做了大量的研究，我们选择其中一个研究结果计算VE： （2-5）其中，和K都是常数，前者和时间有关，后者与细胞特性有关,根据文献采用下面的公式来计算： （2-6）2.4.3外电场下HH模型为了更好的描述在外加电场情况下模型的适用范围，我们基于HH模型用膜电位的变化来表示极低频电场在微观下的作用。为了描述精确我们在考虑VE，的情况下对电场的作用重新建立数学模型。于是改进后的细胞模型如下: （2-7）改写为： （2-8
40、）其中： 为电场在膜电容上产生的感生电流。 图三 外电场下HH神经元的电路模型图三是该HH模型所对应的电路模型。我们从图中可以看到，极低频的外电场的引入，它的作用其实实质上就是给细胞膜增添一个电动势，尽管如此，但是HH模型整体结构却没有发生改变，只是改变了各个离子的反电动势，因为改变了离子的反电动势，电容上也随之产生来一个感应电流。从一般情况下的控制系统的角度来看，加入的外电场的参数对于本来的系统来说是一种干扰，可以算作是扰动，以此来分析系统受到扰动后系统的动态特性。2.5神经元的放电节律我们可以把放电节律变化过程分为下面几个状态：(1) 兴奋，膜电位处于静息电位或者波形处于小幅阈下震荡阶段。
41、由于神经元是可兴奋的，如果逐渐增加刺激，细胞膜的膜电位会在某个时刻突然增大，但这种增大并不会毫无限制，达到顶峰之后会放电，然后就会重新回到静息状态，整个过程循环下来就会产生一个动作电位，如图四 a)所示。 图四 三种不同神经元动作电位序列为了便于观察外界刺激对放电序列的影响，我们将引入以下衡量的标准：峰值(Amplitude, AMP)，是指神经元动作电位的极大值的数值。峰峰间期(interspike interval, ISI)，两个相邻的动作电位触发会经过一定的时间，这个时间就叫做峰峰间期，如果这个峰是周期性的，那他们之间的峰峰间期就不会不同，想要研究神经元信息的编码过程，峰峰间期研究至关
42、重要。周期簇放电的峰峰间期分为两种，时间较短的峰峰间期是一个簇内的两个相邻的动作电位的时间差，而时间较长的峰峰间期是两个相邻的簇之间的时间差，我们称之为簇簇间期(interburst interval, IBI)。平均放电率(Firing Rate, FR)，刺激周期内动作电位的个数与刺激周期数的比。假设刺激周期内动作电位的个数为p，外界刺激周期数为q，则平均放电率的等式为： (2-9)对于刺激为周期性刺激，并且响应也是周期性放电的情况下，FR是一个恒值p/q，此时系统放电存在锁相，锁相比为p:q；系统混沌放电时FR是无规则的，系统不存在锁相现象；神经纤维不放电时锁相比可视为0:1。2.6本章
43、总结1、 介绍了生物神经元。2、 并介绍了有关动作电位和静息电位的概念产生及特点，提供了理论基础。3、 阐述了HH经典神经元模型和外电场作用下细胞模型的动力学方程组成及参数选择，通过HH经典神经元模型和外电场作用下细胞模型的结合，得到了外电场下HH模型。4、 介绍了神经元放电节律的几种典型模式和衡量标准。第三章 MATLAB基本知识3.1 MATLAB组成及功能特点3.1.1 MATLAB平台的组成MATLAB不仅是一门编程语言,还是一个集成的软件平台,包含以下几个主要部分：（1） MATLAB语言MATLAB是一种比较高级的编程类语言,它提供了多种不同的数据类型、丰富的运算符和程序控制语句可
44、以供给相关的用户进行使用。所以即使用户有不同的需求，也没有关系，MATLAB适用性很强，不同行业的人都可以用它来完成自己的任务。（2） MATLAB集成工作环境MATLAB集成工作环境有很多部分组成，主要有程序编辑器、变量查看器、系统仿真器和帮助系统等。在这种情况下，用户就可以根据自己的需要编写，调试，运行程序，直到最终得到自己满意的成果。（3） MATLAB数学函数库MATLAB内部自带多种用于计算的函数库，这其中不仅有我们大众所熟知的数学函数，还有众多专业领域的专业数值计算模型，在很多情况下，用户不用自己编写，而只需要根据自己的需要调用就行了，大大节省了工作时间，提高了工作效率。（4） M
45、ATLAB图形系统MATLAB具有图形操作界面，可以将一些数据以三维或者二维图形的形式显示出来，还可以对窗口进行修改，十分方便用户的使用。（5） Simulink交互式仿真环境通过交互式的仿真环境 Simulink,用户可以通过采用图形化的数学模型,完成对各类不同系统的模型进行建立和系统的仿真,仿真结果也能够以直观的图形方式进行显示。Simulink可以接受用户的键盘和鼠标进行输入,也可以通过程序的语句来实现不同数据的交换,应用也很方便和灵活。（6） MATLAB编译器通过编译器,可以将用MATLAB语言编写的程序编译成脱离MATLAB环境的C源代码、动态连接库或者可以独立运行的可执行文件。（
46、7） 应用程序接口API这是 MATLAB的应用程序接口,它提供了 MATLAB和C、 Fortran、VB、VC等多种语言之间的接口程序库,使用户可以通过在这些语言的程序里调用 MATLAB相关的程序。（8） MATLAB工具箱MATLAB自带的工具箱功能齐全，应有尽有。工具箱由各个领域的顶级专家进行编写, 免去了用户自己编辑的麻烦，用户如果需要直接调用就可以了，这样就大大提高了效率，节省了时间和成本。用户可以直接利用基础程序进行更高层次和深度的研究。例如:控制领域的工具箱就有Systen Identification(系统辨识工具箱)、Optimization(最优化工具箱)等。（9） Notebook工具Notebook可以让用户在Word环境中进行使用 MATLAB所提供的各种有用的资源,为用户营造出一种容文字处理、科学计算、工程设