Plotting sun path diagram problem

2021 8 8
1 个回答
更新时间：2021 8 8
20 次查看（30 天）

0 个投票

I've this code for sun path, when I run it it gives me: Matrix dimensions must agree.
How can I solve it?
clc
clear all
number_of_the_yearly_day = 1:1:365;
T=1:1:24;
declination=23.45*sin((2*pi*(number_of_the_yearly_day-80))/(365));
latitude=31.963158;
Zenith_angle=latitude-declination;
tilt_angle=Zenith_angle;
altitude=1-Zenith_angle;
w=15*(12-T);
X=asind(sin(declination).*sin(latitude)+cos(latitude).*cos(declination).*cos(w));
Y=acosd(sin(altitude).*sin(latitude)-sin(declination)./(cos(altitude).*cos(latitude)));
plot (X,Y)
ylim([-90 90]);
xlabel('Solar Azimuth (deg)');
ylabel('Solar Elevation (deg)');
title('Solar Azimuth and Elevation Angle');

请先登录，再回答此问题。

回答（1 个）

KSSV
KSSV 2021-8-8

0 个投票

clc; clear all;
clc
clear all
number_of_the_yearly_day = 1:1:365;
T=1:1:24;
declination=23.45*sin((2*pi*(number_of_the_yearly_day-80))/(365));
latitude=31.963158;
Zenith_angle=latitude-declination;
tilt_angle=Zenith_angle;
altitude=1-Zenith_angle;
w=15*(12-T);
[declination,w] = meshgrid(declination,w) ;
X=asind(sin(declination).*sin(latitude)+cos(latitude).*cos(declination).*cos(w));
Y=acosd(sin(altitude).*sin(latitude)-sin(declination)./(cos(altitude).*cos(latitude)));
plot (X,Y)
ylim([-90 90]);
xlabel('Solar Azimuth (deg)');
ylabel('Solar Elevation (deg)');
title('Solar Azimuth and Elevation Angle');

2 个评论
显示 无 隐藏 无

ibrahim alzoubi
ibrahim alzoubi 2021-8-8
how can I make cos(w), the materix [declination,w] ?
KSSV
KSSV 2021-8-8
See for meshgrid.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Aerospace Blockset 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by