2
Mensimulasikan proses Gaussian (Ornstein Uhlenbeck) dengan fungsi kovarians yang membusuk secara eksponensial
Saya mencoba untuk menghasilkan banyak undian (yaitu, realisasi) dari proses Gaussian ei(t)ei(t)e_i(t), 1≤t≤T1≤t≤T1\leq t \leq T dengan rata-rata 0 dan fungsi kovarian γ(s,t)=exp(−|t−s|)γ(s,t)=exp(−|t−s|)\gamma(s,t)=\exp(-|t-s|). Apakah ada cara yang efisien untuk melakukan ini yang tidak akan melibatkan perhitungan akar kuadrat dari a T×TT×TT \times Tmatriks kovarians? Atau adakah yang bisa merekomendasikan Rpaket …