AbsoluteCorrelationFunction[data,hspec]
从 data 估计延迟 hspec 处的绝对相关函数.
AbsoluteCorrelationFunction[proc,hspec]
表示对于随机过程 proc 延迟 hspec 处的绝对相关函数.
AbsoluteCorrelationFunction[proc,s,t]
表示随机过程 proc 在时间 s 和 t 处的绝对相关函数.
AbsoluteCorrelationFunction
AbsoluteCorrelationFunction[data,hspec]
从 data 估计延迟 hspec 处的绝对相关函数.
AbsoluteCorrelationFunction[proc,hspec]
表示对于随机过程 proc 延迟 hspec 处的绝对相关函数.
AbsoluteCorrelationFunction[proc,s,t]
表示随机过程 proc 在时间 s 和 t 处的绝对相关函数.
更多信息
- AbsoluteCorrelationFunction 也称为自相关函数.
- 对于 hspec 可以给出下列指定:
-
τ 在时间或者延迟 τ {τmax} 单位间隔从 0 到 τmax {τmin,τmax} 单位间隔从 τmin 到 τmax {τmin,τmax,dτ} 从 τmin 到 τmax,步长为 dτ {{τ1,τ2,…}} 使用显式 {τ1,τ2,…} - AbsoluteCorrelationFunction[{x1,…,xn},h] 等价于
. - 当 data 是包含路径集合的 TemporalData 时,输出表示所有路径的平均.
- AbsoluteCorrelationFunction(其中的过程 proc 在时间 t 处的值为 x[t])由下面给出:
-
Expectation[x[s] x[t]] 对于标量值过程 Expectation[x[s]⊗x[t]] 对于向量值过程 - 符号 ⊗ 表示 KroneckerProduct.
- AbsoluteCorrelationFunction[proc,h] 只有当 proc 是弱平稳过程,并且等价于 AbsoluteCorrelationFunction[proc,0,h] 时才被定义.
- 过程 proc 可以是任意随机过程,例如 ARMAProcess 和 WienerProcess.
范例
打开所有单元 关闭所有单元基本范例 (4)
AbsoluteCorrelationFunction[{2, 3, 4, 3}, 2]data = RandomFunction[ARProcess[{.2, .3, .4}, 1], {1, 10 ^ 4}];ListPlot[AbsoluteCorrelationFunction[data, {30}], Filling -> Axis, PlotRange -> All]AbsoluteCorrelationFunction[BinomialProcess[p], s, t]DiscretePlot3D[Evaluate[% /. {p -> 2 / 3}], {s, 0, 10}, {t, 0, 10}, ExtentSize -> 1 / 2, ColorFunction -> "Rainbow"]AbsoluteCorrelationFunction[WienerProcess[μ, σ], s, t]Plot3D[Evaluate[% /. {μ -> 1, σ -> 2}], {s, 0, 10}, {t, 0, 10}, ColorFunction -> "Rainbow"]范围 (13)
经验估值 (7)
AbsoluteCorrelationFunction[Range[10], 5]AbsoluteCorrelationFunction[Range[10], {9}]AbsoluteCorrelationFunction[Range[10], {1, 9, 2}]ts = TemporalData[TimeSeries, {{{-2.9904395734335427, -1.5547517217370468, -2.007044987914736,
-0.7114933373063236, -1.5920036118254968, 0.23183536513859127, 0.056774681349179124,
0.32394041932443046, 0.10382798659711348, -0.684046765274114, - ... 1687325555,
3.410933732368106, 2.5200924558212163, 2.008558601261935, 0.43151331467691206,
0.6249727955196763}}, {{1, 1000, 1}}, 1, {"Continuous", 1}, {"Discrete", 1}, 1,
{ValueDimensions -> 1, ResamplingMethod -> None}}, False, 10.1];abscorr = AbsoluteCorrelationFunction[ts, {100}]ListPlot[abscorr, Filling -> Axis]data = RandomFunction[ARProcess[{.8}, 1], {0, 500}, 10]abscorr = AbsoluteCorrelationFunction[data, {12}];ListPlot[abscorr, Filling -> 0]proc = MAProcess[{.4, .3, .5, .6, .3}, 1.];data = RandomFunction[proc, {0, 100}];ListPlot[AbsoluteCorrelationFunction[#, {10}]& /@ {proc, data}, Filling -> {1 -> {2}}, PlotStyle -> PointSize[Medium], PlotLegends -> {"process absolute correlation", "sample absolute correlation"}]proc = ARProcess[{{{.3, .1}, {.6, .3}}}, {{1, .3}, {.3, .6}}];
data = RandomFunction[proc, {100}];
cov = AbsoluteCorrelationFunction[data, {-6, 6}];Table[ListPlot[cov["Values"][[All, i, j]], Filling -> 0, PlotLabel -> {i, j}], {i, 1, 2}, {j, 1, 2}]随机过程 (6)
AbsoluteCorrelationFunction[ARProcess[{a}, σ^2], s, t]DiscretePlot3D[Evaluate[% /. {a -> 2 / 3, σ -> 1}], {s, 0, 10}, {t, 0, 10}, ExtentSize -> 1 / 2, ColorFunction -> "Rainbow"]AbsoluteCorrelationFunction[ARProcess[{a}, σ^2], h]DiscretePlot[Evaluate[% /. {a -> 2 / 3, σ -> 1}], {h, 0, 10}, ExtentSize -> 1 / 2, ColorFunction -> "Rainbow", AxesLabel -> {h, None}]AbsoluteCorrelationFunction[OrnsteinUhlenbeckProcess[μ, σ, θ], s, t]Plot3D[Evaluate[% /. {μ -> 1, σ -> 2, θ -> 1 / 3}], {s, 0, 10}, {t, 0, 10}, ColorFunction -> "Rainbow"]AbsoluteCorrelationFunction[OrnsteinUhlenbeckProcess[μ, σ, θ], h]Plot[Evaluate[% /. {μ -> 1, σ -> 2, θ -> 1 / 3}], {h, 0, 10}, ColorFunction -> "Rainbow", Filling -> Axis, AxesLabel -> {h, None}]AbsoluteCorrelationFunction[BinomialProcess[p], s, t]DiscretePlot3D[Evaluate[% /. p -> 1 / 3], {s, 0, 10}, {t, 0, 10}, ExtentSize -> 1 / 2, ColorFunction -> "Rainbow"]antidiagonals = AbsoluteCorrelationFunction[BinomialProcess[1 / 3], t, t + h]DiscretePlot[Evaluate@Table[antidiagonals, {t, 1, 4}], {h, 0, 6}, ExtentSize -> 1 / 2, AxesLabel -> {h, None}, PlotLegends -> (Row[{"t =", #}]& /@ Range[1, 4])]AbsoluteCorrelationFunction[PoissonProcess[μ], s, t]Plot3D[Evaluate[% /. μ -> 3], {s, 0, 10}, {t, 0, 10}, ColorFunction -> "Rainbow"]antidiagonals = AbsoluteCorrelationFunction[PoissonProcess[3], t, t + h]Plot[Evaluate@Table[antidiagonals, {t, 1, 4}], {h, 0, 6}, AxesLabel -> {h, None}, PlotLegends -> (Row[{"t =", #}]& /@ Range[1, 4])]proc1 = ARProcess[{.8}, 1];
proc2 = MAProcess[{-.4, -.7, .6, -1}, 1];
proc3 = SARMAProcess[{}, {}, {3, {.6}, {.1}}, 1];
proc4 = FARIMAProcess[{}, 1 / 3, {}, 1];DiscretePlot[AbsoluteCorrelationFunction[#, h], {h, 0, 15}, ExtentSize -> 1 / 2, PlotRange -> All, PlotLabel -> Head[#]]& /@ {proc1, proc2, proc3, proc4}向量 ARProcess 的绝对互相关图线:
proc = ARProcess[{{{.3, .1}, {.6, .3}}}, {{1, .3}, {.3, .6}}];
cov = AbsoluteCorrelationFunction[proc, h];Table[DiscretePlot[cov[[i, j]], {h, -6, 6}, ExtentSize -> 1 / 2, PlotLabel -> {i, j}], {i, 1, 2}, {j, 1, 2}]应用 (2)
判断下列数据是否可以使用 MAProcess 或者 ARProcess 进行最好的模拟:
ListLinePlot[data = TemporalData[TimeSeries, {{{-1.4735257611469446, -0.9158187257309729, -1.9060205147226934,
-1.6174582169404574, -2.633459084244463, -2.104839842098871, -2.8701276893108894,
-2.4118183939601607, -0.5125633448622677, -1.2356329399850385, 1. ... 26394036266876, -0.16380797769527589, 0.3559909972777121,
1.5343047104249412, 0.10763337590714506}}, {{0, 100, 1}}, 1, {"Continuous", 1},
{"Discrete", 1}, 1, {ResamplingMethod -> {"Interpolation", InterpolationOrder -> 1}}}, False,
10.1]]candidates = {MAProcess[{.5, .4}, 1], ARProcess[{.5, .4}, 1]};SeedRandom[2];
ListLinePlot[RandomFunction[#, {0, 100}]]& /@ candidatesListPlot[AbsoluteCorrelationFunction[data, {10}], Filling -> 0]ARProcess 明显是比 MAProcess 更好的候选模型:
ListPlot[AbsoluteCorrelationFunction[#, {10}], PlotLabel -> Head[#], Filling -> 0]& /@ candidatesproc = OrnsteinUhlenbeckProcess[0, σ, θ];WeakStationarity[proc]abscorr[h_] = AbsoluteCorrelationFunction[proc, h]int[t_] = 1 / (2 * t)Integrate[(1 - Abs[u] / (2 * t))abscorr[u], {u, -2 * t, 2 * t}, Assumptions -> t > 0]Limit[int[t], t -> ∞, Assumptions -> σ > 0 && θ > 0]属性和关系 (13)
proc = ARMAProcess[1, {5 / 6, -1 / 6}, {2 / 3}, 1];
paths = RandomFunction[proc, {0, n = 10}, 10 ^ 4];AbsoluteCorrelationFunction[paths, {10}]sacf = TimeSeriesMapThread[(#2(n + 1) / (n + 1 - #1))&, %]acf = AbsoluteCorrelationFunction[proc, {10}]ListPlot[{sacf, acf}, Filling -> 0, PlotLegends -> {"sample absolute correlation", "process absolute correlation"}]列表的绝对相关函数可以利用 AbsoluteCorrelation 计算:
data = Range[10];
n = Length[data];AbsoluteCorrelationFunction[data, {n - 1}] / (Range[n, 1, -1] / n)Table[ AbsoluteCorrelation[Drop[data, i], Drop[RotateRight[data, i], i]], {i, 0, n - 1}]% - %%AbsoluteCorrelationFunction 是绝对相关矩阵的非对角线项:
𝒫 = WienerProcess[μ, σ];AbsoluteCorrelation[𝒫[{s, t}], 1, 2]AbsoluteCorrelationFunction[𝒫, s, t]Simplify[%% - %, 0 < s < t]延迟 0 处的样本绝对相关函数估计第二个 Moment:
data = Range[10];
AbsoluteCorrelationFunction[data, 0]Moment[data, 2]% - %%样本绝对相关函数与 CovarianceFunction 有关:
data = Range[10];
n = Length[data];AbsoluteCorrelationFunction[data - Mean[data], {n - 1}]CovarianceFunction[data, {n - 1}]% - %%样本绝对相关函数与 CorrelationFunction 有关:
data = Range[10];
n = Length[data];AbsoluteCorrelationFunction[data - Mean[data], {n - 1}]% / First[%]CorrelationFunction[data, {n - 1}]% - %%使用 Expectation 计算绝对相关函数:
𝒫 = WienerProcess[μ, σ];Expectation[x[s] x[t], x𝒫, Assumptions -> s > 0 && t > 0]AbsoluteCorrelationFunction[𝒫, s, t]Simplify[% - %%, s < t]绝对相关函数
与 Moment 函数相关:
𝒫 = GeometricBrownianMotionProcess[μ, σ, θ];AbsoluteCorrelationFunction[𝒫, t, t]Moment[𝒫[t], 2]% - %%//Simplify绝对相关函数
与 CovarianceFunction
相关:
𝒫 = WienerProcess[μ, σ];CovarianceFunction[𝒫, s, t]AbsoluteCorrelationFunction[𝒫, s, t] - Mean[𝒫[s]] Mean[𝒫[t]]% - %%当过程均值为零,绝对相关函数等于 CovarianceFunction:
AbsoluteCorrelationFunction[WienerProcess[], s, t]CovarianceFunction[WienerProcess[], s, t]% - %%绝对相关函数对于 ToInvertibleTimeSeries 是不变的:
proc = MAProcess[{2, 4}, 1];TimeSeriesInvertibility[proc]iproc = ToInvertibleTimeSeries@procAbsoluteCorrelationFunction[#, s, t]& /@ {proc, iproc}data = TemporalData[EventSeries, {{{-0.5997948200136711, -0.6119370022135471, -2.181570597470531,
-3.0852015062413294, -2.5597554304014736, -2.724339030053185, -2.1990400702572104,
-2.2335014895051946, -1.0488409388108924, -0.1005536745721749, 2 ... 63267114671374, 22.486986202712714, 21.766625417091756,
22.44341047189002, 22.079788501970878, 23.102713685532592}}, {{0, 100, 1}}, 1,
{"Continuous", 1}, {"Discrete", 1}, 1, {ResamplingMethod -> None, ValueDimensions -> 1}}, False,
10.1];Mean[data]newdata = Standardize[data, Mean, 1&];ListPlot[AbsoluteCorrelationFunction[#, {20}], Filling -> 0]& /@ {data, newdata}PowerSpectralDensity 对于零均值过程是绝对相关函数的变换:
proc = ARProcess[{a}, σ^2];
cf = AbsoluteCorrelationFunction[proc, h]在合适的参数下使用 FourierSequenceTransform:
ft = FourierSequenceTransform[cf, h, z]PowerSpectralDensity[proc, z]FullSimplify[% - ft]可能存在的问题 (1)
AbsoluteCorrelationFunction 输出可能包含 DifferenceRoot:
AbsoluteCorrelationFunction[ARProcess[{a, b}, σ^2], t] //Simplify[#, t > 0]&使用 FunctionExpand 恢复显式幂:
%//FunctionExpand文本
Wolfram Research (2012),AbsoluteCorrelationFunction,Wolfram 语言函数,https://reference.wolfram.com/language/ref/AbsoluteCorrelationFunction.html.
CMS
Wolfram 语言. 2012. "AbsoluteCorrelationFunction." Wolfram 语言与系统参考资料中心. Wolfram Research. https://reference.wolfram.com/language/ref/AbsoluteCorrelationFunction.html.
APA
Wolfram 语言. (2012). AbsoluteCorrelationFunction. Wolfram 语言与系统参考资料中心. 追溯自 https://reference.wolfram.com/language/ref/AbsoluteCorrelationFunction.html 年
BibTeX
@misc{reference.wolfram_2026_absolutecorrelationfunction, author="Wolfram Research", title="{AbsoluteCorrelationFunction}", year="2012", howpublished="\url{https://reference.wolfram.com/language/ref/AbsoluteCorrelationFunction.html}", note=[Accessed: 16-August-2026]}
BibLaTeX
@online{reference.wolfram_2026_absolutecorrelationfunction, organization={Wolfram Research}, title={AbsoluteCorrelationFunction}, year={2012}, url={https://reference.wolfram.com/language/ref/AbsoluteCorrelationFunction.html}, note=[Accessed: 16-August-2026]}