地球物理反演方法作业及答案.docxVIP

Homework 1: based on Appendix A 1. Norm of the vector（矢量的范数） 在欧式空间下，矢量可以表示为： a=(a1+a2+?+an) 与三维空间中矢量长度的定义类似，我们引入了n维欧式空间中矢量的范数，表示为： a=a12+a22+?+an2 当矢量a是由实数组成时，a的范数a始终是一个非负的值。显然，a的范数即为a与其自身的内积的平方根。 2. Operator; Linear operator（算子；线性算子） 在欧式空间中，对任意一个矢量a，我们通过某种规律的作用使其对应到同在欧式空间中的某个矢量a’，即： a=A(a) 我们就将这种规律称之为算子A。 对欧式空间中的任意两个矢量a1和a2，以及任意两个标量α1和α2，如果有 Aα1a1+α2a2=α1Aa1+α2A(a2) 成立，那么我们就称算子A是一个线性算子。 3. Give an example of operator behavior in geophysics（举出一个地球物理中的算子例子） 地震走时层析是一种常用的地震反演手段，其基本原理就是利用地震射线的路径与走时来建立观测时间与地下速度参数模型的关系。 假设地下模型由矩形网格组成，每个网格内的慢度Si为一个定值。那么，当震源设定好之后（假设震源和检波器都设置在模型边界），单个检波器接收到的走时T与速度之间存在如下关系： T=i=0nSiLi 其中，Li表示每个对应网格内地震射线的路径长度。 那么，对所有检波器接收到的走时数据，我们可以建立如下方程： T=LS 其中 T=T1+T2+?+TkT S是m*n行的列向量，代表模型依次每个网格的慢度。L是一个k行m*n列的矩阵，k是检波点的个数（即射线数量）。 因此，在求取地下模型速度信息的反演过程中，确定和构建矩阵L就成为了关键的因素。在这一地球物理反演问题中的算子A，就是矩阵L与慢度向量S的线性相乘。同时，在这里算子A还是一个线性算子。 4. Norm of operator（算子的范数） 某个算子的范数定义为满足下式： a≤Ma 中的所有M中最小的那个值，即： A=min?{M0,Aa≤Ma} 5. Norm of functional（泛函的范数） 泛函是算子的一种特殊形式，因此泛函的范数的定义与算子的定义类似。 某个泛函的范数定义为满足下式： f(a)≤Ma 中的所有M中最小的那个值，即： f=min?{M0,f(a)≤Ma} 6. Inner product（内积） 在有限维下，内积定义为两个矢量的一种运算： A?B=a1×b1+a2×b2+?+an+bn 推广到一般的线性空间中，内积是指在复数域上的一个线性空间V备有一个正定、对称以及共轭双线性形式。内积运算满足以下性质： 对称性：f,g=(g,f) 线性：f+g,h=f,h+(g,h) 线性：αf,g=α(f,g) 正定性：f,f0并且当且仅当f=0时，f,f=0 7. Hilbert space（希尔伯特空间） 在一个复矢量空间H上给定内积定义后，可以导出任意一个矢量的范数： x=(x,x) 如果空间H对于这个范数来说是完备的，则此空间称为希尔伯特空间。 希尔伯特空间拥有最丰富的几何性质。它是欧式空间的一个推广。在地球物理学中，我们可以把地球物理数据或模型看作是希尔伯特空间中的元素。比如说，我们可以利用相应数据之间的距离来计算观测数据和模拟数据之间的精度。换句话说，我们可以利用希尔伯特空间几何性质中所有的方法和捷径来帮助我们解决地球物理反演问题。 8. L1norm; L2norm; L2Cnorm; Lpnorm; L∞norm（L1范数； L2范数； L2C范数； Lp范数； L∞范数） L1范数：fL1=abfdx L2范数：f,gL2=abfgdx L2C范数：f,gL2C=abfg*dx LP范数：fLP=(abfPdx)1/P L∞范数：f∞=maxa≤x≤bf Homework 2: based on Appendix BCD 1. Inverse operator（逆算子） 对方程 Ax=y 来说，如果其解是唯一的，那么对解空间也即数据空间中的任意一个y’，在模型空间中都存在一个模型x’与之对应，这种y’与x’的对应关系就称为算子A的逆算子A-1 x=A-1y 2. Explain the procedure of Gram-Schmidt orthogonalization process and the role it plays.（阐述Gram-Schmidt正交化方法的流程以及它的作用） 在某些情况下，求解方程组 i=1nΓjiαi=(d0,dj) 是很困难的。显然，如果n个数据{d1,d2,d3,?,dn}是正交的，那么这