首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 671 毫秒
1.
This paper proposes an original method of constructing a criterion matrix for the optimal design of control networks by means of the contraction of the eigenvalues and the rotation of the eigenvectors of a covariance matrix. The rotation of the eigenvectors is carried out by two different methods, that is through independent rotations of the pairs of the eigenvector matrix components and through a procrustean transfomation of the same matrix. The required criterion is that the “essential eigenvector components” must be as orthogonal as possible to the predicted direction of deformation. A Second Order Design problem is then resolved, that is the optimization of the precision of the observations of a local free distance network to be constructed for the study of recent crustal movements in the seismogenetic area of Friuli (Italy). The results obtained show a high degree of correspondence for the two proposed methods of rotation if limited rotations of the eigenvectors are considered.  相似文献   

2.
工业测量中一种二次曲面的拟合方法   总被引:4,自引:0,他引:4  
针对工业测量中二次曲面的拟合问题,阐述了二次多项式的标准化、二次项系数矩阵的特征值和特征向量求取、测量坐标与标准坐标之间的转换关系、特征向量与旋转角的关系等问题,提出了按特征值、旋转角和平移量为参数拟合二次曲面的拟舍方法。算例表明了模型的准确性。  相似文献   

3.
The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (YE Y = (XE X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E Y and E X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.  相似文献   

4.
The mixed spherical map projections of equiareal, cylindric type are based upon the Lambert projection and the sinusoidal Sanson–Flamsteed projection. These cylindric and pseudo-cylindric map projections of the sphere are generalized to the ellipsoid of revolution (biaxial ellipsoid). They are used in consequence by two lemmas to generate a horizontal and a vertical weighted mean of equiareal cylindric map projections of the ellipsoid of revolution. Its left–right deformation analysis via further results leads to the left–right principal stretches/eigenvalues and left–right eigenvectors/eigenspace, as well as the maximal left–right angular distortion for these new mixed cylindric map projections of ellipsoidal type. Detailed illustrations document the cartographic synergy of mixed cylindric map projections. Received: 23 April 1996 / Accepted: 19 April 1997  相似文献   

5.
Y. Yuan  J. Ou 《Journal of Geodesy》2001,75(7-8):438-447
 Ionospheric variation may be considered as a stationary time series under quiet conditions. However, the disturbance of a stationary random process from stationarity results in the bias of corresponding samples from the stationary observations, and in the change of statistical model parameters of the process. From a general mathematical aspect, a new method is presented for monitoring ionospheric variations, based on the characteristic of time-series observation of GPS, and an investigation of the statistical properties of the estimated auto-covariance of the random ionospheric delay when changing the number of samples in the time series is carried out. A preliminary scheme for monitoring ionospheric delays is proposed. Received: 18 August 2000 / Accepted: 12 April 2001  相似文献   

6.
 A comparison was made between two methods for gravity field recovery from orbit perturbations that can be derived from global positioning system satellite-to-satellite tracking observations of the future European gravity field mission GOCE (Gravity Field and Steady-State Ocean Circulation Explorer). The first method is based on the analytical linear orbit perturbation theory that leads under certain conditions to a block-diagonal normal matrix for the gravity unknowns, significantly reducing the required computation time. The second method makes use of numerical integration to derive the observation equations, leading to a full set of normal equations requiring powerful computer facilities. Simulations were carried out for gravity field recovery experiments up to spherical harmonic degree and order 80 from 10 days of observation. It was found that the first method leads to large approximation errors as soon as the maximum degree surpasses the first resonance orders and great care has to be taken with modeling resonance orbit perturbations, thereby loosing the block-diagonal structure. The second method proved to be successful, provided a proper division of the data period into orbital arcs that are not too long. Received: 28 April 2000 / Accepted: 6 November 2000  相似文献   

7.
A least squares solution of sequential array algebra observation equations is derived using spectral decomposition of the normal equation matrix in terms of array algebra. A new direct solution for computation of eigenvectors is derived using the theory of general matrix inverses.  相似文献   

8.
Although the assumption of independence among interaction flows frequently is engaged in spatial interaction modeling, in many circumstances it leads to misspecified models and incorrect inferences. An informed approach is to explicitly incorporate an assumed relationship structure among the interaction flows, and to explicitly model the network autocorrelation. This paper illustrates such an approach in the context of U.S. interstate migration flows. Behavioral assumptions, similar to those of the intervening opportunities or the competing destinations concepts, exemplify how to specify network flows that are related to particular origin–destination combinations. The stepwise incorporation of eigenvectors, which are extracted from a network link matrix, captures the network autocorrelation in a Poisson regression model specification context. Spatial autocorrelation in Poisson regression is measured by the test statistic of Jacqmin-Gadda et al. (Stat Med 16(11):1283–1297, 1997). Results show that estimated regression parameters in the spatial filtering interaction model become more intuitively interpretable.
Yongwan ChunEmail:
  相似文献   

9.
 As either the spatial resolution or the spatial scale for a geographic landscape increases, both latent spatial dependence and spatial heterogeneity also will tend to increase. In addition, the amount of georeferenced data that results becomes massively large. These features of high spatial resolution hyperspectral data present several impediments to conducting a spatial statistical analysis of such data. Foremost is the requirement of popular spatial autoregressive models to compute eigenvalues for a row-standardized geographic weights matrix that depicts the geographic configuration of an image's pixels. A second drawback arises from a need to account for increased spatial heterogeneity. And a third concern stems from the usefulness of marrying geostatistical and spatial autoregressive models in order to employ their combined power in a spatial analysis. Research reported in this paper addresses all three of these topics, proposing successful ways to prevent them from hindering a spatial statistical analysis. For illustrative purposes, the proposed techniques are employed in a spatial analysis of a high spatial resolution hyperspectral image collected during research on riparian habitats in the Yellowstone ecosystem. Received: 25 February 2001 / Accepted: 2 August 2001  相似文献   

10.
A proper perturbation theory of a mathematical model and the quantities derived by means of least-squares adjustments is indispensable if the results have to be interpreted in a wider context. The sensitivity of some characteristic results of least-squares adjustments such as the estimated values of the parameters and their variance–covariance matrix due to imminent uncertainties of the stochastic model is discussed in detail. Linearizations are used with rigorous error measures and interval mathematics. Numerical examples conclude the investigations. Received: 27 December 1997 / Accepted: 19 April 1999  相似文献   

11.
在利用函数映射计算模型间对应关系时,提出了一种校准三维几何模型之间基矩阵的新方法,将模型间对应关系的构建转化为由模型特征函数构建的基矩阵之间的校准运算。首先计算三维模型的Laplace算子,获得模型的特征值和特征向量,并利用所得到的特征向量构建基矩阵;其次,提出了协方差的最小值校准算法,以计算模型基矩阵之间的校准矩阵 S ,并用矩阵 S 对两个模型的函数基进行校准;最后,计算校准模型所有点的高斯曲率来采样源模型尖端特征点,并在校准后的目标模型上遍历所有点,以寻求最优对应点来构建等距变换(或近似等距变换)的三维模型间的对应关系。通过计算采样点与最优对应点的测地错误,以衡量所提算法的匹配准确率。实验结果表明,与已有算法相比,本算法可以较为准确地构建两个或多个模型间的对应关系,同时也克服了模型自身对称性影响对应关系计算的问题。  相似文献   

12.
S. Han 《Journal of Geodesy》1997,71(6):351-361
An integrated method for the instantaneous ambiguity resolution using dual-frequency precise pseudo-range and carrier-phase observations is suggested in this paper. The algorithm combines the search procedures in the coordinate domain, the observation domain and the estimated ambiguity domain (and therefore benefits from the integration of their most positive elements). A three-step procedure is then proposed to enhance the reliability of the ambiguity resolution by: (1) improving the stochastic model for the double-differenced functional model in real time; (2) refining the criteria which distinguish the integer ambiguity set that generates the minimum quadratic form of residuals from that corresponding to the second minimum one; and (3) developing a fault detection and adaptation procedure. Three test scenarios were considered, one static baseline (11.3 km) and two kinematic experiments (baseline lengths from 5.2 to 13.7 km). These showed that the mean computation time for one epoch is less than 0.1 s, and that the success rate reaches 98.4% (compared to just 68.4% using standard ratio tests). Received: 5 June 1996; Accepted: 16 January 1997  相似文献   

13.
In satellite data analysis, one big advantage of analytical orbit integration, which cannot be overestimated, is missed in the numerical integration approach: spectral analysis or the lumped coefficient concept may be used not only to design efficient algorithms but overall for much better insight into the force-field determination problem. The lumped coefficient concept, considered from a practical point of view, consists of the separation of the observation equation matrix A=BT into the product of two matrices. The matrix T is a very sparse matrix separating into small block-diagonal matrices connecting the harmonic coefficients with the lumped coefficients. The lumped coefficients are nothing other than the amplitudes of trigonometric functions depending on three angular orbital variables; therefore, the matrix N=B T B will become for a sufficient length of a data set a diagonal dominant matrix, in the case of an unlimited data string length a strictly diagonal one. Using an analytical solution of high order, the non-linear observation equations for low–low SST range data can be transformed into a form to allow the application of the lumped concept. They are presented here for a second-order solution together with an outline of how to proceed with data analysis in the spectral domain in such a case. The dynamic model presented here provides not only a practical algorithm for the parameter determination but also a simple method for an investigation of some fundamental questions, such as the determination of the range of the subset of geopotential coefficients which can be properly determined by means of SST techniques or the definition of an optimal orbital configuration for particular SST missions. Numerical results have already been obtained and will be published elsewhere. Received: 15 January 1999 / Accepted: 30 November 1999  相似文献   

14.
The ellipsoidal Stokes problem is one of the basic boundary-value problems for the Laplace equation which arises in physical geodesy. Up to now, geodecists have treated this and related problems with high-order series expansions of spherical and spheroidal (ellipsoidal) harmonics. In view of increasing computational power and modern numerical techniques, boundary element methods have become more and more popular in the last decade. This article demonstrates and investigates the nullfield method for a class of Robin boundary-value problems. The ellipsoidal Stokes problem belongs to this class. An integral equation formulation is achieved, and existence and uniqueness conditions are attained in view of the Fredholm alternative. Explicit expressions for the eigenvalues and eigenfunctions for the boundary integral operator are provided. Received: 22 October 1996 / Accepted: 4 August 1997  相似文献   

15.
 Harmonic coefficients of the 2nd degree are separated into the invariant quantitative (the 2nd-degree variance) and the qualitative (the standardized harmonic coefficients) characteristics of the behavior of the potential V 2(t). On this basis the evolution of the Earth's dynamical figure is described as a solution of the time-dependent eigenvalues–eigenvectors problem in the canonical form. Such a canonical quadratic form is defined only by temporal variations of the harmonic coefficients and always remains finite, even within an infinite time interval. An additional condition for the correction or the determination of temporal variations of the 2nd degree is obtained. Temporal variations of the fully normalized sectorial harmonic coefficients are estimated in addition to ˙Cˉ 20, ˙Cˉ 21, and ˙Sˉ 21 of the EGM96 gravity model. In addition, a non-linear hyperbolic model for 2m (t), 2m (t) is constructed. The trigonometric form of the hyperbolic model leads to the consideration of the potential V 2(ψ) instead of V 2(t) within the closed interval −π/2≤ψ≤+π/2. Thus, it is possible to evaluate the global trend of V 2(t), the Earth's principal axes and the differences of the moments of inertia within the whole infinite time interval. Received: 25 September 1998 / Accepted: 28 June 2000  相似文献   

16.
A reliable and accurate gradiometer calibration is essential for the scientific return of the gravity field and steady-state ocean circulation explorer (GOCE) mission. This paper describes a new method for external calibration of the GOCE gradiometer accelerations. A global gravity field model in combination with star sensor quaternions is used to compute reference differential accelerations, which may be used to estimate various combinations of gradiometer scale factors, internal gradiometer misalignments and misalignments between star sensor and gradiometer. In many aspects, the new method is complementary to the GOCE in-flight calibration. In contrast to the in-flight calibration, which requires a satellite-shaking phase, the new method uses data from the nominal measurement phases. The results of a simulation study show that gradiometer scale factors can be estimated on a weekly basis with accuracies better than 2 × 10−3 for the ultrasensitive and 10−2 for the less sensitive axes, which is compatible with the requirements of the gravity gradient error. Based on a 58-day data set, scale factors are found that can reduce the errors of the in-flight-calibrated measurements. The elements of the complete inverse calibration matrix, representing both the internal gradiometer misalignments and scale factors, can be estimated with accuracies in general better than 10−3.  相似文献   

17.
Let there be given a twodimensional symmetric rank two tensor of random type (examples:strain, stress) which is either directly observed or indirectly estimated from observations by an adjustment procedure. Under the assumption of normalityof tensor components we compute the joint probability density functionas well as the marginal probability density functionsof its eigenspectra (eigenvalues) and eigendirections (orientation parameters). Due to the nonlinearity of the relation between eigenspectra-eigendirections and the random tensor components, via the inverse nonlinear error propagationbiases and aliases of their first and centralized second moments (mean value, variance-covariance) are expressed in terms of Jacobianand Hessianmatrices. The joint probability density function and the first and second moments thus form the fundamental of hypothesis testing and qualify control of eigenspectra (eigenvalues, principal components) and eigendirections (orientation parameters, eigenvectors, principial direction) of a twodimensional, symmetric rank two random tensor.  相似文献   

18.
To apply the least squares method for the interpolation of harmonic functions is a common practice in Geodesy. Since the method of least squares can be applied only to overdetermined problem, the interpolation problem which is always under-determined, is often reduced to an overdetermined form by truncating a series of spherical harmonics. When the data points are the knots of a regular grid it is easy to see that the estimated harmonic coefficients converge to the correct theoretical values, but when the observation density is not constant a significant bias is introduced. The result is obtained by assuming that the number of observations tends to infinity with points sampled from a given distribution. Under the same conditions it is shown that quadrature and “collocation-like” formulas displays a statistically consistent behaviour.  相似文献   

19.
It is shown that also in a rank deficient Gauss-Markov model higher weights of the observations automatically improve the precision of the estimated parameters as long as they are computed in thesame datum. However, the amount of improvement in terms of the trace of the dispersion matrix isminimum for the so-called “free datum” which corresponds to the pseudo-inverse normal equations matrix. This behaviour together with its consequences is discussed by an example with special emphasis on geodetic networks for deformation analysis.  相似文献   

20.
The spline interpolation technique is applied to estimate locally the radial component of a planetary gravity field from residual acceleration data along a special direction, the direction of observation from the Earth. After the presentation of the theoretical framework, the method is tested on synthetic and real data in the case of Venus. It is shown that this spline technique can be used succesfully to build local models of radial gravity fields at the planet surface. Received: 13 March 1997 / Accepted: 17 November 1998  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号