lyapunov.m

基于Matlab对约瑟夫森结（Josephson Junction）RCSJ模型的交直流I-V特性及非线性混沌现象进行数值模拟。通过计算机数值模拟得到该模型的非线性微分方程数值解

% the largest lyapunov exponents of josephson
%
% Author's email: [e-mail]
%
clear
clc
global Bc ia w id; 
%
% The Largest Lyapunov Index ia
%
Bc=0.50; 
w=0.66; 
id=0.00;
Z=[];
for ia=linspace(0,0.01,80)
    Z=[Z ia+i*jose_ly(ia,300)];
end
figure
plot(Z,'-') 
title('The Largest Lyapunov Index'),xlabel('ia'),ylabel('lyapunov')
grid on
grid minor
%
% The Largest Lyapunov Index Bc
%
ia=1.212; 
w=0.66; 
id=0.00;
Z=[];
for Bc=linspace(0,0.01,80)
    Z=[Z Bc+i*jose_ly(Bc,300)];
end
figure
plot(Z,'-') 
title('The Largest Lyapunov Index'),xlabel('Bc'),ylabel('lyapunov')
grid on
grid minor
%
% The Largest Lyapunov Index id
%
Bc=0.50; 
w=0.66; 
ia=1.212;
Z=[];
for id=linspace(0,0.01,80)
    Z=[Z id+i*jose_ly(id,300)];
end
figure
plot(Z,'-') 
title('The Largest Lyapunov Index'),xlabel('id'),ylabel('lyapunov')
grid on
grid minor