# sky plot for the satellite

113 views (last 30 days)
vimal kumar chawda on 18 May 2019
Answered: Joga Setiawan on 1 Aug 2020
Hello Schlor
I am having azimuth, zenit and elevation. So how do I plot skyplot for the satellite .
I have code as below are as
polarplot(azi,ele)
It should be only visible only 0° to 90°.
Can anyone help me?
Thank you

Sulaymon Eshkabilov on 18 May 2019
Can you be a little bit more explcit with your problem statement by providing some more info and data (if possible) of your variables ...?
vimal kumar chawda on 20 May 2019
Hallop
I have uploaded the variable of all but please consider azi and ele.
Thank you
vimal kumar chawda on 20 May 2019
Here we are seeing plot as a top view but I need 90 at center. That means its skyplot.
So how can I bring 90° at center and other corner of the image is direction like North =N ; e,w,s.
figure
pax = polaraxes;
theta = 0:0.01:2*pi;
rho = sin(2*theta).*cos(2*theta);
polarplot(theta,rho)
pax.FontSize = 12;
Thank you

vimal kumar chawda on 23 May 2019
I got the result please find the below:
you can use this and found out how skyplot are works . In normal case we are getting polarplotresult top view. In this we get from aside. One can plot for different structure. Hope you enjoy.
pax=polaraxes;
pax.RDir = 'reverse';

Joga Setiawan on 1 Aug 2020
This is an example to get the skyplot of an antenna tracking a LEO satellite, shown in the subplot(223)
t_AZ_EL=[0 41 0.4;
30 43.6 2.4;
60 45.6 3.9;
90 48.1 5.6;
120 51.1 7.4;
150 54.7 9.5;
180 58.5 11.5;
210 62.9 13.4;
240 68 15.4;
270 73.6 16.9;
300 79.7 17.3;
330 86.7 20.3;
360 94.3 21.3;
390 102.1 21.7;
420 110.6 21.6;
450 117.9 21;
480 125.7 19.7;
510 131.9 18.3;
540 137.9 16.5;
570 143.2 14.7;
600 148.3 12.5;
630 151.7 10.9;
660 155.1 9.1;
690 158.4 7.1;
720 161.1 5.4;
750 163.8 3.4;
780 166 1.7;
810 167.3 0.7]
t=t_AZ_EL(:,1);
AZ=t_AZ_EL(:,2);
EL=t_AZ_EL(:,3);
subplot(221)
plot(t,AZ)
grid
xlabel('time (s)')
ylabel('Azimuth (deg)')
subplot(222)
plot(t,EL)
grid
xlabel('time (s)')
ylabel('Elevation (deg)')
subplot(223)
polarplot(AZ*pi/180,EL)
ax = gca;
ax.ThetaZeroLocation = 'top';
ax.RDir = 'reverse';
rlim([0 90])
rtickangle(20)
thetaticks([0 45 90 135 180 225 270 315])
thetaticklabels({'North','45','East','135','South','225','West','315'})

R2019a

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!