دانلود رایگان کد متلب تبدیل مختصات بیضوی به مختصات دکارتی

%programmer Midya Rostami

%this program Convert the

%Elipsoied Cordinate to Decart Cordinate

 clc 

clear

 format long

r= input('Enter the Radious in meter=');

 x= input('Enter the x in degree='); 

 y= input('Enter the y in degree='); 

 z= input('Enter the z in degree='); 

 lambda= atand(y/x); 

 phi= atand(z/(sqrt((x^2)+(y^2)))); 

 Xc=r*cosd(phi)*cosd(lambda);

Yc=r*cosd(phi)*sind(lambda);

 Zc=r*sind(phi); 

 disp('------------------'); 

 disp('--Calculations:---'); 

 disp('------------------'); 

 disp(Xc);

 disp(Yc);

 disp(Zc);

disp('------------------');

منبع

دانلود رایگان کد متلب تبدیل مختصات کارتزین به مختصات منحنی الخط

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);

منبع

دانلود رایگان کد متلب تبدیل مختصات کارتزین به مختصات منحنی الخط

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);

منبع