cladn.m

本软件使用matlab预言编写

function [ N ] = CladN( Lamda,KJ,d )
%UNTITLED1 Summary of this function goes here
%  Detailed explanation goes here
%SetH(1,0.3);
u0=0.001;
u1=10;
u2=(u1-u0)/50;
u_tmp=u0:u2:u1;
%[KJ,d]=LoadH();
R=KJ.*sqrt(sqrt(3)./(2.*pi));   %外径
r=d/2;                          %内径
n=length(Lamda);
for q=1:n
    X=FSM(Lamda(q),u_tmp,KJ,d);
    m=length(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(FSM(Lamda(q),ux1,KJ,d).*FSM(Lamda(q),(ux1+ux2)./2,KJ,d)<=0)
            ux2=(ux1+ux2)./2;
        else
            ux1=(ux1+ux2)./2;
        end
    end
    k0=2.*pi./Lamda(q);
    u=(ux1+ux2)./2;
    ne(q)=sqrt(SiN(Lamda(q)).^2-(u./k0).^2);
    %plot(u_tmp,X,u,FSM(Lamda(q),u),'p'),hold on
end
%axis([u0,u1,-10,10])
%hold off
%grid on
N=ne;