快速高精度的频率空间域声波数值模拟方法研究

频率空间域地震波正演优势明显:无累计误差;多炮反演时,只需几个频率就能得到比拟时间域全带宽反演的结果;容易实现时间域难以精确实现的吸收衰减数值模拟.频率域波动方程被离散化而形成稀疏、非对称、不定、复系数、大规模的系数矩阵,需要求解大规模的线性方程组,最大计算瓶颈是需要海量的计算机内存、导致计算量庞大.采用嵌套剖分法存储...     

利用VC++封装MATLAB 6.5 C++数学库探讨测量平差软件的开发

根据测量数据处理的特点：常常要面对大量的矩阵运算和海量的数据处理,尤其是矩阵的相乘和求逆。通过对vc＋＋6．0和MATLAB6．5特点的分析,探讨了利用Vc＋＋调用MATLAB6．5C／C＋＋数学库快速开发测量平差软件的必要性和可行性,简要地总结和分析了两者联合编程的实现方法。较详细地介绍了作者用VC＋＋调用MATLAB6．5C／C＋＋数学库开发测量平差软件的过程,解决了高程网间接平差的计算。     

基于计算机视觉的地形3D重建

构造3D地形模型有很多种方法．本文提出一种用数码相机获得的未标定照片获取三维地形模型的简单方法．这种方法仅需要从两幅或多幅照片中选取对应点,根据这些对应点,求出基本矩阵,摄像机矩阵,进而进行三维地形表面的估计,实验结果及算法分析表明,本文提出的特征匹配算法精确度较高,实现地形重建方法切实可行,生成的模型可视效果较好。     

数码相机检验中出现的问题及解决方法

数码相机的检验是近景摄影测量中经常碰到的问题。最常用的主要还是基于共线方程的后方交会法。但由于普通数码相机存在较大的镜头畸变,内方位元素也未知,所以仅用单片进行后方交会求取有关参数时,法方程可能会出现病态现象,从而影响最后的检验结果。文中对此现象进行了讨论,并给出了一些实用的解决方法。     

黄河流域湿地景观时空演变格局分析

湿地景观的时空演变及其驱动因素研究是湿地生态恢复与保护的重要参考,本文利用1980—2015年7期黄河流域土地利用数据构建黄河流域湿地矢量数据集,基于景观格局指数对黄河流域湿地的时空演变特征进行定量和定性分析;采用转移矩阵对每两个时期湿地资源的转移类型和数量进行计算与分析。研究结果表明,①1980—2015年,黄河流域湿地率为3.5%,其中自然湿地率约为2.0%,即黄河流域的湿地类型以自然湿地为主。②斑块面积结果显示,黄河流域湿地总面积减少了312.74km²,表现为人工湿地增加,自然湿地减少;自然湿地中沼泽湿地呈增加状态,主要是在2005—2010年增加了979.22km²,滩地呈减少状态,减少了1218.19km²,主要发生在1990—1995年和2005—2010年。斑块密度结果显示,随着河渠和滩地受人为活动影响程度的加大,两者的斑块破碎化程度逐渐加大。最大斑块指数结果显示,水田是黄河流域的主导湿地类型,沼泽和湖泊是自然湿地中的主导湿地类型。③黄河流域湿地类型与非湿地之间转换是黄河流域景观转换的主要组成部分,其中沼泽与滩地是面积发生变化的主要湿地类型。     

直角坐标系的欧拉旋转变换及动力学方程

讨论了在直角坐标系中12种欧拉旋转变换中的6种坐标变换方式,根据矩阵的运算性质,推导了坐标变换的欧拉旋转矩阵和动力学方程,并且用M atlab的矩阵运算功能对所有公式的推导过程进行了验证。     

ω循环型边界条件

通过利用预条件共轭梯度法对对称正定Toeplitz 矩阵系统进行分析，重点介绍了一种新的嵌入式预条件矩阵构造方法，证明了以前的预条件矩阵构造方法大都是这种方法的特例.提出了ω循环型边界条件，并将其与普通循环型及螺旋型边界条件作了分析、比较后得到了一种新的边界条件即混合型边界条件.     

二维视电阻率断面的快速最小二乘反演

快速最小二乘反演是以平滑限定的最小二乘方法为基础,是对二维视电阻率断面进行反演的一种方法。反演过程不需要提供初始模型,在首次迭代时使用一均匀介质地下模型作为初始模型,该模型的视电阻率偏导数值可以用解析法得到。在后面的迭代中,使用了拟牛顿法去修改每一次迭代的偏导数矩阵,避免了偏导数矩阵的直接计算,从而减少了计算时间和存储空间。同时运用牛顿矩阵校正技术解最小二乘方程组也减少了大量的计算时间。总之,该方     

井间电阻率层析成象的平滑度约束反演

井间电阻率层析成象是一种探测地下浅部精细结构的物探方法，它主要用于解决工程地质问题。本文采用２．５－Ｄ有限元法，针对井间电阻率层析成象中的ＡＭ观测系统进行了反演成象计算。文中首先计算了灵敏度矩阵，然后给出了一种电阻率层析成象反演算法—平滑度约束反演，理论模型的计算和实际资料的处理，都证明了该算法的有效性     

罗德里格矩阵在共线方程严密解法中的应用

在共线方程严密解法中,人们习用的是以方位元素X_s,Y_s,Z_s,φ,ω,κ为未知数的严密公式。本文则以罗德里格旋转矩阵的三个独立元素b,a,c取代角方位元素φ,ω,κ,并推导了相应的误差方程式。理论表明,新方法在组成旋转矩阵和计算误差方程系数中只需作加、减、乘、除等运算而不需解求三角函数,因而节省了计算工作量。