在 MATLAB 中生成异方差级数？

Question

我正在尝试在 MATLAB 中生成异方差误差项。

我找到了一个很好的链接来帮助here，我遇到的问题是在 MATLAB 中复制它。

这是我的尝试：

n = [(1:100) (1:100)]';
sigma2 = n.^(1.3);
t = size(n,1);
for i = 1:200
    eps(i) = normrnd(0, sqrt(sigma2(1)));
end
eps = eps'
h = archtest(eps)

但是，异质测试表明我仍然没有异方差数据，任何人都可以看到我哪里出错了。

Answer 1

看起来您将 sigma2 固定在循环内的第一个值中。

替换

eps(i) = normrnd(0, sqrt(sigma2(1)));

与

eps(i) = normrnd(0, sqrt(sigma2(i)));