極座標の時系列を扱っており、予測にカルマン フィルターを適用しています。時系列は衛星軌道に関連しています。
ただし、分散の予測と推定は極座標 [r,theta] で表されます。
関数を使用して予測をデカルト座標に変換する方法を知っています
f(r,theta) <- [r*cos(theta),r*sin(theta)].
しかし、それは線形演算子ではないため、分散を処理する方法がわかりません。
変換を手伝っていただける場合は、データを順番に提供します。
Radius Angle
[1,] "39805.9613778309" "1.46134492279737"
[2,] "39805.9613778309" "1.48689546833425"
[3,] "39805.9613778309" "1.51244601387112"
[4,] "39805.9613778309" "1.537996559408"
[5,] "39805.9613778309" "1.56354710494488"
[6,] "39805.9613778309" "1.58909765048176"
そして、最初の予測の分散行列は次のとおりです。
radius theta
[1,] 5132782 0.000000000
[2,] 0 0.001646994
最初の予測のためにデカルト座標でこの行列を取得する方法を知りたいです。ありがとう!