北工大太阳黑子活动周期

《北工大太阳黑子活动周期》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、信号处理与Matlab太阳黑子活动周期的分析专业:通信工程姓名：高培元学号：12024215指导老师：孙中华、张延华完成时间：2014/12/20摘要：太阳黑子是人类在研究太阳活动中最常见、最重要的现象。太阳黑子的爆发也直接影响到人们的生活方面。根据数字信号处理对数据处理的方法，采MATLAB对太阳黑子数进行处理得出太阳黑子周期约为11年。结果表明该方法对太阳黑子周期研究是直观且有意义的。关键字：太阳黑子数；活动周期；Matlab。一、引言：太阳黑子是在太阳的光球层上发生的一种太阳活动，是太阳活动中最基本、最明

2、显的。一般认为，太阳黑子实际上是太阳表面一种炽热气体的巨大漩涡，温度大约为3000-4500℃。因为其温度比太阳的光球层表面温度要低1000到2000摄氏度，所以看上去像一些深暗色的斑点。太阳黑子虽然颜色较深，但是在观测情况下，与太阳耀斑同样清晰显眼。太阳黑子很少单独活动，通常是成群出现。黑子的活动周期约为11年,天文学家把太阳黑子最多的年份称为太阳活动峰年，太阳黑子最少的年份称为太阳活动谷年。其活跃时会对地球的磁场产生影响，主要是使地球南北极和赤道的大气环流作经向流动，从而造成恶劣天气，使气候转冷。严重时会对

3、各类电子产品和电器造成损害。现今世界公认的最早的太阳黑子记录，是载于《汉书·五行志》中的公元前28年三月出现的太阳黑子：“河平元年……三月己未，日出黄，有黑气大如钱，居日中央。”而后在1840年代德国的一位业余天文学家发现了太阳黑子10-11年的周期变化规律。通过长期的观测，人们还发现太阳黑子在日面上的活动随时间变化的纬度分布也有规律性。二、实验方法在该实验中，对Wolfer数序列做FFT变换后得到Y（长度为n），只取其前个数据的功率谱密度的估计值。原因是时域为离散的实序列的傅立叶变换对应于具有周期性且偶对称的

4、频域特性，因此Y的前个数据已经包含了Wolfer数的全部信息。根据DFT的频域单位与DTFT的频域单位的表达关系式以及与f对应关系，可以看出与f呈线性关系；同样地，因为Y的前个数据已经包含了Wolfer数的全部信息，只取前个数据分析功率-频率图时，对应的横坐标时应取。三、流程图四、MATLAB实验内容，结果，分析(1)读取数据并绘制Wolfer图，观察太阳黑子的活动现象。1、程序如下：sunspot=xlsread('C:UsersAdministratorDesktop1月数据.xls');%读取数据

5、year=sunspot(:,1);%读取年份信息wolfer=sunspot(:,3);%读取黑子活动数据figure%新建图像plot(year,wolfer)%画出时域图xlabel('年份');ylabel('黑子个数');title('时域图')%标注横纵坐标和标题pause截图：数据采用的是比利时皇家天文台中自1750至今的一月份的黑子平均数。2、为使结果更加清晰，画出部分时间段(如实验中采取1800-1850年)的Wolfer图进行观察。程序如下：Figureplot(year(51:101),w

6、olfer(51:101),'b.-');%用1800-1850年数据做图xlabel('年份');ylabel('黑子个数');title('1800-1850年')%标注横纵坐标和标题。截图：(1)对已经得到的Wolfer数应用FFT变换分析它的变化规律。在复平面上绘制出傅立叶系数的分布图。程序如下：G=fft(wolfer);%对全部数据做FFTG(1)=[];%由于算法问题，Y(1)为所有fft数值之和，因此需舍弃第一点pause%观察Wolfer数在FFT分析后得到的在复平面中的分布图。程序如下fig

7、ureplot(G,'ro')%在复平面做图，空心点title('FFT复平面图');xlabel('实轴');%标注横坐标为’实部’ylabel('虚轴');%标注纵坐标为’虚部’pause截图：由图可见，黑子频率沿虚轴严格对称，具有周期性。(1)功率与频率的关系曲线则被定义为周期图。用周期图法对Wolfer数进行功率谱密度分析。程序如下：figuren=length(G);%取FFT结果长度为npower=abs(G(1:n/2)).^2;%取前n/2个数据求其功率nyquist=1/2;%取最大频率为0.

8、5freq=(1:n/2)/(n/2)*nyquist;%将FFT转换为DTFTplot(freq,power)%画周期图xlabel('周期/年份');title('周期图')pause截图：为使结果更加清晰，取与1800-1850年相对应的功率和频率的分量（51-101）作它的周期图。程序如下：figureplot(freq(51:101),power(51:101))%为观察方便