Spectral estimates for saddle point matrices arising in weak constraint four‐dimensional variational data assimilation

We consider the large sparse symmetric linear systems of equations that arise in the solution of weak constraint four‐dimensional variational data assimilation, a method of high interest for numerical weather prediction. These systems can be written as saddle point systems with a 3 × 3 block structure but block eliminations can be performed to reduce them to saddle point systems with a 2 × 2 block structure, or further to symmetric positive definite systems. In this article, we analyse how sensitive the spectra of these matrices are to the number of observations of the underlying dynamical system. We also obtain bounds on the eigenvalues of the matrices. Numerical experiments are used to confirm the theoretical analysis and bounds.     

基因芯片荧光图象采集与分析

描述了基因芯片荧光图象的光学共聚焦成象,大范围高速扫描与控制,数据采集的原理和方法.并讨论了基因芯片荧光图象的分析算法,即模板定位与阈值分割相结合的半自动算法,并将其计算结果与完全采用人工分割算法的结果进行了比较.     

Efficiency and effectiveness of Spanish football teams: a three-stage-DEA approach

In this paper, we apply a three-stage-DEA model to the Spanish Professional Football League, which means separating the teams’ economic behaviour into three components: operating efficiency—of the offence and defence—athletic or operating effectiveness, and social effectiveness. The results showed that the technical inefficiency of the defence is greater than that of the offence, both being caused by aspects linked to the poor management of players’ abilities and by the football team’s size. Teams showed a favourable evolution of their offensive and defensive efficiency during the 2004/2005 season and to a lesser extent in the season before. The point system assigned by the professional football league regulations evaluates the teams’ athletic effectiveness, but we detected that the teams with the most experience perform athletically in a more effective manner. Their social effectiveness is strongly related to the level of play in itself and to factors linked to their PFL ranking: participation in international competitions for important football teams; or the struggle of minor football teams to stay in the first division.     

Topics in data assimilation: Stochastic processes

Stochastic models with varying degrees of complexity are increasingly widespread in the oceanic and atmospheric sciences. One application is data assimilation, i.e., the combination of model output with observations to form the best picture of the system under study. For any given quantity to be estimated, the relative weights of the model and the data will be adjusted according to estimated model and data error statistics, so implementation of any data assimilation scheme will require some assumption about errors, which are considered to be random. For dynamical models, some assumption about the evolution of errors will be needed. Stochastic models are also applied in studies of predictability.

The formal theory of stochastic processes was well developed in the last half of the twentieth century. One consequence of this theory is that methods of simulation of deterministic processes cannot be applied to random processes without some modification. In some cases the rules of ordinary calculus must be modified.

The formal theory was developed in terms of mathematical formalism that may be unfamiliar to many oceanic and atmospheric scientists. The purpose of this article is to provide an informal introduction to the relevant theory, and to point out those situations in which that theory must be applied in order to model random processes correctly.     

长江堤防工程地质信息处理模型及技术方法

长江中下游的主要防洪工程是长江干堤 ,长江堤防线路长、工程浩大、堤基工程地质水文地质条件复杂 ,涉及大量的已有地质数据和资料。今后在整治加固堤防的工程地质勘察工作中 ,还需随时进行数据的收集、存储、检索、分析、统计和维护 ,随时调用这些数据进行综合处理、编制报告和图件 ;在防汛期间可以针对堤防险工险段的险情 ,随时查询和调用有关的资料、数据或者图件 ,提供采取抢险措施决策的地质依据。这既提高了地质数据的利用价值和利用率 ,又保证了信息的准确性、完整性和共享性。因此 ,建立合理的长江堤防工程地质信息处理模型及其成熟的技术方法具有重要的意义。     

Zero weights and non-zero slacks: Different solutions to the same problem

This paper re-assesses three independently developed approaches that are aimed at solving the problem of zero-weights or non-zero slacks in Data Envelopment Analysis (DEA). The methods are weights restricted, non-radial and extended facet DEA models. Weights restricted DEA models are dual to envelopment DEA models with restrictions on the dual variables (DEA weights) aimed at avoiding zero values for those weights; non-radial DEA models are envelopment models which avoid non-zero slacks in the input-output constraints. Finally, extended facet DEA models recognize that only projections on facets of full dimension correspond to well defined rates of substitution/transformation between all inputs/outputs which in turn correspond to non-zero weights in the multiplier version of the DEA model. We demonstrate how these methods are equivalent, not only in their aim but also in the solutions they yield. In addition, we show that the aforementioned methods modify the production frontier by extending existing facets or creating unobserved facets. Further we propose a new approach that uses weight restrictions to extend existing facets. This approach has some advantages in computational terms, because extended facet models normally make use of mixed integer programming models, which are computationally demanding.     

Error estimation,adaptivity and data transfer in enriched plasticity continua to analysis of shear band localization

In this paper, an adaptive FE analysis is presented based on error estimation, adaptive mesh refinement and data transfer for enriched plasticity continua in the modelling of strain localization. As the classical continuum models suffer from pathological mesh-dependence in the strain softening models, the governing equations are regularized by adding rotational degrees-of-freedom to the conventional degrees-of-freedom. Adaptive strategy using element elongation is applied to compute the distribution of required element size using the estimated error distribution. Once a new mesh is generated, state variables and history-dependent variables are mapped from the old finite element mesh to the new one. In order to transfer the history-dependent variables from the old to new mesh, the values of internal variables available at Gauss point are first projected at nodes of old mesh, then the values of the old nodes are transferred to the nodes of new mesh and finally, the values at Gauss points of new elements are determined with respect to nodal values of the new mesh. Finally, the efficiency of the proposed model and computational algorithms is demonstrated by several numerical examples.     

Reducing structural changes in van Emde Boas' data structure to the lower bound for the dynamic predecessor problem

We consider the problem of maintaining a dynamic ordered set of n integers in a universe U under the operations of insertion, deletion and predecessor queries. The computation model used is a unit-cost RAM, with a word length of w bits, and the universe size is |U|=2^w. We present a data structure that uses O(|U|/log|U|+n) space, performs all the operations in O(loglog|U|) time and needs O(loglog|U|/logloglog|U|) structural changes per update operation. The data structure is a simplified version of the van Emde Boas' tree introducing, in its construction and functioning, new concepts, which help to keep the important information for searching along the path of the tree, in a more compact and organized way.     

Determining Modes for Continuous Data Assimilation in 2D Turbulence

We study the number of determining modes necessary for continuous data assimilation in the two-dimensional incompressible Navier–Stokes equations. Our focus is on how the spatial structure of the body forcing affects the rate of continuous data assimilation and the number of determining modes. We treat this problem analytically by proving a convergence result depending on the H ^–1 norm of f and computationally by considering a family of forcing functions with identical Grashof numbers that are supported on different annuli in Fourier space. The rate of continuous data assimilation and the number of determining modes is shown to depend strongly on the length scales present in the forcing.     

试题质量评价研究

教学过程的规范化管理是保证教学质量的重要手段。考试是对教和学的效果的检验，是教学过程管理规范化必须要解决的问题，一直受到教学管理部门的重视，采取了多项措施解决有关问题，其中一个重要的手段是采用试题库。在采用试题库时，必须要有高质量的试题作保证。本探讨基于KDD技术的试题质量评价方法和指标，并通过实例来说明其有效性。