作者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