A Variant Constrained Genetic Algorithm for Solving Conditional Nonlinear Optimal Perturbations简

A variant constrained genetic algorithm （VCGA） for effective tracking of conditional nonlinear optimal perturbations （CNOPs） is presented. Compared with traditional constraint handling methods, the treatment of the constraint condition in VCGA is relatively easy to implement. Moreover, it does not require adjustments to indefinite pararneters. Using a hybrid crossover operator and the newly developed multi-ply mutation operator, VCGA improves the performance of GAs. To demonstrate the capability of VCGA to catch CNOPS in non-smooth cases, a partial differential equation, which has ＂on off＂ switches in its forcing term, is employed as the nonlinear model. To search global CNOPs of the nonlinear model, numerical experiments using VCGA, the traditional gradient descent algorithm based on the adjoint method （ADJ）, and a GA using tournament selection operation and the niching technique （GA-DEB） were performed. The results with various initial reference states showed that, in smooth cases, all three optimization methods are able to catch global CNOPs. Nevertheless, in non-smooth situations, a large proportion of CNOPs captured by the ADJ are local. Compared with ADJ, the performance of GA-DEB shows considerable improvement, but it is far below VCGA. Further, the impacts of population sizes on both VCGA and GA-DEB were investigated. The results were used to estimate the computation time of ~CGA and GA-DEB in obtaining CNOPs. The computational costs for VCGA, GA-DEB and ADJ to catch CNOPs of the nonlinear model are also compared.     

条件非线性最优扰动法在大气与海洋目标观测研究中的应用

对近年来用条件非线性最优扰动法研究大气与海洋目标观测问题的部分工作进行了总结,主要涉及厄尔尼诺-南方涛动（ENSO）事件、黑潮路径变异事件以及阻塞事件。通过研究这些事件发生的最优前期征兆（OPR）和最快增长初始误差（OGE）,发现这些事件的最优前期征兆和最快增长初始误差分别具有空间的高度相似性及其伴随的局地性特征。理想回报试验表明,如果在ENSO事件和黑潮路径变异事件的最快增长初始误差和最优前期征兆所确定的扰动大值区减小初始场误差,上述事件的预报技巧会大幅度提高;最优前期征兆和最快增长初始误差的空间相似性使得在同一敏感区域增加额外观测,不仅有助于捕捉上述异常事件的前期信号,还可以有效减小初始误差,从而提高对该事件的预报技巧。阻塞事件爆发的最优前期征兆和最快增长初始误差的空间相似性和局地性特征在其目标观测研究中的应用,应该是深入研究的课题。     

A systematic sensitivity analysis of constrained cellular automata model for urban growth simulation based on different transition rules

Cellular automata (CA) have emerged as a primary tool for urban growth modeling due to its simplicity, transparency, and ease of implementation. Sensitivity analysis is an important component in CA modeling for a better understanding of errors or uncertainties and their propagation. Most studies on sensitivity analyses in urban CA modeling focus on specific component such as neighborhood configuration or stochastic perturbation. However, sensitivity analysis of transition rules, which is one of the core components in CA models, has not been systematically done. This article proposes a systematic sensitivity analysis of major operational components in urban CA modeling using a stepwise comparison approach. After obtaining transition rules, three stages (i.e. static calibration of transition rules, dynamic evolution with varied time steps, and incorporation with stochastic perturbation) are designed to facilitate a comprehensive analysis. This scheme implemented with a case study in Guangzhou City (China) reveals that gaps in performance from static calibration with different transition rules can be reduced when dynamic evolution is considered. Moreover, the degree of stochastic perturbation is closely related to obtain urban morphology. However, a more realistic (i.e. fragmented) urban landscape is achieved at the cost of decreasing pixel-based accuracy in this study. Thus, a trade-off between pixel-based and pattern-based comparisons should be balanced in practical urban modeling. Finally, experimental results illustrate that models for transition rules extraction with good quality can do an assistance for urban modeling through reducing errors and uncertainty range. Additionally, ensemble methods can feasibly improve the performance of CA models when coupled with nonparametric models (i.e. classification and regression tree).     

Approximate soil-structure interaction analysis by a perturbation approach: The case of stiff soils

An approximate solution of the classical eigenvalue problem governing the vibrations of a structure on an elastic soil is derived through the application of a perturbation analysis. For stiff soils, the full solution is obtained as the sum of the solution for a rigid-soil and small perturbing terms related to the inverse of the soil shear modulus. The procedure leads to approximate analytical expressions for the system frequencies, modal damping ratios and participation factors for all system modes that generalize those presented by other authors for the fundamental mode. The resulting approximate expressions for the system modal properties are validated by comparison with the corresponding quantities obtained by numerical solution of the eigenvalue problem for a nine-story building. The accuracy of the proposed approach and of the classical normal mode approach is assessed through comparison with the exact frequency response of the test structure.     

Effect of stratification on the frequency of bounded Rossby modes over a non-flat bottom

