大三的时候做的一个东西,当时在上海天文台实习了15天,在老师要求下写了下面一大段代码,现在想起来,都不记得这是干的什么了。
人最可悲的,就是很多东西本来学会了,结果没过多久就忘了。
我只记得,当时是美国的一颗卫星,我们登了它的网站,下载了它的观测数据,经过筛选,大概有4万多个点。然后用这些数据导入进来,根据一些计算,求出绝对星等,以及距离,然后绘制南极星空分布图,用的是极坐标。最后证实了,恒星确实是呈网状分布。
据说,恒星的分布就像丝袜一样:相近的恒星喜欢聚成团,然后各个星系团中间是巨大的空洞。
clc;
clear;
format long
Um=0.3;
Uk=0;
Uv=0.7; %基本参数赋值
z=csvread('shixi_gliming.csv',0,6,'G1..G44103');%导入红移z
for i=1:44103
DC(i)=quad('integrand',0,z(i),0.00001); %计算积分
DM(i)=DC(i); %因为k=0
DL(i)=(1+z(i))*DM(i); %计算光度距离
end
dered_r=csvread('shixi_gliming.csv',0,4,'E1..E44103'); %导入视星等r
for i=1:44103
if z(i)<0.1
p0=-0.194858;
p1=0.912570;
else
p0=-0.242659;
p1=1.38731;
end
M(i)=de