- 博客(22)
- 收藏
- 关注
原创 非高斯统计(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
原创 学习笔记(06):21天通关Python(视频课)-案例实操:使用urllib模块读取网络资源及提交请求(上)...
请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课,方可添加图书寄送地址。 本页面,仅为观看视频页面,如需邮寄图书,请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课 疯...
2020-02-03 10:13:08
111
原创 学习笔记(05):21天通关Python(视频课)-网络模块概述与urllib.parse模块讲解
请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课,方可添加图书寄送地址。 本页面,仅为观看视频页面,如需邮寄图书,请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课 疯...
2020-02-03 09:53:13
135
原创 学习笔记(04):21天通关Python(视频课)-案例实操:使用事务控制数据库操作
请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课,方可添加图书寄送地址。 本页面,仅为观看视频页面,如需邮寄图书,请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课 疯...
2020-01-30 15:23:54
75
原创 学习笔记(03):21天通关Python(视频课)-执行查询
请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课,方可添加图书寄送地址。 本页面,仅为观看视频页面,如需邮寄图书,请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课 疯...
2020-01-26 11:18:11
127
原创 学习笔记(02):21天通关Python(视频课)-执行DML语句
请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课,方可添加图书寄送地址。 本页面,仅为观看视频页面,如需邮寄图书,请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课 疯...
2020-01-26 11:00:55
111
原创 学习笔记(01):21天通关Python(视频课)-数据库API与全局变量及核心类基本流程...
请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课,方可添加图书寄送地址。 本页面,仅为观看视频页面,如需邮寄图书,请务必到https://edu.csdn.net/bundled/detail/49?utm_source=banner 下单购买课 疯...
2020-01-26 09:46:38
91
原创 MT2D大地电磁有限元正演程序--单元分析
MT2D有限元单元分析三角单元的形函数假定在三角单元内的电位UUU是线性变化的,即U=NiUi+NjUj+NmUmU=N_iU_i+N_jU_j+N_mU_mU=NiUi+NjUj+NmUm其中,NnN_nNn为形函数,n=i,j,mn=i,j,mn=i,j,mNn=12Δ(an+bnx+cny)N_{n}=\frac{1}{2 \Delta}\left(a_{n}+b...
2019-11-13 21:21:50
1211
原创 大地电磁正演程序MT2D主程序分析
0 相关背景信息见大地电磁二维正演程序--详细介绍(https://blog.csdn.net/spvfly/article/details/92795423)。程序的框架如下图1 主程序MT2D.cpp代码分析01 头文件部分包含四个建好的类#include "mesh2d.h" //网格处理类 包含读入网格数据等功能#include "dofs.h" ...
2019-08-06 07:44:16
1845
原创 大地电磁二维正演程序--详细介绍
0 相关背景1、开发平台vs+intel MKL + Eigen2、正演方法为有限单元法3、程序分为两部分: 网格剖分---使用Triangle程序进行 正演计算---使用任政勇教授开源的MT2D程序(https://sourceforge.net/projects/mt2d/)1 程序大体流程 ...
2019-06-18 17:29:36
4193
原创 三角网格剖分工具 Triangle 安装及使用
0 Triangle资源 下载Triangle源码,见http://www.cs.cmu.edu/~quake/triangle.html1 Triangle安装 triangle目录下共有两个主体程序:triangle网格剖分程序,showme可视化程序。可以直接进入triangle目录,make编译。但因为其中showme只能在X Windows下使用,所以在CentOS下编译...
2019-06-17 10:27:27
6330
8
原创 vs 集成 Intel MKL + Eigen
0 Intel MKL 和 Eigen 简介 Intel数学核心函数库(MKL)是一套高度优化、线程安全的数学例程、函数,面向高性能的工程、科学与财务应用。英特尔MKL的集群版本包括ScaLAPACK与分布式内存快速傅立叶转换,并提供了线性代数(BLAS、LAPACK和Sparse Solver)、快速傅立叶转换、矢量数学(Vector Math)与随机号码生成器支持。具...
2019-06-16 22:01:52
1366
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人