论广义逆矩阵阶数增大的递推算法

本文针对加拿大地球物理学家ＥｄｏＮｇｌａｎｄ教授提出的问题，对Ｍ—Ｐ广义逆矩阵的递推算法进行了探分论述探讨的初步结果和对Ｇｒｅｖｉｌｌｅ递推算法的扩充。文中对算法的扩充给出严格的证明。     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

基于“双差法”地震精确定位研究2002年12月14日玉门5.9级地震的发震断层和发震机制

利用双差地震定位方法对2002年12月14日玉门5．9级地震主震和余震序列进行了精确定位，在此基础上讨论了本次地震的发震断层和发震机制。     

区域-时间-长度算法及其在华北中强地震中短期前兆特征研究中的应用

引入并改进了区域-时间-长度算法（region-time-length algorithm），较之原来的定义更加强调了地震破裂尺度函数对VRTL的贡献，同时保证了震中距函数、时间函数及破裂尺度函数对VRTL影响的权系数相等．华北5级以上地震的回顾性预测研究结果表明，改进的RTL算法所得到的异常大多具有短期或中短期特征. 震前VRTL变化型态可分为两类:Ⅰ类具有相对完整的变化形态，大多呈现从0开始的上升——转折下降或下降——转折上升形态；Ⅱ类变化形态不完整，一般从0开始急剧上升或下降，无明显转折，主震发生在VRTL极值附近．VRTL上升表征地震活动水平相对于背景变化率的增高，意味着地震活动增强；VRTL下降则意味着地震活动相对于背景变化率的降低，在一定程度上表征了地震活动的平静．据此给出由RTL算法进行中、短期异常识别并粗略估计发震时间的方法. 该方法Ⅰ,Ⅱ类异常3个月内发震的R值评分分别为0.6及0.3左右， 6个月内的R评分分别为0.7及0.4左右．对特征时间跨度、特征距离、震级下限等对VRTL的影响, 以及其它一些实际应用中的重要问题进行了初步讨论．      

任意圆弧形凸起地形中隧洞对入射平面SH波的影响

借助辅助函数的思想， 采用波函数展开法给出了圆弧形凸起地形中隧洞对入射平面SH波影响问题的一个解析解．数值结果表明，凸起地形中隧洞的存在以及隧洞大小，对凸起地形表面运动，以及凸起地形中的隧洞的动应力集中均具有显著影响．      

3-D rheologic model of earthquake preparation (II): Strain field and its applications

Introduction Based on the elastic theory of the hard inclusion (Dobrovolskii, 1991), we developed an inclusion theory of rheologic medium, and applied the results of bulk-strain field of a rheologic inclusion model to explain the spatial-temporal evolution process of earthquake precursors (SONG, et al, 2000). In the former paper (SONG, et al, 2003), we derived the viscoelastic displacement field of the rheologic inclusion model on the basis of the analytic expression of displacement field o…     

星载被动微波遥感反演降水算法回顾

根据不同的降水反演方法对多种利用SSM／I、TMI反演降水的算法进行归类总结，按下垫面检测和降水反演算法两大步骤进行简要描述，并分经验法、半经验(半物理)法、物理模式法及物理廓线法4类方法对多种算法作了回顾。     

三维旋转式自动进样器的研制

介绍了自行研制的三维旋转式自动进样器的结构和特点，提出了三维模式、气动升降臂、旋转定位的自动进样器解决方案，并成功地用于自动进样器的研制。文章中给出了三维旋转式自动进样器的组成，气动式升降臂控制原理和系统的软件流程。     

基于EM算法和单幅雷达图像阴影的控制点坡度校正

对在我国现有的条件下进行控制点坡度校正的必要性进行了阐述，并分析了EM算法。根据EM算法以及基于区域增长的余弦散射模型建立了控制点坡度校正模型，并用河北省张北地区的雷达影像进行了实验，取得了较高的精度初值。