太阳黑子的AR（n）模型预测

太阳黑子相对数简称为太阳黑子数，它反映太阳活动强弱的变化，对地球环境和人类活动的影响很大。本文从比利时皇家天文台（the Royal Observatory of Belgium）的太阳黑子指数数据中心（the Sunspot Index Data center）的网站获得了1700—2014 每年的太阳黑子数的数据。利用MATLAB软件结合时间序列建模方法对观测值进行了分析和建模，并利用该模型对未来的太阳黑子数进行了预测，模型的预测效果比较理想。

clear;
clc;
load sunspot_1.dat;
x=sunspot_1（:2）;
n=length（x）;
t=linspace（17002014315）;
mn=mean（x）;  
x1=x-mean（x）;
[r1flags1]=xcorr（x）;
[r2flags2]=xcorr（x1）;

%检验数据平稳性和正态性，如下：
if adftest（x）         
    fprintf（‘平稳\n‘）
else
    fprintf（‘非平稳\n‘）
end

if kstest（x）
    fprintf（‘正态\n‘）
else
    fprintf（‘非正态\n‘）
end

figure（1）;
plot（tx‘color‘‘g‘‘marker‘‘none‘）;
axis（[170020140200]）;grid on;   %标注横纵坐标和标题
xlabel（‘年份Time‘）;ylabel（‘黑子数量‘）;title（‘1700-2014年太阳黑子时序图像‘）;
hold on;
line（[1700 2014][mn mn]‘Linestyle‘‘-.‘‘color‘‘r‘‘LineWidth‘2）;


figure（2）;
plot（tx1‘color‘‘g‘‘marker‘‘none‘）;
line（[1700 2014][0 0]‘Linestyle‘‘-.‘‘color‘‘r‘‘LineWidth‘2）;
axis（[17002014-50150]）;grid on;
xlabel（‘年份‘）;ylabel（‘黑子数量‘）;title（‘零均值后太阳黑子时序图像‘）;