townini6

111111111111
R=6371;
mnoznik=mean(2*pi*R*cosd(x))/360;

for i=1:125
    dt(i)=(t(i+1)-t(i))/3600000;
    dx(i)=(x(i+1)-x(i))*111.1;
    dy(i)=(y(i+1)-y(i))*mnoznik;
    dz(i)=(z(i+1)-z(i))*0.001;
end
Vx=zeros(125,1);
Vy=zeros(125,1);
Vz=zeros(125,1);
Vx=dx./dt;
Vy=dy./dt;
Vz=dz./dt;
V=sqrt((Vx).^2 + (Vy).^2 + (Vz).^2);
vsr=mean(V);
for i=1:124
dV(i)=(V(i+1)-V(i));
end
dV=(dV./3.6)';
for n=1:125
    d(n)=sqrt((dx(n))^2+(dy(n))^2+(dz(n))^2);
end
d=d';

dt2=diff(t(2:126)/1000);
A=dV./dt2;
asr=mean(A);


for i=1:125
    a(i)=(sind((x(i+1)-x(i))/2)^2)+cosd(x(i))*cosd(x(i+1))*(sind((y(i+1)-y(i))/2)^2);
    c(i)=2*atan2(sqrt(a(i)),sqrt(1-a(i)));
    d(i)=R*c(i);
end
s1=sum(d);

s2=0;
for i = 1:124
    s2 = s2 + (V(i)+V(i+1))*dt(i)/2;
end


----------------------------------------------------------------------------------

111111
R=6371;
mnoznik=(mean(2*pi*R*cosd(x))/360)*1000;


x2=x*mnoznik;
y2=y*111100;
t2=t/1000;

s(1)=0;
for n=2:126
    s(n)=s(n-1)+sqrt((x2(n)-x2(n-1))^2+(y2(n)-y2(n-1))^2+(z(n)-z(n-1))^2);
end
s=s';
for n=2:126
    ds(n-1)=sqrt((x2(n)-x2(n-1))^2+(y2(n)-y2(n-1))^2+(z(n)-z(n-1))^2);
end
ds=ds';
for n=2:126
    dt(n-1)=t2(n)-t2(n-1);
end
dt=dt';

v2=ds./dt;

v=diff(s)./diff(t2);
vsr=mean(v);

-------------------------------------------------------
TRZECIE:
x = 0:(2*pi)/100:2*pi;
y1 = sin(x);
y2 = cos(x);
y1=awgn(y1,50);
for i=2:length(x)
    dy1(i-1)=y1(i)-y1(i-1);
end
y1dif1=dy1/((2*pi)/100);
dy2(1)=dy1(1);
for i=2:length(x)-1
    dy2(i)=0.5*y1(i+1)-0.5*y1(i-1);
end
y1dif2=dy2/((2*pi)/100);

-----------------------------
333:
K=10; L=100; N=K*L; n=0:N-1; x = 2*pi/L*n;
y1=sin(x);
y1=awgn(y1,50);
y2 = (2*pi/L)*cos(x); % y2 to pochodna y1 względem indeksu czasu n
M=50;
m=-M:1:M;
h = cos(pi*m)./m;
h(M+1)=0; % odp. impulsowa filtru różnicz.
w = blackman(2*M+1)'; h = h .* w; % opcjonalne wygładzenie brzegów
y2e = conv(y1, h); % operacja splotu (filtracji)

for i=2*M:length(y2e)-2*M
    y2efil(i+1-2*M)=y2e(i);
end

for i=M:length(y2)-M
    y2fil(i+1-M)=y2(i);
end

-------------------------------------------------------