主要内容

使用并行计算实施jackknife

这个例子来自钉书匠功能参考页,但并行运行。

从具有平均0和标准偏差5的正态分布生成大小10000的样本数据。

西格玛= 5;RNG('默认')y = normrnd(0,sigma,10000,1);

跑步钉书匠并行以估计方差。为此,使用实例化创建选项结构并设置使用指α.现场为真。

opts = statset('使用指平行',真的);m = jackknife(@ var,y,1,'选项',选择);

将已知的偏见配方与千刀偏见估计进行比较。

n =长度(y);bias = -sigma ^ 2 / n%已知的偏见公式JBIAS =(n-1)*(平均值(m)-var(y,1))%jackknife偏见估计
使用“本地”配置文件启动并行池(Parpool)连接到并行池(工人数:6)。bias = -0.0025 jbias = -0.0025

比较串行和并行计算需要多长时间。

Tic; m = jackknife(@ var,y,1); toc%串行计算
经过时间为1.638026秒。
tic; m = jackknife(@ var,y,1,'选项',选择); TOC%并行计算
经过时间为0.507961秒。

钉书匠不使用随机数,因此每次提供相同的结果,无论是并行还是串行运行。