beta.m

来自「利用ＭＡＴＬＡＢ实现了对测量中的平面网」· M 代码 · 共 28 行

function [bet,k]=beta(i1,i2,i3)             %检验两方向是否构成角  
global ed dd dd1 sd g f pn ni si m2 md pn dir
k=0;
k1=0;
k2=0;
bet=0;
for i=si(i2):si(i2)+ni(i2)-1
	if(f(i)==i1) 
		k1=1;
		dir1=dir(i);
	end
	if(f(i)==i3) 
		k2=1;
		dir2=dir(i);
	end
end
if(k1==1&k2==1)
	k=1;
	bet=dir2-dir1;
	if(bet<0.)
        bet=2.*pi+bet;
	end
	if(bet>pi)
        bet=2.*pi-bet;
	end
end
return