共查询到20条相似文献，搜索用时 497 毫秒
1.
将TOMS测量的北半球不同纬度各纬向平均逐日O。总量作为无规序列,利用时滞方法重构序列的相空间吸引子,计算其系统状态演化过程中形成的整体结构特征.结果表明,各纬度O3层系统吸引子都具有非整数分维结构；此外,在完整的Lpapunov指数谱中,所有指数分量之和均为负值,且最大指数分量总是正值.反映出O3层系统相空间总体上是收缩的,同时提供了存在奇怪吸引子的证据,表征O3层系统的演化是一种耗散混饨运动.这对于进一步应用非线性动力系统理论和方法,深入研究O3层系统的动力学特征具有重要意义. 相似文献
2.
非线性地震道的混沌反演——Ⅱ.关于Lyapunov指数和吸引子 总被引：29，自引：4，他引：25
在本文上篇《非线性地震道的混沌反演--1.理论和数值试验》中讨论了非线性反问题逐次线性化方法及迭代过程的系统特性、递推法和广义线性反演法输出波阻抗序列的特征。本篇参照混沌理论中关于Lyapunov指数、相空间和吸引子的概念,深入研究逐次线性化迭代的状态变化的阶段性及其定量描述.当引用相空间描述系统的演化时,对于随迭代变化的阻尼因子可以发现几种吸引子.最后,总结了非线性地震混沌反演的要点,指出非线性反演迭代的复杂阶段性要川混沌理论来描述,而Lyapunov指数等才是控制迭代过程的实质性参数. 相似文献
3.
全文以地震道反演的数值实验为例,根据混沌理论讨论了非线性地震反演的以下特性:1.对于带误差的地震数据和极其平滑的初始模型,逐次线性化的迭代过程产生的输出序列最终走向无序,这是非线性反演系统本身的特性决定的;2.迭代（非线性反演）系统是由Poin-care映的或系统方程描述的;3.可以根据地震反问题的特点,用多种不同的方法定义相应的Lyapunov指数,它们的数值和组合是非线性反演系统状态的有效指示;4.取决于系统参数的选取,发现在二维相空间有奇异吸引子的对应物.总之,从混沌理论的角度来研究纯数学的非线性地震反问题,可以揭示某些过去鲜为人知的内在规律性.本文为全文的上篇,主要讨论前两个问题,即理论和地震道反演数值试验结果. 相似文献
4.
非线性地震道的混沌反演——Ⅱ.关于Lyapunov指数和吸引子 总被引：3，自引：0，他引：3
在本文上篇《非线性地震道的混沌反演——1.理论和数值试验》中讨论了非线性反问题逐次线性化方法及迭代过程的系统特性、递推法和广义线性反演法输出波阻抗序列的特征。本篇参照混沌理论中关于Lyapunov指数、相空间和吸引子的概念,深入研究逐次线性化迭代的状态变化的阶段性及其定量描述.当引用相空间描述系统的演化时,对于随迭代变化的阻尼因子可以发现几种吸引子.最后,总结了非线性地震混沌反演的要点,指出非线性反演迭代的复杂阶段性要川混沌理论来描述,而Lyapunov指数等才是控制迭代过程的实质性参数. 相似文献
5.
基于主成分的时间域航空电磁数据神经网络反演仿真研究(英文) 总被引：5，自引：0，他引：5
传统上,时间域航空电磁数据通过拟合迭代反演计算得到大地模型,然而,由于航空电磁数据道间的较强相关性,导致病态反演,并引起超定问题;同时电磁数据的相关性使其与模型参数的映射关系复杂,增加了反演的复杂度。采用主成分分析法将航空电磁数据变换为正交的较少数量的主成分,不仅降低了数据道间的相关性,减小了数据量,同时压制了数据的不相关噪声。本文利用人工神经网络(ANN)逼近主成分与大地模型参数间的映射关系,避免了传统反演算法中雅克比矩阵的复杂计算。层状模型的主成分神经网络与数据神经网络的反演结果对比显示,主成分神经网络反演方法网络结构简单,训练步数少,反演结果好,特别是对于含噪数据。准二维模型的主成分ANN、数据ANN以及Zhody方法的反演结果显示了主成分神经网络具有更接近真实模型的反演效果,进一步证明了主成分神经网络反演方法适合海量航空电磁探测数据反演。 相似文献
6.
通过海城地震余震活动序列资料，根据相空间中余震活动的演化特征，余震活动的功率谱和非线性地震活动演化方程的迭代解详细地研究了余震活动的基本特征．结果显示，余震活动除了通常熟悉的平均的衰减和周期性行为之外，还存在同样显著的混沌行为． 相似文献
7.
结构动力反应分析中的一种显式输入反演方法 总被引：1，自引：0，他引：1
本文简要总结，评述了结构动力反应分析中两类输入反演方法（脉动函数法和直接时域法）的优缺点，在此基础上通过应用文献[1]中提出的3阶显式方法，获得了一个关于结构输入反演问题的新解法－显式直接时域法，新方法集合了脉冲函数法和直接时域法的优点，既可按显式求解动力微方程，又适用于非线性系统的动力反演问题，最后，通过算例对该方法进行了数值验证。 相似文献
8.
9.
蕴震系统是在多种环境动力因子作用下，由多个不同性质、不同层次的子系统强烈相互作用所形成的“复杂巨系统”或“复杂的非线性动力学系统”。本文概述了蕴震系统的自组织演化理论，提出并采用蕴震系统方向信息合成法研究了相应的构造微动态特征。 相似文献
10.
可控源音频大地电磁数据的反演方法 总被引：18，自引：7，他引：18
从反演方程、构造目标函数和求解三方面对用于可控源音频大地电磁法(CSAMT)的实用反演方法中的四种进行了描述．水平层状地层CSAMT法资料的直接反演法首次尝试了一维空间的全资料CSAMT反演，效果较好，但该方法尚难应用于2D、3D复杂介质中；奥克姆反演方法既考虑了横向的光滑函数，又考虑了纵向的光滑函数，得到比较光滑的横向、纵向变化的背景电性结果，但有可能把一些小构造光滑掉．快速松驰反演算法和共轭梯度算法由于计算速度快，占内存少而被用于三维反演中，二者相比，快速松驰算法在求解雅可比矩阵时只做一次正演计算，在更新模型时解小型方程组，所以在速度上更胜一筹．在后三种算法中，由于复杂电性结构无解析解，正演计算都采用数值计算．数值计算的可靠性、速度影响着反演算法的有效性，这方面的研究也将是2D、3D复杂电性结构反演的研究方向之一. 相似文献
11.
在大跨度桥梁设计中,抗风是工程师们最关心的问题之一,也是桥梁动力稳定问题的重要研究领域。本文按照广义力和广义位移之间的对应关系,将相空间的弹性动力学的基本方程乘上相应的虚量,然后积分,代数相加,并考虑到体积力和面积力均为伴生力,建立了适用于研究大跨度桥梁流固耦合初值问题相空间的两类变量的广义拟余能原理。应用该原理,解决了一个非保守桥梁流固耦合问题的具体算例。并且讨论了如何应用非保守系统的广义拟变分原理进行大跨度桥梁近似计算的问题。 相似文献
12.
Based on the approximation by polynomial‐fraction, a series of systematic lumped‐parameter models are developed in this paper for efficiently representing the dynamic behaviour of unbounded soil. Concise formulation is first employed to represent the dynamic flexibility function of foundation with a ratio of two polynomials. By defining an appropriate quadratic error function, the optimal coefficients of the polynomials can be directly solved from a system of linear equations. Through performing partial‐fraction expansion on this polynomial‐fraction and designing two basic discrete‐element models corresponding to the partial fractions, systematic lumped‐parameter models can be conveniently established by connecting these basic units in series. Since the systematic lumped‐parameter models are configured without introducing any mass, the foundation input motion can be directly applied to these models for their applications to the analysis of seismic excitation. The effectiveness of these new models is strictly validated by successfully simulating a semi‐infinite bar on an elastic foundation. Subsequently, these models are applied for representing the dynamic stiffness functions for different types of foundation. Comparison of the new models with the other existing lumped‐parameter models is also made to illustrate their advantages in requiring fewer parameters and featuring a more systematic expansion. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
13.
A time domain boundary element in a cylindrical coordinate system is developed for the analysis of wave propagation in a half space. The integral formulation is based on Graffi's dynamic reciprocal theorem and Stokes' fundamental solutions. The field quantities (displacements and tractions) are expressed as products of Fourier series in the tangential direction and linear polynomials in the other spatial directions. Gaussian integration is used to integrate the non-singular parts of the integral equations, whereas the integration of the singular components, which are either of order 1/r or 1/r2, is handled by special numerical schemes. In the time marching aspect, the field quantities are assumed to vary linearly in the temporal direction as well. Examples for wave propagation due to various forms of surface excitations are reported to demonstrate the accuracy of the method. 相似文献
14.
Nonlinear seismic response analysis of earth dams 总被引：1，自引：0，他引：1
The objective of this paper is to propose a general and efficient numerical procedure for analysing the dynamic response of geotechnical structures, which are considered as both nonlinear and two phase systems. In Section 2, the appropriate coupled dynamic field equations for the response of a two-phase soil system are briefly reviewed. The finite element spatial discretization of the field equations is described and time integration for the resulting nonlinear semi-discrete finite element equations is discussed. In Section 3, iterative techniques are examined for the solution of the global nonlinear system of finite element equations. A large amount of computational effort is expended in the iterative phase of the solution and so the iterative procedure used must be both reliable and efficient. The performance of three iterative procedure is examined: Newton Raphson, Modified Newton Raphson and Quasi-Newton methods, including BGFS and Broyden updates. Finally, in Section 4, the elasto-plastic earthquake response analysis of a two phase nonhomogeneous earth dam is presented. Extensive documentation exists1 for the particular problem selected including recorded earthquake motions at the base and crest of the dam. The results of the numerical calculations are compared to the recorded response of the dam. 相似文献
15.
16.
A. S. Cherepantsev 《Izvestiya Physics of the Solid Earth》2007,43(12):1047-1055
The possibility of determining parameters of complex geophysical processes is considered in terms of nonlinear dynamics. In accordance with modern approaches in the theory of nonlinear dynamic systems, the number of independent parameters controlling the behavior of a nonlinear system can be estimated from the available time realization of only one of these parameters. Model calculations showed that the dimension of the phase space of a dynamic system can be estimated from a sample of one variable. Experimental data on variations in the apparent electric resistivity (AER) and the relative vertical movement of the surface (RVMS) in a seismically active region are analyzed and the dimension of the dynamic system determining its behavior in the phase space is estimated. The resulting estimates of the embedding dimension m = 7?8 for AER variations and m = 6 for RVMS variations possibly characterize the complexity of the dynamic system describing the given fields. The method presented in the paper is also applied to the analysis of the degree of connectivity of different dynamic systems and their parameters. By the connectivity we mean the number of independent parameters simultaneously involved in the formation of the dynamic behavior of various physical fields. The model estimates demonstrate the possibility of such an approach. It is shown that variations in the AER in perpendicular directions are described by a general system of dynamic equations, whereas dynamic systems controlling the AER field and the RVMSs are interconnected only partially. The resulting dimension m = 12 estimated for the AER-RVMS system provides an estimate for the number of common controlling parameters: n = (8 + 6)?12 = 2. The methods and results presented in the paper are applicable to the construction of models of complex geophysical processes and the development and the development of new approaches and methods of identification of prognostic characteristics for the behavior of physical fields of various origins. 相似文献
17.
In this paper, the dynamic stress concentration and scattering of SH-waves by bi-naterial structures that possess an interface elliptic cavity are investigated. First, by using the complex function method, the Green‘s function is constructed. This yields the solution of the displacement field for an elastic half space with a semi-elliptic canyon impacted by an anti-plane harmonic line source loading on the horizontal surface. Then, the problem is divided into an upper and lower half space along the horizontal interface, regarded as a harmony model. In order to satisfythe integral continuity condition,the unknown anti-plane forces are applied to the interface. The integral equations with unknown forces can be established through the continuity condition, and after transformation, the algebraic equations are solved numerically. Finally, the distribution of the dynamic stress concentration factor (DSCF) around the elliptic cavity is given and the effect of different parameters on DSCF is discussed. 相似文献
18.
The delay coordinate technique is examined as an indicator of the regime of particle dynamics for the system of single charged particle motion in magnetic reversals. Examples of numerically integrated trajectories in both static (zero electric field) and time dependent (corresponding nonzero induction electric field) simple models for magnetic reversals are considered. In the static case, the dynamics can in principle be directly classified by constructing Poincaré surfaces of section; here we demonstrate that whilst the Poincaré surface contains the relevant information to classify the dynamics, the corresponding delay coordinate plot can provide a far more sensitive indication of the onset of nonregular behaviour. In the case of nonperiodic time dependence considered here Poincaré plots cannot in general be constructed directly. Nevertheless, delay coordinate plots can still reveal details of the phase space portrait of the system, and here are shown to indicate whether segments of stochastic motion exist in a given trajectory. It is anticipated that the delay coordinate plot technique as realized here will be a valuable tool in characterizing the behaviour in large numbers of trajectories that are evolved in time-dependent systems, thereby giving us insight into the evolution of the distribution function as a whole, either in prescribed fields or in self-consistent numerical simulations. 相似文献
19.
A method is proposed to establish regional design hyetographs for facilitating the determination of design hyetographs at ungauged sites. The method is applied to the central area of Taiwan. First, the single‐station design hyetographs at all rain gauges are analysed using principal components analysis and cluster analysis. The principal components analysis shows that there are six dominant factors, and the cluster analysis indicates that the time to peak rainfall has the largest influence on the classification of hyetographs. It also shows that the single‐station hyetographs in the study area can be classified into three clusters. Finally, the homogeneous regions for these three clusters are delineated and the corresponding regional design hyetographs are proposed. Once the homogeneous regions and the regional hyetographs are available, the design hyetograph at the point of interest can be easily determined. The proposed method is expected to be useful for providing the design hyetographs at ungauged sites. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
20.
首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证. 相似文献