重复观测的结果
y不规则的时间序列表达的
k未知参数
x和一个时间参数
t作为
y
我
=
f
(
x
1
,
x
2
,
…
,
x
k
,
t
)
。这个函数是由期刊(正弦)组件的数量,可以叠加单向趋势变化。此后,连接方程可以写成
y
t
=
b
0
+
b
1
t
+
∑
j
=
1
n
(
米
j
+
一个
j
罪
(
ω
j
t
+
ϕ
j
)
)
,在哪里
b
0的值是线性趋势函数在初始时刻,
b
1是一阶的线性系数(线速度),
米
j是振荡的纵轴的
jth谐波,
一个
j,
ω
j,
φ
j表示振幅、角频率和角的阶段
j分别th谐波。
在缺乏修道院信息数量和贡献水平的期刊谐波在整个角色的变化,一步一步迭代过程似乎是有效的分析。
在第一阶段,人们很自然地确定趋势和删除组件使用方程从系列
y
=
b
0
+
b
1
t
。系统将记者校正方程
v
l
=
(
1
t
1
1
t
2
⋯
⋯
1
t
n
]
⋅
(
b
0
b
1
]
- - - - - -
(
- - - - - -
y
1
- - - - - -
y
2
⋯
- - - - - -
y
n
]
=
一个
l
x
l
+
l
l
,在哪里
v
l代表一个剩余随机噪声项,
l指数反映了一个线性变化的一部分
n值的数量在时间序列分析。
这个任务的条件下可以解决
v
l
T
P
v
l
=
最小值
,在哪里
P是一个
n
×
n权重矩阵。
常见的情况,我们可以有一个女修道院逆权重矩阵
问
l
=
P
- - - - - -
1
=
(
问
11
问
12
⋯
问
1
n
问
21
问
22
⋯
问
2
n
⋯
⋯
⋯
问
n
1
问
n
2
⋯
问
n
n
]
。的最小二乘解(
4)将得到
x
l
=
- - - - - -
(
一个
l
T
P
一个
l
)
- - - - - -
1
一个
l
T
P
l
l
。
第二步,把线性建模组件从初始时间序列,减去(
3)(
2),我们将搜索主导谐波组件。
剩余时间序列的模型将会出现
y
=
∑
j
=
1
k
(
米
j
+
一个
j
罪
(
ω
j
t
+
ϕ
j
)
)
。其未知参数将因此寻找每一个谐波分量
y
=
米
j
+
一个
j
罪
(
ω
j
t
+
ϕ
j
)
。分配由粗糙的模型参数
x
0
=
(
米
0
,
一个
0
,
ω
0
,
φ
0
]
T谐波(
9),由第一次用泰勒级数线性化后,我们将得到一个校正方程矩阵
B:
v
=
B
x
+
l
,并代表
B
0
=
(
1
罪
(
ω
o
t
1
+
ϕ
0
)
一个
o
t
1
因为
(
ω
0
t
1
+
ϕ
0
)
一个
o
因为
(
ω
0
t
1
+
ϕ
0
)
1
罪
(
ω
o
t
2
+
ϕ
0
)
一个
o
t
2
因为
(
ω
0
t
2
+
ϕ
0
)
一个
o
因为
(
ω
0
t
2
+
ϕ
0
)
⋯
⋯
⋯
⋯
1
罪
(
ω
o
t
n
+
ϕ
0
)
一个
o
t
n
因为
(
ω
0
t
n
+
ϕ
0
)
一个
o
因为
(
ω
0
t
n
+
ϕ
0
)
]
。免费的校正方程将收到
l
0
=
(
y
1
- - - - - -
(
米
0
+
一个
0
罪
(
ω
0
t
1
+
ϕ
0
)
)
y
2
- - - - - -
(
米
0
+
一个
0
罪
(
ω
0
t
2
+
ϕ
0
)
)
⋯
y
我
- - - - - -
(
米
0
+
一个
0
罪
(
ω
0
t
我
+
ϕ
0
)
)
]
。作为初始模型(
9)是非线性的,我们将寻找未知参数
x数值迭代过程。
迭代过程将被写成
x
r
+
1
=
x
r
+
f
(
x
r
)
,
(
r
=
0
,
1、2
,
…
)
,在哪里
r是一个迭代的数字。
向量函数的估计
f
(
x
r
)将会由
f
(
x
r
)
=
- - - - - -
(
B
r
T
P
B
r
)
- - - - - -
1
B
r
T
P
l
r
。
然后,迭代过程(
14)将表单
x
r
+
1
=
x
r
- - - - - -
(
B
r
T
P
B
r
)
- - - - - -
1
B
r
T
P
l
r
。签署一些忽视的区别
Δ
x目前,我们继续迭代时(
x
r
+
1
- - - - - -
x
r)将小于或等于
Δ
x。
因此,最后一次迭代的结果将被接受为最终估计的模型(
9)参数
x
=
(
米
一个
ω
ϕ
]
T
。为了估计的统计信心导致特性,我们使用单位重量的评估标准
μ
=
v
T
P
v
n
- - - - - -
4
,在哪里
v是一个
n
×
1残差向量(
11)计算使用(
17)。
实验标准
年代谐波参数估计会决定
(
年代
米
年代
一个
年代
ω
年代
ϕ
]
=
μ
问
我
我
=
μ
(
B
r
T
P
B
r
)
我
我
- - - - - -
1
,在哪里
问
我
我
=
(
B
r
T
P
B
r
)
我
我
- - - - - -
1是倒生的权重矩阵的对角元素的参数。
然后,我们将计算归一化参数估计
Z
=
(
米
/
年代
米
一个
/
年代
一个
ω
/
年代
ω
ϕ
/
年代
ϕ
]
。分配一个置信界限
t
≥
Z
问
%我们确定收到了估计的置信水平。根据解决任务,信心的限制
t被分配使用正态分布。在大多数的研究中,我们可以考虑解决参数统计自信的时候
Z
≥
2。
三种类型的太阳能数据用于运动学建模前预测(
9)更新。所有的新预测表明,24日周期最大阶段可能会发生在2012年或稍后的中间。预测时刻的最大的月度平均值下降到2012年7月- 9月。这个区间的不确定性
+1 - 2个月。预计平均振幅一般周期的形状接近23日最大的从数据
7,
5,
6。SC24将持续不长。它可能会短于10年。下一个最小的开始阶段是初步预测在2017年。太阳黑子数的最大月度平均值Ri被预测为130±20。超过2000±350单位的数以百万计的半球最大月度意味着太阳黑子面积指数的预测。预计每月10.7厘米无线电通量最大价值的200±20学院。宇宙射线的变化在莫斯科站也是建模作为额外的数据来源。其逆转级处理的水平略低于上一个太阳周期的形状。宇宙射线下最少是一年左右转向未来预测太阳能数据最大的引用。 The last estimation of the shift between galactic cosmic ray and solar activity variations was 0.6 years [
15]。