josephson.asv

本程序基于约瑟夫森结的RSJ模型

clear all;
global a;%全局变量的设置至关重要，是整个程序的关键
a=1;
BB=[];
for a=[1:0.5:9]
    % 选用二三阶的龙格库塔函数，其中0到5的步进为x的取值范围，pi/6为x取值为0时y的值
    [x,y]=ode23(@fangcheng,[0:0.001:5],[pi/6]);
    b=find(abs(y-pi/3)<1e-3);%寻找y=pi/3的元素标号
    I1=a*5-5*sin(x(b)/4.8*0.144*pi);%根据相关公式计算电流值
    [x,y]=ode23(@fangcheng,[10:0.001:15],[pi/6]);%二次取值
    c=find(abs(y-pi/3)<1e-3);
    I2=a*5-5*sin(x(c)/4.8*0.144*pi);
    s=length(I1);s1=length(I2);
    A=zeros(1,s+s1);
    A(1:s)=I1';
    A((s+1):(s+s1))=I2';
    B=zeros(2,s+s1);%创建数组
    B(1,:)=A;
    V=(a-0.866)*10;
    B(2,:)=V';
    for r=[1:s+s1]%筛选出去一部分数值
        if B(1,r)<5
            B(2,r)=0;
        end
    end
    BC=[];
    for i=1:length(B)%将数组中的列重排
        BC(:,i)=B(:,(length(B)-i+1));
    end
    BB=[BB BC];%将所有结果存入大数组
end
figure(1);
stairs(BB(2,:),BB(1,:));%阶梯画图
title('微波辐照下约瑟夫森结的I-V特性');
xlabel('\fontsize{10}\bfV(mV)\rightarrow')
ylabel('\fontsize{10}\bfI(mA)\rightarrow')