【数据分析】分数阶泄漏积分点火神经元模型Matlab仿真，模拟单个神经元在恒定注入电流下的膜电位变化、放电频率、峰峰间期等电生理特
✅作者简介：热爱科研的Matlab仿真开发者，擅长毕业设计辅导、数学建模、数据处理、建模仿真、程序设计、完整代码获取、论文复现及科研仿真。
🍎 往期回顾关注个人主页：Matlab科研工作室
👇 关注我领取海量matlab电子书和数学建模资料
🍊个人信条：格物致知,完整Matlab代码获取及仿真咨询内容私信。
🔥 内容介绍
一、背景
（一）神经元模型研究的重要性
神经元作为神经系统的基本组成单元，其电生理特性的研究对于理解神经系统的功能和信息处理机制至关重要。神经元模型是研究神经元行为的重要工具，它通过数学方程来描述神经元在不同条件下的电活动，如膜电位变化、放电模式等。深入研究神经元模型不仅有助于揭示神经系统的奥秘，还为神经科学、认知科学以及神经工程等领域的发展提供理论基础。
（二）传统神经元模型的局限性
传统的神经元模型，如经典的泄漏积分点火（Leaky Integrate - and - Fire, LIF）模型，在解释神经元的一些基本电生理现象方面取得了一定成功。LIF 模型假设神经元膜电位以线性方式随时间积分，并在达到阈值时发放一个脉冲，随后膜电位复位。然而，这种模型基于整数阶微积分，过于简化了神经元复杂的生物物理过程。实际神经元的电生理行为受到多种因素的影响，包括细胞膜的电容、离子通道的动力学等，这些过程往往具有非局部和记忆特性，传统整数阶模型难以准确描述。
（三）分数阶微积分引入的意义
分数阶微积分作为整数阶微积分的推广，能够更好地捕捉系统的非局部和记忆特性。将分数阶微积分引入神经元模型，可以更精确地描述神经元电活动中的复杂动态过程。分数阶导数的阶数可以作为一个额外的参数，用于调整模型对神经元电生理行为的刻画，从而使模型更加灵活和准确地反映真实神经元的特性。
二、分数阶泄漏积分点火神经元模型原理
（三）模型特性
- 非局部性与记忆性
：分数阶导数使得膜电位的变化不仅取决于当前时刻的电流和膜电位，还依赖于过去一段时间内的状态。这更符合神经元实际的电生理过程，例如神经元对之前输入信号的记忆会影响其当前的响应。
- 灵活性
：分数阶导数的阶数 α 可以调节模型的动态特性。不同的 α 值对应不同程度的记忆和非局部特性，从而可以模拟出多种不同类型神经元的电生理行为。例如，较小的 α 值可能对应具有更强记忆性的神经元，而较大的 α 值则使模型更接近传统整数阶模型的行为。
三、模拟神经元电生理特性
（一）膜电位变化模拟
在恒定注入电流 I0 下，通过数值求解分数阶微分方程来模拟膜电位 V(t) 的变化。常用的数值方法如 Adams - Bashforth - Moulton 方法可以有效地处理分数阶微分方程。随着时间推移，膜电位在分数阶泄漏和电流注入的共同作用下逐渐变化。当膜电位达到阈值 θ 时，产生放电脉冲，随后复位到 Vreset，形成周期性的膜电位波动。
（二）放电频率计算
放电频率是指单位时间内神经元发放脉冲的次数。通过记录一段时间内神经元的放电次数 N，并除以该时间段 T，可得到放电频率 f=TN。分数阶模型中，由于其非局部和记忆特性，放电频率可能与传统整数阶模型有所不同，且会受到分数阶导数阶数 α 和注入电流强度 I0 的影响。例如，增加注入电流强度通常会提高放电频率；而改变 α 值可能改变膜电位上升和下降的速率，进而影响放电频率。
（三）峰峰间期（ISI）分析
峰峰间期（ISI）是指相邻两个放电脉冲之间的时间间隔。在分数阶泄漏积分点火神经元模型中，ISI 反映了神经元两次放电之间的时间动态。通过分析 ISI 的分布，可以了解神经元放电的规律性和稳定性。不同的 α 值和注入电流条件下，ISI 分布可能呈现出不同的特征。例如，在某些参数设置下，ISI 可能相对稳定，呈现出近似周期性的分布；而在其他情况下，ISI 可能出现较大波动，反映出神经元放电的不规则性。这种对 ISI 的分析有助于深入理解神经元的信息编码和处理机制。
综上所述，分数阶泄漏积分点火神经元模型通过引入分数阶微积分，为研究神经元的电生理特性提供了一个更准确和灵活的框架。通过模拟膜电位变化、放电频率和峰峰间期等特性，该模型能够揭示神经元复杂的动态行为，为神经科学领域的研究提供有力的工具。
⛳️ 运行结果
📣 部分代码
taum=NetProp.TauM;
refrac=NetProp.Refrac/NetProp.dt;
vth=NetProp.vTh;
vpeak=NetProp.vpeak;
v0=NetProp.v0;
vrest=NetProp.vrest;
Namp=NetProp.Noise;
rm=NetProp.Rm;
dt=NetProp.dt;
v=vrest.*ones(length(t),Ncells);
v(1,:) = v0(1,:);
vrest=vrest.*ones(length(t),Ncells);
sp=zeros(length(t),Ncells);
isstillrefrac=zeros(1,Ncells);
rfcounter=zeros(1,Ncells);
Inhib = zeros(length(t),Ncells);