Memristor+Matlab 仿真（三）

CHAOTIC OSCILLATION IN A SIMPLEST NONLINEAR CIRCUIT WITH LOCALLY-ACTIVE MEMRISTOR
关键词：忆阻器，简易非线性电路，混沌震荡


将微分方程组输入Matlab wenjie.m

function dy=wenjie(~,y)

c = 0.13;   %电容130mF
L = 0.05;   %电感20mH

%i1 = (a + b*y(4)^2)*y(1); %我没有用到i1,所以不需要

dy=[ -1/c*((y(3)^2-y(3)-1)*y(1)+y(2)); 
     1/L*y(1);
     1.8*y(3)-3.9*y(3)^3+1.4*y(1)-1.5*y(3)*y(1);
     y(1);
   ];

电路时域波形 wenjie_run.m

clc;
clear;
close all;

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

x0=[0.1;0.1;0.6794];

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

n=1:1000;
figure
subplot(311)
pl