plotfork.m

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

%
% 分岔图(Fork)
% 
% ****** Author:J.F.Peng([e-mail]) ******
%
clear
clc
%
% tic、toc用来计算程序运行时间，分别表示开始和结束计时
%
tic
global Bc ia w id; 
id=0.00; 
w=0.20; 
ia=1.50;
for Bc=0.00:0.002:2
    [T,Y]=ode45('Josephson_Junction',[0,200],[0;0;0]); 
    data=Y(:,2); 
    n=length(data); 
    % 去掉一些前面的点，防止因计算机配置问题而画不出图，不影响结果
    m=round(n/6); 
    a=data(m-2); 
    b=data(m-1); 
    for i=m:n 
        if b>=a&b>=data(i); 
           plot(Bc,b,'k'); 
           hold on; 
        end 
        a=b; 
        b=data(i); 
    end 
end 
%
% 根据需要选择图像标注语句
%
% ia-u 分岔图（fork）
%
%xlabel('Aternative Current ia')
%ylabel('Voltage u=dφ/dt')
%desc={'Bc=0.50,w=0.66,id=0.00'};
%text(0.5,3,desc)
%
% id-u 分岔图 (fork)
%
%xlabel('Directive Current id')
%ylabel('Voltage u=dφ/dt')
%desc={'Bc=1.00,w=0.50,ia=0.65'};
%text(1,3,desc)
%
% Bc-u 分岔图 (fork)
%
xlabel('Damping Coefficient \betac')
ylabel('Voltage u=dφ/dt')
desc={'ia=1.50,w=0.20,id=0.00'};
text(0.75,8,desc)
toc
t=toc