吴斌宋世超王福来.docVIP

一．引言 太阳黑子是在太阳的光球层上发生的一种太阳活动，是太阳活动中最基本、最明显的活动现象黑子一般成群出现在太阳表面，天文学家又将其称为“黑子群”。黑子的形成周期短，形成后几天到几个月就会消失，新的黑子又会产生。太阳黑子是太阳活动的重要标志，其活动存在着明显的周期性，周期平均为1年。黑子群对地球的磁场和电离层会造成干扰往往会对地球的气候、通信、生物产生不同程度的影响。(R.Wolfer)根据在过去的288 年(1700年~1987 年)间每年太阳黑子出现的数量和大小的观测数据推算出11 年的周期规律。我们现在可以利用Matlab强大的数据处理与仿真功能，对太阳黑子数进行功率谱密度分析从而可以得到对太阳黑子活动周期的结论。 二．实验原理 在该实验中，我先对太阳黑子数进行均值处理，然后用原始数据减去均值再进行自相关处理，使其图形向下平移得到较好的图像。利用自相关函数来对它的规律进行统计，因为自相关处理后其周期是保持不变的。对自相关后的太阳黑子数（amount）序列做FFT变换后得到Y（长度为n），只取其前 个数据的功率谱密度的估计值 。原因是时域为离散的实序列的傅立叶变换对应于具有周期性且偶对称的频域特性，因此Y的前 个数据已经包含了天阳黑子数的全部信息。根据DFT的频域单位 与DTFT的频域单位 的表达关系式 以及 与f对应关系 ，可以看出 与f呈线性关系 ；同样地，因为Y的前 个数据已经包含了太阳黑子个数据分析功率-频率图时，对应的横坐标时应取 。（注： 为采样频率） 三．方案组成及流程图 图一流程图 四．实现步骤及内容 (1) 插入并读取数据并绘制太阳黑子时域图，观察太阳黑子的活动现象。程序如下：figure plot(years,amount)%画时域图（年份，数量） xlabel(Years); ylabel( Sunspot Data ); title(Sunspot Data) %标注横纵坐标和标题 pause (2) 对太阳黑子原始数据求均值并相减后进行自相关处理，并观察图形。程序如下： M=mean（size）；%求太阳黑子原始数据的均值 Z=size-M; C=xcorr(Z,unbiased);%自相关处理 plot(C) pause (3) 对已经得到的自相关数应用FFT技术（傅里叶变换）分析它的变化规律。程序如下： figure Y=fft(C);%傅里叶变换 Y(1)=[];%由于算法问题，Y(1)为所有fft数值之和，因此需舍弃第一个点 plot(Y,ro)%在复平面做图，空心点 title(在复平面的傅里叶系数); xlabel(Real Axis); %标注横坐标为’实部’ ylabel(Imaginary Axis); %标注纵坐标为’虚部’ pause (4) 功率与频率的关系曲线则被定义为周期图。用周期图法对自相关数进行功率谱密度分析。程序如下： figure n=length(Y);%取FFT结果长度为n power = abs(Y(1:n/2)).^2; %取前n/2个数据求其功率 nyquist = 1/2; %取最大频率为0.5 freq = (1:n/2)/(n/2)*nyquist; %将FFT转换为DTFT（离散时间傅里叶变换plot(freq,power) %画周期图 xlabel(cycles/year);title(Periodogram) pause (5) 确定出太阳黑子的活动周期。为清楚起见，画出功率与周期（频率的倒数）的关系曲线图。程序如下： figure period=1./freq; %将横坐标单位转换为年/周期 plot(period,power);%画周期图 axis([0 100 0 2e+7]); ylabel(Power); xlabel(Period (Years/Cycle));%单位为年/周期 pause (6) 在功率与周期关系曲线图中标出功率的最高点，该位置对应的周期则为太阳黑子活动周期。 hold on; index=find(power==max(power));%找到频率最大点，该点横坐标即为太阳黑子周期 mainPeriodStr=num2str(period(index)); plot(period(index),power(index),r., MarkerSize,25);%用实心点指出该点 text(period(index)+2,power(index),[周期 = ,mainPeriodStr]);%文字标注该点 holdoff; 五 实验结果分析与总结 实验得出的太阳黑子的活动周期为11.04年,说明太阳黑子看似无规律的现象，还是有其的宏观规律的，这对我们太阳活动及其周围的影响有巨