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