徕卡激光扫平系统在节水农业中的应用与发展

一、引言农业以土而立、以肥而兴、以水而旺。水是最短缺的农业重要资源之一,也是制约农业可持续发展的关键因素。我国水资源严重紧缺,农业用水矛盾日渐凸显。在我国现有灌溉面积9亿亩,其中有效灌溉面积8.3亿亩,林果牧草灌溉面积7000万亩,在这些灌溉面积中,95%以上采用地面灌溉,很     

求解半变异函数的常用方法与新方法研究

半变异函数的获取既有人为主观求取法,又有运用数学工具客观求取法,但是后者是今后的主要趋势。笔者在前人工作的基础上,归纳出多项式回归法、线性规划法和目标规划法三种求解半变异函数的成熟方法,以及非线性回归理论直接求解法、进化计算和积分面积三种新兴方法,并总结出各种方法的优缺点。最后概括出计算权重的两种方法。     

“准拉格朗日法”和“欧拉法”的数学一致性与三次插值函数算法时间积分方案

从预报方程组通式和欧拉算符出发,用泰勒级数展开,推导出时空间微商余项为二阶、四阶完全预报方程组。可以证明,同阶时空间微商余项的准拉格朗日法和向前差分欧拉法时间积分方案,具有相同的物理意义和数学一致性,相比之下,传统中央差欧拉法时间积分方案是"简单格式"与"增大不可预测计算误差"并存。进而讨论用"三次插值函数"实现二阶时空间微商余项准拉格朗日法、或同阶向前差分欧拉法(可"二选一"),它们应该分别替代"双线性插值"准拉格朗日法和传统时空间中央差欧拉法,因前二者时空间微商余项及计算精度高于后二者。所以,三次插值函数算法可将准拉格朗日法和欧拉法时间积分方案、以及CFL判据统一起来。由于三次插值函数具有对原函数"变量场"的数学定律"收敛性"和二阶可导"最优性",且一次"三次插值函数"运算,即具对网格变量场二阶可导拟合"等价性":不仅拟合变量场斜率、还拟合其曲率和挠率。并因周期"三次插值函数",可作全球变量场"三次插值函数"二阶可导拟合,实现全球"三次"数值模式,并且可按变量场曲率判断,作变量场局域或单点平滑,保持"三次"模式时间积分的稳定性。     

重力波特征线显式数值解法试验

利用CIP(Constrained Interpolation Profile)特征线半拉格朗日解法对平面一维和二维浅水波方程的重力波和普通经纬度网格下的球面浅水波模式进行了数值求解试验,给出了CIP特征线法的数值结果和对应大时间步长积分的对比。结果显示,一维CIP特征线法对重力波的位相、振幅和传播有很好的描述,对于不同时间步长都显示出计算稳定的特性。空间分离的二维特征线法是多维问题计算中简单有效的方法,在平面和球面二维问题的数值结果中都得到了CIP特征线法稳定的计算和良好的波动传播结果,在浅水波方程中可以克服半隐式算法在高分辨率时的大量运算,并有望成为处理重力波的替代算法。     

0679香山局地大暴雨的中小尺度天气分析

利用Doppler雷达、地面自动站、Profile垂直风廓线及GPS水汽分布等多种新型探测资料, 对2006年7月9日夜间发生在北京西郊香山附近的局地大暴雨天气的影响系统和γ-中尺度强降雨落区形成的动力机理进行了精细分析。对雷达等本地多种探测资料的精细研究表明:地形辐合回波带是造成这次过程的主要影响系统。地形辐合回波带上中气旋回波块的滚动更迭是大暴雨落区形成的直接原因; 近地面辐合对大暴雨落区强降雨的发生具有重要作用。大暴雨落区形成阶段近地面3种辐合同时存在:平原东南风与山区偏北风风向切变辐合、平原东南风在山脉阻挡作用下的抬升辐合、大暴雨落区中心的γ-中尺度气旋性辐合。研究还表明:山前近地面地形辐合扰动, 向上传播, 引发边界层扰动的动力过程是香山大暴雨落区形成的主要动力源, 而来自东南方向近地面层的暖湿平流为大暴雨提供了有效的水汽和能量。     

