| 看板: 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 |