資源簡介
太陽黑子相對數簡稱為太陽黑子數,它反映太陽活動強弱的變化,對地球環境和人類活動的影響很大。本文從比利時皇家天文臺(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(‘零均值后太陽黑子時序圖像‘);
評論
共有 條評論