Campos electromagnéticos Transformación de Coordenadas por medio de Software

Alejandro Zabala Camacho

Código: 223655

Campos electromagnéticos

Transformación de Coordenadas por medio de Software

De coordenadas rectangulares a cilíndricas y esféricas

Px=input("ingrese coordenada x del punto");

Py=input("ingrese coordenada y del punto");

Pz=input("ingrese coordenada z del punto");

AX=input("ingrese coordenada x del vector");

AY=input("ingrese coordenada y del vector");

AZ=input("ingrese coordenada z del vector");

if Px^2+Py^2==0

   phi=0;

elseif Py>0

    phi=(acos(Px/(sqrt(Px^2+Py^2))));

else

     phi=((2*%pi)-(acos(Px/(sqrt(Px^2+Py^2)))));

end

ro=sqrt(Px^2+Py^2);

ap=(AX*cos(phi))+(AY*sin(phi));

ao=-(AX*sin(phi))+(AY*cos(phi));

az=AZ;

disp('coordenadas cilindricas del punto');

printf('ro=%f Fi=%f z=%f] ',ro,phi*(180/%pi),Pz);

disp('componentes del vector en coordenadas cilindricas');

printf('[%fap %fafi %faz] ',ap,ao,az);

if  Px^2+Py^2+Pz^2==0 then

  theta=0;

else

    theta=(acos(Pz/(sqrt(Px^2+Py^2+Pz^2))));  

end

r=sqrt(Px^2+Py^2+Pz^2);

ar=(AX*sin(theta)*cos(phi))+(AY*sin(theta)*sin(phi))+(AZ*cos(theta));

atheta=(AX*cos(theta)*cos(phi))+(AY*cos(theta)*sin(phi))-(AZ*sin(theta));

disp('coordenadas esfericas del punto');

printf('r=%f Theta=%f Fi=%f] ',r,theta*(180/%pi),phi*(180/%pi));

disp('componentes del vector en coordenadas esfericas');

printf('[%far %fatheta %fafi] ',ar,atheta,ao);

exec('C:\Users\Alejo\Documents\coor.sce', -1)

Ingrese coordenada x del punto3

Ingrese coordenada y del punto4

Ingrese coordenada z del punto5

Ingrese coordenada x del vector3

Ingrese coordenada y del vector4

Ingrese coordenada z del vector5

 Coordenadas cilíndricas del punto  

ro=5.000000 Fi=53.130102 z=5.000000]  

 Componentes del vector en coordenadas cilíndricas  

[5.000000ap 0.000000afi 5.000000az]  

 Coordenadas esféricas del punto  

r=7.071068 Theta=45.000000 Fi=53.130102]  

 Componentes del vector en coordenadas esféricas  

[7.071068ar -0.000000atheta 0.000000afi]

De coordenadas cilíndricas a rectangulares y esféricas

Pro=input("Ingrese coordenada rho del punto");

Pfi=input("Ingrese coordenada phi del punto en grados");

Pz=input("Ingrese coordenada z del punto");

RO=input("Ingrese coordenada rho del vector");

FI=input("Ingrese coordenada phi del vector");

Z=input("Ingrese coordenada z del vector");

Px=Pro*cos(Pfi*%pi/180);

Py=Pro*sin(Pfi*%pi/180);

Pz=Pz;

Ax=(RO*cos(Pfi*%pi/180))-(FI*sin(Pfi*%pi/180));

Ay=(RO*sin(Pfi*%pi/180))+(FI*cos(Pfi*%pi/180));

Az=Z;

disp('Coordenadas rectangulares del punto');

printf('x=%f y=%f z=%f ',Px,Py,Pz);

disp('Componentes del vector en coordenadas rectangulares');

printf('[%fax %fay %faz] ',Ax,Ay,Az);

if  Px^2+Py^2+Pz^2==0 then

  theta=0;

else

    theta=(acos(Pz/(sqrt(Px^2+Py^2+Pz^2))));  

end

r=sqrt(Px^2+Py^2+Pz^2);

ar=(Ax*sin(theta)*cos(Pfi*%pi/180))+(Ay*sin(theta)*sin(Pfi*%pi/180))+(Az*cos(theta));

atheta=(Ax*cos(theta)*cos(Pfi*%pi/180))+(Ay*cos(theta)*sin(Pfi*%pi/180))-(Az*sin(theta));

ao=-(AX*sin(phi))+(AY*cos(phi));

disp('Coordenadas esféricas del punto');

printf('r=%f Theta=%f Fi=%f ',r,theta*(180/%pi),Pfi);

disp('Componentes del vector en coordenadas esféricas');

printf('[%far %fatheta %fafi] ',ar,atheta,ao);

Ingrese coordenada rho del punto5

...