关于控制迭代过程的结束条件

迭代法属于逐次逼近法,应用很广。它是求解函数方程、方程组和微分方程等的一种基本而重要的数值方法。在具体求解时,需要给出控制迭代过程结束的条件。现有的某些计算方法(或数值分析)的教材和书,在论述求解函数方程的计算方案、流程图(或框图)及算法语言程序时,常常将此条件取为相邻两次迭代值xk、xk-1的差小于给定的准确度本文指出这个条件是不完全的,讨论了该条件的正确表达式,由误差估计式给出了控制迭代过程结束条件的五种形式。     

计算机程序在估价报告中的应用程序

土地估价内业中涉及大量计算及相似的文字描述，如何减轻劳动量，加快评估速度及减少笔误是摆在每个评估师面前的一个课题。针对这一问题，笔者利用常用的文字处理程序(Word)、电子表格程序(Excel)开发了小程序，涉及到评估常见的方法：基准地价系数修正法、市场比较法、土地成本逼近法及评估费的计算。其中土地成本逼近法部分撰写了论文《浅谈在土地估价中采用Excel计算地价、并形成Word文本——以成本逼近法中土地取得费为例》。     

复外差法处理绝对重力仪数字干涉信号

使用复外差方法分析绝对重力仪中自由下落过程产生的激光干涉信号。对采样的扫频信号进行一系列复外差处理,迭代地缩窄其频带至接近零频的位置,估算落体初始速度及加速度,并用非线性拟合得到重力加速度值g。仿真结果表明,该方法的拟合精度可达10^－3 μGal,满足精确测量要求。     

使用黏性逼近方法推广一个平衡及不动点问题

为了研究寻找一般平衡问题的公共解集中共有元素的迭代方法,及无穷非扩张映像族的不动点问题.通过在Hilbert空间中构造非扩张序列Ti得出了一种关于Wn算子的新黏性逼近迭代算法,在一定条件下,得出了其强收敛定理.     

厄瓜多尔俯冲区震间闭锁与粘弹性变形研究北大核心CSCD

以厄瓜多尔俯冲区为例,建立区域三维纯弹性和粘弹性模型,并基于有限元数值模拟方法分析其震间变形过程,同时以十余年的震间GPS速度场为约束,估算不同模型下的最佳震间闭锁深度。结果表明,在厄瓜多尔俯冲区,基于弹性模型的最佳震间闭锁深度约为30~35 km,而基于粘弹性模型的最佳闭锁深度为25 km。同时,相比于弹性模型,粘弹性模型能够更好地解释震间地表变形,说明该区域深部粘弹性圈层在俯冲板块的推挤作用下会引起较为明显的上地壳变形。     

一种有效的速度分析方法

如何得到高分辨率的叠加速度谱是高分辫率地震资料处理的关键.在常规速度分析基础上,介绍了自助速度分析法,它可以很大程度地提高叠加速度谱分辨率.理论及实例证明,该方法具有分辨率高、稳定性强的特点,提高了速度估算精度,特别在中浅层、大偏移距或薄层的情况下,也有很高的分辨率.得到的速度为后续步骤的工作顺利进行提供了保障.     

基于MapReduce的空间敏感性分析并行算法设计

近年来,随着遥感空间数据广泛应用于生态系统,推动了区域尺度生态遥感参数模型的发展。敏感性分析对识别模型关键参数,降低模型不确定性和完善模型具有重要作用。区域尺度的生态遥感参数模型,在进行模型参数敏感性分析时,由于涉及到空间数据的复杂运算,单机环境无法满足快速分析的要求。为了提高生态遥感参数模型空间敏感性分析效率,本文以青藏高原为研究区域,利用植被光合模型VPM(Vegetation Photosynthesis Model)和开源云计算平台Hadoop,设计和实现了基于Sobol′的生态遥感参数模型空间敏感性分析并行算法,并在实验室集群环境下进行算法分析,验证了算法的有效性和适用性。该算法的核心是利用MapReduce并行编程技术,对空间敏感性分析中的地图抽样和模型迭代过程进行任务分割,将分割后的子任务分配至不同的计算节点进行并行计算。实验表明,本文提出的并行策略,能有效缩短地图抽样和模型迭代计算时间,相比于单机算法,并行算法的运行速度提高了14倍左右。     

