he11_effn.m

本软件使用matlab预言编写

function [ N ] = HE11_EffN( Lamda,KJ,d,f )
%UNTITLED1 Summary of this function goes here
%  Detailed explanation goes here
u0=0.001;
%u1=10;
%u2=(u1-u0)/50;
%u_tmp=u0:u2:u1;
%[KJ,d]=LoadH();
%R=KJ.*sqrt(sqrt(3)./(2.*pi));   %外径
%f=0.5;
%r=f.*KJ;                          %内径
n=length(Lamda);
for q=1:n
    n2=CladN(Lamda(q),KJ,d);
    n1=SiN(Lamda(q));
    %n1=n1;
    %n2=1.447;
    %n1=1.45;
    k0=2.*pi./Lamda(q);
    v=k0.*sqrt(n1.^2-n2.^2);
    u2=(v-u0)./100000;
    u1=v-u2;
    u_tmp=u0:u2:u1;
    
    
    X=HE11_tmp(Lamda(q),u_tmp,KJ,d,f,n1,n2);
    m=length(X);
    %plot(u_tmp,X)
    for p=1:m-1
        if((X(p)<0)&(X(p+1)>0))
            u=(u_tmp(p)+u_tmp(p+1))/2;
            ux1=u_tmp(p);
            ux2=u_tmp(p+1);
            break;
        end
    end
    while(abs(ux1-ux2)>0.00001)
        if(HE11_tmp(Lamda(q),ux1,KJ,d,f,n1,n2).*HE11_tmp(Lamda(q),(ux1+ux2)./2,KJ,d,f,n1,n2)<=0)
            ux2=(ux1+ux2)./2;
        else
            ux1=(ux1+ux2)./2;
        end
    end
    %k0=2.*pi./Lamda(q);
    u=(ux1+ux2)./2;
    ne(q)=sqrt(n1.^2-(u./k0).^2);

end

N=ne;