fsm_eh.m

本软件使用matlab预言编写

function [ EH ] = FSM_EH( Lamda,u,KJ,d )
%UNTITLED1 Summary of this function goes here
%  Detailed explanation goes here
n1=1;
n2=SiN(Lamda);
ntmp=(n2.^2)./(n1.^2);
%SetH(2.3,1);
%[KJ,d]=LoadH();
R=KJ.*sqrt(sqrt(3)./(2.*pi));   %外径
r=d/2;                          %内径
k0=2.*pi./Lamda;
v=k0.*sqrt(n2.^2-n1.^2);
%w=k0.*sqrt(n.^2-n1.^2);
w=sqrt(v.^2-u.^2);
%u=k0.*sqrt(n2.^2-n.^2);

g1=besselj(0,r.*u).*bessely(1,R.*u)-bessely(0,r.*u).*besselj(1,R.*u);
g2=besselj(1,r.*u).*bessely(1,R.*u)-bessely(1,r.*u).*besselj(1,R.*u);
g=g1./(u.*r.*g2)-1./(u.*r).^2;

f=(1./r.^4).*(1./u.^2+1./w.^2).*((n2./u).^2+(n1./w).^2);

s1=besseli(2,w.*r)./besseli(1,w.*r)+1./(w.*r);
s2=w.*r.*(1+ntmp).*g./2;
s3=w.*r.*sqrt(((1-ntmp).*g./2).^2+f./n1.^2);
EH=s1+s2+s3;