Elman神经网络在区域速度场建模中的应用

鉴于欧拉矢量参数在局部区域适应性较差,提出一种基于Elman神经网络的速度场逼近方法。首先利用已有欧拉矢量参数估计站点速度,将剩余残差作为Elman神经网络拟合量进行逼近;然后将Elman神经网络估计结果与欧拉矢量计算速度相叠加,获得区域速度场模型。利用山东区域速度场数据进行验证,结果表明,该方法在一定程度上能够削弱系统误差影响,提高计算精度。     

数字高程模型精度评定的基本理论

在剖析随机误差与逼近误差区别的基础上 ,指出数字高程模型 (DEM)内插误差属于逼近误差 ,精度评价方法不应继续延用“中误差”概念 ,而应采用“逼近误差”,给出了常用 DEM内插方法的误差评定模型 ,分析了等高线间距对 DEM内插精度影响的作用机制 ,最后解释了 DEM精度分析中用“中误差”概念无法解释的几个问题。     

基于虚拟植物的幼龄杉木生长模拟

 为估算树木在不同生长环境下的生长量及评价树木与环境之间的相互关系,本文以幼龄杉木为对象,从个体尺度上,根据形态结构生长发育规律,定义了幼龄杉木生长的基本单元,耦合用于表达树木形态发育的L-系统和生理生态模型,形成了L-系统的杉木结构-功能模型。在LSTree软件基础上,集成光合作用、生物量分配等模型,扩展了冠层光合有效辐射(PAR)的空间分布模拟,计算每个叶片上分布的PAR值,为光合作用模型提供核心参数。以福州3-4年生杉木的生长盛期为例,计算不同时刻净光合速率,从而计算各生长阶段的生物量,再把生物量分配到各器官。在生物量驱动下,树木形态结构发生变化,通过迭代,实现生物量驱动下的树木动态生长模拟。净光合速率估算结果与已有相关研究对比,比较合理。这种方法可定量分析、解释植物与环境之间关系。     

A simple lumped mass model to describe velocity of granular flows in a large flume

This paper presents a lumped mass model to describe the run-out and velocity of a series of large flume tests,which was carried out to investigate some propagation mechanisms involved in rapid,dry,dense granular flows and energy transformation when the flows encountered obstacles and reoriented their movement directions.Comparisons between predicted and measured results show that the trend of predicted velocities was basically matched with that of measured ones.Careful scrutiny of test videos reveals that subsequent particles with a higher velocity collided with slowed fronts to make them accelerate. However,this simple model cannot reflect collisions between particles because it treated released materials as a rigid block.Thus,the predicted velocity was somewhat lower than the measured velocity in most cases.When the flow changed its direction due to the variation in slope inclination,the model predicted a decrease in velocity.The predicted decrease in velocity was less than the measured one within a reasonable range of 10% or less.For some cases in which a convexity was introduced,the model also predicted the same trend of velocities as measured in the tests.The velocity increased greatly after the materials took a ballistic trajectory from the vertex of the convexity,and reduced dramatically when they finally made contact with the base of the lower slope.The difference between prediced and measured decrease in velocity was estimated to be about 5% due to the landing.Therefore,the simple lumped mass model based on the energy approach could roughly predict the run-out and velocity of granular flows,although it neglected internal deformation,intergranular collision and friction.     

����Vondrak�ֲ�̽��Ŀ���EKF�㷨��GPS�����е�Ӧ��

???????Vondrak?????????????EKF??,???????GPS??????λ??????????????EKF??????????????????Vondrak??????????????????μ????EKF???????????????д???????λ????????????????EKF???????????????????????д???????????????EKF????????????????д?????????????????????????????????????п????????????????????Ч???     

�������������������ص�FFT�����㷨�Ƚ�

???????????????????????????????Possion???????????????????????FFT????????????????????????????????λ???EGM2008??????????????????????????????????????????о???????????????????????????￡???Possion?????Ч????á?     

卫星重力梯度数据的向下延拓