大气CO215 μm吸收带对温度依赖关系的研究

利用最新版本的大气分子吸收光谱资料HITRAN 2004,以CO215 μm吸收带为例,通过精确的逐线积分模式,详细讨论了温度对谱线半宽度、线强以及吸收系数的影响。结果发现：吸收系数对温度的依赖关系主要受线强对温度依赖关系的影响,受半宽度对温度依赖关系的影响较小;在吸收带的中心区域,温度对线强和吸收系数的影响较弱,而在带翼影响较强。     

星载GPS地球静止轨道卫星自主定轨的新算法

针对地球静止轨道卫星的特点,提出了一种自主定轨的新算法--积分滤波算法,即利用Kalman滤波进行动力学模型数值积分与GPS定轨的有效融合.讨论了该算法的基本过程及其中Kalman滤波的数学模型和重要参数.最后给出了仿真实验过程和结果,证明了该算法的可行性.     

Lippmann-Schwinger积分方程广义超松弛迭代解法及其收敛特性

为克服Born级数只适用于弱散射和短程传播的弱点,从光学散射理论中引入了针对强光学散射问题的Lippmann-Schwinger方程广义超松弛迭代解法,并对这种方法在反射地震条件下的算法表现、收敛特点、计算效率以及计算精度等问题进行了分析和讨论.理论分析和数值计算结果均表明:(1)与散射级数重整化相比,超松弛法具有更坚实的数学基础;(2)超松弛法可以有效地克服Born级数的弱点,得到与有限差分相当的数值模拟结果;(3)在反射地震条件下LS方程超松弛法表现优异,完全可以解决反射地震数值模拟中的强散射问题.     

散射波积分方程的Adomian分解解法

在背景模型基础上,求解模型扰动后的地震波散射场,这是目前地震反演中的一个关键步骤.本文将计算数学中求解非线性积分方程的Adomian分解方法,应用到求解标量波散射场的Lippmann-Schwinger积分方程和Ricatti积分方程中,分别得到了散射场的Born序列解和Rytov序列解.通过一维和二维数值算例说明:在满足一定的条件下,散射场的这两种序列解稳定收敛,与传统的Born和Rytov近似解相比,引入散射序列中的高阶项可以更精确地描述地震波散射场.     

瞬变电磁扩散场到虚拟波场的精细积分变换算法

现有的瞬变电磁解释与反演方法主要对电阻率参数进行解释,较难给出准确的地质构造信息.为了能够较多获得地质构造信息,可将瞬变电磁扩散场信号转换为虚拟波场信号,提高瞬变电磁法的分辨率.本文采用精细积分法,实现了瞬变电磁扩散场到虚拟波场的转换.将一个高度病态的线性方程组转换成求积分的过程,大大降低了解决病态问题的难度.积分步长以2的指数增加,在保证精度的前提下,具有极高的计算效率;设置合理的终止迭代条件,进一步提高本文方法的适用性.对典型地电模型的计算,表明本文方法具有较高的精度;减小两个波峰之间的距离,证明了本文方法具有较好的分辨能力;对含有噪声信号的扩散场进行波场反变换,证明了本文方法具有较好的抗噪性.最后计算三维模型和实测数据的虚拟波场,可以很好地反映出三维地质体的界面信息,证明本文方法稳定、可靠,可以对地下地质目标体的界面进行有效识别.本文方法相比预条件正则化共轭梯度法(PRCG),在计算效率方面提升了 4倍左右;本文方法的相对误差在5％以内,而PRCG的最大相对误差可达40％～50％,精度最高可提升10倍;在分辨率与稳定性方面,也表现出较好的效果.