一种基于Partial EIV模型的多项式拟合解法

针对多项式拟合模型系数矩阵中部分元素是某一自变量的函数的特点,根据Partial EIV模型的解算思想,将系数矩阵中自变量的函数作为随机元素提取,顾及泰勒展开的二阶项,由协方差传播律计算自变量函数的协因数阵进行平差解算。实验结果表明,系数矩阵的元素不再是单独的自变量时,使用该算法可以得到与已有非线性总体最小二乘方法相近的参数结果,从构造随机向量权阵的角度提供了一种新的解算方法。     

介质密度反演偏导矩阵的精确计算方法

实现反演偏导矩阵的计算是基于导数最优化反演方法的关键,然而目前的地震反演几乎都是基于Zoeppritz方程近似实现的,使计算精度和适应范围受到限制.本文利用Zoeppritz方程建立了反射系数对地层介质密度比偏导方程,导出了Zoeppritz方程矩阵元对介质密度比的导数.通过求解偏导方程获得了反射系数对介质密度比偏导数的精确计算(考虑了速度中含介质密度的问题).利用数值算例分析了反射系数对介质密度比偏导数的变化特点.本文采用直接解法求解偏导矩阵方程组,获得了快的计算速度和高的计算精度,为实现地层介质密度反演(包括大角度反演)提供了偏导矩阵的计算方法.     

基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案

作为高精度波形反演或逆时偏移的重要组成部分,地震波数值模拟对计算速度和效率提出了更高要求.GPU通用计算技术的产生及其内在数据并行性,为高效地震波数值模拟应用和研究得以有效开展奠定了基础.本文借助交错网格的Fourier伪谱微分矩阵算子和GPU上高效矩阵乘法,实现了复杂介质中地震波模拟的高效算法.数值试验表明,优化后的GPU计算相比CPU单核计算在大规模二维地震波场计算中获得至少100x以上的加速比.这对我们快速分析目标反射层在地震剖面中同相轴位置,制定优化采集方案具有重要意义.     

地形变应变张量矩阵的不变量分析

在给出正交曲线坐标系的有关位移向量及其全微分、位移梯度矩阵、应变张量矩阵的普适表达式的基础上,又给出了任意两种正交曲线坐标系下的应变张量矩阵的普适转换表达式,并指出:由于该变换矩阵为正交矩阵,故应变张量矩阵为相似矩阵。并对应变张量矩阵的几何物理性质进行了分析,指出任何一种正交曲线坐标系的应变张量矩阵都具有唯一不变的主应变特征多项式,由该矩阵的主应变特征值方程皆可求得地壳质点处的主应变及其主方向,由主方向单位向量又可把该矩阵化为以主应变为对角元素的对角矩阵,该矩阵及其对角矩阵的迹皆为该质点处的体应变,该矩阵的行列式等于该质点处3个主应变的乘积,这些几何物理量皆为该质点处的地应变不变量。     

具有随机状态转移矩阵的Kalman滤波

讨论状态转移矩阵为随机性矩阵时的Kalman滤波方法,并给出滤波的具体算法及分析滤波的收敛特性和参数的取值范围。通过实际算例证实了该滤波方法的有效性。     

基于MRPs无陀螺卫星姿态确定的抗差估计算法

针对精确航天器姿态问题,采用修正罗德里格参数(MRPs)表示姿态,用动力学方程进行角速率的传播,分别基于等价协方差和观测残差统计特性的协方差调整方法,实现UKF算法的抗差效果,并通过仿真验证该方法的有效性。仿真结果表明：在含有粗差的情况下,两种抗差UKF算法均能得到稳定的状态估计值,而第二种方法计算效率更高。     

���α���ת����̽��

????????α???????????????????????????????????????????????λ???????????????????????????λ??????????????t???????1/2????????????????t?????????????????????????????????????????????????????????????????????λ?????????λ???????????????????????????????????????????????????????λ??????????????????????????????????????????λ??????????????λ??????????????????????????????λ?????????????????????????λ?????????????????????????????????????????????????????????Ч?????????????     

地基沉降预测模型的正则化算法

通过分析地基沉降预测模型,指出最小二乘的病态性会导致模型参数求解失败。应用正则化理论,基于矩阵求逆理论,提出了一种沉降预测模型参数的正则化无偏估计算法,说明了新算法的无偏性和方差最小性。在一定条件下,证明了新算法中正则参数的存在性,并给出了正则参数的计算公式。结合文献和工程实例进行的分析表明,新算法降低了矩阵条件数,减轻矩阵病态程度,能有效求得地基沉降预测模型参数。     

基于Markov模型的安宁市土地利用预测

根据安宁市2003～2009年土地利用更新调查数据,利用Markov模型以3年为间隔预测安宁市2012～2021年内的土地利用结构变化。结果表明:2012～2021年安宁市的耕地、园地、林地等农用地大量转变为建设用地;为了保障基本社会生活需求,耕地和园地的农作物生产用地向建设用地转变较少,林地向建设用地转变较多。     

多肢剪力墙自振频率传递矩阵法简化计算

基于连续体系模型与传递矩阵法的结合,提出了一种简化方法来计算多肢剪力墙的动力特性。整个结构体系等效模拟为一个夹层结构梁;建立了相应的基本微分方程;通过求解微分方程得到了楼层间变形函数;利用边界条件及层间变形函数导出了层间传递矩阵;据此计算出了体系的振型与自振频率。算例分析表明本文方法计算结果与文献给出的解是比较吻合的。