% démonstration du calcul de la position par GPS
% exemple à deux dimensions.

t=0:pi/200:2*pi;
r=100;	% rayon de la place
p1=[r*cos(pi/4) r*sin(pi/4)];
p2=[r*cos(-pi/4) r*sin(-pi/4)];
p3=[r*cos(3*pi/2) r*sin(3*pi/2)];
r1=96;
r2=142.5;
r3=156;
subplot(1,2,1)
plot(r*cos(t),r*sin(t),...
	r1*cos(t)+p1(1),r1*sin(t)+p1(2),...
	r2*cos(t)+p2(1),r2*sin(t)+p2(2))
hold on
plot(p1(1),p1(2),'or',p2(1),p2(2),'or')
text(p1(1),p1(2),'P1')
text(p2(1),p2(2),'P2')
text(-50,-130,'Deux messagers')

axis off
axis square
axis([-r r -r r])
subplot(1,2,2)
plot(r*cos(t),r*sin(t),...
	r1*cos(t)+p1(1),r1*sin(t)+p1(2),...
	r2*cos(t)+p2(1),r2*sin(t)+p2(2),...
	r3*cos(t)+p3(1),r3*sin(t)+p3(2))
	hold on
plot(p1(1),p1(2),'or',p2(1),p2(2),'or',p3(1),p3(2),'or')
text(p1(1),p1(2),'P1')
text(p2(1),p2(2),'P2')
text(p3(1),p3(2),'P3')
text(-50,-130,'Trois messagers')
axis off
axis square
axis([-r r -r r])
set(gcf,'color','w')
hold off

