Message Boards Message Boards


Curious about the behaviour of KarhunenLoeveDecomposition

Posted 3 months ago
0 Replies
0 Total Likes

Hello. I have some questions on how to implement propperly the KL expansion in mathematica using KarhunenLoeveDecomposition.

I have constructed the matrix realizationMat. it contains 101 measurements of 90 Gaussian stochastic processes. I implement the expansion as

VA1 = realizationMat[[1]];
VA2 = realizationMat[[2]];
KLVariablesB = KarhunenLoeveDecomposition[{VA1, VA2...}]

However, I have problems understanding the output. In the documentation it says:

rows of the transformation matrix m are the eigenvectors of the covariance matrix formed from the arrays ai.

yet, if I plot rows of m as ListPlot[{KLVariablesB[[2, 1]]}, Joined -> True] I obtain what can be seen in figure 1Figure 1.

Another quote from the documentation is

The transformed arrays bi are uncorrelated, are given in order of decreasing variance, and have the same total variance as ai.

As far as I know, this means that these are the random variables defined in the theory. However, if I plot them, These don't look like random variables, as can be seen in figure 2: enter image description here. These look like eigenvectors but aren't Gaussian eigenvectors, being my data Gaussian.

Can someone please tell me if there is something wrong with my data or why the eigenvectors don't look like eigenvectos and the random variables don't look like random variables?

Best regards. jaime.

E. G. I will attach the code I am using in case someone wants to take a look.

Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
or Discard

Group Abstract Group Abstract