On the inversion of geodetic integrals defined over the sphere using 1-D FFT     

R. V. García  C. A. Alejo 《Journal of Geodesy》2005,79(6-7):331-340

An iterative method is presented which performs inversion of integrals defined over the sphere. The method is based on one-dimensional fast Fourier transform (1-D FFT) inversion and is implemented with the projected Landweber technique, which is used to solve constrained least-squares problems reducing the associated 1-D cyclic-convolution error. The results obtained are as precise as the direct matrix inversion approach, but with better computational efficiency. A case study uses the inversion of Hotine’s integral to obtain gravity disturbances from geoid undulations. Numerical convergence is also analyzed and comparisons with respect to the direct matrix inversion method using conjugate gradient (CG) iteration are presented. Like the CG method, the number of iterations needed to get the optimum (i.e., small) error decreases as the measurement noise increases. Nevertheless, for discrete data given over a whole parallel band, the method can be applied directly without implementing the projected Landweber method, since no cyclic convolution error exists.  相似文献   

一种新的共线方程线性化方法及应用     

李长武  耿则勋 《测绘学院学报》2005,22(2):122-124

分析了共线条件方程两种不同形式之间的关系，提出了一种新的共线条件方程线性化方法。该方法不仅使线性化过程通俗易懂，而且线性化后各改正数的系数简单整齐。文中还将这种新的线性化共线条件方程用于空间后方交会实验，取得了与常规方法完全一致甚至更好的结果，从而说明该方法可以用于摄影测量的研究与生产。  相似文献   

无伸缩缝桥梁的动力特性计算与试验研究   总被引：2，自引：0，他引：2  

彭大文  洪锦祥  郭爱民  方德铭 《地震工程与工程振动》2005,25(2):72-76

福建省永春县上坂大桥是目前我国修建桥长最大，采用整体式桥台的无伸缩缝桥梁。本文通过该桥的动力试验研究了该桥在天然脉动荷载作用下的自振特性，采用结构分析软件ANSYS建立全桥的空间模型，计算该桥的自振特性，与实测结果进行了分析比较。实桥试验还测试了无障碍行车试验和跳车试验这2种车辆激振作用下桥梁的强迫振动特性，并把实测的荷载冲击系数与新桥规的计算结果进行比较。研究结果有助于对该类桥梁动力特性的认识，有利于进一步开展无伸缩缝桥梁的抗震性能研究。  相似文献   

结构-地基动力相互作用时域数值计算模型研究   总被引：3，自引：1，他引：3  

李建波  林皋  陈健云  胡志强 《地震工程与工程振动》2005,25(2):169-176

完全在时域数值分析范畴内，本文对结构一地基动力相互作用系统运动方程进行了深入推导，给出了统一的基本表达形式，简化可获得目前应用的各类主要时域模型。其中，基于结构一地基交界面的位移协调条件及作用力与反作用力平衡方程，推导了动力相互作用分析子结构法模式的时域基本方程。从结构的一般动力平衡方程出发，推导了相互作用直接法模式的基本平衡方程，并基于有限元模型对两类模式间的等价性条件进行了理论验证。  相似文献   

岩石物理的理论模拟和数值实验新方法   总被引：12，自引：13，他引：12       下载免费PDF全文

陶果  岳文正  谢然红  朱益华 《地球物理学进展》2005,20(1):4-11

当代物理实验科学和计算机科学技术的巨大进步为岩石物理的理论和实验研究开辟了新的领域.在细胞自动机基础上发展起来的格子气自动机和格子玻尔滋曼方法就是这种技术进步的产物.本文讨论了格子气自动机的发展及其应用于研究多孔岩石微观孔隙结构对宏观导电特性影响的新方法. 我们开发的模拟多孔岩石导电特性的数值实验新方法，提供了一种可深入到从微观孔隙结构这个尺度上来考察具有复杂孔隙结构的，孔隙度小于10%的低孔、低渗岩石的宏观导电特性及其影响因素.研究表明:Archie 系数m,n并不是传统上的意义，研究各种不同条件下F-、I-Sw关系，结果显示非阿尔奇现象与孔隙介质中孔隙、骨架和流体的混合导电网络有关，得出新的饱和度计算公式.同时还研究了泥质含量和分布、电阻率各向异性对宏观导电特性的影响，提出适用于电阻率各向异性地层的F-、I-Sw的关系.本文最后探讨了应用格子气自动机和格子玻尔滋曼方法进行岩心渗流和核磁共振数值实验的新方法及其发展前景.  相似文献   

