Some remarks on GNSS integer ambiguity validation methods
|
|
- Pierce Manning
- 6 years ago
- Views:
Transcription
1 Some remrks on GNSS integer mbiguity vlidtion methods T. Li nd J. Wng* Ambiguity resolution is n indispensble step in fst nd high precision Globl Nvigtion Stellite System (GNSS) bsed positioning. In generl, mbiguity resolution consists of three steps. The first step is to estimte the mbiguities using lest-squres estimtion process, from which the so clled flot solution or rel vlued solution is obtined. Then in the second step, the flot solution is used to serch for the integer mbiguities. Once integer mbiguities re resolved, the lst step is to pply the integer mbiguities into the models so tht firly ccurte fixed solution cn be generted. Owing to the importnce of the integer mbiguities, one indispensble procedure tht needs to be implemented in the second step is integer mbiguity vlidtion. Over the pst decdes, considerble work hs been concentrted on this procedure nd vrious pproches hve been proposed, such s R-rtio test, F-rtio test, W-rtio test, integer perture estimtor, etc. However, their performnces re controversil. Therefore, in this contribution, n overview of the existing mbiguity vlidtion methods is firstly presented, nd then some numericl nlysis is crried out to evlute their performnces. Keywords: GNSS, Integer lest-squres, Ambiguity vlidtion, Integer perture, Rtio test Introduction In Globl Nvigtion Stellite System (GNSS) positioning, crrier phse mesurements re more precise thn the pseudornges. However, one troublesome problem with the crrier phse mesurements is tht ech crrier phse contins n mbiguity in the number of wvelengths, which hs n intrinsic nture of being integer. Consequently, the most importnt issue in precise GNSS positioning is tht integer mbiguities hve to be resolved correctly, nd once integer mbiguities re obtined, the precision of crrier phses cn be fully exploited in positioning. The whole process of mbiguity resolution includes three steps. Initilly, with n pproprite mthemticl model, the flot solution (or rel vlued solution) cn be estimted by lest-squres or Klmn filter. Then, serching strtegy will be crried out round the flot solution so s to find the integer cndidtes [3]. As there could be severl sets of integer mbiguity cndidtes, vlidtion procedure is required, nd the most likely cndidte should be vlidted nd seprted from the others bsed on sttisticl theory or other methods. To fix the rel vlued mbiguities to the integer cndidtes, in this pper, LAMBDA [12], [13], [14], ws dopted. The lst step is to estimte the unknown prmeters with the correct integer mbiguities. School of Surveying nd Sptil Informtion Systems, University of New South Wles, NSW 2052, Sydney, Austrli *Corresponding uthor, emil jinling.wng@unsw.edu.u For the mbiguity vlidtion step, severl procedures hve been developed, with the qudrtic form of residuls ssocited with the most likely set of integer mbiguities, nd the second qudrtic form of residuls ssocited with the second most likely set of integer mbiguities, such s F-rtio [5], R-rtio [4], W-rtio [24], difference test [23], projector test [7]. With given criticl vlue, the best cndidte will be identified. In such mbiguity vlidtion procedures, the fixed mbiguities re treted s constnt vlues. In nother pproch to mbiguity vlidtion, integer perture (IA) estimtor [18] hs been developed. In [28], the IA estimtor ws considered s frmework for ll other clssicl theory of hypothesis testing methods, nd the geometries of different vlidtion methods re then reflected in the perture pull-in region. Following the definition, mbiguity vlidtion is crried out by determining the criticl vlues bsed on the pre-defined filrte, which sounds more resonble. However, s commented in [9], the theory of the IA estimtor ssumes the most likely set of cndidtes re correct (or true integer vlues) in the designing stge. Through numericl tests, it is found tht determining the criticl vlues ccording to the IA estimtor is not tht strightforwrd, especilly in rel pplictions, which shows tht mbiguity vlidtion by the IA estimtor hs some limittions. In this contribution, the performnces of different mbiguity vlidtion methods hve been studied nd nlysed. This pper firstly strts with introducing the mthemticl models, nd then in the following section, n overview of these vlidtion methods s well s some discussion is presented. After tht, ccording to numericl nlysis, the performnce of different mbiguity ß 2012 Survey Review Ltd Received 5 December 2011; ccepted 5 December Survey Review 2012 VOL 44 NO 326 DOI / Y
2 vlidtion methods nd the limittion of IA theory re extensively nlysed nd discussed. Some conclusions nd suggestions re given in the finl section. Mthemticl models The initil GNSS observtion mthemticl models re constructed with unknown prmeters such s coordintes, integer mbiguities, ionosphere dely, troposphere dely, etc. Ionosphere dely nd troposphere dely certinly hve impcts on the ccurcy, nd in order to eliminte those unnecessry fctors, double differenced functionl models re employed s follows D+w~ 1 l D+rzD+Nze w (1) D+P~D+rze P (2) where D= is the double differencing opertor between stellites nd receivers, w nd P re crrier phse mesurements nd code mesurements respectively, l is the crrier phse wvelength, r is the geometric distnce between stellites nd receivers, N is the integer mbiguities in cycles nd e represents the mesurement errors. With n pproximte rover position given, the bove models re linerised s l~axzv (3) where ~ D+w{ 1 T l D+r,D+P{D+r T,, v~ e w,e P x~ ðx r,d+nþ T, A is the design mtrix of both coordintes A xr nd mbiguities A, r* is the pproximte distnce, nd x r is the coordintes. Thus, the stochstic model for the double differenced mesurements cn be expressed s Dy ð Þ~s 2 0 Q~s2 0 P{1 (4) where D is the covrince mtrix, s 2 0 is priori vrince fctor, nd Q nd P re the cofctor mtrix nd weight mtrix of the mesurements, respectively. By pplying the clssicl lest-squres pproch, whose criterion is v T Pv~min, the unknown prmeters nd their covrince re uniquely estimted s [25] T~ x ~ x r A T {1A PA T Pl Q~ x A T PA " # {1~ Q x Qr x r Q x r Q where represents the rel vlued solution of integer mbiguities. Then the posteriori vrince is 2 TPv s0 ~v ~ V 0 f f with f is the degree of freedom. An overview of existing mbiguity vlidtion methods The unknown prmeters estimted from norml lest-squres procedure re the so clled the flot solution. Considering the integer nture of the mbiguities, n integer lest-squres (ILS) problem needs to be tken into ccount, e.g. [8], [12], [16]. In this pper, (5) (6) LAMBDA is used to serch the integer mbiguities within the constructed hyper-ellipsoid. The remining problem is tht the best integer cndidte should be selected from ll the combintions of integer cndidtes, nd trditionlly, it is ccomplished by investigting the qudrtic forms of residuls. Another proposed theory is the integer perture estimtion, which hs hybrid nture of yielding integer outcomes s well s noninteger outcomes. This pproch llows users to control the pre-defined fil-rte by themselves, yet its ppliction in prctice requires more numericl tests nd nlysis. Bsed on the estimtion from lest-squres, mbiguity vlidtion will strt from, 2 s 0 nd Q, nd we ssume tht the flot solution of the mbiguity is normlly distributed s N,Q ð Þ. The qudrtic form of residuls for flot solution is define s V 0, while the qudrtic form of the fixed solution residuls is V, while TQ R~ {1 {^ {^ mens the geometric distnce between V 0 nd V [25]. We ssume V 1 nd V 2 re the minimum nd second minimum qudrtic forms of the residuls in the mbiguity fixed solutions, which correspond to ă 1, ă 2 nd R 1, R 2, then mbiguity vlidtion methods cn be generlly divided into two min groups [9]. Sttisticl tests bsed on the best nd second best mbiguity cndidtes Trditionlly, mbiguity vlidtion methods consider both the best nd second best mbiguity cndidtes nd the reltionship between them re nlysed from sttistic point of view. Most commonly used methods re the F- rtio, R-rtio, W-rtio, difference test nd projector test, nd their detiled definitions re expressed s follows: Rtio tests The first pproch proposed for mbiguity vlidtion is the F-rtio test [5], [27] (here inverse version is used for nlysis), which is the rtio between V 1 nd V 2 F~ V 1 ~ V 0zR 1 (7) V 2 V 0 zr 2 Normlly, the test sttistic for F-rtio is ssumed to hve n F-distribution, with degrees of freedom for V 1 nd V 2 respectively. Then, criticl vlue c could be obtined by specifying significnce level. However, ccording to, e.g. [24] nd [15], the numertor nd denomintor re not independent, nd F-rtio does not follow n F-distribution. Sometimes, there re lrge discrepncies between the obtined results nd their true vlues, which clerly show tht F-rtio is tht not relible, even though, certin empiricl vlues, e.g. 1?5 (inversed 0?67), 2 (inversed 0?5), 3 (inversed 0?3), perform stisfctorily [2], [10]. An lterntive populr test similr to the F-rtio test is the rtio between R 1 nd R 2, known s the R-rtio [4], [27]. The inverse of R-rtio is constructed s R~ R 1 ƒc (8) R 2 As described in the bove formul, the distribution of R- rtio is lso unknown. Hence, it would be incorrect to determine the performnce of the rtio test on the bsis of the distributionl results s provided by the clssicl theory of hypothesis testing [22], [28]. Theoreticlly, the Survey Review 2012 VOL 44 NO
3 chosen of criticl vlue is groundless by n empiricl vlue. In the originl work of [4], the criticl vlues rnges from 5 to 10 were given. However, in prctice, mny reserchers vlidte the mbiguities bsed on the R-rtio, with criticl vlues like 2 (inversed 0?5), 2?5 (inversed 0?4), see e.g. [6], [11], [22]. In [24], the discrimintion procedure is constructed either by compring the likelihood of two integer cndidtes or by rtificilly nesting the two compred models by nesting prmeter, nd then the W-rtio is defined s follows d W~ c (9) 1=2 ½vrðdÞŠ where d~v 2 {V 1, vrðdþ~r 2 Q d (10) where r 2 could be decided by users either from the priori vrince s 2 0 or from the posteriori vrince 2 s 0. By pplying the vrince covrince propgtion lw, the vrince for TQ {1 d could be Q d ~4 ^ 2{^1 ^ 2{^1. Assuming tht W nd W s re the two rtios corresponding to priori vrince s 2 0 nd posteriori vrince 2 s 0,they re supposed to hve truncted stndrd norml distribution nd truncted student t-distribution respectively, from which the criticl vlues could be esily obtined. Under the ssumption tht the fixed mbiguities re deterministic quntities, the W-rtios cn provide rigorous confidence level for the vlidtion test. Difference test In [23], the vlidtion pproch is to nlyse the difference between R 1 nd R 2. First, globl model test is crried out. When it pssed, the mbiguity could be vlidted by compring the difference test with its criticl vlue s follows R 2 {R 1 c (11) where c is non-negtive sclr, which is user defined tolernce level. This test ccepts the best solution if the flot solution is much closer to the best solution thn to the second best solution. Criticl vlues of 15 nd 12 were suggested in [23] nd [7]. Still, the determintion of the criticl vlue, which is lso empiricl, is controversil, s the distribution of the difference test is unknown. A looser setting of criticl vlues increses the performnce, but lso increses the number of wrongly ccepted cndidtes (type I errors). On the other hnd, the difference sttistic is dependent on the priori vrince fctor given in eqution (4). Projector test The projector test is bstrcted from [1], nd proposed in [7], lso deducted in [24], the null hypothesis is tht there is no outlier so ^1 is ccepted, while the lterntive hypothesis is tht there is n outlier in the direction of (^ 2{^1); consequently, we hve y~a ^ 1zA ^ 2{^1 cza xr x rzv With quntity c nd its vrince s TQ ^2{^ 1 {1 ^{^1 c^~ TQ, s^2c ~ 4V 1{c 2 Qd (12) {1 ðn{4þq d ^ 2{^1 ^ 2{^1 Then the definition of the test sttistic is c =s ~t(n24) c (in this contribution, only bseline components re considered) if ^ 2{^ 1 is selected. If ^2 is selected, the test sttistic is then non-centrl t-distribution. The quntity hs been proved to be lwys smller thn or equl to 0?5, so tht the definition nd the distribution re ctully not strict s well. Besides, s discussed in [24], the derivtion is unfortuntely not rigorous becuse the non-centrl prmeter is formulted with the estimted vrince insted of its known vlue. So fr, there hs been no criterion identifying which of the mbiguity vlidtion methods proposed is the best, since the criticl vlues re chosen differently for different distributions. Their performnces should be evluted nd tested through vrious mens, including theoreticl nlysis nd evlution with the ground truth for the integer mbiguities. Integer perture theory On the bsis of ILS estimtor, IA theory ws firstly introduced by [18], nd the IA estimtor - is developed s -~ X ( zv z z 1{ X h i ) v z (13) z [ Z n z [ Z n With the indictor function v z (x) defined s v z ðxþ~ 1 if x[v z (14) 0 otherwise where the V z re the perture pull-in regions, nd their union V5R n is the perture spce, which lso hs property of trnsltionl invrint. With the bove definition, three outcomes cn be distinguished s [V Success: correct integer estimtion [V=V Filure: incorrect integer estimtion [/ V Undecided: mbiguity not fixed to n integer The corresponding probbilities of success (P s ), filure (P f ) nd undecided (P u ) re given by ð P s ~Pð-~ Þ~ ðxþdx ð P f ~ V=V f P u ~1{P s {P f V f ðxþdx (15) In the cse of GNSS model, f ðxþ represents the probbility density function of the flot mbiguities, nd is usully ssumed to be normlly distributed. The IA estimtor llows the user to choose predefined fil-rte, nd then determine the criticl vlue ccordingly. With the determintion of the criticl vlue, the user is cpble of compring different vlidtion methods by clculting the corresponding success rtes. Ellipsoidl integer perture A rther strightforwrd wy of determining the IA pull-in region is the ellipsoidl integer perture (EIA) [19], which cn be constructed s E z ~E 0 zz, E 0 ~S 0 \C e,0, Vz[Z n (16) With n S 0 being the lest-squres o pull-in region nd C e,0 ~ x[r n 2 jjjxjj Q ƒc 2, n origin centred ellipsoidl 232 Survey Review 2012 VOL 44 NO 326
4 region of which the size is controlled by the perture prmeter (or criticl vlue, in the sequel, criticl vlue is used insted) e with kxk 2 Q ~x T Q {1 x The structure of the EIA pull-in region is considered to hve Chi-squre distribution, nd thus the EIA probbilities of success, filure, nd undecided re given s P f ~ X P x 2 ðn,l z Þƒc 2 z [ Z n = f0g P s ~P x 2 ðn,0þƒc 2 P u ~1{P f {P s (17) where x 2 ðn,l z Þ represents rndom vrible hving non-centrl Chi-squre distribution, nd n nd l z re the degree of freedom nd the non-centrl prmeter l z ~z T Q {1 z respectively. Optiml integer perture Optiml integer perture (OIA), which ims t the mximistion of the success rte, is proposed by [20], [28] s mx P s subject to : P f ~b with the pull-in region V 0 5S 0, nd b is the pre-defined fil-rte. The relistion of OIA could be considered s the rtio between the PDF of the mbiguity residul nd the PDF of the mbiguities s P f exp { ð1=2þe^zz 2 e^ e^ z [ Z ~ n Q ƒc (18) f e^ exp { ð1=2þe^ 2 Q with e^ ~ {^. According to the definition of OIA, it hs the dvntge of obtining the highest success rte. Wheres its ppliction in relity is quite time consuming, this cn be known from the PDF of the mbiguity residul. In [17], n exct formul of f e^ ðxþ is derived, nd in order to clculte f e^ ðxþ, summtion over infinite integers should be crried out. This is, however, impossible to operte in prctice. As consequence, severl finite integers re serched within n ellipsoidl spce, nd pplied [29]. In the next section, we will discuss the following two questions through numericl nlysis: (1) how to compre different vlidtion methods; nd (2) how to determine criticl vlues from given fil-rte for different vlidtion methods. Numericl nlysis According to the bove mentioned sttements, different mbiguity vlidtion pproches re hrd to compre directly, s the different chosen criticl vlues. However, under the frmework of the IA estimtion, this is possible if pre-defined fil-rte is provided. In this section, the comprison of different vlidtion methods is firstly studied by simultion, nd then the limittions of IA theory re discussed. In ddition, the performnce 1 Determintion of criticl vlues with given fil-rte of both IA nd non-ia bsed methods re nlysed in rel pplictions. Comprison of different vlidtion methods With the sme pre-defined fil-rte given, one is ble to select the criticl vlue for different vlidtion methods, nd then determine their corresponding success rtes. By compring their success-rtes, we re ble to figure out which method my hve the best performnce. Detiled procedures re shown s Fig. 1. The geometry informtion of n existing GNSS model is entirely reflected in Q, so Monte Crlo simultions were used to investigte the performnce of different vlidtion methods with the informtion of Q given. By following Fig. 1, rndom genertor ws used to simulte number of smples (e.g., ) of the estimted flot mbiguities. Here, bivrite norml distribution with the following vrince covrince mtrix is considered first [28] Q~ ½0 : 1021, {0 : 0364; {0 : 0364, 0 : 1100Š; Survey Review 2012 VOL 44 NO
5 The reltionships between the fil-rte nd the criticl vlue for ech of the mbiguity vlidtion methods re lwys on one to one bsis, which implies tht the filrte chnges with different chosen criticl vlues for ech vlidtion method. The corresponding success rtes obtined with given fil-rte re shown in Tble 1. It is noted tht the OIA, which hs the highest success rte, performs the best, no mtter how big the fil-rte is. In this cse, the R-rtio, W -rtio nd difference test, which re ll treted s seprted IA method, hve performnce quite close to tht of the optiml IA estimtor. When the fil-rte increses from 0?005 to 0?02, these three tests hve quite similr success rte, but it is impossible to judge which one is the best mong these three. The reson is tht the smple size of simultion will hve some impct on the results. Such impcts re lrge enough to cuse slight differences in the success rtes. As for the projector test s n IA method, its success rte is smller thn the bove three. For the EIA, it is esy to mnipulte. However, its performnce is not tht good. In [9], combintion of R-rtio nd overlpped EIA ws proposed. To void the conservtiveness of the EIA, the EIA is llowed to be overlpped to increse the success rte, nd for the purpose of ensuring tht the best mbiguity cndidte is sttisticlly better thn the others, the R-rtio is combined with the overlpped EIA. Theoreticlly, however, the definition of the EIA is not overlpped, nd the overlpped EIA increses the fil-rte s well. Aprt from this, the combintion of the EIA nd the R-rtio mkes the success rte decrese, which is not the purpose of the overlpped EIA. F-rtio s n IA method performs slightly worse thn the other vlidtion methods. The reson could be tht, s shown in eqution (7), nd lso commented in [28], both flot solution residuls nd mbiguity residuls re considered, so tht the dimension of perture spce is different. Tble 2 presents four-dimensionl rel dt exmple, with the vrince-covrince mtrix s follows: Q~ ½1:1834, {1:3348, 0:2814, {0:6537; {1:3348, 1:6088, {0:2619, 0:8123; 0:2814, {0:2619, 0:0975, {0:1139; {0:6537, 0:8123, {0:1139, 0:4177Š Tble 1 Two-dimensionl cse for success rtes nd criticl vlues with the pre-defined-fil rte P f FTIA RTIA WTIA DTIA PTIA EIA OIA c P s c P s c P s For the fil-rte of 0?001, the OIA still yields the highest success rte, nd the difference test, in this cse, hs much better performnce thn the others. When the filrte is 0?01, the W-rtio test is more preferble thn ny other non-ia vlidtion methods. Menwhile, in both cses, the EIA hs the worst performnce. It should be mentioned tht for four-dimensionl simultion, it consumes hours of computtion for ech method when the user choose the Mtlb build-in function fzero ( good tool to find zero of function ner roughly given criticl vlue) to run the simultion. The numericl results revel three issues: (1) the OIA in these two cses performs the best in the simultions, but its efficiency is the worst mong ll the vlidtion methods, which indictes tht this method is hrdly pplicble in rel time GNSS positioning; (2) The higher the dimension, the more time the simultions will consume; (3) using the function fzero to serch for the criticl vlue, the pre-defined fil-rte should be resonbly determined by the user, since for certin stellite geometry, there is n upper bound for the fil-rte. Limittions of the IA methods: pre-defined filrte versus the criticl vlue The integer perture theory llows the user to determine the criticl vlues bsed on the user preferred fil-rte. Theoreticlly, once the stellite geometry is obtined, the reltionship between the pre-defined fil-rte nd the criticl vlue hs lredy been determined, which mens tht when fil-rte is pre-defined, the criticl vlue will be uniquely fixed under the current stellite geometry. However, the problem is tht their reltionship cnnot be mthemticlly represented; only simultions re fesible. Thus the wy of clculting criticl vlue with pre-defined fil-rte will be entirely dependent on simultions. In order to generte relible results, the number of smples should be s lrge s possible. As suggested in [19], [28], desirble smple size N should be lrger thn In order to show the influence of the smple size N, two simultions were crried out to compre with the vlues given in [28], note tht the sme Q 5[0?0865, 20?0364; 20?0364, 0?0847] ws utilised. First, the smple size ws incresed. As listed in Tble 3, we utilised the sme criticl vlue, nd the Tble 3 Comprison of R-rtio s n IA estimtor N c P f P s Ref Tble 2 Four-dimensionl cse for success-rtes nd criticl vlues with the pre-defined fil-rtes P f FTIA RTIA WTIA DTIA PTIA EIA OIA c P s c P s Survey Review 2012 VOL 44 NO 326
6 2 Fil-rte chnges with both smple size nd number of simultion for R-rtio simulted results fluctuted ner the given vlue. Even if the smple size is constnt, s shown in Tble 4, the filrte nd success rte re not stble. The smple size N undoubtedly influences the fil-rte with the criticl vlue known (or criticl vlue with the fil-rte given). Hence, the smple size N, to some extent, determines the relibility of Integer Aperture theory. As shown in Figs. 2 4, when the smple size increses from 5000 to , the fil-rte nd the success rte re becoming smoother. This tells us tht the lrger the smple size, the more relible results we cn hve. However, s routine, the simulted results my be problemtic if the ccurcy of the criticl vlue is importnt. According to numericl nlysis, smple size lrger thn will hve more relible solution. However, with powerful PC, such simultions (given the criticl vlue to clculte the fil-rte) for one solution will tke over 1?5 s nd thus hrdly be usble for rel time opertions. Limittions of the IA methods in rel pplictions However, the problem is tht incresing the smple size ggrvtes the computtion burden, which is troublesome problem in kinemtic positioning. The wy of determining the criticl vlue more rpidly nd relibly bsed on the pre-defined fil-rte is perhps nother mjor issue. A fesible wy could be using look-up tble (here only for R-rtio), which ws creted in [21], 3 Fil-rte chnges with both smple size nd number of simultion for W -rtio [22]. The look-up tble is generted bsed on the lower bound of ILS fil-rte (expressed by the fil-rte of integer bootstrpping), nd then pply liner interpoltion to find the corresponding criticl vlue. In order to nlyse the performnce of the look-up tble, rel dt were collected t smpling rte of 1 s on 12 September 2009, in Sydney, Austrli, nd then processed by different session lengths. With single frequency crrier phse mesurements nd n elevtion ngle of 15u, five stellites were used to execute the double difference so tht four mbiguities need to be resolved. As shown in Tble 5, it is obvious tht, by determining the criticl vlue c from the look-up tble, when the session length is longer thn 3 min, the resolved mbiguities were ccepted s undoubtedly correct. However, the truth is tht correct mbiguities (given by processing the whole dtset) could only be obtined when the session length is longer thn 6 min. This tells us tht using look-up tble to determine the criticl vlue is not verstile, nd fixed Tble 4 Ten simultions of R-rtio with the smple size N c P f P s Ref Fil-rte chnges with both smple size nd number of simultion for difference test Survey Review 2012 VOL 44 NO
7 5 W-rtio (W nd W s ) vlue chnges with the epoch number 6 R-rtio nd F-rtio vlues chnge with epoch number look-up tble sometimes does not properly reflect the chnging of stellite geometry, even though it is more efficient to operte the simultions. In Tble 6, the differences between the look-up tble nd simultion re given. With the criticl vlues from the look-up tble, the simulted fil-rtes re not quite close to 0?001 nd 0?01 respectively, with notble error of 10 20%, nd in the cse of the sme fil-rtes given, the simulted criticl vlues lso hve some discrepncies with the look-up tble results. Note tht the Mtlb build-in function fzero ws utilised here to find the corresponding criticl vlues more efficiently nd rigorously (compred with the repets of the simultion). However, n even more noticeble phenomenon is neither using look-up tble nor simultion gives the correct decision in this cse (type I error). Performnces of both the IA nd non-ia bsed methods in rel ppliction Another rel dtset ws collected on 9 June 2010, in Sydney, Austrli. A session length of 10 min (smpling rte: 1 s) ws used here to briefly show the performnce of both trditionlly used methods to determine the criticl vlue nd using simultions bsed on the IA theory. Only single frequency code nd crrier phse were utilised, nd fter double differencing, there re Tble 5 Ambiguity vlidtion by using the look-up tbles to determine the criticl vlues 7 Difference nd projector tests vlues chnge with epoch number Tble 6 An insight into the determintion of criticl vlues from simultions nd look-up tble min c c P f P f P f P f c c min c c P f P f P f P f c c Integer mbiguity P f P f Session length/min sv17 sv3 sv18 sv21 RTIA c Criticl vlue c 01 P s Criticl vlue c 02 P s Vlidtion results* CR CR WA WA WA CA CA CA *CR5correctly rejected, WA5wrongly ccepted, CA5correctly ccepted. 236 Survey Review 2012 VOL 44 NO 326
8 seven mbiguities to be estimted with the degree of freedom s 4, nd the dt were processed on n epoch by epoch bsis. Figures 5 7 show the rtio vlues between the best solution nd the second best solution. As shown in Tble 7, F-rtio dopts commonly used criticl vlue 0?67, which yields 560 ccepted epochs out of the totl 600 epochs. The R-rtio ccepts 506 epochs with n empiricl criticl vlue of 0?5. For the W -rtio nd W s -rtio here, with the truncted norml distribution nd the truncted Student s t-distribution respectively, by significnce level, the corresponding criticl vlue cn be generted. Five hundred nd sixty-eight epochs were ccepted for W with confidence level of 0?99. For the W s, confidence level of 0?99 with degree of freedom s 4 llows 438 epochs ccepted. For the difference test, if n empiricl vlue of 15 ws chosen, it is too conservtive to operte nd the ccepted number of epoch is zero. The empiricl vlue suggested in [7] is lso too conservtive in this cse. When specifying criticl vlue of 0?87 for the projector test [28], 412 epochs were ccepted. However, mong the 600 best mbiguity combintions, the totl numbers of correct mbiguities re 597, nd for the W-rtio test, nd the projector test, the three sets of wrong mbiguities were correctly rejected, wheres for the F-rtio nd R-rtio tests, those three epochs hd been wrongly ccepted. In the cse of running simultions to determine the criticl vlue, for ech epoch, with given fil-rte, the geometry of flot solution ws used to obtin the criticl vlues. Owing to the short observtion time spn, the stellite geometry chnges slightly, nd so does the criticl vlues determined from the simultions. Menwhile, s mentioned before, there re hevy computtion burdens for higher dimension cses (more thn 2 h per simultion); consequently, the pproximte criticl vlues were used, nd the performnce of the IA theory in rel pplictions ws shown in Tble 8. The results re, however, not preferble compring the results from the non-ia s listed in Tble 7 (except the WTIA, which performed the best in both pproches). When fil-rte of 0?01 is given, the corresponding criticl vlues cn be simulted, nd then the IA bsed criticl vlues were pplied. For the EIA, once the criticl vlue obtined, we could compre its rel fil-rte with the predefined fil-rte. For the OIA, eqution (18) ws implemented. Besides the hevy computtion burden, it cn be seen tht in this cse, the OIA does not mke the optiml decision, nd the non-ia methods outperform the IA bsed method for most of the time in this rel ppliction. Another issue is tht the success rtes re extremely low, much smller thn the ILS success rte (59?3%), nd the probbility of undecided solution is quite high. This phenomenon is not consistent with the rel performnce, which cn be judged from the number of the ccepted solutions. Moreover, the IA estimtor committed type I error s did some conventionl methods. Obviously, in this rel ppliction, the success rte does not reflect the true correct rte properly, nd the IA estimtor does not perform resonbly. Concluding remrks To sum up, this contribution presents n overview of current mbiguity vlidtion pproches which hve been numericlly compred nd nlysed. Trditionlly, the comprison mong different mbiguity vlidtion methods is difficult s the determintion of criticl vlues vry with different methods. Within the frmework of the IA theory, the users re ble to control the fil-rte of mbiguity vlidtion, nd different vlidtion methods could be compred. The two-dimensionl simultion indictes tht OIA performs the best s the mximistion of success rte, yet the computtion burden is extremely hevy. While for the other vlidtion methods, R-rtio, W-rtio nd difference test re preferred. However, higher dimensionl cse gives slightly different results, nd this shows tht in prcticl pplictions, such simultions cnnot show which mbiguity vlidtion method is the best in terms of mbiguity success rte, becuse the mbiguity success rte is not the correct rte of mbiguities, which hs lso been shown in the prcticl exmple discussed in Tble 8. Therefore, for the IA bsed methods, there re some issues tht need to be emphsised: (1) even if we cn determine the criticl vlue from given fil-rte by using look-up tble or simultion, the relibility of both the look-up tble nd the simultion needs to be pid more ttention, s they re lwys ffected by the smple size. In order to obtin relible results, the smple size should be s lrge s possible, which is, however, contrdictory to the computtionl efficiency. Numericl results hve shown tht simultions re Tble 7 The performnce of the non-ia bsed methods Non-IA bsed methods F-rtio R-rtio W -rtio W s -rtio D-test P-test Criticl vlues No. of ccepted solutions Performnce in three specil epochs* WA WA CR CR CR CR *The best solutions from those three epochs re not correct. Tble 8 The performnce of the IA bsed methods IA bsed methods FTIA RTIA WTIA DTIA PTIA EIA OIA Criticl vlues (P f ) No. of ccepted solutions P s P u Performnce in three specil epochs* WA WA CR CR CR CR WA *The best solutions from those three epochs re not correct. Survey Review 2012 VOL 44 NO
9 extremely time consuming nd not stble, nd re ccordingly not pplicble in prctice; (2) it should be noted tht the success rte defined under the frmework of IA theory is not equl to correct rte, which is unknown during mbiguity resolution. The rel GPS dtset hs been used to evlute both the conventionl nd IA bsed methods with the ground truth. In terms of mbiguity correct rte, the W -rtio hs outperformed ll the other vlidtion methods. Moreover, it should be stressed tht the current mbiguity vlidtion methods endevour to discriminte the best integer mbiguity from the others, wheres the correctness of the resolved best integer mbiguity is only vlidted indirectly. Therefore, their performnces in vlidting the correctness of the mbiguities should be evluted by using the ground truth in vrious ppliction scenrios. Acknowledgements The first uthor is PhD student sponsored by Chinese Scholrship Council (CSC) for his studies t the University of New South Wles. References 1. Brd, W., A Testing Procedure for Use in Geodetic Networks. Publictions on Geodesy, 2/5, Netherlnds Geodetic Commission, Delft. 2. Chen, Y., An Approch to Vlidte the Resolved Ambiguities in GPS Rpid Positioning. Proceedings of the Interntionl Symposium on Kinemtic Systems in Geodesy, Geomtics nd Nvigtion, 3 6 June, Bnff, Alt, Cnd. 3. Counselmn, C. C. nd Abbot, R. I., Method of Resolving Rdio Phse Ambiguity in Stellite Orbit Determintion. Journl Geophysicl Reserch, 94(B6): Euler, H. J. nd Schffrin B., On Mesure for the Discernbility between Different Ambiguity Solutions in the Sttic- Kinemtic GPS Mode. IAG Symposi no 107, Kinemtic Systems in Geodesy, Surveying, nd Remote Sensing, Springer, Berlin/ Heidelberg/New York: Frei, E. nd Beutler, G., Rpid Sttic Positioning Bsed on the Fst Ambiguity Resolution Approch FARA: Theory nd First Results. Mnuscript Geodetic, 15(4): Hn, S. nd Rizos, C., Vlidtion nd Rejection Criteri for Integer Lest-Squres Estimtion. Survey Review, 33(260): Hn, S., Qulity Control Issues Relting to Instntneous Ambiguity Resolution for Rel-Time GPS Kinemtic Positioning. Journl of Geodesy, 71(6): Hssibi, A. nd Boyd, S., Integer Prmeter Estimtion in Liner Models with Applictions to GPS. IEEE Trnsctions on Signl Processing, 46(11): Ji, S. Y., Chen, W., Ding, X. L., Chen, Y. Q., Zho, C. M. nd Hu, C. W., Ambiguity Vlidtion with Combined Rtio Test nd Ellipsoidl Integer Aperture Estimtor. Journl of Geodesy, 84(8): Leick, A., GPS Stellite Surveying, New York, John Wiley & Sons, 3rd edition. 11. Tksu, T. nd Ysud, A., Development of the Low-Cost RTK-GPS Receiver With n Open Source Progrm Pckge RTKLIB. Proceedings of the Interntionl Symposium on GPS/ GNSS, 4 6 November, Jeju, Kore. 12. Teunissen, P. J. G., Lest-Squres Estimtion of the Integer GPS Ambiguities. Invited lecture, Section IV, Theory nd Methodology, Proceedings of the IAG Generl Meeting, 8 15 August, Beijing, Chin. 13. Teunissen, P. J. G., A New Method for Fst Crrier Phse Ambiguity Estimtion. Proceedings of the IEEE Position Loction nd Nvigtion Symposium, April, Ls Vegs, NV: Teunissen, P. J. G., The Lest-Squres Ambiguity Decorreltion Adjustment: A Method for Fst GPS Integer Ambiguity Estimtion. Journl of Geodesy, 70(1 2): Teunissen, P. J. G., Success Probbility of Integer GPS Ambiguity Rounding nd Bootstrpping. Journl of Geodesy, 72(10): Teunissen, P. J. G., An Optimlity Property of the Integer Lest-Squres Estimtor. Journl of Geodesy, 73(11): Teunissen, P. J. G., The Prmeter Distributions of the Integer GPS Model. Journl of Geodesy, 76(1): Teunissen, P. J. G., Integer Aperture GNSS Ambiguity Resolution. Artificil Stellites, 38(3): Teunissen, P. J. G., 2003b. A Crrier Phse Ambiguity Estimtor with Esy-to-Evlute Fil-Rte. Artificil Stellites, 38(3): Teunissen, P. J. G., Optiml Integer Aperture Estimtion. Artificil Stellite, to be published. 21. Teunissen, P. J. G., On GNSS Ambiguity Acceptnce Tests. Proceedings of the Interntionl Globl Nvigtion Stellite Systems Society IGNSS Symposium 2007, 4 6 December, Sydney, NSW, Austrli. 22. Teunissen, P. J. G., The GNSS Ambiguity Rtio-Test Revisited: A Better Wy of Using It. Survey Review, 41(312): Tiberius, C. C. J. M. nd de Jonge, P. J., Fst Positioning Using the LAMBDA Method. Proceedings of DSNS-95, April, Bergen, Norwy, no Wng, J., Stewrt, M.P. nd Tskiri, M., A Discrimintion Test Procedure for Ambiguity Resolution On-the-fly. Journl of Geodesy, 72(11), Wng, J., Stewrt, M. P. nd Tskiri, M., A Comprtive Study of the Integer Ambiguity Vlidtion Procedures. Erth, Plnets & Spce, 52(10): Wng, J., Stochstic Modeling for RTK GPS/Glonss Positioning. Journl of the US Institute of Nvigtion, 46(4): Verhgen, S., Integer Ambiguity Vlidtion: An Open Problem? GPS Solutions, 8(1): Verhgen, S., The GNSS Integer Ambiguities: Estimtion nd Vlidtion. PhD thesis, Publictions on Geodesy, 58, Netherlnd Geodetic Commission, Delft. 29. Verhgen, S., On the Probbility Density Function of the GNSS Ambiguity Residuls. GPS Solutions, 10(1): Survey Review 2012 VOL 44 NO 326
The steps of the hypothesis test
ttisticl Methods I (EXT 7005) Pge 78 Mosquito species Time of dy A B C Mid morning 0.0088 5.4900 5.5000 Mid Afternoon.3400 0.0300 0.8700 Dusk 0.600 5.400 3.000 The Chi squre test sttistic is the sum of
More informationData Assimilation. Alan O Neill Data Assimilation Research Centre University of Reading
Dt Assimiltion Aln O Neill Dt Assimiltion Reserch Centre University of Reding Contents Motivtion Univrite sclr dt ssimiltion Multivrite vector dt ssimiltion Optiml Interpoltion BLUE 3d-Vritionl Method
More informationTests for the Ratio of Two Poisson Rates
Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson
More information13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationHybrid Group Acceptance Sampling Plan Based on Size Biased Lomax Model
Mthemtics nd Sttistics 2(3): 137-141, 2014 DOI: 10.13189/ms.2014.020305 http://www.hrpub.org Hybrid Group Acceptnce Smpling Pln Bsed on Size Bised Lomx Model R. Subb Ro 1,*, A. Ng Durgmmb 2, R.R.L. Kntm
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationAcceptance Sampling by Attributes
Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationNon-Linear & Logistic Regression
Non-Liner & Logistic Regression If the sttistics re boring, then you've got the wrong numbers. Edwrd R. Tufte (Sttistics Professor, Yle University) Regression Anlyses When do we use these? PART 1: find
More informationStudent Activity 3: Single Factor ANOVA
MATH 40 Student Activity 3: Single Fctor ANOVA Some Bsic Concepts In designed experiment, two or more tretments, or combintions of tretments, is pplied to experimentl units The number of tretments, whether
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationCBE 291b - Computation And Optimization For Engineers
The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
More informationPredict Global Earth Temperature using Linier Regression
Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationEstimation of Binomial Distribution in the Light of Future Data
British Journl of Mthemtics & Computer Science 102: 1-7, 2015, Article no.bjmcs.19191 ISSN: 2231-0851 SCIENCEDOMAIN interntionl www.sciencedomin.org Estimtion of Binomil Distribution in the Light of Future
More informationDriving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d
Interntionl Industril Informtics nd Computer Engineering Conference (IIICEC 15) Driving Cycle Construction of City Rod for Hybrid Bus Bsed on Mrkov Process Deng Pn1,, Fengchun Sun1,b*, Hongwen He1, c,
More informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationMonte Carlo method in solving numerical integration and differential equation
Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationCredibility Hypothesis Testing of Fuzzy Triangular Distributions
666663 Journl of Uncertin Systems Vol.9, No., pp.6-74, 5 Online t: www.jus.org.uk Credibility Hypothesis Testing of Fuzzy Tringulr Distributions S. Smpth, B. Rmy Received April 3; Revised 4 April 4 Abstrct
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationCHM Physical Chemistry I Chapter 1 - Supplementary Material
CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion
More information1 Online Learning and Regret Minimization
2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in
More informationTesting categorized bivariate normality with two-stage. polychoric correlation estimates
Testing ctegorized bivrite normlity with two-stge polychoric correltion estimtes Albert Mydeu-Olivres Dept. of Psychology University of Brcelon Address correspondence to: Albert Mydeu-Olivres. Fculty of
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationContinuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht
More informationDiscrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17
EECS 70 Discrete Mthemtics nd Proility Theory Spring 2013 Annt Shi Lecture 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion,
More information4.4 Areas, Integrals and Antiderivatives
. res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationNumerical Integration
Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl
More informationPhysics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018
Physics 201 Lb 3: Mesurement of Erth s locl grvittionl field I Dt Acquisition nd Preliminry Anlysis Dr. Timothy C. Blck Summer I, 2018 Theoreticl Discussion Grvity is one of the four known fundmentl forces.
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
More information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
More informationA signalling model of school grades: centralized versus decentralized examinations
A signlling model of school grdes: centrlized versus decentrlized exmintions Mri De Pol nd Vincenzo Scopp Diprtimento di Economi e Sttistic, Università dell Clbri m.depol@unicl.it; v.scopp@unicl.it 1 The
More informationDuality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.
Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we
More informationChapter 5 : Continuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll
More information19 Optimal behavior: Game theory
Intro. to Artificil Intelligence: Dle Schuurmns, Relu Ptrscu 1 19 Optiml behvior: Gme theory Adversril stte dynmics hve to ccount for worst cse Compute policy π : S A tht mximizes minimum rewrd Let S (,
More information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS
ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:
More informationProbability Distributions for Gradient Directions in Uncertain 3D Scalar Fields
Technicl Report 7.8. Technische Universität München Probbility Distributions for Grdient Directions in Uncertin 3D Sclr Fields Tobis Pfffelmoser, Mihel Mihi, nd Rüdiger Westermnn Computer Grphics nd Visuliztion
More informationSection 11.5 Estimation of difference of two proportions
ection.5 Estimtion of difference of two proportions As seen in estimtion of difference of two mens for nonnorml popultion bsed on lrge smple sizes, one cn use CLT in the pproximtion of the distribution
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More informationLecture 3 Gaussian Probability Distribution
Introduction Lecture 3 Gussin Probbility Distribution Gussin probbility distribution is perhps the most used distribution in ll of science. lso clled bell shped curve or norml distribution Unlike the binomil
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More informationOrdinary differential equations
Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties
More informationDiscrete Mathematics and Probability Theory Summer 2014 James Cook Note 17
CS 70 Discrete Mthemtics nd Proility Theory Summer 2014 Jmes Cook Note 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion, y tking
More informationFor the percentage of full time students at RCC the symbols would be:
Mth 17/171 Chpter 7- ypothesis Testing with One Smple This chpter is s simple s the previous one, except it is more interesting In this chpter we will test clims concerning the sme prmeters tht we worked
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More information1B40 Practical Skills
B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need
More informationMIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:
1 2 MIXED MODELS (Sections 17.7 17.8) Exmple: Suppose tht in the fiber breking strength exmple, the four mchines used were the only ones of interest, but the interest ws over wide rnge of opertors, nd
More informationChapters 4 & 5 Integrals & Applications
Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions
More informationChapter 2 Fundamental Concepts
Chpter 2 Fundmentl Concepts This chpter describes the fundmentl concepts in the theory of time series models In prticulr we introduce the concepts of stochstic process, men nd covrince function, sttionry
More informationp-adic Egyptian Fractions
p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction
More informationDiscrete Least-squares Approximations
Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve
More informationUNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction
Lesson : Logrithmic Functions s Inverses Prerequisite Skills This lesson requires the use of the following skills: determining the dependent nd independent vribles in n exponentil function bsed on dt from
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
More informationMATH 144: Business Calculus Final Review
MATH 144: Business Clculus Finl Review 1 Skills 1. Clculte severl limits. 2. Find verticl nd horizontl symptotes for given rtionl function. 3. Clculte derivtive by definition. 4. Clculte severl derivtives
More informationMath 426: Probability Final Exam Practice
Mth 46: Probbility Finl Exm Prctice. Computtionl problems 4. Let T k (n) denote the number of prtitions of the set {,..., n} into k nonempty subsets, where k n. Argue tht T k (n) kt k (n ) + T k (n ) by
More informationNew data structures to reduce data size and search time
New dt structures to reduce dt size nd serch time Tsuneo Kuwbr Deprtment of Informtion Sciences, Fculty of Science, Kngw University, Hirtsuk-shi, Jpn FIT2018 1D-1, No2, pp1-4 Copyright (c)2018 by The Institute
More informationImproved Results on Stability of Time-delay Systems using Wirtinger-based Inequality
Preprints of the 9th World Congress he Interntionl Federtion of Automtic Control Improved Results on Stbility of ime-dely Systems using Wirtinger-bsed Inequlity e H. Lee Ju H. Prk H.Y. Jung O.M. Kwon S.M.
More informationHere we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system.
Section 24 Nonsingulr Liner Systems Here we study squre liner systems nd properties of their coefficient mtrices s they relte to the solution set of the liner system Let A be n n Then we know from previous
More informationLecture 21: Order statistics
Lecture : Order sttistics Suppose we hve N mesurements of sclr, x i =, N Tke ll mesurements nd sort them into scending order x x x 3 x N Define the mesured running integrl S N (x) = 0 for x < x = i/n for
More informationA New Statistic Feature of the Short-Time Amplitude Spectrum Values for Human s Unvoiced Pronunciation
Xiodong Zhung A ew Sttistic Feture of the Short-Time Amplitude Spectrum Vlues for Humn s Unvoiced Pronuncition IAODOG ZHUAG 1 1. Qingdo University, Electronics & Informtion College, Qingdo, 6671 CHIA Abstrct:
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationSolution for Assignment 1 : Intro to Probability and Statistics, PAC learning
Solution for Assignment 1 : Intro to Probbility nd Sttistics, PAC lerning 10-701/15-781: Mchine Lerning (Fll 004) Due: Sept. 30th 004, Thursdy, Strt of clss Question 1. Bsic Probbility ( 18 pts) 1.1 (
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk bout solving systems of liner equtions. These re problems tht give couple of equtions with couple of unknowns, like: 6 2 3 7 4
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationReinforcement learning II
CS 1675 Introduction to Mchine Lerning Lecture 26 Reinforcement lerning II Milos Huskrecht milos@cs.pitt.edu 5329 Sennott Squre Reinforcement lerning Bsics: Input x Lerner Output Reinforcement r Critic
More informationTime Truncated Two Stage Group Sampling Plan For Various Distributions
Time Truncted Two Stge Group Smpling Pln For Vrious Distributions Dr. A. R. Sudmni Rmswmy, S.Jysri Associte Professor, Deprtment of Mthemtics, Avinshilingm University, Coimbtore Assistnt professor, Deprtment
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationSection 14.3 Arc Length and Curvature
Section 4.3 Arc Length nd Curvture Clculus on Curves in Spce In this section, we ly the foundtions for describing the movement of n object in spce.. Vector Function Bsics In Clc, formul for rc length in
More informationPART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point.
PART MULTIPLE CHOICE Circle the pproprite response to ech of the questions below. Ech question hs vlue of point.. If in sequence the second level difference is constnt, thn the sequence is:. rithmetic
More informationRiemann is the Mann! (But Lebesgue may besgue to differ.)
Riemnn is the Mnn! (But Lebesgue my besgue to differ.) Leo Livshits My 2, 2008 1 For finite intervls in R We hve seen in clss tht every continuous function f : [, b] R hs the property tht for every ɛ >
More informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
More informationThe Predom module. Predom calculates and plots isothermal 1-, 2- and 3-metal predominance area diagrams. Predom accesses only compound databases.
Section 1 Section 2 The module clcultes nd plots isotherml 1-, 2- nd 3-metl predominnce re digrms. ccesses only compound dtbses. Tble of Contents Tble of Contents Opening the module Section 3 Stoichiometric
More informationUnit #9 : Definite Integral Properties; Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More informationFig. 1. Open-Loop and Closed-Loop Systems with Plant Variations
ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses
More informationAdvanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004
Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More informationNatural examples of rings are the ring of integers, a ring of polynomials in one variable, the ring
More generlly, we define ring to be non-empty set R hving two binry opertions (we ll think of these s ddition nd multipliction) which is n Abelin group under + (we ll denote the dditive identity by 0),
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationGenetic Programming. Outline. Evolutionary Strategies. Evolutionary strategies Genetic programming Summary
Outline Genetic Progrmming Evolutionry strtegies Genetic progrmming Summry Bsed on the mteril provided y Professor Michel Negnevitsky Evolutionry Strtegies An pproch simulting nturl evolution ws proposed
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 7.. Introduction to Integrtion. 7. Integrl Clculus As ws the cse with the chpter on differentil
More informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
More informationImproper Integrals. Type I Improper Integrals How do we evaluate an integral such as
Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
More informationChapter 3 Polynomials
Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling
More informationConstruction and Selection of Single Sampling Quick Switching Variables System for given Control Limits Involving Minimum Sum of Risks
Construction nd Selection of Single Smpling Quick Switching Vribles System for given Control Limits Involving Minimum Sum of Risks Dr. D. SENHILKUMAR *1 R. GANESAN B. ESHA RAFFIE 1 Associte Professor,
More informationAn approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95
An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech
More informationDOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES
DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m
More information12 TRANSFORMING BIVARIATE DENSITY FUNCTIONS
1 TRANSFORMING BIVARIATE DENSITY FUNCTIONS Hving seen how to trnsform the probbility density functions ssocited with single rndom vrible, the next logicl step is to see how to trnsform bivrite probbility
More information