از سری کد ها و مجموعه ی کاربردی زبان متلب در علم ژئوماتیک در ادامه کد ها را ببینید:
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 |
clc clear format long a=input('enter the a parameter Vaule='); f=input('enter the Scale factor ='); %f=1/298.25722221 b=a*(1-f); e=2*f-f^2; X=input('Enter X='); Y=input('Enter Y='); Z=input('Enter Z='); P=sqrt(X^2+Y^2); lan=2*atend((-X+sqrt(X^2+Y^2))/Y); N0=a; h0=((X^2+Y^2+Z^2)^.5)-(a*b)^.5; phi0=atand(Z/P*(1-(e^2*N0)/(N0+h0)))^(-1); N=1;h=1;phi=1; while(abs(h-h0)>.00001)&&(abs(phi-phi0)>0000001) N=a/(1-(e^2)*(sin(phi0))^2)^.5; h=P/cos(phi0)-N0; phi=atand(z/p*(1-e^2*N0/(N0+h0))^-1); N0=N;h0=h;phi0=phi; end phi=abs([fix((fix(phi)-phi)*60)+round(((ix((ix(phi)-phi)*60))-((fix(phi)-phi)*60))*60)]); phi=abs([fix((fix(lan)-lan)*60)+round(((ix((ix(lan)-phi)*lan))-((fix(lan)-lan)*60))*60)]); disp(phi); disp(lan); |
