太阳黑子活动周期的分析借鉴.pdfVIP

一、 引言 太阳黑子是人们最早发现也是人们最熟悉的一种太阳表面活动。 因为太阳内 部磁场发生变化， 太阳黑子的数量并不是固定的， 它会随着时间的变化而上下波 动，每隔一定时间会达到一个最高点，这段时间就被称之为一个太阳黑子周期。 太阳黑子的活动呈现周期性变化是由施瓦贝首次发现的。沃尔夫 (R.Wolfer) 继 而推算出 11年的周期规律。 实际上， 太阳黑子的活动不仅呈 11年的周期变化， 还 有海耳在研究太阳黑子磁场分布时发现的 22年周期；格莱斯堡等人发现的 80年周 期以及蒙德极小期等。 由于太阳黑子的活动规律极其复杂， 时至今日科学家们仍 在努力研究其内在的规律和特性。 事实上，对太阳黑子活动规律的研究不仅具有 理论意义， 而且具有直接的应用需求。 太阳黑子的活动呈现周期性变化的， 沃尔 夫 ( R.Wolfer) 根据在过去的 288 年( 1700年~1987 年) 间每年太阳黑子出现的数 量和大小的观测数据推算出 11 年的周期规律。我们利用 Matlab 强大的数据处理 与仿真功能， 对Wolfer 数进行功率谱密度分析从而可以得到对太阳黑子活动周期 的结论。 二、 实验原理 n （长度为n） 在该实验中， 对Wolfer数序列做 FFT变换后得到 Y ，只取其前 2 个数据的功率谱密度的估计值 |Y |2 。原因是时域为离散的实序列的傅立叶变 n 换对应于具有周期性且偶对称的频域特性，因此 Y 的前 个数据已经包含了 2 Wolfer数的全部信息。 根据 DFT 的频域单位 k 与DTFT 的频域单位 的表达关 2 f 系式 k 以及 与f对应关系 =2 ，可以看出 k 与f 呈线性关系 N f s f n f s k ；同样地，因为 Y 的前 个数据已经包含了 Wolfer数的全部信息， N 2 n fs 只取前 个数据分析功率 -频率图时，对应的横坐标时应取 f k 。（注： 2 N