netzbild.m

利用ＭＡＴＬＡＢ实现了对测量中的平面网

function kk=netzbild
global ed dd pn m1 m2 m3 x y e d sid g f dir s t az c ai bi fi
kk=0;
title1   = '误差椭圆比例  ';
prompt011  = {'请输入误差椭圆比例'};
lines= 1;
d11= {'100'};
c12 = inputdlg(prompt011,title1,lines,d11);%输入误差椭圆比例对话框
scale=str2num(str2mat(c12));
hold off
title(' 控制网网图 ');
xlabel('Y坐标');
ylabel('X坐标');
plot(y(1:ed),x(1:ed),'k^')
hold on;
plot(y((ed+1):end),x((ed+1):end),'k.')
for i=1:ed+dd%写点名
    text(y(i)+3,x(i)+3,num2str(pn(i)));
end
for i=1:m1%画边长连线
    plot([y(e(i)) y(d(i))],[x(e(i)) x(d(i))],'r');
end
for i=1:m2%画方向连线
    plot([y(g(i)) y(f(i))],[x(g(i)) x(f(i))],'b');
    direction=rad_dms(dir(i));
end
for i=1:m3%画方位角连线
    plot([y(s(i)) y(t(i))],[x(s(i)) x(t(i))],'y');
end
for i=1:dd%画误差椭圆
    fe(i)=fi(i)*pi/180.;
plot(y(i+ed)+(ai(i)*cos(0:0.01:2*pi)*sin(fe(i))+bi(i)*sin(0:0.01:2*pi)*cos(fe(i)))*scale,x(i+ed)+(ai(i)*cos(0:0.01:2*pi)*cos(fe(i))-bi(i)*sin(0:0.01:2*pi)*sin(fe(i)))*scale,'r');
end
axis equal;%使得X,Y比例一致
zoom on;%利用鼠标左键放大图形，右键缩小图形
return