josephson_junction.m

来自「基于Matlab对约瑟夫森结（Josephson Junction）RCSJ模型」· M 代码 · 共 19 行

%
function dy=Josephson_Junction(t,y)
% 在程序中参数t代表文中的τ
% Bc 为阻尼参数, w 为归一化的外加交流电流圆频率
% Bc、ia、id、w、τ均为无量纲化时引入的变量，其中ia、id为各电流无量纲化时引入的
% y(1)= φ  约瑟夫森结中心处波函数ψ1、ψ2的相位差
% y(2)=dφ/dt 即归一化电压u
% 以下为约瑟夫森结RCSJ模型的无量纲化微分方程
% dφ/dt=u
% du/dt=-Bc*u-sin(φ)+id+ia*sin(w*t)
%
% ****** Author:J.F.Peng([e-mail]) ******
%
global Bc ia w id;
dy=zeros(3,1);
dy(1)=y(2);
dy(2)=id+ia*sin(w*t)-sin(y(1))-Bc*y(2);
dy(3)=w;