time=(0:10:80)'; %matematiksel metodu kullanabilmek adına temsili bir zaman aralığı belirledim
latitude = [41.324729, 41.323208, 41.321382 , 41.323086 , 41.323512 , 41.326006 , 41.328622, 41.329656, 41.332090 ]'; %enlem verileri
longtitude = [26.810825, 26.814471, 26.817654 , 26.823382 , 26.828486 , 26.831727 , 26.836426, 26.840557, 26.844203 ]'; %boylam verileri
data = table(time,latitude,longtitude); %enlem ve boylam verilerinden data table oluşturulur.
subplot(2,2,1);plot(time,latitude,'o'); xlabel('Zaman , sn'); ylabel('Enlem , DD'); title('Enlem-Zaman '); %enlem ve boylam verileri grafiklerini çizer
subplot(2,2,2);plot(time,longtitude,'o'); xlabel('Zaman , sn'); ylabel('Boylam , DD'); title('Boylam-Zaman ');
p_latitude=polyfit(data.time, data.latitude,5);%enlem için 5.dereceden polinomik denklem oluşturur.
p_=polyfit(data.time, data.longtitude,5); %boylam için 5.dereceden polinomik denklem oluşturur.
p_latitude_2=polyfit(data.time, data.latitude,3); %enlem için 3.dereceden polinomik denklem oluşturur.
p_longtitude_2=polyfit(data.time, data.longtitude,3); %boylam için 3.dereceden polinomik denklem oluşturur.
%Burada 5.dereceden ve 3.dereceden denklemlerin oluşturulma nedeni bu iki
%yaklaşımın ayrı ayrı nasıl davrandığını gözlemlemektir.
t=1;
i=10;
while i<105
predict_latitude1(t)=polyval(p_latitude,i); %oluşturulan 5.dereceden enlem denklemine data table'daki değerler girilir
predict_longtitude1(t)=polyval(p_boylam,i); %oluşturulan 5.dereceden boylam denklemine data table'daki değerler girilir
predict_latitude2(t)= polyval(p_latitude_2,i) %oluşturulan 3.dereceden enlem denklemine data table'daki değerler girilir
predict_longtitude(t)= polyval(p_longtitude_2,i) %oluşturulan 3.dereceden boylam denklemine data table'daki değerler girilir
i=i+5;
t=t+1;
end
% disp(enlemTahmini);
% disp(boylamTahmini);
t=1; i=10;
subplot(2,2,3);
geoscatter(latitude,longtitude,'p','y'); title("Beşinci Dereceden Polinomik Yaklaşım");
while i<105
hold on
geoscatter(predict_latitude1(t),predict_longtitude1(t),'*'); %5.dereceden polinomik yaklaşımı haritada çizer
hold off
t=t+1;
i=i+5;
end
t=1; i=10;
subplot(2,2,4);
geoscatter(latitude,longtitude,'p','y'); title("Üçüncü Dereceden Polinomik Yaklaşım");
while i<105
hold on
geoscatter(predict_latitude2(t),predict_longtitude(t),'*'); %3.dereceden polinomik yaklaşımı haritada çizer
hold off
t=t+1;
i=i+5;
end