【发布时间】:2016-11-28 11:15:42
【问题描述】:
以下是我老师给我的问题,
- 生成一个包含 N = 1000 个随机变量的独立观测值的序列,分布:(c) 参数 λ = 1 的指数,通过 倒置法。
- 呈现以图形方式获得的序列(e 点生成的序列除外),例如(一) 一。在坐标中绘图(编号 obs., obs 的值) ii.在坐标中绘图 (obs no n, obs. no n + i) 对于 i = 1, 2, 3。 iii.绘制一些所谓的协方差函数 价值观。即和平均值:
我写了如下代码,
(*****************************************************************)
(*Task 01(c) and 02(a)*)
(*****************************************************************)
n = 1000;
taskC = Table[-Log[RandomReal[]], {n}];
ListPlot[taskC, AxesLabel->{"No. obs", "value of the obs"}]
i = 1;
ListPlot[Table[
{taskC[[k]], taskC[[k+i]]},
{k, 1, n-i,1}],
AxesLabel->{"obs.no.n", "obs.no.n+1"}]
i++;
ListPlot[Table[
{taskC[[k]], taskC[[k+i]]},
{k, 1, n-i,1}],
AxesLabel-> {"obs.no.n", "obs.no.n+2"}]
i++;
ListPlot[Table[
{taskC[[k]], taskC[[k+i]]},
{k,1,n-i,1}],
AxesLabel->{"obs.no.n", "obs.no.n+3"}]
avg = (1/n)*Sum[taskC[[i]], {i,n}];
ListPlot[Table[1/(n-tau) * Sum[(taskC[[i]]-avg)*(taskC[[i+tau]] - avg), n], {tau, 1,100}],
Joined->True,
AxesLabel->"Covariance Function"]
他发表了评论,
协方差函数的图应该从 0-shift 开始。笔记 对于大于 0 的班次,您正在估计之间的协方差 独立观察值为零,而对于 0 班次,您是 估计观察的方差很大。因此对比 这两种情况之间的差异清楚地表明,观察 是不相关的。
我做错了什么?
如何更正我的代码?
【问题讨论】:
标签: random wolfram-mathematica stochastic-process