This work attempts to express and analyze the challenges, induced by stratification, affecting the Rossby-topographic eigenmodes of a closed domain with a general uneven bottom of arbitrary shape filled with a uniform fluid in the unperturbed configuration. The modified eigenmodes have been computed analytically: stratification is introduced in the mathematical form of a perturbation of a homogeneous fluid over a non-flat bottom. The eigenmodes lose their barotropic character and differences appear in the dynamical fields (velocity and pressure) from upper to lower layer, as expected. Expressions for the baroclinic and ageostrophic velocity components due to the perturbation are given. The analysis is carried out in the frame of linear shallow water approximation. All terms have been retained apart from nonlinear advection in the governing equations. We find that the frequencies of the eigenmodes change; an analytical expression of frequency correction as a function of layer density difference and interface depth is found. Initial results for some elementary geometrical settings with a waveguide bottom are determined and expressed in a concise, easily readable closed form. The results obtained in the shallow water approximation are expanded in series with respect to the Rossby number. Next, they are compared with the frequency correction obtained in an alternative framework in which the quasi-geostrophic approximation is used, and a purely baroclinic perturbation is imposed from the outset as the result of the introduction of stratification in the otherwise homogeneous fluid. In this scenario, reduced gravity and the ratio of upper to lower layer depth are, in turn, used as the expansion parameters in lieu of the Rossby number.     

基于 Collocation 积分的 GPS 轨道仿真

在对G PS仿真的理论实施基础上,为方便求取卫星任意历元的位置、速度、加速度和加加速度,介绍了轨道仿真中涉及的摄动力模型,详细讨论了轨道仿真中初值的确定及轨道外推的方法,选用Collocation积分方法实现了轨道仿真,利用附加确定性参数的方法提高了模型精度,满足了1米精度的定轨要求。     

A High Order Perturbation Analysis of the Sitnikov Problem

The Sitnikov problem is one of the most simple cases of the elliptic restricted three body system. A massless body oscillates along a line (z) perpendicular to a plane (x,y) in which two equally massive bodies, called primary masses, perform Keplerian orbits around their common barycentre with a given eccentricity e. The crossing point of the line of motion of the third mass with the plane is equal to the centre of gravity of the entire system. In spite of its simple geometrical structure, the system is nonlinear and explicitly time dependent. It is globally non integrable and therefore represents an interesting application for advanced perturbative methods. In the present work a high order perturbation approach to the problem was performed, by using symbolic algorithms written in Mathematica. Floquet theory was used to derive solutions of the linearized equation up to 17th order in e. In this way precise analytical expressions for the stability of the system were obtained. Then, applying the Courant and Snyder transformation to the nonlinear equation, algebraic solutions of seventh order in z and e were derived using the method of Poincaré–Lindstedt. The enormous amount of necessary computations were performed by extensive use of symbolic programming. We developed automated and highly modularized algorithms in order to master the problem of ordering an increasing number of algebraic terms originating from high order perturbation theory.     

总体最小二乘拟合问题求解方法的比较研究

目前对总体最小二乘求解方法的研究,出现了奇异值分解的总体最小二乘法、顾及自变量和因变量误差的总体最小二乘法及正交总体最小二乘法.在模型推导的基础上,本文对3种总体最小二乘法在直线和平面拟合中求解的参数及其精度进行了分析,通过与最小二乘法的比较表明,总体最小二乘法得到的拟合结果更加稳健,且以正交总体最小二乘法的拟合结果为最优.     

日本南部黑潮路径发生弯曲的最优前期征兆及其发展机制

基于正压出入流模式, 利用条件非线性最优扰动(CNOP)方法研究初始异常的位置与模态对日本南部黑潮路径变异的影响。以模式模拟出的黑潮平直路径的平衡态作为参考态, 计算CNOP, 考察该扰动随时间的发展, 并与随机扰动的发展进行对比。结果表明, CNOP 能够导致黑潮弯曲路径发生, 随机扰动则不能。因此, CNOP 可以作为导致日本南部黑潮路径发生弯曲的一种最优前期征兆。通过分析CNOP 和随机扰动的发展过程, 可以得出: (1) CNOP 使黑潮发展成弯曲路径的过程是一个气旋涡向下游传播并增长的过程。(2) 气旋涡的向东传播都是非线性项的作用, 也就是涡度平流造成的。(3) CNOP和随机扰动发展过程中所产生的气旋涡均会传播到下游区域, 但是CNOP 产生的气旋涡能够增强, 最终导致弯曲路径发生, 而随机扰动产生的气旋涡则会减弱, 并不能导致弯曲路径发生。分析发现, 在CNOP 实验中, 非线性作用使气旋涡增大; 但在随机扰动实验中, 非线性作用使气旋涡减弱, 所以非线性作用对日本南部黑潮路径发生弯曲有重要影响。(4) 底摩擦效应对日本南部黑潮路径变异影响较小。本文揭示的黑潮路径发生弯曲的最优前期征兆及其非线性发展机制, 对提高黑潮路径变异的预报技巧具有重要意义。     

考虑地球扁率摄动影响的初轨计算方法

在二体问题意义下的短弧定轨,Laplace型方法是最主要最典型的一种初轨计算方法。若测角资料达到10^-4-10^-5精度(相当于2″—20″之间),那么要使定轨精度达到与其相应的程度,地球非球形引力位中的扁率项摄动应该考虑,在此前提下,同样可以采用相应的Laplace型定轨方法。即给出这种严格包含扁率摄动的初轨计算方法的原理和具体计算过程以及计算实例,除采用多资料定轨方法外,这种方法也是提高初轨计算精度的一种途径,它同样可用于多资料的情况,这种方法对于大扁率主天体(即中心天体)的卫星定轨将更有实用价值。