how to plot a subject class

Hello

I want to run same code 5 times with different satellite. How can I do it?

each satellite has different parameters are as below :-

function gps_(a,e,eta,omega_1,omega_2)=(26500,0.01,398600.44,0,0)

function gnss_(a,e,eta,omega_1,omega_2)=(42164.140,0,63,0,0)

function galileo_(a,e,eta,omega_1,omega_2)=(29994,0,56,0,0)

function champ_(a,e,eta,omega_1,omega_2)=(6838,0.004,87,0,0)

function molniya_(a,e,eta,omega_1,omega_2)=(26554,0.7,65,245,270)

so above 5 satellite how can I repeat 5 parameters again and again? Left hand side in brackets notation of the 5 parameters and on right hand side the value of the 5 parameters. Thank you :)

k=1;
    for i=0:10:2*T
    M= n *(i-0);                 %mean anomly
    E0=M;
    err=1;
    while(err>10^-12)
    E= M+e*sin(E0);
    err=abs(E-E0);         % I want to understand this step
    E0=E;                  % eccentric anomly
    end
    v = 2*atan((sqrt(1+e/1-e)) *(tan(E/2)));     %atan2 is not working   % true anomly


    r = a*(1-e*cos(E));                          %distance of the satellite

    r_b= [r*cos(v);r*sin(v);0]                  %position vector in the orbit system
    r_b1(k,:)=[a*(cos(E)-e);a*(sqrt(1-(e*e)))*sin(E);0];   %position vector in the orbit system


       v_b(k,:)=[-a*n*sin(E)*a/r;a*n*sqrt(1-(e*e))*cos(E)*a/r]    %velocity vector in the orbit system
    fi=deg2rad(0)
    
    eta=deg2rad(55)
    omega1=deg2rad(0)
    omega2=deg2rad(0)
     r_3= [cos(-fi) sin(-fi) 0;-sin(-fi) cos(-fi) 0;0 0 1]
     r_1=[1 0 0;0 cos(-eta) sin(-eta);0 -sin(-eta) cos(-eta)]
     r_i= r_3*r_1*r_3*r_b
     r_i1(k,:)=transpose(r_i)

        w_a=((2*pi)/86164)
       phi=w_a*(i-0)
       r_phi= [cos(phi) sin(phi) 0;-sin(phi) cos(phi) 0;0 0 1]
       r_e=r_phi*r_i   
       r_e1(k,:)=r_phi*r_i 
       z_lat(k,:)= r_e1(k,3)
       y_lat(k,:)=r_e1(k,2)
       x_lat(k,:)=r_e1(k,1)
      lat(k,:)= atan2((z_lat(k,:)),(sqrt((x_lat(k,:)*x_lat(k,:))+(y_lat(k,:)*y_lat(k,:)))))
      long(k,:)= atan2(y_lat(k,:),x_lat(k,:))  
         
   
   r_w=[4075.53022;931.78130;4801.61819]
   lat_w(k,:)=atan2(r_w(3,1),(sqrt(r_w(1,1)*r_w(1,1)+r_w(2,1)*r_w(2,1))))
   long_w(k,:)=atan2(r_w(2,1),r_w(1,1))
   A  = [-sin(lat_w(k,:))*cos(long_w(k,:)) -sin(long_w(k,:)) cos(lat_w(k,:))*cos(long_w(k,:));-sin(lat_w(k,:))*sin(long_w(k,:)) cos(long_w(k,:)) cos(lat_w(k,:))*sin(long_w(k,:));cos(lat_w(k,:)) 0 sin(lat_w(k,:))]
   r_top = transpose(A)*(r_e-r_w)
   r_top1(k,:)=transpose(A)*(r_e-r_w)

   
   S=sqrt(r_top(1,1).^2+r_top(2,1).^2+r_top(3,1).^2)
   S_1(k,:)=sqrt(r_top1(k,1).^2+r_top1(k,2).^2+r_top1(k,3).^2)
   
zenit(k,:)=acos(r_top(3,1)/S);
Elevation(k,:)=(pi/2) -zenit(k,:)
Azi(k,:)= atan2(-r_top1(k,2),r_top1(k,1))
 
   k=k+1;
    end
    
   fig1= plot3(r_e1(:,1),r_e1(:,2),r_e1(:,3))   
    fig2=plot3(r_i1(:,1),r_i1(:,2),r_i1(:,3))
    
    pax=polaraxes;
    pax.ThetaDir = 'clockwise';
  fig3=polarplot(Azi,rad2deg(Elevation))
  pax.RDir = 'reverse';
  
  
    %fig4=polarplot(zenit,rad2deg(Elevation))
    fig5=plot(long_deg,lat_deg)
    hold on
    plot(coastlon,coastlat)
    hold off
   mean_r_e=mean(r_e)
  closest= mean(r_e-r_w) 

please guide me ?

Thank you

Kind regards,

Vims Rocz

[Geoff Hayes]

vimal - how do the parameters a,e,eta,omega_1,omega_2 correspond to your code? Is omega_1 the same as omega1?

Thie first thing you may want to do is change your code from a script to a function so that you can pass in the parameters to the function for each of your different satellites. For example,

function mySatelliteCode(a,e,eta,omega_1,omega_2)

and save this to a file named mySatelliteCode.m. You would then update your code so that it uses the input parameters instead of what you already have. You would then call this function (from the command line or another function or script) and pass in these parameters for the satellite of interest

>> mySatelliteCode(26500,0.01,398600.44,0,0)

[vimal kumar chawda]

yes, Omega_1 is same as omega1 sorry for typing errror.

function mySatelliteCode(a,e,eta,omega_1,omega_2)

I have to make my file name as mySatelliteCode. Then below changes but how it will take gps then gnss then galileo then champ then molniya. Can you please give one small example so that I can impliment over here.

like as below:-

mySatelliteCode(26500,0.01,398600.44,0,0)
mySatelliteCode(42164.140,0,63,0,0)
mySatelliteCode(29994,0,56,0,0)
mySatelliteCode(6838,0.004,87,0,0)
mySatelliteCode(26554,0.7,65,245,270)
code

but i want to differentiate by mentioning name of the satellite as well?

you can find in highlighted blue section and 5 satellite parameters as top. Currently i am having 5 different file for each satellite but want to unite it.

Thank you very much for the response.

[Geoff Hayes]

Then below changes but how it will take gps then gnss then galileo then champ then molniya

You call the mySatelliteCode function five times, once for each satellite.

but i want to differentiate by mentioning name of the satellite as well?

I'm not sure what you mean...where do you want the satellite name to be mentioned? In a figure? Just pass the name of the satellite as an additional parameter to your function. Then reference this parameter in your code.

[vimal kumar chawda]

like calling five times means

sorry sir but I didnot get you. can you please provide me an example so that i can do? how to call?

2- where do you want the satellite name to be mentioned? In a figure? Just pass the name of the satellite as an additional parameter to your function. Then reference this parameter in your code.

mySatelliteCode_gps(26500,0.01,398600.44,0,0)
mySatelliteCode_gnss(42164.140,0,63,0,0)
mySatelliteCode_galileo(29994,0,56,0,0)
mySatelliteCode_champ(6838,0.004,87,0,0)
mySatelliteCode_molniya(26554,0.7,65,245,270)

how to call this satellite parameters 5 times ?

Here I have to mention the parameters which i have mentioned in the code are as function mySatelliteCode(a,e,eta,omega_1,omega_2) so how one can do both at the same time?

thank you sir and sorry for late response I was having an exam.