在利用空域法恢复地球重力场时,向下延拓是卫星重力梯度数据预处理必不可少的步骤。将航空重力数据处理中的球内Dirichlet法、泊松积分迭代法、谱方法引入卫星重力梯度数据的向下延拓中,建立了相应的数学模型,解决了传统的球内Dirichlet问题存在的数值矛盾,利用模拟的卫星重力梯度数据对3种方法的向下延拓效果进行了分析和讨论。结果表明：当延拓距离为5 km时,谱方法所获得的延拓结果精度最高,其次为球内Dirichlet法,泊松积分迭代法精度最差;当延拓距离是250 km时,泊松积分迭代法的精度最好,其次为谱方法,球内Dirichlet法的精度最差。     

����Landweber�������ȿ���������С���������㷨�о�

???????Landweber???????x ^?? _k+1=x ^?? _k-f??(x ^?? _k) ^*(f(x ^?? _k)-y ^??) ???з??????????????????????????????????????????????????1/????? |(B??(x _k)B(x _k))|???????????μ???????????????С?????????Landweber???????, x ^?? _k+1=x ^?? _k-??(B??(x ^?? _k)B(x ^?? _k)) ^*B??(x ^?? _k)(f(x ^?? _k)-y ^??) ??????????????????????????????????????????????????????????????????ó?????[WTBZ]Landweber??????????????????????????????????????????С????????????     

利用地震矩张量反演鲜水河断裂带现今运动学特征

本文探讨了利用地震矩反演断裂形变带运动学参数的基本理论和方法，将其初步应用于鲜水河断裂形变带变形分析和运动机制的研究。结果表明，鲜水河断裂带呈现出走向拉伸、倾向压缩的形变格局，由地震矩反演的断裂带剪切形变速率（１０．９ｍｍ／ａ）与用地质学估算方法（１７ｍｍ／ａ）和现今地壳形变测量（８ｍｍ／ａ）的结果相当。同时，反演出的应变主方向能解释鲜水河断裂现今活动分段性特征以及多种滑动方式共存的现状，从而证明该方法是目前研究区域运动学问题切实有效的手段之一。     

病态总体最小二乘问题的谱修正迭代法

针对病态问题岭估计法及正则化方法存在破坏方程的等量关系、解是有偏估计的问题,在最小二乘谱修正迭代法的基础上,提出病态总体最小二乘问题的谱修正迭代法,并推导了总体最小二乘谱修正迭代法及其改进算法的具体公式。通过算例验证和阐明了TLS谱修正迭代法在病态问题中的有效性和易受迭代初值影响的缺点。     

Density Perturbation and Energy Flux of Internal Waves from Velocity Data

Internal waves play a crucial role in ocean mixing, and density perturbation and energy flux are essential quantities to investigate the generation and propagation of internal waves. This paper presents a methodology for calculating density perturbation and energy flux of internal waves only using a velocity field that is based on linearized equations for internal waves. The method was tested by numerical simulations of internal waves generated by tidal flowing over a Gaussian topography in a stratified fluid. The density perturbations and energy fluxes determined using our method that only used velocity data agreed with density perturbations and energy fluxes determined by the equation of state based on temperature data. The mean relative error (MRE) and root mean square error (RMSE) between the two methods were lower than 5% and 10% respectively. In addition, an experiment was performed to exam our method using the velocity field measured by Particle Image Velocimetry (PIV), and the setup of the experiment is consistent with the numerical model. The results of the experiments calculated by the methods using PIV data were also generally equal to those of the numerical model.     

独立变量频率法计算P—T大气稳定度频率及其误差分析与订正

本文把确定大气稳定度的因子:太阳高度角、云量和风速看作相互独立的变量,并用韦伯分布拟合风频,通过计算三个独立变量的频率的方法计算大气稳定度频率。成都地区实例计算表明,该法与传统的P-T法统计结果十分吻合,但却具有所需资料少,计算工作量小等明显优点,通过对该法的误差进行分析,提出合理的订正方法,得出精确的订正结果。并提出更精确的三参数风频拟合设想。通过与传统P-T法比较,指出了传统P-T法风速分级与实测风速整数取值之间的矛盾,给出了切实可行的解决方法。     

VARIATIONAL ITERATION SOLVING METHOD FOR SEA-AIR OSCILLATOR MODEL OF INTERDECADAL CLIMATE FLUCTUATIONS

1IN T R O D U C T IO NInteractionsbetween the ocean and atmosphere con-tributeto climatefluctuationosver a broad spectrum oftime scales.Studiesof those interactionhsave farfo-cused on El Ni o-SouthernOscillation(ENSO)phe-nomenon thathas a period of3to4yea…