r4 - 31 Mar 2009 - 02:02:18 - YanWeiYou are here: TWiki >  Documentation Web  >  AskedQuestions > ChangeSolverPPPFromStaticToKinematic

Question

Hi,all. I want to process spaceborne gps data(champ grace etc.) with SolverPPP?. As example 8 and example 9 are all process static obs data with a static model, I have to make some changes on SolverPPP?. I make a copy of SolverPPP? and change class name to SolverLEO?.After going through the code, I found lots of work to do,and wish you can give me some help or suggestions. SolverPPP? use a static model state vector (ztrop_wet,px py pz,cdt,[N...]), and I want to use a new kinematic state vector (px,py,pz,vx,vy,vz,cdt,[N...]),[px py px] and [vx vy vz] stand for position and velocity.

1 I try to change the default equation of the solver.*defaultEqDef. If I add [vx vy vz] to enum TypeID?::ValueType{}. Is the follow changes ok?

enum TypeID::ValueType
{
    ...
    vx
    vy
    vz
    Last
    Placeholder = Last+1000
}

*2 change the kalman matrixs

phiMatrix=
{
    Ip3      Icdt3           o           o
    o         Iv3            o           o
    o          o           Ocdt1         o
    o          o             o          Inm
}
qMatrix
{
    Qp3        o            o            o
    o          Qv3          o            o
    o          o   
}

Dear Mr. Wei:

In order to do kinematic PPP-like processing with SolverPPP? you do not have to change SolverPPP?. You just have to change the default stochastic model being used, with instructions like:

WhiteNoiseModel newCoordinatesModel(100.0);
pppSolver.setCoordinatesModel(&newCoordinatesModel);

I understand you want to add new velocity state variables to the filter, so you should modify TypeID? to include them. Just don't forget to change both TypeID?.hpp and TypeID?.cpp in a coherent way.

However, I must warn you that your problem will NOT come from the filter: A LEO moves very fast, so you'll have several issues to deal with:

  • If you want to use precise ephemeris, you must implement a way to also use 30 seconds precise clocks.

  • You must carefully tune the cycle slip objects and arc marker to fit your applications.

  • You'll have a lot of short arcs due to the high satellite speed. This will impair proper convergence of ambiguities.

  • Adding new variables to your filter will only worsen the ambiguities problems.

I suggest you to start with a very simple filter, and then add more variables to your state vector in a very careful way.

I also suggest you to check this thesis written by one of my colleagues:

http://www.tdx.cat/TDX-0312109-123420/

Best regards,

Dago

Dear Dago:

Thanks for your suggestions and doctoral thesis,I'll have a try.

Best regards,

Wei

-- DagobertoSalazar - 30 Mar 2009

Edit | WYSIWYG | Attach | Printable | Raw View | Backlinks: Web, All Webs | History: r4 < r3 < r2 < r1 | More topic actions
 
WEBLOGOALT
This site is powered by the TWiki collaboration platformCopyright amp;© by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback