Facebook
From Sharp Marmoset, 9 Years ago, written in Plain Text.
Embed
Download Paste or View Raw
Hits: 793
  1. function [fib,lb,Ab] = kivioj(fi,l,A)
  2. a=6378137;
  3. e2=0.00669438002290;
  4. s=25000;
  5. A=66*pi/180;
  6. fi=0.90805;
  7. l=0.32100;
  8. n=round(s/1000);
  9. ds=s/n
  10. for i=1:n
  11.     M(i)=a*(1-e2)/sqrt((1-e2*(sin(fi(i))^2))^3);
  12.     N(i)=a/sqrt(1-e2*(sin(fi(i)))^2);
  13.     dfi(i)=cos(A(i))*ds/M(i);
  14.     dA(i)=sin(A(i))*tan(fi(i))*ds/N(i);
  15.     fim(i)=fi(i)+(1/2)*dfi(i);
  16.     Am(i)=A(i)+(1/2)*dA(i);
  17.     Mm(i)=a*(1-e2)/sqrt((1-e2*(sin(fim(i))^2))^3);
  18.    
  19.     Nm(i)=a/sqrt(1-e2*(sin(fim(i)))^2);
  20.     dfim(i)=ds*cos(Am(i))/Mm(i);
  21.     dlm(i)=ds*sin(Am(i))/Nm(i)*cos(fim(i));
  22.     dAm(i)=ds*sin(Am(i))*tan(fim(i))/Nm(i);
  23.     fi(i+1)=fi(i)+dfim(i);
  24.     l(i+1)=l(i)+dlm(i);
  25.     A(i+1)=A(i)+dAm(i);
  26. end
  27. fib=fi(n)
  28. lb=l(n)
  29. Ab=A(n)
  30. end
  31.