Сгенерируйте истинный маленький круг, loxomdromic маленький круг и ограничивающий случай большого круга.
Отобразите оси карты с ортогональной проекцией.
figure; axesm ortho; gridm on; framem on setm(gca,'Origin', [45 30 30], 'MLineLimit', [75 -75],... 'MLineException',[0 90 180 270])
Задайте три центральных точки на сфере.
A = [45 90]; B = [0 60]; C = [0 30];
Создайте три маленьких круга вокруг этих трех центральных точек. Первым, sca
, является истинный маленький круг. Вторым, scb
, является loxomdromic маленький круг. Третьим, scc
, является большой круг как ограничивающий случай маленького круга.
sca = scircle1(A(1), A(2), 20);
scb = scircle2(B(1), B(2), 0, 150);
scc = scircle1('rh',C(1), C(2), 20);
Отобразите точки и их соответствующие маленькие круги с различными цветами. Маркируйте маленькие круги.
plotm(A(1), A(2),'ro','MarkerFaceColor','r') plotm(B(1), B(2),'bo','MarkerFaceColor','b') plotm(C(1), C(2),'mo','MarkerFaceColor','m') plotm(sca(:,1), sca(:,2),'r') plotm(scb(:,1), scb(:,2),'b--') plotm(scc(:,1), scc(:,2),'m') textm(50,0,'Normal Small Circle') textm(46,6,'(20\circ from point A)') textm(4.5,-10,'Loxodromic Small Circle') textm(4,-6,'(20\circ from point C') textm(-2,-4,'in rhumb line sense)') textm(40,-60,'Great Circle as Small Circle') textm(45,-50,'(90\circ from point B)')