估计从UUV航行参数的混合坐标系下的EKF算法

在主从式UUV 协作系统中，由于定位和导航的需要，要求尽快估计出从UUV 的航行参数，但通常所用的递推最小二乘（RLS）算法，其初始方位测量对滤波结果影响大且存在收敛速度慢、计算精度低的缺点，难以满足应用需求，而推广卡尔曼滤波（EKF）算法能较好地克服上述问题。在直角坐标系下（CEKF），方位信息与距离信息相互耦合导致初始振荡剧烈，改为混合坐标系（MEKF）后问题得到了极大的改善。最后，通过仿真及现场试验验证了此改进方法的有效性。

基于光谱分析的主动式距离图像获取方法

Rainbow三维摄像机是一种基于光谱分析的快速三维信息获取方法。该方法利用连续变化的彩色光谱照射景物 ,彩色CCD摄像机摄取的景物图像将呈现有规律的颜色变化 ,而且不同的颜色 (波长 )构成了不同的空间颜色面。通过标定这些颜色面和摄像机成象模型 ,即可计算出图像中各点的三维坐标。

基于光谱分析的主动式距离图象获取方法

Ｒａｉｎｂｏｗ 三维摄像机是一种基于光谱分析的快速三维信息获取方法。该方法利用连续变化的彩色光谱照射景物，彩色ＣＣＤ 摄像机摄取的景物图象将呈现有规律的颜色变化，而且不同的颜色（波长）构成了不同的空间颜色面。通过标定这些颜色面和摄像机成象模型，即可计算出图象中各点的三维坐标。该文重点讨论实现该方法的标定技术和颜色分类技术，最后给出实验结果。

并联机器人机床运动学正、逆解

基于一般ＳＴＥＷＡＲＴ机构研制的并联机器人机床是新一代智能化金属切削加工机床．然而，机床的运动学位置正、逆解呈强非线性，求解困难．出于机床精度的需要，本研究的模型样机在结构上采用了滚珠丝杠传动，因此又带来了关节运动耦合，导致机床运动学位置正、逆解求解更加复杂．利用运动学等效的原则，引入整机等效串联机构及分支等效串联机构，以等效广义坐标为中间变量建立机床运动学正、逆解求解迭代算法．仿真与控制实验表明，该算法具有收索速度快便于实际应用等特点。

基于运动耦合结构的并联机器人的运动学正逆解

根据运动学等效的原则，在并联机器人中引入等效串联机器人及分支等效串联机器人，以等效广义坐标为中间变量建立并联机器人运动学正道解求解算法。该算法能有效处理结构带来的运动耦合，并且规划的软件具有自动生成迭代初始点、避免多解性以及便于实际应用等特点，从而为并联机器人的结构设计与创新提供了理论支持。

基于动力学约束的机械臂最优关节轨迹规划

本文提出了基于机械臂关节驱动力矩约束方程规划其关节最优运动轨迹的一种有效方法.该方法运用矩阵范数理论简化机械臂的动力学约束方程;在机械臂的关节空间内采用归一化的无因次量运用非线性规划法优化其运动轨迹.将所规划的无因次量轨迹方程作为机械臂产生实际运动轨迹的发生器,通过给定机械臂各运动段的起始和终止关节坐标,由系统的动力学约束方程计算出整个运动段所允许的最短运行时间,即生成所期望的运动轨迹.本文的轨迹规划方法计算效率高,可用于在线轨迹规划,文中通过算例证实了该方法的实用性.

一种高智能化的通用自动胶片判读仪

自动通用胶片判读仪是一种高智能化的精密光学测试设备,采用了计算机控制飞点扫描技术、精密光学测量技术、图象跟踪测量与信息处理技术。飞点管分辨率达4096×4096象元,通过光学系统胶片上获得6.55μm高分辨率,飞点扫描方式灵活多样且可随意控制,通用于目前我国靶场所有的电影经纬仪和高速摄影机35mm胶片的数据判读。具有自动判读和半自动判读两种工作模式,自动判读的速度为5帧/秒,自动判读的精度为σ=±0.011mm,半自动判读的精度为σ=±0.009mm,测量数据可以记录、打印和显示。

多关节机器人工作空间的分析与评价方法

本文运用 H-D 变换的基本原理,结合极坐标变换导出了产生n自由度多关节机器人工作空间的递推算法,当给定了机器人的结构尺寸,即可将机器人工作空间在一特定平面内的边界图形用计算机打出并计算出机器人工作空间容积。在本文的另一部分介绍了两种计算机器人工作空间的性能指标,最后用几个机器人的结构参数进行计算和讨论。

景物的三维信息获取及其运动分析

