作者aqwsderfgtyh (aqwsderfgtyhjuik)
看板sky
标题[讨论] 我需要求救
时间Mon May 2 17:08:42 2011
关於太阳轨迹角度的观测,希望大家可以教我要怎麽打成程式
或是把它变成XYZ 谢谢
我用matlab算出三个角度
fi(入射角)
alpha(高度角)
as(方位角)
自己写出的程式如下
L=input('纬度');
hs=90
n=1
for n=1:1:365
for hs=90:-15:-90
dels=23.45*sind(360*(284+n)/365);
sin_alpha=cosd(L)*cosd(hs)*cosd(dels)+sind(L)*sind(dels);
alpha=asind(sin_alpha);
fi=acosd(sin_alpha);
as=acosd(secd(alpha)*(cosd(L)*sind(dels)-cosd(dels)*sind(L)*cosd(hs)));
fprintf(1, 'n %4.0f hs %4.0f alpha %4.4f fi %4.4f as %4.4f \n', n, hs , alpha , fi , as \n);
end
end
已经跑成功了,接下来只要抓出我需要的日子,调整n的大小就可以画图了
但是我不会画...
请问有人能帮忙吗? 谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 118.168.182.222
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 118.168.182.222
※ 编辑: aqwsderfgtyh 来自: 118.168.182.222 (05/02 17:09)
1F：推 kyphosis:要不要去matlab板问? 05/02 21:22
2F：推 knml:画图应该很简单吧,都跑成功了画图也没问题 05/02 22:26
3F：推 ckpt:这数值有几个问题 1.方位角在上午和下午时 是相同的 05/02 22:56
4F：推 ckpt:2.没有加入公转轨道偏心率 我一直画不出来8字轨迹 05/02 22:58
5F：→ ckpt:然後要画图 在两个end前一行 加入 plot(n,as);hold on; 即可 05/02 22:59
6F：→ ckpt:1.的解决方式 可能可以设两回圈 跑 hs = 0~90 和 -90~0 05/02 23:01
7F：→ ckpt:再调整 as=acosd plot(变数1,变数2) 2D 可自己放 05/02 23:02
8F：→ ckpt:想画3D 也是用 plot3(x1,x2,x3);hold on; 但问题2最好弄一下 05/02 23:04