Generate AR1 spatial data
This script generates artificial spatial data using a first order spatial autoregressive process (AR1)
The process is
X(i,j) = PHI*(X(i,j-1)+X(i,j+1)+X(i-1,j)+X(i+1,j) + error
To generate a 10x10, periodic spatial pattern, with normal error distribution, and a variance of SIGMA, the command is
spatialPatternAR1([10,10],PHI,SIGMA)
or equivalently,
spatialPatternAR1([10,10],PHI,SIGMA,1)
the same pattern, but with a rigid boundary is generated with
spatialPatternAR1([10,10],PHI,SIGMA,0)
Different error distributions can be used. Using the current script this can be done by extracting the generator matrix, M
M = spatialPatternAR1([10,10],PHI,SIGMA,1,1)
M is a 100x100 matrix.
Then given a random vector, error (size = 100x1) with indepedent elements sampled from the required error distribution, the spatial pattern is given by
x = reshape(M * error,[10,10])
引用格式
Jon Yearsley (2024). Generate AR1 spatial data (https://www.mathworks.com/matlabcentral/fileexchange/5099-generate-ar1-spatial-data), MATLAB Central File Exchange. 检索来源 .
MATLAB 版本兼容性
平台兼容性
Windows macOS Linux类别
标签
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!版本 | 已发布 | 发行说明 | |
---|---|---|---|
1.0.0.0 | Update for BSD License |