采用模糊熵函数对图象象素分类作出整体最优分类评价,实现了区域分割.利用矩及其函数做为各区域的特征表达,构成以区域为基元的符号特征集并描述图象内容。根据立体图象对间的几何关系,解出各区域(基元)的相对三维坐标。与象索匹配相比较,它可以获得较高精度的三维信息和可描述的景物信息.通过获取不同时空的各区域(基元)三维信息,确定了它们的空间运动状态。联系这些状态,构造出景物中物体间的空间关系和近似模型,实现了对景物的3-D识别和描述。

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

电离层参量M(3000)F2和hmF2的经验正交函数方法建模研究

The ionospheric parameter M(3000)F2 (the so-called transmission factor or the propagation factor) is important not only in practical applications such as frequency planning for radio-communication but also in ionospheric modeling. This parameter is strongly anti-correlated with the ionospheric F2-layer peak height hmF2，a parameter often used as a key anchor point in some widely used empirical models of the ionospheric electron density profile (e.g., in IRI and NeQuick models). Since hmF2 is not easy to obtain from measurements and M(3000)F2 can be routinely scaled from ionograms recorded by ionosonde/digisonde stations distributed globally and its data has been accumulated for a long history, usually the value of hmF2 is calculated from M(3000)F2 using the empirical formula connecting them. In practice, CCIR M(3000)F2 model is widely used to obtain M(3000)F2 value. However, recently some authors found that the CCIR M(3000)F2 model has remarkable discrepancies with the measured M(3000)F2, especially in low-latitude and equatorial regions. For this reason, the International Reference Ionosphere (IRI) research community proposes to improve or update the currently used CCIR M(3000)F2 model. Any efforts toward the improvement and updating of the current M(3000)F2 model or newly development of a global hmF2 model are encouraged. In this dissertation, an effort is made to construct the empirical models of M(3000)F2 and hmF2 based on the empirical orthogonal function (EOF) analysis combined with regression analysis method. The main results are as follows: 1. A single station model is constructed using monthly median hourly values of M(3000)F2 data observed at Wuhan Ionospheric Observatory during the years of 1957–1991 and compared with the IRI model. The result shows that EOF method is possible to use only a few orders of EOF components to represent most of the variance of the original data set. It is a powerful method for ionospheric modeling. 2. Using the values of M(3000)F2 observed by ionosondes distributed globally, data at grids uniformly distributed globally were obtained by using the Kriging interpolation method. Then the gridded data were decomposed into EOF components using two different coordinates: (1) geographical longitude and latitude; (2) modified dip (Modip) and local time. Based on the EOF decompositions of the gridded data under these two coordinates systems, two types of the global M(3000)F2 model are constructed. Statistical analysis showed that the two types of the constructed M(3000)F2 model have better agreement with the observational M(3000)F2 than the M(3000)F2 model currently used by IRI. The constructed models can represent the global variations of M(3000)F2 better. 3. The hmF2 data used to construct the hmF2 model were converted from the observed M(3000)F2 based on the empirical formula connecting them. We also constructed two types of the global hmF2 model using the similar method of modeling M(3000)F2. Statistical analysis showed that the prediction of our models is more accurate than the model of IRI. This demonstrated that using EOF analysis method to construct global model of hmF2 directly is feasible. The results in this thesis indicate that the modeling technique based on EOF expansion combined with regression analysis is very promising when used to construct the global models of M(3000)F2 and hmF2. It is worthwhile to investigate further and has the potential to be used to the global modeling of other ionospheric parameters.

中高层大气对低层大气声重波响应的分析

Acoustic Gravity waves (AGW) play an important role in balancing the atmospheric energy and momentum budget. Propagation of gravity wave in the atmosphere is one of the important factors of changing middle and upper atmosphere and ionosphere. The purpose of this dissertation is to study the propagation of gravity wave in a compression atmosphere whit means of numerical simulation and to analyze the response of middle and upper atmosphere to pulse disturbance from lower atmosphere. This work begins with the establishment of 2-D fully nonlinear compressible atmospheric dynamic model in polar coordinate, which is used ton numerically study gravity wave propagation. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. We also simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model and analyze the data we obtained by using Fourier Transform (FT), Short-time Fourier Transform (STFT) and Empirical Mode Decomposition (EMD) method which is an important part of Hilbert-Huang Transform (HHT). The research content is summarized in the following: 1. By using a two-dimensional full-implicit-continuous-Eulerian (FICE) scheme and taking the atmospheric basic motion equations as the governing equations, a numerical model for nonlinear propagation of acoustic gravity wave disturbance in two-dimensional polar coordinates is solved. 2. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. Results of numerical simulation show that the acoustic gravity wave packets propagate steadily upward and keep its shape well after several periods. 3. We simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model, and obtain the distribution of a certain physical quantity in time and space from earth’s surface to 300km above. The results reveal that the response of ionosphere occurs at a large horizontal distance from the source and the disturbance becomes greater with increasing of height. The situation when the direction of the background wind is opposite to or the same as the direction of disturbed velocity of gravity-wave is studied. The results show that gravity wave propagating against winds is easier than those propagating along winds and the background wind can accelerate gravity wave propagation. Just upon the source, an acoustic wave component with period of 6 min can be found. These images of simulation are similar to observations of the total electron content (TEC) disturbances caused by the great Sumatra-Andaman earthquake on December 26 in 2004. 4. Using the EMD method the disturbed velocity data of a certain physical quantity in time and space can be decomposed into a series of intrinsic mode function (IMF) and a trend mode respectively. The results of EMD reveal impact of the gravity wave frequency under the background winds.

