Jack O'Brien edited untitled.tex  over 7 years ago

Commit id: 42f5eaab3e17f151a3442001284907f990c1aab8

deletions | additions      

       

  Where $G_{i}$ is the {\em observation matrix} which in our case will simply represent whether or not the value of the light curve was observed, and $W_{i}$ is a random zero mean vector following a set of known covariance matricies, $R_{i}$, representing the observation error.  The operation of the Kalman filter can be broken down into a prediction step and an update step. The prediction step involves finding the best linear predictor of our system, $\^{X}$, $\hat{X_{i}}$,  which is where we think our system should be according to our CARMA model (how is this done?). The update step compares our prediction to our true value and attempts to correct for the differences in order to find the best linear estimator, $X^{\dagger}$ $X^{\dagger}_{i}$  for the state of the system (how is this done). (talk about how the likelihood is then calculated) Estimation of these vectors, and their corresponding covariance matricies matricies, $\Sigma_{i}$,$\Sigma^{\dagger}_{i}$  is accomplished through the following recursive algorithm: %AGN light curves only give us a picture of the effects that the underlying physical processes have on the AGN itself, but they don't directly tell us what those physical processes are or how they evolve. Futher more, observation error makes it even more difficult to determine their structure. To understand these processes, we need a method of inspecting the state of the system based on our external observations that can also take into account our observation error. The Kalman filter is a digital filter that attempts to reduce observation error and intrinsic system noise by making predictions about future values of the time-series based on an assumed model using Markov chains. Using a Kalman filter together with our CARMA model, we can accurately model theses underlying physical processes and begin analyzing the true structure of the driving system. In order to effectively use a Kalman filter, we must transform our observations, $y(t)$, into state space, $x(t)$.