植被—侵蚀状态图在典型流域的应用   总被引：11，自引：1，他引：11  

王兆印  郭彦彪  李昌志  王费新 《地球科学进展》2005,20(2):149-157

运用植被—侵蚀动力学耦合方程组，研究了流域植被与侵蚀在人类活动特别是植树造林和工程治理影响下的演变。动力学方程组中的4个参数由实测植被侵蚀资料算出，利用这些参数发展出植被—侵蚀状态图作为分析工具，分析在停止人类干预后植被和土壤侵蚀的发展趋势。流域的植被—侵蚀状态有4种，即向着侵蚀模数增大和植被覆盖度降低方向发展的A区，向着侵蚀模数降低和植被完善方向发展的C区，植被和侵蚀模数都增长的B区，植被和侵蚀模数都减少的D区。将耦合方程组和植被—侵蚀状态图应用于黄土高原王家沟、安家沟流域，云南小江流域和北京西山地区，结果表明动力学方程组准确地模拟了植被和滑动侵蚀模数的演变过程。对植被侵蚀状态图的分析表明，在干旱寒冷的黄土高原增加植被对于减少侵蚀的作用显著，治理应该首先增加植被使其进入D区，再结合控制侵蚀进入C区，治理后仍应加强管理和防止人为破坏；在雨量较丰、气温较高的小江流域，减少侵蚀对改善植被的作用大，治理措施首选控制侵蚀；北京西山地区介于两者之间。  相似文献   

运用半离散中心迎风格式计算二维浅水方程的研究          下载免费PDF全文

陈建忠  史忠科 《水科学进展》2005,16(6):853-857

以三阶中心加权本质无振荡重构为基础,采用一维一维进行计算的方法,给出了求解二维浅水方程的高分辨率三阶半离散中心迎风格式.引入的重构方法既提高了格式的精度,又保证格式是无振荡的.时间的离散用最优的三阶SSP(Strong Stability Preserving)Runge-Kutta方法.源项的离散用辛普森公式.计算方法保持了中心差分格式简单的优点,即不需用黎曼解算器和进行特征分解过程.数值模拟结果与其它方法所得结果一致,表明了方法的有效性和稳定性.  相似文献   

第二类非线性Fredholm型积分方程数值解   总被引：1，自引：0，他引：1  

刘唐伟应正卫  吴志强 《华东地质学院学报》2005,28(3):294-296

配置法研究了地球物理中常见的第二类非线性 Fredholm 型积分方程的数值解法,将第二类非线性 Fredholm 型积分方程转化为非线性代数方程组进行求解,采用高斯数值积分公式,给出了数值计算的具体实例.利用Matlab软件的符号运算功能编程计算,克服了非线性方程难于变成求解的困难,数值例子表明该方法编程简便有效.对非线性积分方程和非线性代数方程组的求解都有重要价值.  相似文献   

Accuracy and efficiency considerations for wide-angle wavefield extrapolators and scattering operators     

C. J. Thomson 《Geophysical Journal International》2005,163(1):308-323

  相似文献   

Full three-dimensional tomography: a comparison between the scattering-integral and adjoint-wavefield methods   总被引：2，自引：0，他引：2  

Po Chen  Thomas H. Jordan  Li Zhao 《Geophysical Journal International》2007,170(1):175-181

This paper analyses the computational issues of full 3-D tomography, in which the starting model as well as the model perturbation is 3-D and the sensitivity (Fréchet) kernels are calculated using the full physics of 3-D wave propagation. We compare two formulations of the structural inverse problem: the adjoint-wavefield (AW) method, which back-propagates the data from the receivers to image structure, and the scattering-integral (SI) method, which sets up the inverse problem by calculating and storing the Fréchet kernels for each data functional. The two inverse methods are closely related, but which one is more efficient depends on the overall problem geometry, particularly on the ratio of sources to receivers, as well as trade-offs in computational resources, such as the relative costs of compute cycles to data storage. We find that the SI method is computationally more efficient than the AW method in regional waveform tomography using large sets of natural sources, although it requires more storage.  相似文献