太阳黑子活动周期的分析

一、 引言

太阳黑子是人们最早发现也是人们最熟悉的一种太阳表面活动。因为太阳内部磁场发生变化，太阳黑子的数量并不是固定的，它会随着时间的变化而上下波动，每隔一定时间会达到一个最高点，这段时间就被称之为一个太阳黑子周期。太阳黑子的活动呈现周期性变化是由施瓦贝首次发现的。沃尔夫 (R.Wolfer)继而推算出11年的周期规律。实际上，太阳黑子的活动不仅呈11年的周期变化，还有海耳在研究太阳黑子磁场分布时发现的22年周期；格莱斯堡等人发现的80年周期以及蒙德极小期等。由于太阳黑子的活动规律极其复杂，时至今日科学家们仍在努力研究其内在的规律和特性。事实上，对太阳黑子活动规律的研究不仅具有理论意义，而且具有直接的应用需求。太阳黑子的活动呈现周期性变化的，沃尔夫(R.Wolfer)根据在过去的288 年(1700年~1987 年)间每年太阳黑子出现的数量和大小的观测数据推算出11 年的周期规律。我们利用Matlab 强大的数据处理与仿真功能，对Wolfer 数进行功率谱密度分析从而可以得到对太阳黑子活动周期的结论。

二、 实验原理

在该实验中，对Wolfer 数序列做FFT 变换后得到Y （长度为n ），只取其前

n

2

个数据的功率谱密度的估计值2|Y |。原因是时域为离散的实序列的傅立叶变换对应于具有周期性且偶对称的频域特性，因此Y 的前

n

2

个数据已经包含了Wolfer 数的全部信息。根据DFT 的频域单位k 与DTFT 的频域单位ϖ的表达关系式2k N

π

ϖ=

⋅以及ϖ与f 对应关系=2s f f ϖπ⋅，可以看出k 与f 呈线性关系

s f f k N

=

⋅；同样地，因为Y 的前n

2个数据已经包含了Wolfer 数的全部信息，

只取前

n

2

个数据分析功率-频率图时，对应的横坐标时应取2

s f f k N =⋅。（注：

s f 为采样频率）

三、方案组成与流程图

1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111 1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111

四、实现步骤及内容

(1) 读取数据并绘制Wolfer 图，观察太阳黑子的活动现象。程序如下： load sunspot_sun.dat %读取数据 year=sunspot(:,1); %读取年份信息

wolfer=sunspot(:,2); %读取黑子活动数据 figure%新建图像

plot(year,wolfer) %画出时域图

xlabel('Years'); ylabel(' Sunspot Data '); title('Sunspot Data') %标注横纵坐标和标题

pause

1700

1750

1800

18501900

1950

2000

020406080100120140160180200Years

S u n s p o t D a t a

Sunspot Data

为清楚起见，画出部分时间段(如实验中采取开始50年)的Wolfer 图进行观察。