上海天文台实习的一个项目-根据卫星数据绘制南极星空分布图

       大三的时候做的一个东西，当时在上海天文台实习了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