spvfly-CSDN博客

原创非高斯统计（L2范数）的地球物理反问题

最小二乘法和L2L_2L2范数在反问题求解中应用广泛，同时这种情形可以通过假设数据与先验模型参数遵从高斯统计来合理解释。但这个假设并不总是恰当的，有些数据集服从其他分布。双边指数概率密度函数就是一个简单的变通选择。假设具有相同的均值⟨d⟩\left \langle {d}\right \rangle⟨d⟩和方差σ2\sigma^2σ2，对高斯和指数分布两者进行比较：两者的概率密度分布如下图所示，红色曲线对应指数分布，蓝色对应高斯分布，可以看到指数分布有更长的拖尾。那么也就是说，对于远离⟨d⟩\l

2020-07-14 12:28:48 563

原创地球物理反问题的非唯一性和局部平均

0对于线性的反问题Gm=d\mathbf{Gm=d}Gm=d，有三种不同的方法求解：长度理论广义逆概率密度上述三种方法从不同的角度出发，但其推导的最终结果是相同的。而且，对于反问题可以从向量空间的角度来进行考虑，这一观点对于上述三种方法均是适用的，用向量空间可以更容易和形象的去探索非唯一性的本质。1 非唯一性和零向量的关系假设，反问题Gm=d\mathbf{Gm=d}Gm=d的解不是唯一的（非唯一性），那么必然至少存在两个不同的解，记为m(1),m(2)(m(1)≠m(2))\mathbf

2020-07-14 09:59:09 464

原创有限单元：单元+自然坐标+插值函数

有限元某些边界问题可转变为变分问题（即泛函极值问题）。泛函的一般形式是函数的积分，有限元法是求解泛函极值的一种方法。在有限元法的计算过程中，会遇到大量积分计算，用一般的方法计算积分，常常很麻烦。为了简化这些计算，引入自然坐标和等参单元。自然坐标有限元法将区域剖分为单元，采用单元上的局部坐标，有利于积分计算。自然坐标是一种局部坐标，用一组不超过1的无量纲数来确定单元上点的位置。长度坐标（一维自然坐标）定义...

2020-07-12 17:22:04 2098 1

原创地球物理的向量空间

0 模型空间和数据空间从向量空间观点来看，将数据和模型参数写为向量d⃗\mathbf{\vec{d}}d和m⃗\mathbf{\vec{m}}m（以后简写为d\mathbf{{d}}d和m\mathbf{{m}}m），那么可将向量空间记为S(d)\mathbf{S(d)}S(d)和S(m)\mathbf{S(m)}S(m)。如下图所示，某一d\mathbf{{d}}d和m\mathbf{{m}}m均可表示为空间中的一个向量。这样的话：线性方程d=Gm\mathbf{d=Gm}d=Gm可以解释为向量从

2020-07-11 22:21:53 435

原创线性高斯反问题的解--最大似然法

0 一组测量的均值假设，在一个实验中对一模型参数m1m_1m1在有噪声的情况下测量了N次，每次都采集到一个独立数据did_idi。那么，如果随机变量是高斯型的，那么它们的联合概率密度函数能够以方差σ2\sigma^2σ2和均值μ\muμ的形式来描述，即p(d)=σ−N(2π)−N/2exp⁡[−12σ−2∑i=1N[di−μ]2]p(\mathbf{d})=\sigma^{-N}(2 \pi)^{-N / 2} \exp \left[-\frac{1}{2} \sigma^{-2} \sum_{i

2020-07-08 16:57:05 650

原创线性高斯反问题--广义逆

0 解与算符对于线性反问题Gm=d\mathbf{Gm=d}Gm=d的方法，长度方法基于分析解的两个属性：预测误差和解的简单程度。而广义逆的方法，则希望通过研究算符矩阵来获得更多关于反问题的属性。假设求解反问题Gm=d\mathbf{Gm = d}Gm=d，存在一个广义逆可求解模型估计，mest=G−gd\mathbf{m^{est}=G^{-g}d}mest=G−gd，这里的G−gG^{-g}G−g为广义逆。通过广义逆，可以求取：1）数据分辨率矩阵2）模型分辨率矩阵3）单位协方差矩阵1 分辨

2020-07-04 16:43:14 740

原创线性高斯反问题的解--长度方法1

3.1 长度首先，以一个简单的线性回归问题（一条直线拟合数据的问题），来直观的体会“长度”。a)对数据（z,d）的最小二乘直线拟合；b)ei=diobs−dipree_i = d_{i}^{obs} - d_{i}^{pre}ei=diobs−dipre那么，最佳的拟合直线所具有模型参数（截距和斜率）将使总误差E最小，即E=∑i=1nei2=eTe=minE = \sum\limits_{i=1}^{n}e_{i}^{2}=\bf {e^Te} = minE=i=1∑nei2=eTe=mi

2020-06-24 22:09:32 407

原创有限元法求变分问题示例

v[y]=∫01(12y′2+y)dx取极值v[y]=\int_{0}^{1}\left(\frac{1}{2} y^{\prime 2}+y\right) d x \quad 取极值v[y]=∫01(21y′2+y)dx取极值y(0)=0,y(1)=1y(0)=0, y(1)=1y(0)=0,y(1)=1第一步：区域剖分将自变量x的取值区域[0,1]，以步长h，剖分为n份，等分点x0=0,x1,x2,...,xn=1x_0=0,x_1,x_2,...,x_n=1x0=0,x1,x2,..

2020-05-09 16:02:01 1191

原创学习笔记(08):21天通关Python（仅视频课）-使用shell调试

本页面购买不发书！！！仅为视频课购买！！！请务必到https://edu.csdn.net/bundled/detail/49下单购买课+书。本页面，仅为观看视频页面，如需一并购买图书，请务必到https://edu.csdn.net/bundled/detail/49下单购买课程+图书！！！疯狂Python精讲课程覆盖《疯狂Py...

2020-04-04 11:31:43 116

原创学习笔记(07):21天通关Python（仅视频课）-Scrapy下载和安装及创建项目

本页面购买不发书！！！仅为视频课购买！！！请务必到https://edu.csdn.net/bundled/detail/49下单购买课+书。本页面，仅为观看视频页面，如需一并购买图书，请务必到https://edu.csdn.net/bundled/detail/49下单购买课程+图书！！！疯狂Python精讲课程覆盖《疯狂Py...

2020-04-04 10:51:31 122

原创地球物理正反演回顾

地球物理勘探，是指通过测量地球及其周边所存在的人工/天然物理场，研究与地下介质相关物理场的时空变化规律，进而解释推断地下介质分布、地下构造发育等地质相关的真实情况（状态）。从上述的概念来看，大致的逻辑规律应该是正过程：地球--》物理场--》观测数值，反过程：观测数值--》物理场--》地球。所以，地球物理学是一门交叉学科，即地质+物理+数学+....。对于特定的一个问题，通常的方法是采集--&g...

2020-03-01 11:46:18 4845