Memristor+Matlab 仿真（二）

忆阻器的定义

Any 2-terminal device exhibiting a pinched hysteresis loop which always passes through the origin in the voltage-current plane when driven by any periodic input current source, or voltage source, with zero DC component is called a memristor. - L.O.Chua

myddefun.m

function dy = myddefun(t,y)

vm=sin(0.3*t);

dy=[
   1.8*y(1)-3.9*y(1)^3+1.4*vm-1.5*y(1)*vm
       ];

wenjie.m

clear;clc
RelTol=1e-5;		
AbsTol=1e-6;		
options=odeset('RelTol',RelTol,'AbsTol',AbsTol);

x0=0.1;

t0=0:0.01:100;
[t,x]=ode45('myddefun',t0,x0,options);

n=1:10000;
vm=sin(0.3*t);
I1 = vm .*(x.^2-x-1); % current  

figure
subplot(311)
plot(n, x(1:10000,1),'k');
%plot(sol.x,sol.y)

subplot(312)
plot(vm, I1,'k');

ω=0.3时

Q:1.为什么ω>0.3时时是混乱的图像
A:1.混乱是因为忘记保存了两边Vm不一样

Lissajous图形
定义：两个沿着互相垂直方向的正弦振动的合成的轨迹。
作用：测量出两个信号的频率比与相位差。