lambda.ps

Lamda算法：GPS快速整周模糊度算法的核心

lineto closepath clip}if/showpage{}N/erasepage{}N/copypage{}N newpath}N
/@endspecial{count ocount sub{pop}repeat countdictstack dcount sub{end}
repeat grestore SpecialSave restore end}N/@defspecial{SDict begin}N
/@fedspecial{end}B/li{lineto}B/rl{rlineto}B/rc{rcurveto}B/np{/SaveX
currentpoint/SaveY X N 1 setlinecap newpath}N/st{stroke SaveX SaveY
moveto}N/fil{fill SaveX SaveY moveto}N/ellipse{/endangle X/startangle X
/yrad X/xrad X/savematrix matrix currentmatrix N TR xrad yrad scale 0 0
1 startangle endangle arc savematrix setmatrix}N end

%%EndProcSet
TeXDict begin 39158280 55380996 1000 600 600 (lambda.dvi)
@start /Fa 171[41 6[58 77[{TeXBase1Encoding ReEncodeFont}2
61.4362 /Times-Bold rf
%DVIPSBitmapFont: Fc cmsy10 6 1
/Fc 1 107 df<12C0B3B3AD02317AA40E>106 D E
%EndDVIPSBitmapFont
/Fd 139[14 1[19 3[25 4[14 25 3[25 26[17 73[{
TeXBase1Encoding ReEncodeFont}7 49.8132 /Times-Italic
rf
%DVIPSBitmapFont: Fe cmex10 10 1
/Fe 1 113 df<1B301B781BF8A2F201F0A2F203E0A2F207C0A2F20F80A2F21F00A21A3E
A262A262A24F5AA24F5AA24F5AA262190FA24FC7FCA2193EA261A261A24E5AA24E5AA24E
5AA24E5AA24EC8FCA2183EA260131001305E13F800014C5A1203D80FFC4B5A121DD838FE
4B5A12F0D8407F4B5A12004DC9FC6D7E173E6D7E5F6D7E5FA26D6C495AA26D6C495AA26D
6C5C1607A26D6C495AA2027F49CAFCA291383F803EA25EEC1FC05EEC0FE0EDE1F0EC07F1
EDF3E0A26EB45AA26E5BA26E90CBFCA25D157E157C15384D64788353>112
D E
%EndDVIPSBitmapFont
/Fg 133[29 33 33 48 33 33 18 26 22 33 33 33 33 52 18
33 1[18 33 33 22 29 33 29 33 29 11[48 41 37 44 2[48 48
59 41 2[22 48 48 37 41 48 44 44 48 12[33 33 4[17 22 17
44[{TeXBase1Encoding ReEncodeFont}47 66.4176 /Times-Roman
rf /Fh 171[30 6[44 26[25 25 49[{TeXBase1Encoding ReEncodeFont}4
49.8132 /Times-Roman rf /Fi 134[42 2[42 46 28 32 37 46
46 42 46 69 23 2[23 46 42 28 37 46 37 46 42 11[60 55
46 60 4[78 55 2[32 5[60 1[60 7[42 42 42 42 42 42 42 42
42 42 23 47[{TeXBase1Encoding ReEncodeFont}41 83.022
/Times-Bold rf /Fj 171[61 6[86 77[{TeXBase1Encoding ReEncodeFont}2
90.9091 /Times-Bold rf
%DVIPSBitmapFont: Fk cmsy10 7.39999 3
/Fk 3 107 df<B712F8A325037B9030>0 D<133013381330A400E0131C00F0133C007C13
F8381E31E0380733803801FE00EA0078A2EA01FE38073380381E31E0387C30F800F0133C
00E0131C00001300A413381330161A7C9C1F>3 D<12E0B3B3B3A7033D79AD11>106
D E
%EndDVIPSBitmapFont
%DVIPSBitmapFont: Fl cmmi10 7.39999 1
/Fl 1 60 df<1238127C12FEA212FF127F123B1203A41206A2120CA21218123812701220
08137B8611>59 D E
%EndDVIPSBitmapFont
/Fm 133[24 5[17 5[31 4[17 6[31 31 18[41 4[20 6[38 66[{
TeXBase1Encoding ReEncodeFont}9 61.4362 /Times-Italic
rf
%DVIPSBitmapFont: Fn cmr10 7.39999 1
/Fn 1 62 df<B81280A3CBFCA9B81280A3290F7D9630>61 D E
%EndDVIPSBitmapFont
%DVIPSBitmapFont: Fp cmr10 10 6
/Fp 6 94 df<146014E0EB01C0EB0380EB0700130E131E5B5BA25B485AA2485AA212075B
120F90C7FCA25A121EA2123EA35AA65AB2127CA67EA3121EA2121F7EA27F12077F1203A2
6C7EA26C7E1378A27F7F130E7FEB0380EB01C0EB00E01460135278BD20>40
D<12C07E12707E7E7E120F6C7E6C7EA26C7E6C7EA21378A2137C133C133E131EA2131F7F
A21480A3EB07C0A6EB03E0B2EB07C0A6EB0F80A31400A25B131EA2133E133C137C1378A2
5BA2485A485AA2485A48C7FC120E5A5A5A5A5A13527CBD20>I<15301578B3A6007FB812
F8B912FCA26C17F8C80078C8FCB3A6153036367BAF41>43 D<007FB812F8B912FCA26C17
F8CCFCAE007FB812F8B912FCA26C17F836167B9F41>61 D<EAFFF8A4EAF000B3B3B3B3A3
EAFFF8A40D5378BD17>91 D<EAFFF8A4EA0078B3B3B3B3A3EAFFF8A40D537FBD17>93
D E
%EndDVIPSBitmapFont
%DVIPSBitmapFont: Fq cmsy10 10 8
/Fq 8 108 df<007FB81280B912C0A26C17803204799641>0 D<0060150600F8150F6C15
1F007E153F6C157E6C6C14FC6C6CEB01F86C6CEB03F06C6CEB07E06C6CEB0FC06C6CEB1F
80017EEB3F006D137E6D6C5A90380FC1F8903807E3F0903803F7E06DB45A6D5B6EC7FCA2
4A7E497F903803F7E0903807E3F090380FC1F890381F80FC90383F007E017E7F49EB1F80
4848EB0FC04848EB07E04848EB03F04848EB01F84848EB00FC48C8127E007E153F48151F
48150F00601506282874A841>2 D<EB0380497EA7397803803C00FC147E00FE14FE397F
8383FC393FC387F8390FE38FE03903FBBF803900FFFE00EB3FF8EB0FE0A2EB3FF8EBFFFE
3903FBBF80390FE38FE0393FC387F8397F8383FC39FE0380FE00FC147E0078143C390007
C000A76D5A1F247BA62A>I<EF0180EF07C0171FEF7F80933801FF00EE07FCEE1FF0EE7F
C04B48C7FCED07FCED1FF0ED7FC04A48C8FCEC07FCEC1FF0EC7FC04948C9FCEB07FCEB1F
F0EB7FC04848CAFCEA07FCEA3FF0EA7FC048CBFC5AEA7F80EA3FE0EA0FF8EA03FEC66C7E
EB3FE0EB0FF8EB03FE903800FF80EC3FE0EC0FF8EC03FE913800FF80ED3FE0ED0FF8ED03
FE923800FF80EE3FE0EE0FF8EE03FE933800FF80EF3FC0170FEF038094C7FCAE007FB812
80B912C0A26C1780324479B441>20 D<181EA4181F84A285180785727EA2727E727E8519
7E85F11F80F10FC0F107F0007FBA12FCBCFCA26C19FCCCEA07F0F10FC0F11F80F13F0019
7E61614E5A4E5AA24E5A61180F96C7FCA260181EA4482C7BAA53>33
D<EC01F8140FEC3F80ECFC00495A495A495AA2130F5CB3A7131F5C133F49C7FC13FEEA03
F8EA7FE048C8FCEA7FE0EA03F8EA00FE137F6D7E131F80130FB3A7801307A26D7E6D7E6D
7EEC3F80EC0FF814011D537ABD2A>102 D<12FCEAFFC0EA07F0EA01FCEA007E7F80131F
80130FB3A7801307806D7E6D7EEB007EEC1FF0EC07F8EC1FF0EC7E00495A495A495A5C13
0F5CB3A7131F5C133F91C7FC137E485AEA07F0EAFFC000FCC8FC1D537ABD2A>I<007013
1C00F0131EB3B3B3B3A80070131C175277BD2A>107 D E
%EndDVIPSBitmapFont
/Fr 134[50 3[55 33 39 44 55 55 50 55 83 28 2[28 55 50
33 44 1[44 55 50 13[55 2[61 78 5[39 4[72 12[50 50 50
50 50 50 50 2[25 33 45[{TeXBase1Encoding ReEncodeFont}33
99.6264 /Times-Bold rf /Fs 133[42 42 42 42 1[42 42 42
42 42 42 42 42 42 42 42 1[42 42 42 42 42 42 42 42 42
3[42 1[42 42 42 1[42 42 42 1[42 42 42 42 42 1[42 42 42
42 42 4[42 42 2[42 2[42 1[42 42 1[42 5[42 42 1[42 42
42 42 42 1[42 42 42 39[{.85 ExtendFont TeXBase1Encoding ReEncodeFont}58
83.022 /Courier rf /Ft 171[38 6[55 77[{
.167 SlantFont TeXBase1Encoding ReEncodeFont}2 61.4362
/Times-Roman rf /Fu 139[23 1[28 1[42 42 42 65 23 2[23
42 1[28 37 42 1[42 37 12[51 6[74 51 7[60 1[55 60 7[42
42 42 42 42 42 42 42 42 42 2[28 21 44[{
.167 SlantFont TeXBase1Encoding ReEncodeFont}32 83.022
/Times-Roman rf /Fv 119[20 51[38 6[55 26[31 31 32[20
16[{TeXBase1Encoding ReEncodeFont}6 61.4362 /Times-Roman
rf
%DVIPSBitmapFont: Fw cmmi10 10 2
/Fw 2 60 df<121C127FEAFF80A5EA7F00121C0909798817>58 D<121C127FEAFF80A213
C0A3127F121C1200A412011380A2120313005A1206120E5A5A5A12600A19798817>I
E
%EndDVIPSBitmapFont
/Fx 87[28 45[32 37 1[55 37 42 23 32 32 42 42 42 42 60
23 37 23 23 42 42 23 37 42 37 42 42 6[46 1[51 69 51 60
46 42 51 60 51 60 55 69 46 55 37 28 60 60 51 51 60 55
51 51 6[28 42 42 42 42 42 42 42 42 42 42 1[21 28 21 2[28
28 28 36[42 2[{TeXBase1Encoding ReEncodeFont}69 83.022
/Times-Italic rf /Fy 26[59 113[50 2[46 12[46 27[50 71[{}5
83.022 /Symbol rf /Fz 105[42 1[37 37 10[28 13[37 42 42
60 42 42 23 32 28 42 42 42 42 65 23 42 23 23 42 42 28
37 42 37 42 37 3[28 1[28 51 1[60 78 60 60 51 46 55 1[46
60 60 74 51 60 32 28 60 60 46 51 60 55 55 60 3[47 1[23
23 42 42 42 42 42 42 42 42 42 42 23 21 28 21 47 1[28
28 28 1[69 3[28 16[28 12[46 46 2[{TeXBase1Encoding ReEncodeFont}82
83.022 /Times-Roman rf /FA 134[60 2[60 66 40 47 53 66
66 60 66 100 33 2[33 66 60 40 53 66 53 66 60 11[86 80
66 86 4[113 80 2[47 5[86 1[86 7[60 60 60 60 60 60 60
60 60 1[33 47[{TeXBase1Encoding ReEncodeFont}40 119.552
/Times-Bold rf /FB 136[51 1[51 51 51 51 2[51 51 51 51
2[51 51 51 51 51 51 51 51 51 16[51 5[51 9[51 13[51 51
2[51 46[{.85 ExtendFont TeXBase1Encoding ReEncodeFont}24
99.6264 /Courier rf /FD 139[21 29 25 2[37 37 4[21 37
2[33 1[33 1[33 10[54 8[66 26[37 37 37 1[19 1[19 44[{
TeXBase1Encoding ReEncodeFont}17 74.7198 /Times-Roman
rf /FE 134[50 2[50 1[28 39 33 2[50 50 78 28 2[28 50 50
33 44 50 44 1[44 11[72 61 3[55 1[72 89 2[39 2[72 1[61
72 66 8[28 7[50 50 50 3[25 44[{TeXBase1Encoding ReEncodeFont}32
99.6264 /Times-Roman rf /FF 139[40 3[72 72 72 112 40
2[40 72 2[64 72 1[72 64 12[88 6[128 88 2[48 4[104 1[96
104 6[40 12[48 45[{TeXBase1Encoding ReEncodeFont}21 143.462
/Times-Roman rf end
%%EndProlog
%%BeginSetup
%%Feature: *Resolution 600dpi
TeXDict begin
%%PaperSize: A4

%%EndSetup
%%Page: 1 1
1 0 bop 405 637 a FF(The)35 b(LAMBD)-6 b(A-Method:)40
b(Matlab)2266 585 y FE(TM)2466 637 y FF(Implementation)1444
820 y FD(V)-8 b(ersion)19 b(2.1,)f(12th)i(March)g(2001)1626
1160 y FE(Peter)25 b(Joosten)1105 1277 y(Mathematical)f(Geodesy)h(and)f
(Positioning)1190 1393 y(Ci)n(vil)f(Engineering)h(and)h(Geosciences)913
1509 y(Delft)g(Uni)n(v)o(ersity)e(of)i(T)-7 b(echnology)h(,)23
b(The)i(Netherlands)1132 1625 y(Email:)k FB([e-mail])
1547 1919 y FE(12th)24 b(March)h(2001)0 2382 y FA(1)119
b(Intr)n(oduction)0 2667 y Fz(One)16 b(of)f(the)h(k)o(e)o(ys)f(to)h
(precise)f(positioning)f(using)i(GPS)g(is)h(to)f(resolv)o(e)e(the)i
(double)e(dif)n(ferenced)g(carrier)g(phase)i(ambiguities.)22
b(One)0 2767 y(of)h(the)h(methods)f(to)g(do)h(this)g(is)h(the)e(LAMBD)m
(A-method,)f(de)n(v)o(eloped)f(at)j(the)g(Delft)g(Uni)n(v)o(ersity)e
(of)i(T)-6 b(echnology)h(.)32 b(LAMBD)m(A)0 2866 y(stands)20
b(for)g(\223Least-squares)f(AMBiguity)h(Decorrelation)e
(Adjustment\224.)24 b(When)c(using)f(the)h(LAMBD)m(A-method,)e(the)i
(inte)o(ger)0 2966 y(least)k(squares)f(ambiguity)e(estimates)j(are)f
(computed)e(in)i(tw)o(o)h(steps.)34 b(First)24 b(the)f(ambiguities)f
(are)h(decorrelated,)f(by)h(means)f(of)0 3066 y(the)27
b Fx(Z)t Fz(-transformation.)44 b(Then)27 b(the)g(inte)o(ger)g
(minimization)f(problem)f(is)k(solv)o(ed)d(by)h(a)h(discrete)f(search)g
(o)o(v)o(er)f(an)i(ellipsoidal)0 3165 y(re)o(gion,)18
b(the)j(ambiguity)d(search)i(ellipsoid.)0 3365 y(The)i(method)f(has)i
(been)e(introduced)f(in)j([7)o(])g(by)f(P)-9 b(.J.G.)22
b(T)-6 b(eunissen.)31 b(The)22 b(method)f(w)o(as)i(\002rst)g
(implemented)e(in)h(F)o(ortran-77,)f(by)0 3464 y(P)-9
b(.J.)23 b(de)g(Jonge)g(and)f(C.C.J.M.)i(T)m(iberius.)33
b(A)24 b(detailed)e(description)g(of)g(this)i(implementation)d(is)j(gi)
n(v)o(en)e(in)h([3)o(].)34 b(The)23 b(code)f(w)o(as)0
3564 y(translated)e(into)f(Matlab)729 3534 y Fv(TM)872
3564 y Fz(by)h(K.)g(Borre)g(\(see)h([6)o(]\))0 3763 y(When)e(a)h
(research)e(on)h(the)h(stochastic)f(properties)f(of)h(the)h(estimated)f
(inte)o(ger)f(ambiguities)g(w)o(as)j(started)e(in)g(1998,)f(initially)i
(this)0 3863 y(Matlab)236 3833 y Fv(TM)379 3863 y Fz(implementation)f
(w)o(as)k(used.)30 b(Ho)n(we)n(v)o(er)20 b(it)j(w)o(as)f(found)f(that)g
(some)h(impro)o(v)o(ements)d(were)j(necessary)-5 b(.)29
b(One)21 b(of)h(the)0 3962 y(biggest)e(dra)o(wbacks)f(of)i(the)g
(earlier)f(v)o(ersion)g(w)o(as)i(the)f(f)o(act)g(that)g(it)g(w)o(as)h
(not)f(possible)f(to)h(\002nd)g(more)f(than)g(tw)o(o)h(candidates.)26
b(As)0 4062 y(a)21 b(result)f(of)g(it,)g(it)h(w)o(as)g(decided)e(to)i
(create)f(a)g(ne)n(w)g(v)o(ersion,)f(based)h(on)f(both)h(the)g
(description)f(in)41 b([3)o(])20 b(and)g([6)o(].)0 4261
y(Compared)f(to)h(the)g(earlier)g(v)o(ersion,)f(this)h(ne)n(w)g(v)o
(ersion)f(of)n(fers:)138 4527 y(-)42 b(Impro)o(v)o(ed)17
b(readability)h(of)i(the)h(softw)o(are.)138 4693 y(-)42
b(A)17 b(more)g(modular)e(approach,)g(enabling)h(one)h(to)g(perform)e
(for)i(e)o(xample)e(only)i(the)g(decorrelation-step,)d(only)j(the)g
(search-)208 4793 y(step,)j(or)g(both)f(steps.)138 4959
y(-)42 b(It)18 b(is)h(no)n(w)e(possible)h(to)g(\002nd)g(a)g(number)e
(of)i(candidates)f(dif)n(ferent)f(from)h(2,)i(which)e(can)h(be)g
(useful)f(for)h(research)f(purposes.)138 5125 y(-)42
b(Sometimes)18 b(intermediate)g(output)f(can)i(be)g(v)o(ery)f(helpful)g
(for)g(research)g(purposes,)g(sometimes)g(it)i(can)f(be)g(rather)f
(anno)o(y-)208 5224 y(ing.)24 b(It)d(is)g(no)n(w)e(possible)h(to)h
(switch)f(on)g(or)g(of)n(f)f(this)i(output.)138 5390
y(-)42 b(In)19 b(the)g(original)g(code,)g(not)g(much)f(of)i(Matlab)1551
5360 y Fv(TM)1693 5390 y Fz(')-5 b(s)21 b(v)o(ectorized)d(computation)f
(possibilities)j(are)f(used.)25 b(Using)19 b(this)h(can)208
5490 y(impro)o(v)o(e)e(both)h(readability)g(of)h(the)g(code)f(and)h
(computational)e(speed.)1869 5871 y(1)p eop
%%Page: 2 2
2 1 bop 0 -47 a Fu(The)20 b(LAMBD)m(A-Method,)e(Matlab)1070
-77 y Ft(TM)1212 -47 y Fu(implementation)2002 b(2)p 0
-14 3780 4 v 138 272 a Fz(-)42 b(The)18 b(original)f(code)h(is)i(a)f
(rather)e(e)o(xtensi)n(v)o(e)h(set)h(of)f(subroutines.)23
b(Readability)18 b(w)o(as)i(impro)o(v)o(ed)15 b(by)k(reducing)d(the)j
(number)208 372 y(of)g(subroutines.)138 525 y(-)42 b(The)18
b(dimension)g(of)g(the)h(problem)f(is)i(deri)n(v)o(ed)d(from)h(the)h
(input)f(using)g(the)h Fs(size)p Fz(-function.)24 b(This)19
b(decreases)g(the)g(number)208 625 y(of)g(ar)o(guments)g(to)h
(functions,)e(and)i(therefore)e(the)j(probability)d(of)i(mistak)o(es.)0
999 y FA(2)119 b(Ambiguity)31 b(r)n(esolution)f(in)h(brief)0
1301 y Fr(2.1)99 b(P)o(arameter)26 b(estimation)0 1556
y Fz(GPS)32 b(ambiguity)d(resolution)h(is)i(the)f(process)f(of)h
(resolving)f(the)h(unkno)n(wn)d(c)o(ycle)j(ambiguities)f(of)h(the)g
(double)e(dif)n(ference)0 1656 y(carrier)e(phase)h(data)g(as)h(inte)o
(gers.)49 b(The)28 b(GPS)h(models)f(on)f(which)h(ambiguity)f
(resolution)g(is)i(based,)h(can)e(all)g(be)h(cast)g(in)f(the)0
1755 y(follo)n(wing)19 b(conceptual)f(frame)h(of)h(linear\(ized\))f
(observ)n(ation)f(equations:)236 1996 y Fx(E)6 b Fq(f)p
Fx(y)p Fq(g)17 b Fp(=)h Fx(Aa)12 b Fp(+)g Fx(Bb)g Fp(+)g
Fx(e)2768 b Fz(\(1\))0 2237 y(Where)24 b Fx(y)g Fz(is)i(the)e(gi)n(v)o
(en)f(GPS)h(data)g(v)o(ector)m(,)g Fx(a)g Fz(and)g Fx(b)g
Fz(are)g(the)g(unkno)n(wn)e(parameter)g(v)o(ectors)h(of)h(order)f
Fx(n)h Fz(and)30 b Fx(p)24 b Fz(respecti)n(v)o(ely)-5
b(,)0 2336 y(and)22 b(where)h Fx(e)g Fz(is)h(the)f(noise)g(v)o(ector)f
(of)g(order)g Fx(m)p Fz(.)34 b(The)22 b(matrices)h Fx(A)g
Fz(and)g Fx(B)g Fz(are)g(the)g(corresponding)d(design)i(matrices)h(of)f
(order)0 2436 y Fx(m)12 b Fq(\002)g Fx(n)21 b Fz(and)g
Fx(m)12 b Fq(\002)18 b Fx(p)j Fz(respecti)n(v)o(ely)-5
b(.)27 b(The)21 b(data)g(v)o(ector)f Fx(y)i Fz(will)h(usually)d
(consist)i(of)f(the)h(\223observ)o(ed)d(minus)i(computed\224)e(single-)
i(or)0 2536 y(dual-frequenc)o(y)16 b(DD)21 b(phase)e(and/or)g(pseudo)g
(range)g(\(code\))g(observ)n(ations,)f(accumulated)g(o)o(v)o(er)h(all)h
(observ)n(ation)e(epochs.)24 b(The)0 2635 y(entries)d(of)g(v)o(ector)f
Fx(a)i Fz(are)f(then)g(the)g(DD)h(carrier)e(phase)h(ambiguities,)g(e)o
(xpressed)f(in)h(units)g(of)g(c)o(ycles)h(rather)e(than)h(range.)27
b(The)o(y)0 2735 y(are)d(kno)n(wn)f(to)h(be)g(inte)o(gers.)36
b(The)24 b(entries)g(of)g(v)o(ector)f Fx(b)i Fz(will)f(consist)h(of)f
(the)g(remaining)e(unkno)n(wn)g(parameters,)i(such)g(as)h(for)0
2834 y(instance)19 b(baseline)f(components)f(\(coordinates\))g(and)h
(possibly)h(atmospheric)e(delay)i(parameters)f(\(troposphere,)e
(ionosphere\).)0 2934 y(The)26 b(procedure)e(which)i(is)h(usually)f
(follo)n(wed)f(for)g(solving)h(the)g(abo)o(v)o(e)f(model)g(can)h(be)h
(di)n(vided)e(into)h(three)g(steps)h(\(for)e(more)0 3034