An Aggregate Constraint Method for Inequality-constrained Least Squares Problems

The inequality-constrained least squares (ICLS) problem can be solved by the simplex algorithm of quadratic programming. The ICLS problem may also be reformulated as a Bayesian problem and solved by using the Bayesian principle. This paper proposes using the aggregate constraint method of non-linear programming to solve the ICLS problem by converting many inequality constraints into one equality constraint, which is a basic augmented Lagrangean algorithm for deriving the solution to equality-constrained non-linear programming problems. Since the new approach finds the active constraints, we can derive the approximate algorithm-dependent statistical properties of the solution. As a result, some conclusions about the superiority of the estimator can be approximately made. Two simulated examples are given to show how to compute the approximate statistical properties and to show that the reasonable inequality constraints can improve the results of geodetic network with an ill-conditioned normal matrix.     

中纬全站仪GSI数据生成CASS9.1绘图坐标文件编程设计

利用全站仪测量碎部点是目前1∶500数字地籍碎部测量主要的作业方法。为了节省外业设站、定向时间,提高作业效率,可通过中纬全站仪采集测站点到碎部点的距离与角度等数据(下文称GSI数据),再用Visual Basic 6.0编写程序,根据测站点坐标及GSI数据生成CASS9.1绘图坐标文件。相关测量实例证明了上述编程可达到预期效果。     

基于动恣规划的整体影像匹配

影像匹配是数字摄影测量工作的基础。考虑到多条核线上的匹配点，本文提出了一种基于动态规划的整体影像匹配方法，构造了一种新的立体匹配模型，并将该模型转化成为简单的一维动态规划模型，采用对角线的方法划分阶段，用视差较和相关系数来构造阶段的代价函数，从而得到了基于动态规划的整体最优匹配结果。     

锥束CT图像重建算法的快速实现

本文基于锥束CT滤波反投影重建的FDK算法,通过两种算法改进并结合基于共享内存的OpenMP并行技术和代码优化,实现了锥束CT图像的快速重建。基于锥束CT实际投影数据的重建结果表明,图像重建速度得到了较大的提高,断层图像重建质量与FDK原型算法相当。     

An application of Karmarkar's interior-point linear programming algorithm for multi-reservoir operations optimization

Optimization of multi-reservoir systems operations is typically a very large scale optimization problem. The following are the three types of optimization problems solved using linear programming (LP): (i) deterministic optimization for multiple periods involving fine stage intervals, for example, from an hour to a week (ii) implicit stochastic optimization using multiple years of inflow data, and (iii) explicit stochastic optimization using probability distributions of inflow data. Until recently, the revised simplex method has been the most efficient solution method available for solving large scale LP problems. In this paper, we show that an implementation of the Karmarkar's interior-point LP algorithm with a newly developed stopping criterion solves optimization problems of large multi-reservoir operations more efficiently than the simplex method. For example, using a Micro VAX II minicomputer, a 40 year, monthly stage, two-reservoir system optimization problem is solved 7.8 times faster than the advanced simplex code in MINOS 5.0. The advantage of this method is expected to be greater as the size of the problem grows from two reservoirs to multiples of reservoirs. This paper presents the details of the implementation and testing and in addition, some other features of the Karmarkar's algorithm which makes it a valuable optimization tool are illuminated.     

Matlab与高级编程语言在科学计算中的联合运用

在数字地震观测、地震预报及地震工程研究中,Matlab作为一套功能强大的科学计算及数学分析软件,得到广泛的应用.本文介绍如何将Matlab高效的数值计算和分析能力与高级编程语言结合的方法.     

DAS24-3B数据采集系统采集控制模块的工作原理

DAS24-3B是中国科学院地质与地球物理研究所研制的宽频带地震数据采集系统。采集控制模块是该系统的重要组成部分。本文阐述了控制模块的工作原理，介绍了DAS24-3B系统的基本组成和采集器工作过程，并给出该系统如何与笔记本端应用程序界面通讯和联络的操作建议，使该系统能够在工作中发挥更大的作用。     

Optimization and capacity expansion of a water distribution system

This paper develops an iterative procedure for capacity expansion studies for water distribution systems. We propose a methodology to analyze an existing water distribution system and identify the potential bottlenecks in the system. Based on the results, capacity expansion alternatives are proposed and evaluated for improving the efficiency of water supply. The methodology includes a network flow based optimization model, four evaluation indices, and a series of evaluation steps. We first use a directed graph to configure the water distribution system into a network. The network flow based model optimizes the water distribution in the system so that different expansion alternatives can be evaluated on a comparable basis. This model lends itself to linear programming (LP) and can be easily solved by a standard LP code. The results from the evaluation tool help to identify the bottlenecks in the water distribution system and provide capacity expansion alternatives. A useful complementary tool for decision making is composed of a series of evaluation steps with the bottleneck findings, capacity expansion alternatives, and the evaluation of results. We apply the proposed methodology to the Tou-Qian River Basin, located in the northern region of Taiwan, to demonstrate its applicability in optimization and capacity expansion studies.     

Optimizing layout of pumping well in irrigation district for groundwater sustainable use in northwest China

Optimizing layout of pumping well plays a vital role in curbing the groundwater level decline. A novel optimization model is presented in this study. First, the optimal well number is obtained by taking into account factors of local economy and environment based on nonlinear programming model. Then, the well spatial layout assessment model is attained based on information entropy weight and technique for order preference by similarity to ideal solution (TOPSIS). After that, the relative closeness to positive ideal solution of alternative (c_i) on the rationality of well spatial layout in cultivated land is calculated, and a set of alternatives are ranked according to the descending order of c_i. Finally, the well optimization layout is obtained by combining the optimal well number with well spatial layout assessment result based on the GIS data of pumping wells. As a case study, this method was applied in Yongchang Irrigation District of Shiyang River Basin, the arid region of northwest China. Results show that under the conditions of sustainable use of water resources, the irrigation district needed 724 wells for irrigation, with a decrease of 31.0% when compared with the existing number of wells. The wells with low flow rate and operating efficiency distributed in high density where groundwater is over‐exploitation were recommended to be closed. This well optimization layout method is expected to play a significant role in helping make plans for exploiting groundwater at more sustainable level, curbing the groundwater level decline trend, and improving the local ecological environment. Copyright © 2015 John Wiley & Sons, Ltd.     

大地电磁三维正演并行算法研究

大地电磁法三维正演算法计算量大,用传统的串行程序计算相当耗时。而三维正演是逐个频率按顺序计算的,并行性好,适合并行运算。结合MPI自身的优越性,在深入分析大地电磁三维正演串行程序实现流程的基础上,确定了并行计算的思路,实现了三维正演的并行计算。通过三个理论模型对实现的三维正演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维正演并行程序运行结果正确,效率提高明显。此思路可为解决其它地球物理超大计算量问题所借鉴。