BRIAN中的STDP示例
STDP简单来说就是根据脉冲到达前后神经元的时间差来更新连接这两个神经元之间的突触权重。
普通更新是这样的:
从公式可以直观看出来:权重更新的前提是连接当前突触的前后两个神经元都要发射脉冲,也就是都要被激活,不然的话,整个网络是不存在更新的。
官方文档最后给了一个例子,但是没有说明,这里就简单补充说明一下,最后再验证一下多个脉冲发射的情况。
首先是例子的解释:
start_scope()taupre = taupost = 20*ms
Apre = 0.01
Apost = -Apre*taupre/taupost*1.05
tmax = 50*ms
N = 100# Presynaptic neurons G spike at times from 0 to tmax
# Postsynaptic neurons G spike at times from tmax to 0
# So difference in spike times will vary from -tmax to +tmax
G = NeuronGroup(N, 'tspike:second', threshold='t>tspike', refractory=100*ms)
H = NeuronGroup(N, 'tspike:second', threshold='t>tspike', refractory=100*ms)
G.tspike = 'i*tmax/(N-1)'
H.tspike = '(N-1-i)*tmax/(N-1)'S = Synapses(G, H,'''w : 1dapre/dt = -apre/taupre : 1 (event-driven)dapost/dt = -apost/taupost : 1 (event-driven)''',on_pre='''apre += Aprew = w+apost''',on_post='''apost += Apostw = w+apre''')
S.connect(j='i')run(tmax+1*ms)plot((H.tspike-G.tspike)/ms, S.w)
xlabel(r'$\Delta t$ (ms)')
ylabel(r'$\Delta w$')
axhline(0, ls='-', c='k');
这个例子给的比较直观,只要按照官方的说明文档看过前面的,基本能够理解:
首先,定义了一些相关的基本参数,随后定义了两组神经元,每组神经元的个数是100.
根据连接方式,可以看出来,前后两组神经元连接是一对一的。可以画图看看。
def visualise_connectivity(S):Ns = len(S.source)Nt = len(S.target)figure(figsize=(10, 4))subplot(121)plot(zeros(Ns), arange(Ns), 'ok', ms=10)plot(ones(Nt), arange(Nt), 'ok', ms=10)for i, j in zip(S.i, S.j):plot([0, 1], [i, j], '-k')xticks([0, 1], ['Source', 'Target'])ylabel('Neuron index')xlim(-0.1, 1.1)ylim(-1, max(Ns, Nt))subplot(122)plot(S.i, S.j, 'ok')xlim(-1, Ns)ylim(-1, Nt)xlabel('Source neuron index')ylabel('Target neuron index')
visualise_connectivity(S)
运行上面可视化之后可以看到:
显示有点密集,到那时意思是没有问题的,就是一对一的连接。
这里将每个神经元的发射方式定义为超过固定的时间,也就是说第一个神经元超过0就发放,后面的时间根据表达式G.tspike = 'i*tmax/(N-1)'
来计算。
后面神经元发放时间刚好是相反的。
这样一来,按照规定的时间运行,就可以根据定义的规则对权重进行更新,就有了最后的可视化结果:
横坐标是时间差,纵坐标是更新的权重。【正是因为有了每个神经元只发放一次的约束,所以采能这样直接画出来,不然的话,如果发放多次,需要用多次更新的累加和来更新权重】
如果将运行时间缩短,也就是修改代码run(tmax/2+3*ms)
那么出来的图像就会发生变化:
可以看到,只有中间时间差更新了权重,而其他突触还是初始值0.
为什么是中间的呢?可以这么直观理解一下:
第一组发射的时间假设是0,2,4,6,8
第二组发射的时间就是8,6,4,2,0
如果我们只运行四个步骤,第一组0和第二组0代号的神经元就不会更新。严格来说,当前第四步也不会更新,因为神经元的发射条件是大于当前时间。这也就是为什么给的官方示例代码运行的时候多加了1ms
run(tmax+1*ms)
如果单个神经元发射多次呢?官方没有给代码,我这里给一个根据官方代码改写的示例:
start_scope()taupre = taupost = 20*ms
Apre = 0.01
Apost = -Apre*taupre/taupost*1.05
tmax = 50*ms
N = 1# Presynaptic neurons G spike at times from 0 to tmax
# Postsynaptic neurons G spike at times from tmax to 0
# So difference in spike times will vary from -tmax to +tmax
indices_g = [0,0,0,1]
times_g = [0,10,40,40]*ms
G = SpikeGeneratorGroup(2, indices_g, times_g)
indices_h = [0]
times_h = [30]*ms
H = SpikeGeneratorGroup(1, indices_h, times_h)
# G.tspike = [0,40]*ms
# H.tspike = [30]*msS = Synapses(G, H,'''w : 1dapre/dt = -apre/taupre : 1 (event-driven)dapost/dt = -apost/taupost : 1 (event-driven)''',on_pre='''apre += Aprew = w+apost''',on_post='''apost += Apostw = w+apre''')
S.connect()g_mon = SpikeMonitor(G)run(tmax+1*ms)plot(S.w)
ylabel(r'$\Delta w$')
axhline(0, ls='-', c='k');
为了让同一个神经元在不同的时刻产生脉冲,这里修改原始代码,使用了SpikeGeneratorGroup
,通过参数indices_g = [0,0,0,1] times_g = [0,10,40,40]*ms
让代号为0的神经元在0ms,10ms,40ms,发射三次脉冲。神经元的连接很简单,如下图所示
权重更新会是什么结果呢?
因为只有两个突触,权重看图不如直接看数值:
可以看到,两个权重都变成了负数,因为权重的初始值是0,第二组的单个神经元发射时间是30ms,而第一组中的两个神经在40ms的时候都发射过脉冲,根据赫布规则,也就是STDP,权重应该减弱,但是显然,第一个权重,要大些,这是因为除了40ms的脉冲之外,第一个权重对应的神经元在30ms之前也发射过脉冲,具体的计算过程如下图所示:
分别是两个权重的来源,就是套最开始的公式。
之前对于STDP都是感性的认知,这是第一次用BRIAN来从数值角度计算,所以记录一下,以备未来肯定 会忘记的自己查看。
BRIAN中的STDP示例相关推荐
- python中typing.NamedTuple示例
python中typing.NamedTuple示例 from typing import * # 导入依赖项 class Friend(NamedTuple): """ ...
- python把坐标写入文本_Python实现将数据写入netCDF4中的方法示例
本文实例讲述了Python实现将数据写入netCDF4中的方法.分享给大家供大家参考,具体如下: nc文件为处理气象数据文件.用户可以去https://www.lfd.uci.edu/~gohlke/ ...
- SQL Server:日志备份和差异备份还原中的常见问题示例
日志备份和差异备份还原中的常见问题示例 代码 1 --创建测试 2 CREATE DATABASE db 3 GO 4 5 --正常备份 6 BACKUP DATABASE db TO ...
- python画图代码星星-Python中turtle作图示例
在Python里,海龟不仅可以画简单的黑线,还可以用它画更复杂的几何图形,用不同的颜色,甚至还可以给形状填色. 一.从基本的正方形开始 引入turtle模块并创建Pen对象: >>> ...
- 【Android 逆向】Android 进程注入工具开发 ( 系统调用 | Android NDK 中的系统调用示例 )
文章目录 一.系统调用 二.Android NDK 中的系统调用示例 一.系统调用 在 " 用户层 " , 运行的都是用户应用程序 ; 用户层 下面 是 驱动层 , 驱动层 下面是 ...
- 在 wxWidgets 中播放声音的示例
在 wxWidgets 中播放声音的示例 在 wxWidgets 中播放声音的示例 在 wxWidgets 中播放声音的示例 #include "wx/wxprec.h" #ifn ...
- boost::gil模块数字扩展中的 resize_view() 示例
boost::gil模块数字扩展中的 resize_view(示例 实现功能 C++实现代码 实现功能 boost::gil模块数字扩展中的 resize_view() 示例 C++实现代码 #inc ...
- boost::gil模块数字扩展中的 resample_pixels() 示例
boost::gil模块数字扩展中的 resample_pixels 示例 实现功能 C++实现代码 实现功能 boost::gil模块数字扩展中的 resample_pixels() 示例 C++实 ...
- java中get接口示例_Java即时类| 带示例的get()方法
java中get接口示例 即时类的get()方法 (Instant Class get() method) get() method is available in java.time package ...
最新文章
- Django之BBS博客项目
- 2020年 HackerEarth 调查:Go 语言成为最受欢迎的语言 开发者最关心职业发展路径...
- 利用Jquery处理跨域请求
- 32位算术逻辑运算单元alu设计_依据基本原理构建现代计算机(二)—— ALU
- 用于构建高级媒体应用程序的工具
- 答非所问:产品质量怎样?发现了很多BUG
- 高清电脑壁纸2k图片,4k批量免费下载,壁纸打包下载,美女壁纸
- 二维码图片生成(带文字显示)
- [51Nod]NOIP2018提高组省一冲奖班模测训练(一)题解
- 国内外镜像下载合集(详细最终版)
- linux如何ping多个ip,如何同时Ping多个IP地址并保存记录 | 个人信息技术网
- 响应式织梦模板家政服务保姆月嫂类网站
- AUTOCAD——基线标注
- 在eclipse上配置tomcat,安装eclipse的JST Server Adapters Extensions时报错 Anerroroccurredwhilecollectingitems...
- 如何选择合适的渗压计?
- PHP(Laravel)作为客户端调用GRPC
- uniapp 分享缩略图过大怎么办_新版本微信下,如何设置分享到朋友圈的缩略图?...
- Python安装(3.8、3.9, 3.10)、Python的VS Code配置
- 201621123063叶斓浈 第一周实验报告
- 半导体后道工艺及设备