三峡库区堆积层斜坡稳定性分级与应用研究

Considering the lacking of standard for the classification of Accumulative slopes so far, research working was conducted based on the results of geological investigation, data analysis and experiment carried out in Wanzhou. By mean of statistical method and grey system, the author studied in detail inflationary factors to Accumulative slopes. In order to study the mechanism of Rock-Soil Aggregate (RSA), numerical testing method was used. Coordinates in the two and three dimensional space and its corresponding rock fragments in the sample were generated randomly by VB and Particle flow code. After being built the models of RSA with different rock content, uniaxial and triaxial numerical simulation tests were carried out respectively. In order to study the effect of rainfall in Accumulative slopes, in situ infiltration testing had been conducted on site in Wanzhou, Three Gorges Area. Relationship between the infiltration rate and amount of precipitation has been obtained. Eleven factors are considered in the classification of Accumulative slopes in this paper.(1)On the basis of four basic factors and four inducing factors, sum-and-difference method for the classification system has been built. (2)After weight of factors being determined by analytic hierarchy process and membership function of Accumulative slopes stability being built in virtue of fuzzy mathematics, AHP-FM model of Accumulative slopes stability has been completed. In the end of this paper, having been applied on stability of Accumulative slopes in Three Gorge area and compared with result by limit equilibrium, classification system has good effect.

电离层胶片频高图数字化转换与分析研究

The ionogram acquired with the ionospheric vertical sounding method is the oldest data in the history of ionospheric research. Using of modern microelectronics and computer technology to digitalize, analyse and preserve the huge amount of historical film ionogram has become more and more important and urgent. This paper introduced the progress of the film ionogram digitalization by using digital image processing technologies to correct and repair film ionogram and convert them in an exchangeable format. An analysis and conversion software, basing on this method, has been developed for the film ionogram analysis, and then it introduces the application of this software by combining the SAO Explorer program for Wuhan film ionogram and pseudo-color ionogram in Yamagawa in Japan. It shows that our method is reliable，and the developed software is used friendly and provides a positive solution in digitalization and analysis of huge amount of historical film ionogram. Firstly, we briefly introduce the film ionogram and the process of its digitalization. By observing a amount of film ionogram, we obtain some common characteristics of the digitalized film ionogram following as: (1) the image rotation are caused by scanning; (2) the vertical axis of a large number of film ionogram exist more or less tilt and bending ; (3) coordinates of the film ionogram appear the non-uniformity phenomena result from the instability of driving motor rotation and the error of altitudinal cursor orientation. Moreover, based on the characteristics of the film ionogram and the SAO Explorer software which is widely used for the digital ionogram analysis in the world, a new method has been developed for film ionogram procession. The method contains the image geometric correction and film ionogram format conversion. The image geometric correction includes such as image rotation correction, vertical correction and coordinates scale correction. After geometric correction, the BMP file format images will be converted to the SBF file format images. Then, we also discuss the data format converting methods, which include two methods of the image data mapping basing on the normalization and logarithm, and the method of the preprocessing of the noise filtering and the threshold setting. Combining with SAO Explorer software, we successfully obtain ionospheric parameters and electron profile from the converted SBF file format digital ionograms. Based on the above method, we developed the software for the film ionogram to realize its correction analysis and conversion of the image format, and then give a introduction for its function and operation. Subsequently, the software are applied into the Wuhan film ionogram which separately observed in the high solar activity year and the low in 1980s last century. The results reveal the converted SBF digital ionogram almost preserve the all echo information of the film ionogram. Furthermore, we expressly discuss the application to the Wuhan film ionogram in 1958 in order to validate the applicability and credibility of the software. And it is showed that the important information of the film ionogram are maintained into the SBF digital ionogram. It is represented that there is credibility for conversion of the software when it applied in the older film ionogram. In sum, this software could apply to the digitalization and analysis of huge amount of historical film ionogram. Last, we extended the function of the software by bring some new conversion method and used it to apply to the pseudo-color ionogram of yamagawa in Japan. The results show that the converted ionogram information basically maintain the importantly ionogram information and the error of scaling of converted SBF file format image is almost acceptable, though there is no preprocessing for the original ionogram. Hence, we could extend the applicable range of the software and apply it to all kinds of simulative ionogram imaging by improving the method and software.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.