besl.m

来自「Matlab设计FIR,IIR」· M 代码 · 共 13 行

function S=besl(x)


eps=10^(-9);
n=1;s=1;D=1;
S=0;
while D>(eps*s)
    T=x/(2*n);
    n=n+1;
    D=D*T^2;
    s=s+D;
end
S=s;