共查询到19条相似文献,搜索用时 609 毫秒
1.
2.
3.
分析了求解变分不等式Newton方法的半局部收敛性,建立了类似于Kantorovich定理的收敛性结果。该结果不仅为判断Newton方法的收敛性提供了可计算的充分条件,也给出了Newton方法的收敛域以及问题解的存在区域。同时,文章也得到了Newton方法的若干收敛性质,包含收敛阶以及可计算的误差估计式等。 相似文献
4.
ρ(L_(r。w))的界和迭代法收敛性刘光清,辛迪远(长春地质学院基础科学系.长春130026)(四平电视大学.四平136000)本文给出了广义AOR法(广义SSOR,广义SOR,广义G-S及广义西法)迭代矩阵谱半径较理想的上界及其收敛准则,推广了文?.. 相似文献
5.
利用Matlab函数识别沉积物中的米兰柯维奇旋回信号 总被引:1,自引:0,他引:1
针对旋回分析中需要对数据进行的各种处理,利用Matlab中集成的滤波(filter)函数、插值(interp1)函数、求均值(mean)函数、去线性趋势(detrend)函数、求方差(variance)函数及谱分析(pwelch)函数等的组合,编制了识别沉积物中米兰柯维奇旋回信号的程序。选择窗口宽度(window)约等于有用信号中最低频率信号的波长、交叠点数(noverlap)为窗口宽度的50%、傅里叶变换的数据点数(nfft)大于全部采样数据的长度且为2的最小整数次幂时,可获得较好的谱估计结果。其不足之处是不能估计背景的红噪声谱及给出谱峰值的置信水平。 相似文献
6.
井眼轨道的软着陆设计模型的求解可以归结为一个七元非线性方程组的求解问题。前人给出了数值迭代求解算法,然而并没有证明该迭代算法的收敛性,并且该算法是否收敛严重依赖于用户给出的迭代初始值。通过一系列的消元、化简的数学技巧,将七元非线性方程组化简为一元多项式方程,并在此基础上给出了软着陆设计模型的一个新算法。理论分析和实际算例表明,新算法的主要计算工作量是求多项式方程的非负实数根,其他未知数与实数根是简单的函数关系,计算量很小。新算法克服了迭代算法的初值依赖性以及迭代过程可能发散等缺陷,并且在设计模型有多个解的情况下,可以同时求出这些解。 相似文献
7.
大地电磁反演问题通常表述为目标函数最优化,难点是多参数、非线性和不适定性,局部和全局方法都不能实现快速全局优化[4].针对局部线性方法易使解陷入局部极值,严重依赖初始模型,而传统的遗传算法在优化应用中存在局部搜索能力弱、早熟收敛等问题.这里引进一种求解一维大地电磁测深反演问题的实数编码广义遗传算法.该算法利用拟网格法初始种群和综合交叉策略,克服了早熟收敛现象,从而提高了遗传寻优的效率.理论模型反演与其它方法比较,结果说明遗传算法具有不依赖初始模型,不容易陷入局部极小,多点多路径概率搜索,以及隐合并行性等优点. 相似文献
8.
9.
浸润线的计算是尾矿坝安全状况评估的重要环节,对尾矿库安全起到至关重要的作用。为了能够用解析的方法得出浸润线,首先略掉无关紧要的部分,将尾矿坝的渗流分析区域简化为一个矩形,然后在矩形区域中选定两种符合工程实际的边界条件,最后求解渗流微分方程即求出这两种边界条件下的浸润线。求解中涉及到的偏微分方程无论是齐次还是非齐次的,都通过傅里叶级数法得到收敛的无穷级数解。该方法的突出特点就是得到的最终结果是解析解。实际应用时通过与实测浸润线数据拟合确定浸润线求解公式中的常数。通过对某尾矿坝浸润线的计算发现前面得到的公式可以用于实际浸润线的计算。 相似文献
10.
Laplace方程Cauchy问题的一种数值解法 总被引:2,自引:1,他引:2
研究了Laplace方程Cauchy问题的数值求解,该问题是一个典型的病态问题。利用格林(Green)公式将Laplace方程的Cauchy问题转化为Hausdorff矩问题。本文利用一种新方法,即矩问题的积分方程方法,求解矩问题,设计了二维Laplace方程Cauchy问题稳定的算法,给出了近似解的误差估计,并对二维Cauchy问题进行了数值模拟。 相似文献
11.
基于三维非线性有限元的边坡稳定分析方法 总被引:3,自引:0,他引:3
刚体极限平衡法不能反映岩体中实际的应力分布,而基于有限元的强度折减系数法在判断收敛性方面存在一些问题。为了解决这些问题,采用多重网格法,分别建立用于有限元计算的结构网格和用于计算滑面稳定安全系数的滑面网格,可以方便地获得任意滑面或滑块的稳定安全系数,从而将非线性有限元和极限平衡分析结合起来。为了提高计算规模和计算精度,采用有限元并行计算程序TFINE.Pfem进行计算,分析了网格密度对计算结果精度的影响,并应用于锦屏高边坡的稳定分析中。与刚体极限平衡法结果的对比分析表明,由于考虑了计算过程中的非线性应力调整,该方法的计算结果比刚体极限平衡法偏大,而且更符合实际情况。 相似文献
12.
传递函数分析在边坡振动台模型试验的应用探讨 总被引:4,自引:1,他引:3
传递函数分析对获取振动系统的动力特性参数及估算动力响应具有重要意义。边坡地震模拟振动台模型试验可求算两种传函:反映边坡各点的场地地震动特征的绝对加速度传函(基于测点实测加速度时程)和反映坡体本身振动性质的相对加速度传函(基于测点与坡脚即台面间的实测加速度相对时程)。以弹性体振动方程分析和试验实测,探讨了边坡振动台模型试验中动力特性参数计算的传函选择和以传递函数对动力响应的估算,得出结论:(1)不同于机械模态分析(以绝对传函的模和虚部计算动力特性参数具有等价性),对于边坡模型试验,绝对传函实部和模混入了台面加速度影响,以绝对传函虚部或相对传函虚部才能正确计算动力特性参数;(2)传递函数和台面输入波形Fourier谱的逐频相乘得到动力响应估算谱,与测点实测响应Fourier谱吻合较良好,这种估算方法有助于对不同地震波作用下的边坡动力响应作出总体评价,并指导模型试验中地震波形的合理选取。 相似文献
13.
14.
The real spectral analysis of SP anomalies due to 2-D horizontal circular cylinder is carried out using Hartley transform. The Hartley transform is an alternate means of realizing spectral analysis in real domain unlike the Fourier spectral analysis. The Hartley transform yields a straightforward interpretation of SP anomalies caused by horizontal circular cylinder wherein all the parameters are derived independently as a function of frequency. A theoretical example illustrates the procedure. The effect of random noise, up to 10% of white Gaussian noise, on the interpretation scheme was also studied and found to be of negligible. The field example of the ??Sulleymonkey?? anomaly in the Ergoni copper district, Turkey exemplifies the applicability of the proposed method. 相似文献
15.
Binh-Le Ly 《国际地质力学数值与分析法杂志》1978,2(4):405-409
A method is derived to estimate the degree of consolidation of a saturated, non-homogeneous clay layer. In this method, the diffusion equation govering the dissipation of pore water pressure is transformed so that the pore water pressure can be expressed by a Fourier series. The Fourier coefficients are then evaluated along the trajectory of the diffusion process. The Fourier series does not converge point-wise to the pore water pressure; instead, it converges to the pore water pressure in the mean., The Fourier series can be used to estimate the degree of consolidation, because, by Schwartz's inequality, convergence in the mean implies convergence in the degree of consolidation. 相似文献
16.
本文给出了瞬变弹性动力问题的一个边界元法。该法是利用威尔逊——θ法的差分公式,把运动方程化为椭圆型微分方程。根据贝蒂定理和动力点荷载的特解,可获得动力问题的边界积分方程。这个解法是在真实时间域内逐步求解的,不需要使用拉氏变换。一个应力波传播的数值算例证实了该方法使用方便,且解答精度较高。 相似文献
17.
A boundary element method formulation in the real time domain in the transient elastodynamic problems is presented. The equation motion of the linear isotropic elastiic bodare transformed into the elliptic partial differential eqnations by means of the differential formulations of the Wilson-θ method. By utilizing the particular solution corresponding to a dynamic point force in an infinite medium and the Betti’s theorem the bundary integral equations may be obtain. As this approach formulates and solves the problem in real time domain in conjunction with a time-stepping algorithm, it does not require a nimerical inversion.Which transformes the results from the Laplace or Fourier ransformed domain to the time domain. A numerical example was shown in order to demonstrate the applicability of the present method to the practical problems. 相似文献
18.
针对向下延拓的不适定问题,提出了一种新的位场向下延拓方法-导数迭代法。从位场垂向一阶导数的定义出发,将观测面和向上延拓、向下延拓同等高度平面上的位场值近似联系起来,采用迭代法进行逐次逼近,推导出空间域位场向下延拓的导数迭代法的递推公式。考虑到空间域向上延拓积分方程实现的复杂性,对递推公式进行快速傅里叶变换,整理得到波数域中的导数迭代公式,同时证明了该方法的收敛性。模型检验和实例分析均表明迭代法相对直接FFT法具有稳定性强和下延深度大的优点。 相似文献
19.
波动方程角度域共成像道集 总被引:1,自引:0,他引:1
针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础. 相似文献