MATLAB Answers

sky plot for the satellite

113 views (last 30 days)
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

  3 Comments

Sulaymon Eshkabilov
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
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
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.ThetaDir = 'clockwise';
pax.FontSize = 12;
Please find the code below?
Thank you

Sign in to comment.

Accepted Answer

vimal kumar chawda
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.ThetaDir = 'clockwise';
fig3=polarplot(Azi,rad2deg(Elevation))
pax.RDir = 'reverse';

  0 Comments

Sign in to comment.

More Answers (1)

Joga Setiawan
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;
d = ax.ThetaDir;
ax.ThetaDir = 'clockwise';
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'})

  0 Comments

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!