[PP,dd,m,r,nf,Xiter] = ESDPDtracking(PP,m,n,r,nf,degree,movingspeed)
    Test SDP-based Sensor Localization and Tracking System
    Input
    PP: 2xn matrix representing all point on 2D
    n : the number of total points
    m : the number of anchor points; the first m columns of PP
    r : the radio range
    nf: noisy factor
    degree: the limit on the number of edges connected to any sensor point, typicallu set it to 5-10
    movingspeed: the moving distance of every sensor per snapshot, typically 1-3% of the diameter of the region
      
    Output
    PP: 2xn matrix representing all point on 2D, the last m columns are anchors
    dd: the uppertriangle distance matrix from PP(:i) to PP(:,j), i<j.
    Xiter: position estimation of the last snapshot for all n-m unknown points made by the method    

**************************************************************************************************

    [PP,dd,m,r,nf, Xiter] = ESDP(PP,m,n,r,nf,degree)
    or
    [PP,dd,m,r,nf, Xiter] = ESDPDual(PP,m,n,r,nf,degree)
    Test SDP-based Sensor Localization System
    Input
    PP: 2xn matrix representing all point on 2D
    n : the number of total points
    m : the number of anchor points; the first m columns of PP
    r : the radio range
    nf: noisy factor
    degree: the limit on the number of edges connected to any sensor point, typicallu set it to 5-10

    Output
    PP: 2xn matrix representing all point on 2D, the last m columns are anchors
    dd: the uppertriangle distance matrix from PP(:i) to PP(:,j), i<j.
    Xiter: position estimation for all n-m unknown points made by the method

**************************************************************************************************

    [Y_opt,Xiter]=FULLSDPD(PP,m,dd);

    Xiter = ESDPD(PP,m,dd)

    Xiter = ESDPDualD(PP,m,dd)

    Sensor Localization Solvers
    Input 
      PP  : 2xn matrix representing all point true locations on 2D
      where n is the number of total sensor and anchor points
      m   : the number of Anchor points; last m columns of PP; it should be at least 3
      dd : the sparse (and noisy) distance matrix between PP(:,i) and PP(:,j)
    Output
      Xiter: position estimation for all n-m unknown points made by the method
      Y_opt: SDP matrix representing the second moment of estimation 
    
    Developed by Zizhuo Wang, Song Zheng, Boyd and Yinyu Ye for Demonstration;
    and the steepest descent code is due to Biswas, Toh and Ye;
    see paper http://www.stanford.edu/~yyye/relaxationsdp6.pdf
    Edge based decompositionThree test SDP relaxaion models for sensor network localization in 2-D: