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

				%
				% 分岔图(Fork)
				% 
				% ****** Author:J.F.Peng(jingyujiafu@163.com) ******
				%
				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