-
1971年,文献[1]从对称性的角度在理论上预言了忆阻器(memristor)的存在。1976年,文献[2]将忆阻器的概念扩展至忆阻系统(memristive systems)。2008年,文献[3]首次发现了纳米忆阻器件。由此证明了忆阻器的存在,并掀起了忆阻器研究的热潮。忆阻器在非易失性存储器[4]、逻辑门运算[5]、人工神经网络[6]、生物学[7]、忆阻混沌电路及动力学研究[8-10]等多个领域具有丰富的应用前景。
忆阻器描述的是电荷q与磁通量
$\varphi $ 之间的关系[1],若磁通量$\varphi $ 是电荷q的单值函数,则忆阻M(q)是电荷q的单值函数,称之为荷控忆阻器。同样地,若电荷q是磁通量$\varphi $ 的单值函数,则忆导W($\varphi $ )是磁通量$\varphi $ 的单值函数,则称之为磁控忆阻器(或忆导器)。根据荷控与磁控忆阻器的定义,两者的特征分别由
$\varphi $ −q曲线或q−$\varphi $ 曲线决定。基于q−$\varphi $ 三次非线性关系,文献中[9]提出了一种三阶磁控忆阻器模型。基于$\varphi $ −q二次非线性关系,文献[11]提出了一种荷控忆阻器模型。基于q与$\varphi $ 之间分段线性关系,文献[12]提出了荷控和磁控忆阻器模型。滤波器是一种对不同频率范围内的信号具有不同抑制作用的频率选择电路。由电阻和电容组成的无源滤波是最常见的滤波器,而基于忆阻器的滤波器却比较少见。文献[13]使用HP忆阻器模型搭建了高低通滤波器,并进行了仿真分析;文献[14]基于三阶磁控忆阻器搭建了高低通滤波器,并进行了仿真和实际电路分析;文献[15]基于三阶磁控忆阻器搭建了有源低通滤波器。
本文基于q−
$\varphi $ 之间三次多项式关系,提出了一种新型的三阶磁控忆阻器模型。与文献[9-10,14-15]中使用的三阶磁控忆阻器模型不同,本文描述q−$\varphi $ 关系的三次多项式中二次项${\varphi ^2}$ 的系数不为0,模型更具有一般性。同时基于新型三阶磁控忆阻器模型设计了一阶低、高通滤波器。仿真与实际电路分析结果表明:本文设计的新型三阶磁控忆阻器模型具备忆阻器的3个基本特征[16],即忆阻器两端电压与流经的电流呈"8"字形滞回曲线,随着频率的升高,滞回曲线所包围的面积减小,当频率趋近于无穷时,滞回曲线退化成单值函数。 -
如前所述,磁控忆阻器中电荷q是磁通量
$\varphi $ 的单值函数,本文将函数q($\varphi $ )进行泰勒级数展开,其表达式为:$$ q(\varphi ) = \sum\limits_{n = 0}^\infty {{a_n}{\varphi ^n}} $$ (1) 式中,an为
${\varphi ^n}$ 的系数。为了简化问题,本文只将该多项式展开到三阶,有:$$ q = \alpha \varphi + \beta {\varphi ^2} + \gamma {\varphi ^3} $$ (2) 式中,
$\alpha $ 、$\,\beta $ 和$\gamma $ 分别为各次项的系数,其值均为常数。与文献[9-10,14-15]使用的三阶磁控忆阻器模型相比,本文提出的磁控忆阻器模型在描述q−$\varphi $ 关系时多了$\,\beta {\varphi ^2}$ 分量,使之更具有一般性。当$\,\beta $ =0时,此时忆阻器模型与文献[9-10,14-15]提出的忆阻器模型一致。结合
$W(\varphi ) = {\rm{d}}q(\varphi )/{\rm{d}}(\varphi)$ ,可以得出相应的忆导W($\varphi $ )为:$$ W(\varphi ) = \alpha + 2\beta \varphi + 3\gamma {\varphi ^2} $$ (3) 新型三阶磁控忆阻器与无
$\varphi $ 的二次分量的三阶磁控忆阻器的q−$\varphi $ 曲线和W−$\varphi $ 曲线如图1所示,其中$\alpha $ =1 S,$\beta $ =1 S/Wb,$\gamma $ =1 S/Wb2。由电压v、磁通量
$\varphi $ 和时间t三者之间的关系式$v(t) = {\rm{d}}\varphi (t)/{\rm{d}}t$ 以及电压v、电流i以及电导W三者之间的关系式$i(t) = W(\varphi )v(t)$ 可以得出:$$ \left\{ {\begin{aligned} & {i(t) = (\alpha + 2\beta \varphi + 3\gamma {\varphi ^2})v(t)} \\ & {\varphi (t) = {\varphi _0} + \int_0^t {v(\tau ){\rm{d}}\tau } } \end{aligned}} \right. $$ (4) 式中,
${\varphi _0}$ 为磁通量的初始量。可以看出,磁控忆阻器的输出不仅跟此刻的输入有关,还与此前的输入有关,体现了忆阻器的“记忆”特性。另外,为了使磁控忆阻器的伏安特性曲线始终在第一、三象限内[11,16],即忆导值始终大于0,则要求$\alpha $ 、$\,\beta $ 和$\gamma $ 三者之间的关系为:$$ \left\{ {\begin{aligned} & {4{\beta ^2} - 4\alpha 3\gamma < 0} \\ & {3\gamma > 0} \end{aligned}} \right. $$ (5) 为了检验该模型是否正确,本文设置一个电压源输出的电压信号为:
$$ \left\{ {\begin{aligned} & {v(t) = A\cos (\omega t)}&{t \geqslant 0} \\ & {v(t) = 0}&{t < 0} \end{aligned}} \right. $$ (6) 由式(4)与式(6)可得电流i(t)的表达式:
$$ \begin{split} & i(t) = \left(\alpha A + \frac{{3\gamma {A^2}}}{{4{\omega ^2}}}\right)\cos (\omega t) +\\ & \beta \frac{{{A^2}}}{\omega }\sin (2\omega t) - \frac{{3\gamma {A^3}}}{{4{\omega ^2}}}\cos (3\omega t) \end{split} $$ (7) 现给定输入电压源,选取幅值A、频率
$\omega $ ,以及参数$\alpha $ 、$\beta $ 和$\gamma $ 的值如表1所示。根据式(7)在MATLAB软件中仿真可得该三阶磁控忆阻器的二端电压以及对应的电流波形和伏安特性曲线,如图2所示。表 1 三阶磁控忆阻器模型参数值
参数 参数值 A/V 1 $\omega $/rad·s−1 π $\alpha $/S 1 $2\beta $/S·Wb−1 1 $3\gamma $/S·Wb−2 1 从图2b的三阶磁控忆阻器伏安特性曲线可以看出,该滞回曲线在原点相交,位于第一、三象限,且形状呈斜“8”字形,说明该三阶磁控忆阻器模型符合忆阻器的特性。
-
上一节忆阻器的建模过程中对电压的各种运算都是相对于零电位而言的,如果直接根据式(4)搭建忆阻器等效电路,那么该等效电路的一端必然会接地。
为了克服忆阻器等效电路一端接地的局限性,使之成为“浮地”二端口,文献[17-20]提出了多种浮地忆阻器等效电路,然而其中的电路结构较为复杂。为简化电路,本文将相对于零电位的电压转换为二端口之间的“电位差”[14],从而得到如下所示的三阶磁控忆阻器的二端口“浮地”表达式:
$$ \left\{ {\begin{aligned} & {i(t) = (\alpha + 2\beta \varphi + 3\gamma {\varphi ^2})({v_1} - {v_2})} \\ & {\varphi (t) = {\varphi _0} + \int\nolimits_0^t {[{v_1}(\tau ) - {v_2}(\tau )]{\rm d}\tau } } \end{aligned}} \right. $$ (8) 式中,v1和v2分别代表忆阻器等效电路两端的电压。
-
忆阻器等效电路二端口之间的电位差可以利用由运算放大器构成的减法器计算得出。之后便可以将计算得出的电位差进行加法、乘法以及积分运算。具体的电路图实现如图3所示。
该忆阻器等效电路由15个电阻、一个电容、5个AD712运算放大器以及3个AD633乘法器构成,整个等效电路由15 V直流电源驱动。由于本文中所提出的新型三阶磁控忆阻器模型多了
$\,\beta {\varphi ^2}$ 分量,因此该等效电路中使用乘法器U6来计算该分量。运算放大器U1A构成一个电压跟随器,它具有高输入阻抗和低输出阻抗的特点,输入电流极少进入运算放大器,起到与前级电路隔离的作用。其输出电压v3等于输入电压v1。
运算放大器U1B与电阻R1~R4构成一个减法器。由于电阻R1~R4的阻值相等,因此输出电压v4与输入电压v2、v3之间的关系为:
$$ {v_4} = {v_3} - {v_2} $$ (9) 运算放大器U2A与电阻R5~R7以及电容C1构成一个积分器。由于实际运算放大器并不是理想的,因此电阻R7的作用是为了减小误差。电阻R6的作用是为了防止积分漂移所造成的积分饱和或截止现象,同时,为了减小积分误差,电阻R6>>R5。输出电压v5与输入电压v4之间的关系为:
$$ {v_5} = \frac{{ - 1}}{{{R_5}{C_1}}}\int\nolimits_0^t {{v_4}(\tau ){\rm{d}}\tau } = \frac{{ - 1}}{{{R_5}{C_1}}}\varphi (t) $$ (10) U4~U6为AD633乘法器,根据其数据手册知其输入输出之间的关系式为:
$$ W = \frac{{(X1 - X2)(Y1 - Y2)}}{{10V}} + Z $$ (11) 根据式(11)可以直接计算出电压v6与电压v7的表达式:
$$ {v_6} = \frac{{{v_4}{v_5}^2}}{{100}} $$ (12) $$ {v_7} = - \frac{{{v_4}{v_5}}}{{10}} $$ (13) 运算放大器U2B与电阻R13、R14构成一个反相放大器,由于电阻R13=R14,因此输出电压v8与输入电压v2之间的关系为:
$$ {v_8} = - {v_2} $$ (14) 运算放大器U3A与电阻R8~R12构成一个反相加法器。同样地,电阻R12的作用是为了减小计算误差。由于电阻R8~R11的阻值均为10 kΩ,因此输出电压v9与输入电压v6、v7、v8之间的关系为:
$$ {v_9} = - ({v_6} + {v_7} + {v_8}) $$ (15) 最后,根据电阻R15两端的电压即可计算出流经电阻R15的电流。正如前面所提及,几乎没有电流流入运算放大器U1A,因此流经R15的电流即为流经忆阻器的电流。结合式(8)~式(15),即可得出电流i(t)的表达式为:
$$ \begin{split} &\quad\;\; i(t) = \frac{{{v_1} - {v_9}}}{{{R_{15}}}} = \frac{{{v_1} - {v_2}}}{{{R_{15}}}} + \\ & \frac{{{v_1} - {v_2}}}{{10{R_5}{C_1}{R_{15}}}}\varphi (t) + \frac{{{v_1} - {v_2}}}{{100{R_5}^2{C_1}^2{R_{15}}}}{\varphi ^2}(t) \end{split} $$ (16) 对比式(8)与式(16),即可得出重要参数
$\alpha $ 、$\ \beta $ 、$ \gamma $ 的表示式:$$ \left\{ {\begin{split} & {\alpha = \frac{1}{{{R_{15}}}}} \\ & {2\beta = \frac{1}{{10{R_5}{C_1}{R_{15}}}}} \\ & {3\gamma = \frac{1}{{100{R_5}^2{C_1}^2{R_{15}}}}} \end{split}} \right. $$ (17) -
在Cadence(PSpice)软件上对上述三阶磁控忆阻器的伏安特性进行仿真分析。根据式(5)中
$\alpha $ 、$\,\beta $ 和$\gamma $ 三者之间的关系,并经过一些调优,电阻电容的设置如下:C1=200 nF、R5=10 kΩ、R15=10 kΩ。可以计算出$\alpha $ =0.1 mS、$2\beta $ =5 mS/Wb、$3\gamma $ =0.25 S/Wb2。为了测试三阶磁控忆阻器等效电路的伏安特性,使用如图4所示的测试电路。正弦电源V3的两端接入忆阻器等效电路的MEM_PLUS端和MEM_MINUS端。为了体现该忆阻器等效电路的浮地特性,电压源两端均未接地。同时,在正弦电压源端标记一个电流探针,用以测量经过忆阻器的电流大小,将横坐标改为电压源的幅值便可画出该忆阻器等效电路的伏安特性曲线。
将正弦电压源的幅值固定为5 V,分别设置频率为100、200 Hz、10 kHz,其伏安特性曲线如图5所示。从图5中可以看出,该忆阻器模型符合忆阻器的3个特征[16]:当忆阻器两端施加周期信号时,其伏安特性曲线为滞回曲线,且经过原点;随着频率增加,其滞回曲线面积减小;当频率趋近于无穷时,其滞回曲线退化成为单值函数,忆阻器退化成为线性电阻。
将电压源的频率固定为100 Hz,分别设置振幅为3、5、7 V,其伏安特性曲线如图6所示。从图6中可以看出:随着忆阻器等效电路两端电压增加,其滞回曲线所围成的面积也会随之增加。
-
基于磁控忆阻器(忆导)的WC低通滤波器电路图如图7所示。
根据基尔霍夫定律,流经忆导器W与电容C的电流相等,联立式(8)可以列出:
$$ \left\{ {\begin{aligned} & {C\frac{{{\rm{d}}{v_{\rm{out}}}}}{{{\rm{d}}t}} = W(\varphi )({v_{{\rm{in}}}} - {v_{{\rm{out}}}})} \\ & {W(\varphi ) = \alpha + 2\beta \varphi + 3\gamma {\varphi ^2}} \\ & {\varphi (t) = {\varphi _0} + \int\nolimits_0^t {[{v_{{\rm{in}}}}(\tau ) - {v_{{\rm{out}}}}(\tau )]{\rm d}\tau } } \end{aligned}} \right. $$ (18) 由于基于磁控忆阻器的低通滤波器为时变系统[14],不能使用一般的系统传输函数进行分析,因此从定性的角度对该系统进行分析。根据式(7)可以得出,三阶磁控忆阻器在低频信号的驱动下,会产生明显的二次和三次谐波分量。由电容的传输函数
${H_C}({\rm{j}}\omega ) = 1/{\rm{j}}\omega C$ 可知,频率越高,电容C的阻抗越小。因此,WC低通滤波器相较于RC低通滤波器,在其通带内(即低频信号)的输出增益较低。基于图3所示的三阶磁控忆阻器(忆导)W,搭建出低通WC滤波器,并在PSpice软件中进行仿真,与此同时,搭建出一个RC低通滤波器进行比较,如图8所示。
图9为WC低通滤波器与RC低通滤波器的仿真结果。从图中可以看出,WC低通滤波器输出增益较RC低通滤波器低,且相位变化也与RC低通滤波器有所不同。表明了磁控忆阻器由于对低频信号产生高次谐波分量,其具有非线性特性。
-
基于磁控忆阻器(忆导)的9高通滤波器电路图如图10所示。
同样根据基尔霍夫定律,流经忆导器W与电容C的电流相等,联立式(8)可以列出:
$$ \left\{ {\begin{split} & {C\frac{{{\rm{d}}({v_{{\rm{in}}}} - {v_{{\rm{out}}}})}}{{{\rm{d}}t}} = W(\varphi ){v_{{\rm{out}}}}} \\ & {W(\varphi ) = \alpha + 2\beta \varphi + 3\gamma {\varphi ^2}} \\ & {\varphi (t) = {\varphi _0} + \int\nolimits_0^t {{v_{{\rm{out}}}}(\tau ){\rm{d}}\tau } } \end{split}} \right. $$ (19) 根据忆阻器的基本特征:随着频率增加,其滞回曲线面积减小;当频率趋近于无穷时,其滞回曲线退化成为单值函数,忆阻器退化成为线性电阻。因此WC高通滤波器与RC高通滤波器在高频带输出增益几乎一致。
基于图3所示的三阶磁控忆阻器(忆导)W,搭建WC高通滤波器,并在PSpice软件中进行仿真,同时,搭建出一个RC高通滤波器进行比较,如图11所示。
图12为WC高通滤波器与RC高通滤波器的仿真结果。从图中可以看出,WC高通滤波器输出增益与RC高通滤波器输出增益几乎无差异,且两者的输出相位也类似。这表明忆阻器在高频信号的激励下,退化成为线性电阻。
A Novel Cubic Flux-Controlled Memristor Model and Its Filter Study
-
摘要: 该文基于电荷量q与磁通量φ的三次多项式关系,提出一种新型的三阶磁控忆阻器模型,并利用MATLAB进行了仿真验证。利用“浮地”二端口等效电路原理,设计实现三阶磁控忆阻器及其低通与高通滤波器的等效电路,并在PSPICE上进行了仿真分析。最后,搭建出三阶磁控忆阻器及其低通与高通滤波器的实际电路,并利用示波器进行了测试。实验与仿真结果表明:基于该新型三阶磁控忆阻器的低通滤波器较基于电阻的输出增益低;而基于该新型三阶磁控忆阻器的高通滤波器与基于电阻的输出增益基本一致。
-
关键词:
- 三阶磁控忆阻器 /
- 高通滤波器 /
- 低通滤波器 /
- MATLAB与PSPICE仿真
Abstract: In this paper, based on the cubic polynomial relation between charge q and magnetic flux$ \varphi$ , a novel cubic flux-controlled memristor model is proposed and verified by MATLAB simulation. Applying the principle of "floating ground" two port equivalent circuit, the equivalent circuit of the cubic flux-controlled memristor and its low-pass and high-pass filters are designed and implemented, and the simulation is carried out on PSPICE. Finally, the actual circuits of the cubic flux-controlled memristor and its low-pass and high-pass filters are built and tested by oscilloscope. The experimental and simulation results show that the low-pass filter based on the novel flux-controlled memristor has a lower output gain than the resistor-based one, while high-pass filter based on the novel flux-controlled memristor is basically consistent with the output gain based on resistor. -
表 1 三阶磁控忆阻器模型参数值
参数 参数值 A/V 1 $\omega $ /rad·s−1π $\alpha $ /S1 $2\beta $ /S·Wb−11 $3\gamma $ /S·Wb−21 -
[1] CHUA L O. Memristor—The missing circuit element[J]. IEEE Trans Circuit Theory, 1971, 18(5): 507-519. doi: 10.1109/TCT.1971.1083337 [2] CHUA L O, KANG S M. Memristive devices and systems[J]. Proceedings of the IEEE, 1976, 64(2): 209-223. doi: 10.1109/PROC.1976.10092 [3] STRUKOV D B, SNIDER G S, STEWART D R, et al. The missing memristor found[J]. Nature, 2008, 453: 80-83. doi: 10.1038/nature06932 [4] HO Y, HUANG G M, LI P. Dynamical properties and design analysis for nonvolatile memristor memories[J]. IEEE Transactions on Circuits and Systems I, 2011, 58(4): 724-736. doi: 10.1109/TCSI.2010.2078710 [5] PREZIOSO M, RIMINUCCI A, GRAZIOSI P, et al. A single-device universal logic gate based on a magnetically enhanced memristor[J]. Advanced Materials, 2013, 25(4): 534-538. doi: 10.1002/adma.201202031 [6] PERSHIN Y V, VENTRA M D. Experimental demonstration of associative memory with memristive neural networks[J]. Neural Networks the Official Journal of the International Neural Network Society, 2010, 23(7): 881-886. doi: 10.1016/j.neunet.2010.05.001 [7] PERSHIN Y V, la FONTAINE S, di VENTRA M. Memristive model of amoeba learning[J]. Physical Review E, 2009, DOI: 10.1103/PhysRevE.82.019904. [8] MUTHUSWAMY B, CHUA L O. Simplest chaotic circuit[J]. International Journal of Bifurcation & Chaos, 2010, 20(5): 1567-1580. [9] ZHONG G Q. Implementation of Chua’s circuit with a cubic nonlinearity[J]. IEEE Transactions on Circuits & Systems I Fundamental Theory & Applications, 1994, 41(12): 934-941. [10] 包伯成, 胡文, 许建平, 等. 忆阻混沌电路的分析与实现[J]. 物理学报, 2011, 60(12): 63-70. BAO Bo-cheng, HU Wen, XU Jian-ping, et al. Analysis and implementation of memristor chaotic circuit[J]. Journal of Physics, 2011, 60(12): 63-70. [11] DANG B, YU Y, MA H, et al. A novel charge-controlled memristor model and its emulator circuit[C]//The 18th International Conference on Communication Technology (ICCT). [S.l.]: IEEE, 2018: 586-590. [12] ITOH M, CHUA L O. Memristor oscillators[J]. International Journal of Bifurcation & Chaos, 2008, 18(11): 3183-3206. [13] 田晓波. 忆阻器电路特性与应用研究[D]. 合肥: 国防科学技术大学, 2009. TIAN Xiao-bo. Research on circuit characteristics and application of memristor[D]. Hefei: National University of Defense Technology, 2009. [14] 杨凌, 蒲中柱, 黄麟, 等. 有源磁控忆阻器的等效电路及滤波频率特性研究[J]. 电子元件与材料, 2015, 34(7): 5-10. YANG Ling, PU Zhong-zhu, HUANG Lin, et al. Active flux-controlled memristor emulator and the frequency response analysis of memristor based filter[J]. Electronic Components and Materials, 2015, 34(7): 5-10. [15] 欧青立, 赵平荣, 游淼. 基于三次非线性忆阻器的有源低通滤波器[J]. 应用科技, 2014, 41(5): 6-10. OU Qing-li, ZHAO Ping-rong, YOU Miao. Active low-pass filter characteristics based on the cubic nonlinear memristor[J]. Applied Science and Technology, 2014, 41(5): 6-10. [16] ADHIKARI S P, SAH M P, KIM H, et al. Three fingerprints of memristor[J]. IEEE Transactions on Circuits & Systems I Regular Papers, 2013, 60(11): 3008-3021. [17] LIU W, WANG F Q, MA X K. Exponential flux-controlled memristor model and its floating emulator[J]. Chinese Physics B, 2015, 24(11): 581-587. [18] SANCHEZ-LOPEZ C, MENDOZA-LOPEZ J, CARRASCO-AGUILAR M A, et al. A simple floating memristor emulator circuit based on current conveyors[C]//International Conference on Electrical Engineering. [S.l.]: IEEE, 2013: DOI: 10.1109/ICEEE.2013.6676031. [19] SANCHEZ-LOPEZ C, MENDOZA-LOPEZ J, CARRASCO-AGUILAR M A, et al. A floating analog memristor emulator circuit[J]. IEEE Transactions on Circuits and Systems Ⅱ, 2014, 61(5): 1. doi: 10.1109/TCSII.2014.2321638 [20] KIM H, SAH M P, YANG C, et al. Memristor emulator for memristor circuit applications[J]. IEEE Transactions on Circuits & Systems I, 2012, 59(10): 2422-2431.