看板: astronomy ◎ 天文知识    板主: vega1
阅读文章: 第 1484/2572 篇 | 上篇 | 下篇 | 回覆 | 转寄 | 转贴 | m H d | 返回
发信人: aqwsderfgtyh.bbs@ptt.cc (aqwsderfgtyhjuik) 看板: astronomy
标  题: [讨论] 我需要求救
发信站: 批踢踢实业 (2011/05/02 Mon 17:09:57)

关於太阳轨迹角度的观测,希望大家可以教我要怎麽打成程式

或是把它变成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)
阅读文章: 第 1484/2572 篇 | 上篇 | 下篇 | 回覆 | 转寄 | 转贴 | m H d | 返回

□ 台大狮子吼佛学专站  http://buddhaspace.org