tworitdash
-
Posts
2 -
Joined
-
Last visited
Content Type
Profiles
Forums
Events
Posts posted by tworitdash
-
-
I want to develop a 2D random field and its change with time with constant velocity. My process:
1. Define a 2D grid (not the fields yet)
[x,y]with
n×npoints
2. Define 1D time axis
[t]with
ntelements
3. Find the Lagrangian distance
Dbetween the points in space and with time [ with advection velocities in
xand
y, i.e.,
[vx,vy]]. For now
vxand
vyare constants.
D[t,(xi,yi),(xj,yj)]=(xi−vxt−xj)2+(yi−vyt−yj)2−−−−−−−−−−−−−−−−−−−−−−−−−−√Dhas size
[nt×n2×n2]4. Use this distance matrix and the time axis to formulate a covariance cube of size
[nt×n2×n2]ρ(D,t;θ)=e−(D/bs)cse−(t/bt)ct1−θe−(D/bs)cse−(t/bt)ctWhere
[bs,cs]are spatial correlation parameters and
[bt,ct]are time correlation parameters. A non-zero
θmakes this spatio-temporal correlation non-separable.
Output Expected:
A 2D Random field and its variation with respect to time that moves with a velocity of
[vx,vy]. I haven't used a 2D random field yet. I defined the grid only and based on Lagrangian distances of every pixel in space and time I defined a covariance matrix (Also includes the velocity information). I want to reproduce the 2D random fields out of this covariance structure.
So, I have the correlation of all spatial points with all spatial points in space $[n^2 \times n^2]$ and also the change in the spatial correlation with time
[nt×n2×n2]I have no clue how to proceed next. I read ideas in papers to make it a Multivariate autoregressive process (MAR) and find coefficients based on this covariance matrix.
The first confusion is if this AR model is for my random field itself with time and covariance is just a tool to solve that AR process or the AR model is for the 2D covariance matrices with time.
The second confusion is that if it is an AR process of the 2D random field itself with time, then, how this covariance matrix will help me solve for the random field? Because the covariance matrix is a
[nt×n2×n2]matrix.
Is there any other way I can solve this? I have the correlation structure already so in my head I somehow believe it should be simple, but when I think more I am confused.
0
2D space and 1D time evolution of a random field
in Applied Mathematics
Posted
There are some typographical mistakes in the post.
The Lagrangian distance:
$$ D[t, (x_i, y_i), (x_j, y_j)] = \sqrt{(x_i - v_xt - x_j)^2 + (y_i - v_yt - y_j)^2} $$
The correlation function:
$$ \rho(D, t; \theta) = \frac{e^{-(D/b_s)^{c_s}} e^{-(t/b_t)^{c_t}}}{1 - \theta e^{-(D/b_s)^{c_s}} e^{-(t/b_t)^{c_t}}} $$