limitfunc.m

高频地波雷达完整仿真matlab程序 网上找到的

function [sitaL,sitaR]=limitfunc(eta)
 if (abs(eta)>=0)&(abs(eta)<=sqrt(2))
     sitaL=-pi;
     sitaR=pi;
 elseif abs(eta)>sqrt(2)
     sitaR=pi-acos(2/eta^2);
     sitaL=-sitaR;
 end