matlab蔡氏电路

function chua_memristor
clear all
figure（1）
%This program originally computed by S. Sabarathinam 
%based on ‘Memristor based Chaotic Circuits-M. Bharathwaj‘
%don‘t forgot to rate my files...
%To vary the value of ‘r‘ you can get the bifurcation scenario.
       display（‘S. SABARATHINAM>>>please rate my files..‘）;
            r=2000.0; c1=6.8e-9;	c2=68.0e-9;	bind=18.0e-3;
            alpha = -0.667D-03;beta=0.029D-03; xx=bind*c2;tau=3.499e-05;
          cn1=tau/c1*r;cn2=tau/c2*r; cn3=tau/bind;
 %                     
            tspan = 0:0.01:500;  
            x10 = 0.1; x20 = -0.2; x30 = -0.003;
            x40 = 0.05; 

            y0 = [x10; x20; x30 ; x40 ];
           op=odeset（‘absto