rsj.m

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

clear all;
global a;
B=[];
for a=[1:3]
    [x,y]=ode45(@fangcheng,[0:0.001:5],[pi/6]);
    b=find(abs(y-pi/3)<1e-3);%寻找特定y值
    %%因为是单值，所以没有进行二次取值
    V=(a-sin(y(b)))'*10;
    I=zeros(1,length(V));
    I=5*a;
    A=zeros(2,length(V));
    A(1,:)=I;
    A(2,:)=V;
    B=[B A]; %将电压电流值存贮画图
    
end
figure(1);
plot(B(2,:),B(1,:));
hold on
I=[0:0.5:5];
V=zeros(1,11);
figure(1);
plot(V,I);
hold on;%为了在一个图上画所采用的命令
V=[0 1.3351];I=[5,5];
figure(1);
plot(V,I);
hold off;
title('RSJ直流模型图');
xlabel('\fontsize{10}\bfV(mV)\rightarrow')
ylabel('\fontsize{10}\bfI(mA)\rightarrow')