CN105787291B - 一种模拟实现Morris‑Lecar神经元模型的电路 - Google Patents

一种模拟实现Morris‑Lecar神经元模型的电路 Download PDF

Info

Publication number
CN105787291B
CN105787291B CN201610066207.8A CN201610066207A CN105787291B CN 105787291 B CN105787291 B CN 105787291B CN 201610066207 A CN201610066207 A CN 201610066207A CN 105787291 B CN105787291 B CN 105787291B
Authority
CN
China
Prior art keywords
input
circuit unit
operational amplifier
input terminal
resistance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610066207.8A
Other languages
English (en)
Other versions
CN105787291A (zh
Inventor
胡晓宇
刘崇新
倪骏康
沈天时
李傲岸
刘航
李石磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201610066207.8A priority Critical patent/CN105787291B/zh
Publication of CN105787291A publication Critical patent/CN105787291A/zh
Application granted granted Critical
Publication of CN105787291B publication Critical patent/CN105787291B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/06Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons
    • G06N3/063Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using electronic means
    • G06N3/065Analogue means
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B5/00ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Data Mining & Analysis (AREA)
  • Biotechnology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Physiology (AREA)
  • Neurology (AREA)
  • Artificial Intelligence (AREA)
  • Computational Linguistics (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Amplifiers (AREA)

Abstract

本发明公开了一种模拟实现Morris‑Lecar神经元模型的电路,所述电路包括:第一双曲正切函数电路单元、第二双曲正切函数电路单元、双曲余弦函数电路单元、第一比例积分电路单元和第二比例积分电路单元。所述电路由运算放大器、晶体管、加法器、乘法器、电阻等基本的模拟电子器件组成。本发明所实现电路结构简单,电路参数调节灵活方便,电路输出精度高,随着电路参数的变化,本发明所实现的电路可以模拟Morris‑Lecar神经元模型静息态、峰放电(spiking)、簇放电(bursting)等多种放电形式。

Description

一种模拟实现Morris-Lecar神经元模型的电路
技术领域:
本发明属于神经元电路模拟技术领域,特别涉及一种模拟实现Morris-Lecar神经元模型的电路。
背景技术:
神经元是神经系统的基本单位,人类神经系统的神经元数量数以亿计。神经元在不同的外界环境下或受到外界刺激时,能够产生电信号,这些电信号可以在不同的神经元之间传递,神经系统通过神经元的这种电信号接受或者发送功能感受外界环境变化,或者传递大脑发送的信号指令,以实现与外界环境的交互。截止目前,相关学者建立了多种数学模型来描述神经元的动力学行为,如Hodgkin-Huxley模型、FitzHugh-Naguma模型、Hindmarsh-Rose模型、Morris-Lecar模型等。通过分析不同条件下神经元模型的动力学行为变化,可以预测神经元不同环境下的放电行为,已有的研究表明,人类所患帕金森症、癫痫等疾病与神经系统的异常放电有直接关系,因此,对单个神经元、耦合神经元以及神经元网络动力学行为的研究能够为这些疾病的治疗提供潜在的理论支持。
神经元模型一般由多个微分方程组成,随着神经元数量的增加,对神经元系统动力学行为分析的计算负担将大大增加,计算的实时性将不能得到保障,因此,通过电路来模拟神经元的放电行为,成为解决实时性问题的重要途径。由于神经元模型中多包含指数函数、双曲函数等非线性函数,这些函数很难用现有的模拟器件来直接实现,现有的神经元电路模拟通常采用将这些非线性函数线性化,或者采用动力学行为相似的函数来代替这些非线性函数的方法来实现,这也造成了误差较大、精度不高等问题,因此,设计新颖的模拟电路,来精确实现这些非线性函数,进一步实现神经元模型的模拟意义重大。
发明内容:
本发明的目的是针对现有技术存在的问题,提供了一种模拟实现Morris-Lecar神经元模型的电路。
为达到上述目的,本发明采用如下技术方案来实现的:
一种模拟实现Morris-Lecar神经元模型的电路,包括第一双曲正切函数电路单元、第二双曲正切函数电路单元、双曲余弦函数电路单元、第一比例积分电路单元和第二比例积分电路单元;其中,
所述的第一双曲正切函数电路单元和第二双曲正切函数电路单元均由输入函数运算模块、双曲正切函数运算模块、双端输入转单端输出模块和加法运算模块组成;
所述双曲余弦函数电路单元由输入函数运算模块、第一指数函数电路单元、第二指数电路单元和加法运算模块组成;
所述第一比例积分电路单元由乘法运算模块、电流源、电压源和比例积分模块组成,其包含六个输入端,输入端D1A、输入端D1B、输入端D2A、输入端D2B、输入端D3、输入端D4和一个输出端,输出端电压即为神经元模型的膜电位;
所述第二比例积分电路单元由减法运算模块、乘法运算模块、积分运算模块和比例运算模块组成,其包含两个输入端,输入端B21,输入端B22和一个输出端;
第一双曲正切函数电路单元的输入端连接第一比例积分电路单元的输出端,第一双曲正切函数电路单元的输出端连接第一比例积分电路单元的输入端D1A;第二双曲正切函数电路单元的输入端连接第一比例积分电路单元的输出端,第二双曲正切函数电路单元的输出端连接第二比例积分电路单元的输入端B21;第二比例积分电路单元的输出端连接第一比例积分电路单元的输入端D2A;双曲余弦函数电路单元的输入端连接第一比例积分电路单元的输出端,双曲余弦函数电路单元的输出端连接第二比例积分电路单元的输入端B22;第一比例积分电路单元的输入端D1B连接第一比例积分电路单元的输出端;第一比例积分电路单元的输入端D2B连接第一比例积分电路单元的输出端;第一比例积分电路单元的输入端D3连接第一比例积分电路单元的输出端;第一比例积分电路单元的输入端D4连接电流源。
本发明进一步的改进在于,第一双曲正切函数电路单元中,输入函数运算模块由运算放大器、电压源及电阻组成;
双曲正切函数运算模块由双极性晶体管对、电压源、电流源及电阻组成;
双端输入转单端输出模块由运算放大器及电阻组成;
加法运算模块由加法器、电压源及电阻组成;
第一双曲正切函数电路单元的输入端即其输入函数运算模块输入端,其输入经电阻R1连接到运算放大器U1的反相输入端,运算放大器U1的输出经电阻R2反馈连接到运算放大器U1的反相输入端,同时,运算放大器U1的输出端连接到双曲正切函数运算模块的输入端,运算放大器U1的同相输入端连接电压源Vth1的负极,电压源Vth1的正极接地;双曲正切函数运算模块由一个输入端和两个输出端,其输入端即是双极性晶体管Q1的基极,双极型晶体管Q1和Q2的集电极分别经电阻R3和R4连接至电压源V1的正极,电压源V1的负极接地,双极型晶体管Q1和Q2的发射极连接在一起并经电流源I1接地,双极性晶体管Q2的基极接地,双曲正切函数运算模块的两个输出端分别经Q1和Q2的集电极引出;双端输入转单端输出模块的输入端分别连接双曲正切函数运算模块的两个输出端,两个输入信号分别经电阻R7和R5连接到运算放大器U2同相输入端和反相输入端,运算放大器U2的同相输入端同时经电阻R8接地,运算放大器U2的输出端经电阻R6反馈连接至运算放大器U2的反相输入端,双端输入转单端输出模块的输出端连接至加法运算模块的输入端;加法运算模块输入端即为加法器的输入端,加法器的另一输入端连接至电压源V2的正极,电压源V2的负极接地,加法器的输出端即为第一双曲正切函数电路单元的输出端。
本发明进一步的改进在于,第二双曲正切函数电路单元中,输入函数运算模块由运算放大器、电压源及电阻组成;
双曲正切函数运算模块由双极性晶体管对、电压源、电流源及电阻组成;
双端输入转单端输出模块由运算放大器及电阻组成
加法运算模块由加法器、电压源及电阻组成;
第二双曲正切函数电路单元的输入端即其输入函数运算模块输入端,其输入经电阻R9连接到运算放大器U3的反相输入端,运算放大器U3的输出经电阻R10反馈连接到运算放大器U3的反相输入端,同时,运算放大器U3的输出端连接到双曲正切函数运算模块的输入端,运算放大器U3的同相输入端连接电压源Vth2的正极,电压源Vth2的负极接地;双曲正切函数运算模块由一个输入端和两个输出端,其输入端即是双极性晶体管Q3的基极,双极型晶体管Q3和Q4的集电极分别经电阻R3和R4连接至电压源V3的正极,电压源V3的负极接地,双极型晶体管Q3和Q4的发射极连接在一起并经电流源I2接地,双极性晶体管Q4的基极接地,双曲正切函数运算模块的两个输出端分别经Q3和Q4的集电极引出;双端输入转单端输出模块的输入端分别连接双曲正切函数运算模块的两个输出端,两个输入信号分别经电阻R15和R13连接到运算放大器U4同相输入端和反相输入端,运算放大器U4的同相输入端同时经电阻R16接地,运算放大器U4的输出端经电阻R14反馈连接至运算放大器U4的反相输入端,双端输入转单端输出模块的输出端连接至加法运算模块的输入端;加法运算模块输入端即为加法器的输入端,加法器的另一输入端连接至电压源V4的正极,电压源V4的负极接地,加法器的输出端即为第二双曲正切函数电路单元的输出端。
本发明进一步的改进在于,双曲余弦函数电路单元中,输入函数运算模块由运算放大器、电压源及电阻组成;
第一指数函数电路单元由运算放大器、电压源、双极性晶体管及电阻组成;
第二指数函数电路单元由电压反相电路、运算放大器、电压源、双极性晶体管及电阻组成;
加法运算模块由加法器实现;
双曲余弦函数电路单元的输入端即其输入函数运算模块输入端,其输入经电阻R22连接到运算放大器U7的反相输入端,运算放大器U7的输出经电阻R23反馈连接到运算放大器U7的反相输入端,同时,运算放大器U7的输出端连接到双曲正切函数运算模块的输入端,运算放大器U7的同相输入端连接电压源Vch的正极,电压源Vch的负极接地;第一指数函数电路单元的输入端连接输入函数运算模块的输出端,第一指数函数电路单元的输出端连接加法运算模块的一个输入端;第二指数函数电路单元的输入端连接输入函数运算模块的输出端,第二指数函数电路单元的输出端连接加法运算模块的另一个输入端;加法运算模块的两个输入端分别连接第一指数函数电路单元的输出端和第二指数函数电路单元的输出端,加法运算模块的输出端连接第二比例积分电路单元的输入端B22。
本发明进一步的改进在于,双曲余弦函数电路单元中实现双曲余弦函数运算的电路是由两个指数函数电路组成,两个指数函数电路的输入量幅值相等符号相反。
本发明进一步的改进在于,第一指数函数电路单元和第二指数函数电路单元能够实现指数运算且输入值的符号是正值或者是负值。
本发明进一步的改进在于,第一指数函数电路单元的输入端即双极性晶体管Q5的基极,双极性晶体管Q5的集电极连接至运算放大器U8的反相输入端,同时运算放大器U8的反相输入端经电阻R24连接至电压源V5的正极,电压源V5的负极接地,运算放大器U8的同相输入端经电阻R25接地,运算放大器U8的输出端连接至双极性晶体管Q5的发射极,同时双极性晶体管Q5的发射极连接至双极性晶体管Q6的发射极,双极性晶体管Q6的集电极连接至运算放大器U9的反相输入端,运算放大器U9的同相输入端经电阻R28接地,运算放大器U9的输出端经电阻R27连接至运算放大器U9的反相输入端,运算放大器U9的输出端即为第一指数函数电路单元的输出端。
本发明进一步的改进在于,第二指数函数电路单元的输入端即为电压反相电路的输入端,第二指数函数电路单元的输入端连接输入函数运算模块的输出端,电压反相电路的输入信号经电阻R29连接至运算放大器U10的反相输入端,运算放大器U10的同相输出端经电阻R31接地,运算放大器U10的输出端经电阻R30连接至算放大器U10的反相输入端,运算放大器的输出端即为电压反相器的输出端,电压反相器的输出端连接至双极性晶体管Q7的基极,双极性晶体管Q7的集电极连接至运算放大器U11的反相输入端,同时运算放大器U11的反相输入端经电阻R32连接至电压源V6的正极,电压源V6的负极接地,运算放大器U11的同相输入端经电阻R33接地,运算放大器U11的输出端连接至双极性晶体管Q7的发射极,同时双极性晶体管Q7的发射极连接至双极性晶体管Q8的发射极,双极性晶体管Q8的集电极连接至运算放大器U12的反相输入端,运算放大器U12的同相输入端经电阻R36接地,运算放大器U12的输出端经电阻R37连接至运算放大器U12的反相输入端,运算放大器U12的输出端即为第二指数函数电路单元的输出端。
本发明进一步的改进在于,第一比例积分电路单元的输入端D1A连接第一双曲正切函数电路单元的输出端,其输入端D2A连接第二比例积分电路单元的输出端,其输入端D3连接第一比例积分电路单元的输出端,其输入端D1B连接第一比例积分电路单元的输出端,其输入端D2B连接第一比例积分电路单元的输出端;输入端D1B的输入信号连接至电压源V7的正极,电压源V7的负极连接乘法器的一个输入端,乘法器的另一个输入端连接输入端D1A的输入信号,乘法器的输出端经电阻R38连接至运算放大器U13的反相输入端;输入端D2B的输入信号连接至电压源V8的负极,电压源V8的正极连接乘法器的一个输入端,乘法器的另一个输入端连接输入端D2A的输入信号,乘法器的输出端经电阻R39连接至运算放大器U13的反相输入端;输入端D3的输入信号连接至电压源V9负极,电压源V9的正极经电阻R37连接至运算放大器U13的反相输入端;输入端D4连接至电流源I3的负极,电流源I3的正极接地;运算放大器U13的输出端经电容C2连接至运算放大器U13的反相输入端。
相对于现有技术,本发明提出一种模拟实现Morris-Lecar神经元模型的电路,采用普通的模拟电子器件,按照Morris-Lecar神经元模型的数学表达式,巧妙的设计电路使其能够完美的呈现数学表达式中的非线性运算关系,进而使得电路输出能够精确重现Morris-Lecar神经元模型的各种放电模式。
Morris-Lecar神经元模型中的双曲正切函数运算采用所述的第一双曲正切函数电路单元和第二双曲正切函数电路单元来实现,双曲正切函数运算借助于一组对偶双极性晶体管对来实现,晶体管对的对偶结构保证了电路在输入为零时输出也为零,即有效的抑制了零点漂移,极大的提高了双曲正切函数模拟电路实现的精度。
在模拟电路中,双曲余弦函数很难采用模拟器件直接实现,本发明中,Morris-Lecar神经元模型中的双曲余弦函数运算采用所述的双曲余弦函数电路单元来实现,所设计的电路通过将双曲余弦函数运算拆解成两个指数运算之和,利用双极性晶体管中电压电流的指数关系来进一步模拟双曲余弦函数运算,使双曲余弦函数电路单元结构简单,动态调节范围宽,精度较高。
发明所实现电路采用普通的模拟电子器件,造价低效果好,并且电路结构简单,电路参数调节灵活方便。得益于对双曲正切函数电路单元与双曲余弦函数的电路单元的巧妙设计,避免将模型中的非线性函数进行线性化或者化简替代,使得电路输出稳定,精度较高。随着电路参数的变化,本发明所实现的电路可以模拟Morris-Lecar神经元模型静息态、峰放电(spiking)、簇放电(bursting)等多种放电形式,为实现耦合神经元以及神经元网络的实时分析、计算及应用提供了良好的支撑。
综上所述,本发明提出的一种模拟实现Morris-Lecar神经元模型的电路,未对系统中的非线性函数进行线性化或者化简替代,模拟电路是对Morris-Lecar神经元模型的完全实现,保证了电路能够准确呈现Morris-Lecar神经元模型的各种放电行为,为实现耦合神经元以及神经元网络的实时分析、计算及应用提供了良好的支撑。
附图说明:
图1为模拟实现Morris-Lecar神经元模型的电路图;
图2为刺激电流I=30μA时神经元静息态波形图;
图3为刺激电流I=45μA时神经元峰放电波形图;
图4为刺激电流I=100μA时神经元峰放电波形图;
图5为刺激电流i=50(1+sin(2πft))μA,f=1时神经元簇放电波形图。
具体实施方式:
为了详细说明本发明的技术内容、电路结构、所实现的目的及效果,以下结合具体实施方式并配合附图给予详细说明。
Morris-Lecar神经元模型由包含两个微分方程的方程组组成,这一二维微分方程组如下所示:
式中:V和W为系统变量,分别表示神经元膜电位和离子通道门电位,C为神经元膜电容,gCa、gK和gL分别表示钙离子通道、钾离子通道和漏离子通道最大电导,VCa、VK和VL分别表示钙离子通道、钾离子通道和漏离子通道稳态电位,M(V)和W(V)分别表示钙离子通道和钾离子通道打开概率的稳态值,V1、V2、V3和V4为系统参数,I为外界刺激电流。
图1为本发明所实现的Morris-Lecar神经元模型的电路图,其中A区域为所述的第一双曲正切函数电路单元,B1区域为所述的第二双曲正切函数电路单元,B2区域为所述的第二比例积分电路单元,C区域为双曲余弦函数电路单元,C1区域和C2区域包括在C区域内,C1区域为第一指数函数电路单元,C2区域为第二指数函数电路单元,D区域为第一比例积分电路单元。
本发明一种模拟实现Morris-Lecar神经元模型的电路,包括第一双曲正切函数电路单元、第二双曲正切函数电路单元、双曲余弦函数电路单元、第一比例积分电路单元和第二比例积分电路单元。
所述的第一双曲正切函数电路单元和第二双曲正切函数电路单元均由输入函数运算模块、双曲正切函数运算模块、双端输入转单端输出模块和加法运算模块组成;
所述双曲余弦函数电路单元由输入函数运算模块、第一指数函数电路单元、第二指数电路单元和加法运算模块组成;
所述第一比例积分电路单元由乘法运算模块、电流源、电压源和比例积分模块组成,其包含六个输入端,输入端D1A、输入端D1B、输入端D2A、输入端D2B、输入端D3、输入端D4和一个输出端,输出端电压即为神经元模型的膜电位;
所述第二比例积分电路单元由减法运算模块、乘法运算模块、积分运算模块和比例运算模块组成,其包含两个输入端,输入端B21,输入端B22和一个输出端。
具体的连接方式是:第一双曲正切函数电路单元的输入端连接第一比例积分电路单元的输出端,第一双曲正切函数电路单元的输出端连接第一比例积分电路单元的输入端D1A;第二双曲正切函数电路单元的输入端连接第一比例积分电路单元的输出端,第二双曲正切函数电路单元的输出端连接第二比例积分电路单元的输入端B21;第二比例积分电路单元的输出端连接第一比例积分电路单元的输入端D2A;双曲余弦函数电路单元的输入端连接第一比例积分电路单元的输出端,双曲余弦函数电路单元的输出端连接第二比例积分电路单元的输入端B22;第一比例积分电路单元的输入端D1B连接第一比例积分电路单元的输出端;第一比例积分电路单元的输入端D2B连接第一比例积分电路单元的输出端;第一比例积分电路单元的输入端D3连接第一比例积分电路单元的输出端;第一比例积分电路单元的输入端D4连接电流源。
所述的第一双曲正切函数电路单元中,输入函数运算模块由运算放大器、电压源及电阻组成,双曲正切函数运算模块由双极性晶体管对、电压源、电流源及电阻组成,双端输入转单端输出模块由运算放大器及电阻组成,加法运算模块由加法器、电压源及电阻组成;第一双曲正切函数电路单元的输入端即其输入函数运算模块输入端,其输入经电阻R1连接到运算放大器U1的反相输入端,运算放大器U1的输出经电阻R2反馈连接到运算放大器U1的反相输入端,同时,运算放大器U1的输出端连接到双曲正切函数运算模块的输入端,运算放大器U1的同相输入端连接电压源Vth1的负极,电压源Vth1的正极接地;双曲正切函数运算模块由一个输入端和两个输出端,其输入端即是双极性晶体管Q1的基极,双极型晶体管Q1和Q2的集电极分别经电阻R3和R4连接至电压源V1的正极,电压源V1的负极接地,双极型晶体管Q1和Q2的发射极连接在一起并经电流源I1接地,双极性晶体管Q2的基极接地,双曲正切函数运算模块的两个输出端分别经Q1和Q2的集电极引出;双端输入转单端输出模块的输入端分别连接双曲正切函数运算模块的两个输出端,两个输入信号分别经电阻R7和R5连接到运算放大器U2同相输入端和反相输入端,运算放大器U2的同相输入端同时经电阻R8接地,运算放大器U2的输出端经电阻R6反馈连接至运算放大器U2的反相输入端,双端输入转单端输出模块的输出端连接至加法运算模块的输入端;加法运算模块输入端即为加法器的输入端,加法器的另一输入端连接至电压源V2的正极,电压源V2的负极接地,加法器的输出端即为第一双曲正切函数电路单元的输出端。
双极型晶体管Q1和Q2可选用NPN型晶体管Q2N222,运算放大器U1和U2采用型号为μA741,电压源Vth1=0.89mV,V1=12V,V2=0.5V,电流源I1=0.55mA,电阻R2=28.9kΩ,R1=R5=R6=R7=R8=10kΩ,R3=R4=1kΩ。
所述的第二双曲正切函数电路单元中,输入函数运算模块由运算放大器、电压源及电阻组成,双曲正切函数运算模块由双极性晶体管对、电压源、电流源及电阻组成,双端输入转单端输出模块由运算放大器及电阻组成,加法运算模块由加法器、电压源及电阻组成;第二双曲正切函数电路单元的输入端即其输入函数运算模块输入端,其输入经电阻R9连接到运算放大器U3的反相输入端,运算放大器U3的输出经电阻R10反馈连接到运算放大器U3的反相输入端,同时,运算放大器U3的输出端连接到双曲正切函数运算模块的输入端,运算放大器U3的同相输入端连接电压源Vth2的正极,电压源Vth2的负极接地;双曲正切函数运算模块由一个输入端和两个输出端,其输入端即是双极性晶体管Q3的基极,双极型晶体管Q3和Q4的集电极分别经电阻R3和R4连接至电压源V3的正极,电压源V3的负极接地,双极型晶体管Q3和Q4的发射极连接在一起并经电流源I2接地,双极性晶体管Q4的基极接地,双曲正切函数运算模块的两个输出端分别经Q3和Q4的集电极引出;双端输入转单端输出模块的输入端分别连接双曲正切函数运算模块的两个输出端,两个输入信号分别经电阻R15和R13连接到运算放大器U4同相输入端和反相输入端,运算放大器U4的同相输入端同时经电阻R16接地,运算放大器U4的输出端经电阻R14反馈连接至运算放大器U4的反相输入端,双端输入转单端输出模块的输出端连接至加法运算模块的输入端;加法运算模块输入端即为加法器的输入端,加法器的另一输入端连接至电压源V4的正极,电压源V4的负极接地,加法器的输出端即为第二双曲正切函数电路单元的输出端。
双极型晶体管Q3和Q4可选用NPN型晶体管Q2N222,运算放大器U3和U4采用型号为μA741,电压源Vth2=8.9mV,V3=12V,V4=0.5V,电流源I2=0.55mA,电阻R10=30kΩ,R9=R16=R13=R14=R15=10kΩ,R11=R12=1kΩ。
所述的第一双曲正切函数电路单元和第二双曲正切函数电路单元,其优点在于:双曲正切函数运算的实现采用一对双极性晶体管,由于双极性晶体管对的结构对称,可以保证双极性晶体管对在输入信号为零时输出信号也为零,即零点漂移为零,进一步保证了整个电路单元具有较小的级联误差。
所述的第二比例积分电路单元,其具有两个输入端(B21和B22)和一个输出端,第二比例积分电路单元的输入端B21连接至第二指数函数电路单元的输出端,第二比例积分电路的输入端B22连接至双曲余弦函数电路单元的输出端,第二比例积分电路单元的输出端连接至第一比例积分电路单元的输入端D2A,第二比例积分电路单元的输入端B21即为减法运算模块的正相输入端,减法运算模块的输出端连接至乘法运算模块的一个输入端,乘法运算模块的另一个输入端即为第二比例运算电路单元的输入端B22;乘法运算模块的输出端连接至积分运算模块的输入端,输入信号经电阻R17连接至运算放大器U5的反相输入端,运算放大器U5的同相输入端经电阻R18接地,运算放大器U5的输出端经电容C1连接至运算放大器U5的反相输入端,运算放大器U5的输出即为积分运算模块的输出;积分运算模块的输出连接至比例运算模块的输入,输入信号经电阻R19连接至运算放大器U6的反相输入端,运算放大器U6的同相输入端经电阻R21接地,运算放大器U6的输出端经电容R20连接至运算放大器U6的反相输入端,运算放大器U6的输出即为第二比例积分电路单元的输出。
运算放大器采用型号为μA741,乘法器型号为AD633,R17=R18=1kΩ,C1=1μF,R19=150kΩ,R20=R21=10kΩ。
所述的双曲余弦函数电路单元中,输入函数运算模块由运算放大器、电压源及电阻组成,第一指数函数电路单元由运算放大器、电压源、双极性晶体管及电阻组成,第二指数函数电路单元由电压反相电路、运算放大器、电压源、双极性晶体管及电阻组成,加法运算模块由加法器实现;双曲余弦函数电路单元的输入端即其输入函数运算模块输入端,其输入经电阻R22连接到运算放大器U7的反相输入端,运算放大器U7的输出经电阻R23反馈连接到运算放大器U7的反相输入端,同时,运算放大器U7的输出端连接到双曲正切函数运算模块的输入端,运算放大器U7的同相输入端连接电压源Vch的正极,电压源Vch的负极接地;第一指数函数电路单元的输入端连接输入函数运算模块的输出端,第一指数函数电路单元的输出端连接加法运算模块的一个输入端;第二指数函数电路单元的输入端连接输入函数运算模块的输出端,第二指数函数电路单元的输出端连接加法运算模块的另一个输入端;加法运算模块的两个输入端分别连接第一指数函数电路单元的输出端和第二指数函数电路单元的输出端,加法运算模块的输出端连接第二比例积分电路单元的输入端B22。
所述的双曲余弦函数电路单元中,实现双曲余弦函数运算的电路是由两个指数函数电路组成,两个指数函数电路的输入量幅值相等符号相反。
所述的双曲余弦函数电路单元,其优点在于:将双曲余弦函数运算拆解成两个指数运算之和,利用双极性晶体管中电压电流的指数关系来进一步模拟双曲余弦函数运算,使双曲余弦函数电路单元结构简单,动态范围宽,精度较高。
所述的第一指数函数电路单元和第二指数函数电路单元能够实现指数运算且输入值的符号既可以是正值也可以是负值,
所述的第一指数函数电路单元的输入端即双极性晶体管Q5的基极,双极性晶体管Q5的集电极连接至运算放大器U8的反相输入端,同时运算放大器U8的反相输入端经电阻R24连接至电压源V5的正极,电压源V5的负极接地,运算放大器U8的同相输入端经电阻R25接地,运算放大器U8的输出端连接至双极性晶体管Q5的发射极,同时双极性晶体管Q5的发射极连接至双极性晶体管Q6的发射极,双极性晶体管Q6的集电极连接至运算放大器U9的反相输入端,运算放大器U9的同相输入端经电阻R28接地,运算放大器U9的输出端经电阻R27连接至运算放大器U9的反相输入端,运算放大器U9的输出端即为第一指数函数电路单元的输出端。
所述的第二指数函数电路单元,第二指数函数电路单元的输入端即为电压反相电路的输入端,第二指数函数电路单元的输入端连接输入函数运算模块的输出端,电压反相电路的输入信号经电阻R29连接至运算放大器U10的反相输入端,运算放大器U10的同相输出端经电阻R31接地,运算放大器U10的输出端经电阻R30连接至算放大器U10的反相输入端,运算放大器的输出端即为电压反相器的输出端,电压反相器的输出端连接至双极性晶体管Q7的基极,双极性晶体管Q7的集电极连接至运算放大器U11的反相输入端,同时运算放大器U11的反相输入端经电阻R32连接至电压源V6的正极,电压源V6的负极接地,运算放大器U11的同相输入端经电阻R33接地,运算放大器U11的输出端连接至双极性晶体管Q7的发射极,同时双极性晶体管Q7的发射极连接至双极性晶体管Q8的发射极,双极性晶体管Q8的集电极连接至运算放大器U12的反相输入端,运算放大器U12的同相输入端经电阻R36接地,运算放大器U12的输出端经电阻R37连接至运算放大器U12的反相输入端,运算放大器U12的输出端即为第二指数函数电路单元的输出端。
V5=V6=0.5V,R22=10kΩ,R24=R25=R26=R27=R28=10kΩ
R29=R30=R31=R32=R33=R34=R35=R36=10kΩ.R23=7.47kΩ,Vch=5.13mV,双极型晶体管均选用NPN型晶体管Q2N222,运算放大器采用型号为μA741。
所述的第一比例积分电路单元,其输入端D1A连接第一双曲正切函数电路单元的输出端,其输入端D2A连接第二比例积分电路单元的输出端,其输入端D3连接第一比例积分电路单元的输出端,其输入端D1B连接第一比例积分电路单元的输出端,其输入端D2B连接第一比例积分电路单元的输出端;输入端D1B的输入信号连接至电压源V7的正极,电压源V7的负极连接乘法器的一个输入端,乘法器的另一个输入端连接输入端D1A的输入信号,乘法器的输出端经电阻R38连接至运算放大器U13的反相输入端;输入端D2B的输入信号连接至电压源V8的负极,电压源V8的正极连接乘法器的一个输入端,乘法器的另一个输入端连接输入端D2A的输入信号,乘法器的输出端经电阻R39连接至运算放大器U13的反相输入端;输入端D3的输入信号连接至电压源V9负极,电压源V9的正极经电阻R37连接至运算放大器U13的反相输入端;输入端D4连接至电流源I3的负极,电流源I3的正极接地;运算放大器U13的输出端经电容C2连接至运算放大器U13的反相输入端。电流源I3为外界刺激电流,所外界刺激电流大小的不同,神经元将呈现不同的放电行为,如图2、图3、图4和图5所示。
运算放大器采用型号为μA741,
R17=R18=1kΩ,C1=1μF,R19=150kΩ,R20=R21=10kΩ,V7=120mV,V8=80mV,V9=60mV,R37=0.5kΩ,R38=0.25kΩ,R39=0.125kΩ,R40=1kΩ,C2=5μF。
本发明提供的一种模拟实现Morris-Lecar神经元模型的电路,未对系统中的非线性函数进行线性化或者化简替代,模拟电路是对Morris-Lecar神经元模型的完全实现,保证了电路能够准确呈现Morris-Lecar神经元模型的各种放电行为,为实现耦合神经元以及神经元网络的实时分析、计算及应用提供了良好的支撑。

Claims (9)

1.一种模拟实现Morris-Lecar神经元模型的电路,其特征在于,包括第一双曲正切函数电路单元、第二双曲正切函数电路单元、双曲余弦函数电路单元、第一比例积分电路单元和第二比例积分电路单元;其中,
所述的第一双曲正切函数电路单元和第二双曲正切函数电路单元均由输入函数运算模块、双曲正切函数运算模块、双端输入转单端输出模块和加法运算模块组成;
所述双曲余弦函数电路单元由输入函数运算模块、第一指数函数电路单元、第二指数电路单元和加法运算模块组成;
所述第一比例积分电路单元由乘法运算模块、电流源、电压源和比例积分模块组成,其包含六个输入端,输入端D1A、输入端D1B、输入端D2A、输入端D2B、输入端D3、输入端D4和一个输出端,输出端电压即为神经元模型的膜电位;
所述第二比例积分电路单元由减法运算模块、乘法运算模块、积分运算模块和比例运算模块组成,其包含两个输入端,输入端B21,输入端B22和一个输出端;
第一双曲正切函数电路单元的输入端连接第一比例积分电路单元的输出端,第一双曲正切函数电路单元的输出端连接第一比例积分电路单元的输入端D1A;第二双曲正切函数电路单元的输入端连接第一比例积分电路单元的输出端,第二双曲正切函数电路单元的输出端连接第二比例积分电路单元的输入端B21;第二比例积分电路单元的输出端连接第一比例积分电路单元的输入端D2A;双曲余弦函数电路单元的输入端连接第一比例积分电路单元的输出端,双曲余弦函数电路单元的输出端连接第二比例积分电路单元的输入端B22;第一比例积分电路单元的输入端D1B连接第一比例积分电路单元的输出端;第一比例积分电路单元的输入端D2B连接第一比例积分电路单元的输出端;第一比例积分电路单元的输入端D3连接第一比例积分电路单元的输出端;第一比例积分电路单元的输入端D4连接电流源。
2.根据权利要求1所述的一种模拟实现Morris-Lecar神经元模型的电路,其特征在于,第一双曲正切函数电路单元中,输入函数运算模块由运算放大器、电压源及电阻组成;
双曲正切函数运算模块由双极性晶体管对、电压源、电流源及电阻组成;
双端输入转单端输出模块由运算放大器及电阻组成;
加法运算模块由加法器、电压源及电阻组成;
第一双曲正切函数电路单元的输入端即其输入函数运算模块输入端,其输入经电阻R1连接到运算放大器U1的反相输入端,运算放大器U1的输出经电阻R2反馈连接到运算放大器U1的反相输入端,同时,运算放大器U1的输出端连接到双曲正切函数运算模块的输入端,运算放大器U1的同相输入端连接电压源Vth1的负极,电压源Vth1的正极接地;双曲正切函数运算模块由一个输入端和两个输出端,其输入端即是双极性晶体管Q1的基极,双极型晶体管Q1和Q2的集电极分别经电阻R3和R4连接至电压源V1的正极,电压源V1的负极接地,双极型晶体管Q1和Q2的发射极连接在一起并经电流源I1接地,双极性晶体管Q2的基极接地,双曲正切函数运算模块的两个输出端分别经Q1和Q2的集电极引出;双端输入转单端输出模块的输入端分别连接双曲正切函数运算模块的两个输出端,两个输入信号分别经电阻R7和R5连接到运算放大器U2同相输入端和反相输入端,运算放大器U2的同相输入端同时经电阻R8接地,运算放大器U2的输出端经电阻R6反馈连接至运算放大器U2的反相输入端,双端输入转单端输出模块的输出端连接至加法运算模块的输入端;加法运算模块输入端即为加法器的输入端,加法器的另一输入端连接至电压源V2的正极,电压源V2的负极接地,加法器的输出端即为第一双曲正切函数电路单元的输出端。
3.根据权利要求1所述的一种模拟实现Morris-Lecar神经元模型的电路,其特征在于,第二双曲正切函数电路单元中,输入函数运算模块由运算放大器、电压源及电阻组成;
双曲正切函数运算模块由双极性晶体管对、电压源、电流源及电阻组成;
双端输入转单端输出模块由运算放大器及电阻组成
加法运算模块由加法器、电压源及电阻组成;
第二双曲正切函数电路单元的输入端即其输入函数运算模块输入端,其输入经电阻R9连接到运算放大器U3的反相输入端,运算放大器U3的输出经电阻R10反馈连接到运算放大器U3的反相输入端,同时,运算放大器U3的输出端连接到双曲正切函数运算模块的输入端,运算放大器U3的同相输入端连接电压源Vth2的正极,电压源Vth2的负极接地;双曲正切函数运算模块由一个输入端和两个输出端,其输入端即是双极性晶体管Q3的基极,双极型晶体管Q3和Q4的集电极分别经电阻R3和R4连接至电压源V3的正极,电压源V3的负极接地,双极型晶体管Q3和Q4的发射极连接在一起并经电流源I2接地,双极性晶体管Q4的基极接地,双曲正切函数运算模块的两个输出端分别经Q3和Q4的集电极引出;双端输入转单端输出模块的输入端分别连接双曲正切函数运算模块的两个输出端,两个输入信号分别经电阻R15和R13连接到运算放大器U4同相输入端和反相输入端,运算放大器U4的同相输入端同时经电阻R16接地,运算放大器U4的输出端经电阻R14反馈连接至运算放大器U4的反相输入端,双端输入转单端输出模块的输出端连接至加法运算模块的输入端;加法运算模块输入端即为加法器的输入端,加法器的另一输入端连接至电压源V4的正极,电压源V4的负极接地,加法器的输出端即为第二双曲正切函数电路单元的输出端。
4.根据权利要求1所述的一种模拟实现Morris-Lecar神经元模型的电路,其特征在于,双曲余弦函数电路单元中,输入函数运算模块由运算放大器、电压源及电阻组成;
第一指数函数电路单元由运算放大器、电压源、双极性晶体管及电阻组成;
第二指数函数电路单元由电压反相电路、运算放大器、电压源、双极性晶体管及电阻组成;
加法运算模块由加法器实现;
双曲余弦函数电路单元的输入端即其输入函数运算模块输入端,其输入经电阻R22连接到运算放大器U7的反相输入端,运算放大器U7的输出经电阻R23反馈连接到运算放大器U7的反相输入端,同时,运算放大器U7的输出端连接到双曲正切函数运算模块的输入端,运算放大器U7的同相输入端连接电压源Vch的正极,电压源Vch的负极接地;第一指数函数电路单元的输入端连接输入函数运算模块的输出端,第一指数函数电路单元的输出端连接加法运算模块的一个输入端;第二指数函数电路单元的输入端连接输入函数运算模块的输出端,第二指数函数电路单元的输出端连接加法运算模块的另一个输入端;加法运算模块的两个输入端分别连接第一指数函数电路单元的输出端和第二指数函数电路单元的输出端,加法运算模块的输出端连接第二比例积分电路单元的输入端B22。
5.根据权利要求4所述的一种模拟实现Morris-Lecar神经元模型的电路,其特征在于,双曲余弦函数电路单元中实现双曲余弦函数运算的电路是由两个指数函数电路组成,两个指数函数电路的输入量幅值相等符号相反。
6.根据权利要求4所述的一种模拟实现Morris-Lecar神经元模型的电路,其特征在于,第一指数函数电路单元和第二指数函数电路单元能够实现指数运算且输入值的符号是正值或者是负值。
7.根据权利要求4所述的一种模拟实现Morris-Lecar神经元模型的电路,其特征在于,第一指数函数电路单元的输入端即双极性晶体管Q5的基极,双极性晶体管Q5的集电极连接至运算放大器U8的反相输入端,同时运算放大器U8的反相输入端经电阻R24连接至电压源V5的正极,电压源V5的负极接地,运算放大器U8的同相输入端经电阻R25接地,运算放大器U8的输出端连接至双极性晶体管Q5的发射极,同时双极性晶体管Q5的发射极连接至双极性晶体管Q6的发射极,双极性晶体管Q6的集电极连接至运算放大器U9的反相输入端,运算放大器U9的同相输入端经电阻R28接地,运算放大器U9的输出端经电阻R27连接至运算放大器U9的反相输入端,运算放大器U9的输出端即为第一指数函数电路单元的输出端。
8.根据权利要求4所述的一种模拟实现Morris-Lecar神经元模型的电路,其特征在于,第二指数函数电路单元的输入端即为电压反相电路的输入端,第二指数函数电路单元的输入端连接输入函数运算模块的输出端,电压反相电路的输入信号经电阻R29连接至运算放大器U10的反相输入端,运算放大器U10的同相输出端经电阻R31接地,运算放大器U10的输出端经电阻R30连接至算放大器U10的反相输入端,运算放大器的输出端即为电压反相器的输出端,电压反相器的输出端连接至双极性晶体管Q7的基极,双极性晶体管Q7的集电极连接至运算放大器U11的反相输入端,同时运算放大器U11的反相输入端经电阻R32连接至电压源V6的正极,电压源V6的负极接地,运算放大器U11的同相输入端经电阻R33接地,运算放大器U11的输出端连接至双极性晶体管Q7的发射极,同时双极性晶体管Q7的发射极连接至双极性晶体管Q8的发射极,双极性晶体管Q8的集电极连接至运算放大器U12的反相输入端,运算放大器U12的同相输入端经电阻R36接地,运算放大器U12的输出端经电阻R37连接至运算放大器U12的反相输入端,运算放大器U12的输出端即为第二指数函数电路单元的输出端。
9.根据权利要求1所述的一种模拟实现Morris-Lecar神经元模型的电路,其特征在于,第一比例积分电路单元的输入端D1A连接第一双曲正切函数电路单元的输出端,其输入端D2A连接第二比例积分电路单元的输出端,其输入端D3连接第一比例积分电路单元的输出端,其输入端D1B连接第一比例积分电路单元的输出端,其输入端D2B连接第一比例积分电路单元的输出端;输入端D1B的输入信号连接至电压源V7的正极,电压源V7的负极连接乘法器的一个输入端,乘法器的另一个输入端连接输入端D1A的输入信号,乘法器的输出端经电阻R38连接至运算放大器U13的反相输入端;输入端D2B的输入信号连接至电压源V8的负极,电压源V8的正极连接乘法器的一个输入端,乘法器的另一个输入端连接输入端D2A的输入信号,乘法器的输出端经电阻R39连接至运算放大器U13的反相输入端;输入端D3的输入信号连接至电压源V9负极,电压源V9的正极经电阻R37连接至运算放大器U13的反相输入端;输入端D4连接至电流源I3的负极,电流源I3的正极接地;运算放大器U13的输出端经电容C2连接至运算放大器U13的反相输入端。
CN201610066207.8A 2016-01-29 2016-01-29 一种模拟实现Morris‑Lecar神经元模型的电路 Expired - Fee Related CN105787291B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610066207.8A CN105787291B (zh) 2016-01-29 2016-01-29 一种模拟实现Morris‑Lecar神经元模型的电路

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610066207.8A CN105787291B (zh) 2016-01-29 2016-01-29 一种模拟实现Morris‑Lecar神经元模型的电路

Publications (2)

Publication Number Publication Date
CN105787291A CN105787291A (zh) 2016-07-20
CN105787291B true CN105787291B (zh) 2018-04-17

Family

ID=56402614

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610066207.8A Expired - Fee Related CN105787291B (zh) 2016-01-29 2016-01-29 一种模拟实现Morris‑Lecar神经元模型的电路

Country Status (1)

Country Link
CN (1) CN105787291B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3064384B1 (fr) * 2017-03-23 2019-05-03 Stmicroelectronics Sa Bloc refractaire pour dispositif integre de neurone artificiel
CN109447255B (zh) * 2018-11-29 2021-05-18 西北工业大学 一种模拟实现简化Hodgkin-Huxley神经元模型的电路
CN109934338B (zh) * 2019-03-03 2024-03-19 广西师范大学 一种实现神经元模型的硬件电路
CN115062772B (zh) * 2022-06-10 2023-09-26 常州大学 一种简单的rc型神经元簇放电电路
CN114861903B (zh) * 2022-06-15 2023-05-26 兰州交通大学 一种时滞耦合神经元模型的硬件电路

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1194469C (zh) * 2001-12-28 2005-03-23 北京六合万通微电子技术有限公司 低功耗的模拟信号采样保持电路
CN1428934A (zh) * 2001-12-28 2003-07-09 北京六合万通微电子技术有限公司 以神经元电路为基本单元的电压比较器
US20070022070A1 (en) * 2005-03-15 2007-01-25 Wells Richard B Forgetful logic for artificial neural networks
US8996431B2 (en) * 2012-07-25 2015-03-31 Hrl Laboratories, Llc Spike domain neuron circuit with programmable kinetic dynamic, homeostatic plasticity and axonal delays
CN103716014B (zh) * 2013-12-04 2016-04-06 浙江大学城市学院 一种基于神经元mos管的差分型双边沿触发器设计
CN104689473B (zh) * 2015-02-02 2016-12-07 天津大学 基于fpga的电刺激下神经元随机响应及共振实验平台

Also Published As

Publication number Publication date
CN105787291A (zh) 2016-07-20

Similar Documents

Publication Publication Date Title
CN105787291B (zh) 一种模拟实现Morris‑Lecar神经元模型的电路
CN106815636B (zh) 一种基于忆阻器的神经元电路
CN103559328B (zh) 一种忆容器的实现电路及其实现方法
CN110097182B (zh) 用神经激活梯度λ控制的三维Hopfield神经网络模型实现电路
Li et al. On the number of limit cycles of a perturbed cubic polynomial differential center
CN101295454B (zh) 一种无电感蔡氏电路
Herceg et al. Arduino and numerical mathematics
CN206147706U (zh) 神经元电路
CN105207771A (zh) 一种基于单片电流反馈运算放大器的双涡卷隐藏蔡氏混沌吸引子产生电路
Siemon Launching mathematical futures: The key role of multiplicative thinking
Sivakumar et al. Shooting type Laplace–Adomian decomposition algorithm for nonlinear differential equations with boundary conditions at infinity
Steger et al. Teaching battery basics in laboratories: Comparing learning outcomes of hands-on experiments and computer-based simulations
Puchta On the role of mathematics and mathematical knowledge in the invention of Vannevar Bush's early analog computers
Suits Electronics for Physicists: An Introduction
CN206282337U (zh) 脉冲神经电路
Liu et al. Reachability and controllability of linear switched impulsive systems
Kudryashov On solutions of generalized modified Korteweg–de Vries equation of the fifth order with dissipation
Roy et al. 14 Mathematical modelling of semiconductor devices and circuits
CN205407839U (zh) 一种网格多翅膀超混沌信号发生器
Tunç et al. Boundedness of solutions to differential equations of fourth order with oscillatory restoring and forcing terms
CN109447255A (zh) 一种模拟实现简化Hodgkin-Huxley神经元模型的电路
Liu et al. An Intelligent Recognition Method of Circuit Connection in Three-dimensional Virtual Experiments
Hagigat Using MUTISIM software to reinforce use and application of Norton’s theory in electrical circuits
Hagigat Using MUTISIM software to reinforce use and application of Thevenin’s theory in electrical circuits
Sahni et al. Mathematical Modeling, Computational Intelligence Techniques and Renewable Energy: Proceedings of the Third International Conference, MMCITRE 2022

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180417

Termination date: 20220129

CF01 Termination of patent right due to non-payment of annual fee