下载作业帮安装包
扫二维码下载作业帮
1.75亿学生的选择
拉普拉斯方程通解为何二维拉普拉斯方程Uxx+Uyy=0的解U(x,y)一定能写成f(x+iy)+g(x-iy)的形式
这是因为：Uxx+Uyy=0,通解为U(x,y)=f(x+iy)+g(x-iy),你求导就知道为什么了.具体,怎么算.你看看它的通解,是不是跟欧拉公式相似? 给我邮箱,我给你. 给你了,你邮箱那里
为您推荐：
其他类似问题
扫描下载二维码【复变函数课件第3章4调和函数】-学网-中国IT综合门户网站-提供健康,养生,留学,移民,创业,汽车等信息
【复变函数课件第3章4调和函数】
来源：互联网 更新时间： 22:53:47 责任编辑：李志喜字体：
本文由炅斯贡献
§3.4 解析函数与调和函数的关系 复习 二元函数 u( x , y ) u 一阶偏导数： 有两个 一阶偏导数： x , u y
有四个 二阶偏导数： xx , u xy , u 二阶偏导数： 二阶偏导数连续 如果二阶偏导数连续， 则 连续，
u yx , u yy u xy= u yx
偏导数, 偏导数, 且满足 Laplace方程 u xx + u yy Laplace方程 调和函数 则称 u( x , y ) 为区域D内 的调和函数 为区域D
二阶连续 在区域D 具有二阶 定义3.3 定义3.3 如果 二元函数 u( x , y ) 在区域D内具有二阶连续
都是D 的调和函数， 定义3.4 定义3.4 如果 u( x , y ) 与 v ( x , y ) 都是D内 的调和函数， 且满足 u x= v y , u y = ? v x 调和函数， 则称 v ( x , y ) 为 u( x , y ) 的共轭 调和函数，
P54定理 P54定理3.13 定理3.13
在区域D f ( z ) = u( x , y ) + iv ( x , y ) 在区域D内 解析 连续， ?ux , u y, v x , v y 在D内连续，且满足
连续. ? uxx , uxy , uyx , uyy 连续. 调和函数 ? uxx+ u yy= 0, u 是D内 的调和函数
? f ′( z ) = ux + iv x = v y ? iuy
v xx + v yy = 0,
u x= v y , u y = ? v x
解析， 解析，
调和函数 也是D v 也是D内 的调和函数
调和函数。 v 是实部 u 的共轭 调和函数。 在区域D ? f ( z )= u + iv 在区域D内 解析
调和函数 函数， 60页12 设 u( x , y ) 是D内 的调和函数， 60页
?u ?u 解析函数 函数？ , f ( z ) 是否是D内 的解析函数？ ?i f ( z) = 是否是D ?x ?y 是 调和函数 函数， 因为 u( x , y ) 是D内 的调和函数， ?u ?u 连续， 所以(1) 可微， 所以(1) u( x , y ) 的二阶偏导数 连续， , ? 可微， ?x ?y ? ?u ? 2u ? 2u ? ?u (2) 由 + 2 = 0 得到 ( ) = (? ) 2 ?y ?y ?y ?x ?x ?x ?u ? ? ?u ? 2u ? 2u ( ) = ? (? ) = 由 得到 ?x ?y ?y ? x ?x ?y ?y?x ? u ?u f ( z) = 解析函数 函数。 ?i 是D内 的解析函数。 ?x ?y 调和函数 函数， 若 v ( x , y ) 是D内 的调和函数， 同样可证： 同样可证： ?v ?v 解析函数 函数。 是D内 的解析函数。 +i 3 则 f (z) = ?x ?y
调和函数 调和函数 证明： 也是调和 例2 证明： 一对共轭 调和函数 的乘积 也是调和函数 调和函数， 证： 设在区域D内 v ( x , y ) 为 u( x , y ) 的共轭 调和函数， 在区域D 解析函数 则 f ( z )= u( x , y ) + iv ( x , y ) 是区域D内 的解析函数 是区域D 解析函数 [ f ( z )] = u 2 ? v 2 + i 2uv 是区域D内 的解析函数 是区域D
调和函数 函数， 2uv 是调和函数， uv
是调和函数. 调和函数 函数. 证明： 证明：
解析函数 是区域D 60页14.设 f ( z ) = u + iv 是区域D内 的解析函数 60页
也是解析函数， 解析函数 i f ( z ) = i (u ? iv) = v + iu = v ? iu 也是解析函数， 可微， 证： 由 f ( z ) = u + iv 解析得到 u, v 可微， 且 ?v ?u ?v ?u =? = , 可微，
因此 v , ? u 可微， 且 ?x ?x ?y ?y ?v ? ( ? u ) ?v ? ( ? u ) i f ( z ) = v ? iu 解析 =? = ?x ?y 4 ?y ?x ? u 是 v 的共轭 调和函数
例4 设 v ( x , y )= e p x sin y 求p的值，使 的值， 调和函数 并求出解析函数 f ( z ) = 函数， v ( x , y ) 为调和函数，
?v px 解 = pe sin y ?x ?v px = e cos y ?y
? 2v 2 px = p e sin y 2 ?x 2 ? v px = ?e sin y 2 ?y
p = ±1 f ′( z ) = ux + iv x = e x (cos y + i sin y ) p=1 = v y + iv x = e x cos y + ie x sin y z f ′( z ) = e z f ( z) = e + C ′( z )= e ? x cos y ? ie ? x sin y f p = ?1 ?x = e (cos y ? i sin y ) = e ? x?iy f ′( z ) = e ? z f ( z ) = ?e ? z + C
例5 求证 u( x , y ) = ( x ? y )( x + 4 xy + y ) 调和函数 函数， 为调和函数， 求出解析函数 f ( z ) = u + iv 2 3 3 2 证明 u( x , y ) = x +3x y ?3xy ? y u x = 3x 2 +6xy ?3 y 2 u xx = 6x +6 y 2 2 u yy = ?6x ?6 y u y = 3x ?6xy ?3 y 调和函数 函数， u xx + u yy = 0, u( x , y ) 为调和函数， 2 2 2 2 ? iv y f ′( z ) = u x + iux= 3x +6 xy ?3 y ? i 3 x + i 6 xy + i 3 y
= 3 (x ? y +2 xyi ) ?3i (x ? y +2 xyi ) 2 2 = 3 ( x + iy) ?3i ( x + iy ) 2 = 3(1 ? i ) z 3 6 f ( z ) = (1 ? i ) z +C
例6 设 v ( x , y ) = 解
求出解析函数 f ( z ) = u + iv ,
为调和函数， 调和函数， 函数
x ?y vy = 2 2
( x + y )2
f ′( z ) = u x + iv x = v y + iv x
1 x ? y ? i 2 xy ( x ? yi ) = 2 = = 2 2 2 2 (x + y ) ( z?z ) z 1 1 f (2)= 0, C= f ( z )= ? + C 2 z 1 1 f (z) = ? 2 z
例7 设 u( x , y ) = 2( x ? 1) y 为调和函数， 调和函数 函数， 求出解析函数 f ( z )= u + iv , f (2) = ? i 解 u( x , y ) = 2xy ?2 y ux = 2 y uy = 2x ? 2 f ′( z ) = u x + iv x
= ux ? iu y
= 2 y ? i 2 x +2i = ?2i ( x + yi ) +2i
f ( z ) = ? i ( z ? 1) +C
f ′( z ) = ?2i ( z ? 1)
f (2) = ? i +C = ? i f ( z ) = ? i ( z ? 1)
y 调和函数 函数， 为调和函数，x & 0 例8 设 v ( x , y )= arctan x
求出解析函数 解
f ( z ) = u + iv
x 1 1 ?v ? = 2 = 2 y x ?y x + y2 1+ 2 x y 1 y ?v ?( ? 2 )= ? 2 = y2 x + y2 x ?x 1+ 2
1 z = f ′( z ) = ux + iv x = v y + iv x = = x2 + y2 z ? z z f ( z ) = ln z + C
相关文章：
上一篇文章：下一篇文章：
最新添加资讯
24小时热门资讯
Copyright © 2004- All Rights Reserved. 学网 版权所有
京ICP备号-1 京公网安备02号对地球表面地形地貌离散表达的数字高程模型的构建方法
专利名称对地球表面地形地貌离散表达的数字高程模型的构建方法
技术领域本发明属于地球观测与导航技术领域，涉及地理信息系统曲面生成方法，尤其是 数字高程模型的构建方法。
背景技术数字高程模型(Digital Elevation Model，DEM)是对地球表面地形地貌的离散表 达，是GIS进行各种空间分析的数据基础，广泛应用于地貌定量分析、流域分析、土壤侵蚀 分析、遥感图像辅助分类、道路设计、军事演练等领域，DEM作为数字地形模拟的重要成果已 经成为国家空间数据基础设施的基本内容之一，并被纳入数字化空间数据框架进行规模化生产。自1958年美国麻生理工学院Miller教授为公路设计首次提出DEM以来，DEM已 经取得了长足的进展。由于DEM的广泛应用，对DEM精度的研究也受到有关专家的高度重 视。在近30年的发展中，产生了许多关于DEM精度评价指标，如中误差、标准差、误差直方 图、空间相关系数等，并提出了许多改进DEM精度的方法，这主要表现在以下几个方面(1)原始数据精度的评价分析由于DEM是对地球表面地形地貌的一种离散的近似数字表达，因此不可避免的含 有地形表达误差、人为采样误差等。原始数据误差决定着DEM精度，因此在进行各种插值之 前，应分析采样点误差，剔除误差较大点，充分挖掘原始采样数据的有用信息，提高插值精 度。如等高线数据用于DEM构建时，如果仅考虑等高线提供的离散数据点，忽略对应区域山 脊点(线)、山谷点(线)等地形特征点(线)这些隐藏的信息，会导致构建的DEM与原始 地形不符。(2)各种内插算法精度的对比评价各种内插算法都是为不同的应用目的而提出的，因此在进行插值前应根据采样点 的空间分布详细对比分析各种算法的使用条件、最优参数、插值精度等，确定最佳插值方法 或者对原始插值算法改进，提高插值精度。如反距离权重插值(IDW)没有考虑插值点周围 采样点的空间分布、密度，而将模型参数设成定值导致插值结果出现“牛眼”等异常现象。(3)专业数据的综合利用DEM—直是水文分析的重要数据来源，高精度DEM可以提高各种水文分析的精 度。如 Hutchinson 以粗糙度罚函数(roughness penalty)为基本方程[1]Hutchinson M F. Optimising the degree of data smoothing for locally adaptive finite element bivariate smoothing splines. Australian & New Zealand Industrial and Applied Mathematics Journal，4_796，弓| 入了 jftM胃fe (drainage enforcement algorithm)生成的DEM极大提高了流域面积、水系位置、水系长度等计算结果精度。由于 充分利用了流域水系矢量数据，并对计算数据进行了专业处理，因此Hutchinson研发的 ANUDEM已成为世界上流行的建立水文、地貌关系正确的DEM专业化软件。但需要注意的是， 尽管对DEM的改进提高了水文模型的精度，但是改进的DEM同时改变了原始DEM的高程和对应地形参数。现有的空间插值模型获取内插点数值时往往忽略采样点的空间相关性，使得已知 点之间的关系遭到破坏，最终导致较大的误差甚至出现各种错误。在以往的研究中，Evans 虽然认识到坡面、坡向和曲率是反映局部地形规律的重要变量[2]EvanS I S. An Integrated System of Terrain Analysis & Slope Mapping. Dept. of Geography, University of Durham, 1979，但他并没有将它们用于表达数字地面模型。事实上根据微分几何学原理，坡 面、坡向和曲率只是地面剖面线的决定要素，数字地面模型的决定要素是曲面的第一类基 本量和第二类基本量。
为了克服现有技术的不足，且从本质上解决曲面建模的误差问题，本发明的目的 是提供一种能对地球表面地形地貌离散表达的数字高程模型的构建方法。为达成所述目的，本发明提供对地球表面地形地貌离散表达的数字高程模型的构 建方法，该方法的实现步骤包括步骤Sl 利用数据采集设备对地球表面原始数据进行高程采样，获得地球表面多 个采样点的高程值；对多个采样点的高程值进行空间离散化处理，得到多个离散化的地形 曲面网格点；将每个采样点的高程值赋给离其最近的离散化的地形曲面网格点，生成多个 已知采样点高程值的地形曲面网格点；步骤S2 利用插值方法，根据已知采样点高程值的地形曲面网格点对未采样点的 地形曲面网格点进行空间预插值，求算出未采样点的地形曲面网格点高程值；步骤S3 根据所有地形曲面网格点高程值计算每个地形曲面网格点的第一类基 本量E，F，G和第二类基本量L，M, N，其中第一类基本量表示地形曲面上曲线的弧长、地形 曲面上两个方向的夹角和地形曲面域的面积，第二类基本量刻画了地形曲面空间中的弯曲 性；用高斯方程作为基本方程和高程采样点作为约束方程建立地形曲面方程组SUn+1 = Tn, 其中，S为地形曲面方程组的系数矩阵，U为待求的数字高程模型数据，T为计算出的地形曲 面常数项，η为地形曲面方程组迭代次数；步骤S4 利用高斯-赛德尔迭代算法解算地形曲面方程组SUn+1 = Tn，整个解算过 程为迭代计算，每次迭代结束得出数字高程模型的模拟结果；步骤S5 判断数字高程模型模拟结果是否满足精度要求，如果满足要求，执行步 骤6 ；如果不满足要求，执行步骤3 ；步骤S6 停止迭代计算，输出一组有序数值阵列形式表示地面高程实体的数字高 程模型，该数字高程模型用于对地球表面地形地貌进行离散表达。优选地，设所述地形曲面ζ的表达式为z = u(x，y)，其中x，y分别表示地形曲面 点的横、纵坐标；函数u表示地形曲面点横、纵坐标(x，y)和地形曲面高程点ζ的映射关系。优选地，所述第一类基本量E，F，G和第二类基本量L，M，N的表达式为 其中叫，
分别表示函数u对
地形曲面点的横、纵坐标x，y求一阶偏导数；uxy，Uxx，Uyy分别表示函数u对地形曲面点的横、纵坐标x，y求二阶偏导数。优选地，所述地形曲面方程组SUn+1 = Tn为大型稀疏线性方程组，且该方程组系数 矩阵S每行最多有九个非零元素。优选地，所述解算地形曲面方程组所需计算的时间与地形曲面网格点总数呈线性关系。为达成所述目的，本发明的第二方面，是提供一种地球表面地形地貌离散表达的 数字高程模型用作水文模型、地貌模型的数据基础。本发明的有益效果本发明以GS解算地形曲面方程组(该过程简称为HASM-GS) 的计算时间与地形曲面网格点总数呈线性相关，比采用直接法计算时间降低两个数量级。 本发明生成的数字高程模型精度要比传统的方法(SPLINE，IDW和KRIGING)高多倍。本发 明可作为水文模型、地貌模型的数据基础。选择高斯合成曲面作为研究对象如图2示出，将HASM-GS与MATLAB提供的经典方 程组解算方法比较表明(表1)，各种方法在达到相同的精度2. 65e-006时，HASM-GS耗时远 远小于MATLAB提供的各种函数。除去HASM-GS，耗时最少的为二次共轭梯度法(CGS)，该法 耗时约为HASM-GS的12. 409倍；耗时最多的为最小二乘QR分解法(LSQR)，该法耗时约为 HASM-GS的470. 912倍。根据不同的模拟区域验证HASM-GS的计算时间表明如表2和图3 示出，HASM-GS计算时间与地形曲面网格点总数成线性相关，HASM-GS极大的提高了数字高 程模型的生成速度。由表2和图3可见，HASM-GS计算时间与网格数呈非常好的线性相关 关系，其回归方程为t = 2. 036Xl(T5gn+0. 010 (R2 = 1) (1)式(1)中，t为计算时间，gn为网格数目。表1均勻采样条件下各种方法计算效率比较
误差最小各方法耗时与方法中文名称RMSE耗时(S)迭代步数HASM-GS 比值HASM-GSGauss-Seidel 法2.65e-11BICG双共轭梯度法2.65e-219. 936 表2HASM-GS计算时间与模拟区域网格数比较
图1本发明流程图；图2高斯合成曲面；图3HASM-GS计算时间和地形曲面网格点总数的回归曲线；图4 (a)-(d)各种方法的检核点离散图及对应的相关系数；
具体实施例方式下面结合附图详细说明本发明技术方案中所涉及的各个细节问题。本发明以曲面论为理论基础建立了对地球表面地形地貌离散表达的高精度数字 高程模型(HASM)。经典的曲面论是本发明的出发点，也是HASM区别于以往数字高程模型 生成方法的不同之处。基于曲面论基本定律建立的数字高程模型理论更加成熟，使其精度 远远高于其它方法。本发明是以高斯(Gauss)方程作为基本方程，以高程采样点作为约束 条件，将数字高程模型模拟转换为解算地形曲面方程组SUn+1 = Tn，该地形曲面方程组为大 型稀疏线性方程组，且该方程组系数矩阵S每行最多有九个非零元素。求解该大型稀疏线 性方程组即可得到高精度的数字高程模型。请参见图1本发明方法的流程图，具体实施如 下利用数据采集设备(如全站仪、水准仪、全球定位系统(GPS)等)对地球表面原始 数据进行高程采样，采集地球表面多个采样点的高程值，其中，高程数据源主要有地形图、航天航空遥感图像、雷达、野外实测等；将多个采样点的高程值进行空间离散化处理，得到 离散化的地形曲面网格点；将每个采样点的高程值赋给离其最近的地形曲面网格点，这些 网格点即为已知采样点高程值的地形曲面网格点；利用插值算法(如线性插值、三次样条 插值等)借助已知知采样点高程值的地形曲面网格点对未知采样点的地形曲面网格点空 间预插值，求算未采样点的地形曲面网格点的高程值。利用所有地形曲面网格点的高程值 计算所得地形曲面网格点的第一类基本量E，F，G和第二类基本量L，M, N，其中第一类基本 量可以表示地形曲面上曲线的弧长、地形曲面上两个方向的夹角和地形曲面(域)的面积， 第二类基本量刻画了地形曲面空间中的弯曲性。以高斯(Gauss)方程作为基本方程(请 详见公式2)，以高程采样点作为约束方程，建立对地球表面地形地貌离散表达的高精度数 字高程模型地形曲面方程组
，该方程组为大型稀疏线性方程组，且方程组系数 矩阵S每行最多有九个非零元素，U为待求数字高程模型数据，T为计算出的地形曲面常数 项，η为方程组迭代次数；解算方程组称为内迭代，更新T称为外迭代。利用高斯-赛德尔 (Gauss-Seidel)解算地形曲面方程组，SUn+1 = Τη，整个解算过程为迭代计算，每次迭代结束 即可得出高精度数字高程模型的模拟结果。设地形曲面表达式为
其中x，y分别表示地形曲面点的横、纵坐标；函 数u表示地形曲面点横、纵坐标(x，y)和地形曲面高程点ζ的映射关系，则Gauss方程的有
限差分形式为 式中，Uiij表示地形曲面第i行j列点高程值；i、j为自然数，
为地 形曲面的第二类克里斯托费尔符号，Xi, Yj表示地形曲面点的平面坐标；Ω\3Ω表示所有地 形曲面网格点的集合；h表示网格分辨率屯
」、队,」表示地形曲面网格点的第 一类基本量和第二类基本量，它们的有限差分形式为 式中，化,,表示地形曲面网格点的高程值。地形曲面的第二类克里斯托费尔符号
高斯(Gauss)方程的矩阵表达式为 式中，
，L，C分别表示
矩阵的行列数；Α，B分别表示方程(2)左端系数矩阵；Dn，EEn分别表示方程(2)右端向量。在高程采样点处的模拟值等于已知高程采样值，因此高精度数字高程模型(HASM) 数值模拟转换为等式约束的最小二乘问题，
,, 式中，P
&P，Qn均为矩阵集合；s. t.表示等式约束条件；M表示
高程采样点矩阵，H表示采样点高程值矩阵；当高程采样点在(i，j)时， 对于充分大的采样点权重因子λ，数字高程模型求算转换为解算无约束的最小二 乘问题，即HASM表达式为 也就是求解 式中，
其中解算方程组(15)称为内
迭代，更新右端向量Tn称为外迭代。利用Gauss-Seidel (GS)分量形式解算方程组(15)可表示为for k = 0，1,2,...Y1 = Ukfor i = 1,2，…，ηYi+1 = Yi+α
end for iuk+1 = yn+1stopping criteriaend for k式中uk表示数字高程模型第k次模拟值，k为0时，Uk用Utl来表示数字高程模型 初始模拟值；其余各个符号的表达式为两=(1，0，0，…，0)，e2 = (0，1，0，···，())，…，en =(0,0,0, -,1) Wi = -PiA^Pi= (S^yi)-Tei表示单位向量；CIi表示数字高程模型 每次迭代更新项；Pi为临时变量；Sii表示HASM系数矩阵第i行，i列元素；Si表示HASM的 系数矩阵第i行向量Ji表示计算出的地形曲面常数项。公式(16)本发明称为高精度数字 高程模型_高斯赛德尔迭代(HASM-GS)。利用计算机和MATLAB语言，根据曲面论的思想，先计算地形曲面上每个网格点的 第一类基本量E，F，G和第二类基本量L，M, N,然后根据曲面的差分形式建立曲面方程，利用 GS迭代法即(16)求解该方程即可得到高精度数字高程模型。实例选自山东省济南市某测区，该测区面积为3. 7km2，测区平均高程为180. 5m，高 程标准差为60. 5m。为了提高高程采样点的精度，首先在该测区均勻布设GPS控制点321 个，然后利用全站仪采集9756个高程点，高程采样点高程精度为5cm，平面精度为1cm。高 程采样点的平均密度为2637p0intS/km2，对应的高程采样点的间隔(ES)为19.5m。随机选 择9260个点用于模拟DEM，其余496个点用于精度检验。除了 HASM-GS外，传统的插值方法 (SPLINE, KRIGING和IDW)也用于构建DEM，并将其精度与HASM-GS进行比较。精度指标除 了采用中误差(RMSE)外，平均绝对误差(MAE)也用于精度评价，平均绝对误差(MAE)的计 算公式为 式中，EVi表示高程观测值；SVi表示数字高程模型模拟值；N为检核的高程点数目。 各种方法的模拟结果如表3所示。表3各种方法的插值精度比较(比值经典方法模拟误差与HASM-GS比值) 由表3可见，无论是MAE还是RMSE作为精度指标，HASM-GS插值精度要远高于传统 的方法，SPLINE, KRIGING, IDff 的 MAE 分别是 HASM-GS 的 9. 90 和 8. 3383 倍；RMSE 分别是 HASM-GS 的 7. 28 和 8. 6510 倍。为了更准确的描述每个检核点的模拟值与观测值之间的差异，本发明提供了检核 点的离散图请见图4(a)-图4(d)，其中，横轴表示观测值(EV)，纵轴表示模拟值(SV)，因此 当观测值等于模拟值时，检核点落在y = χ这条直线上。由图可见，除了 HASM-GS外，每种
9方法均存在模拟值与观测值相差较大的检核点。利用SPSS软件计算各种方法的模拟值和 观测值的拟合曲线相关系数得，HASM-GS的相关系数最大，表明精度最高(图4a)；而IDW相 关系数最小，表明精度最低(图4b)。图4(a)HASM-GS ；相关系数0. 9999943 ；图4(b) IDW ； 相关系数0. 9995738 ；图 4(c)KRIGING ；相关系数0. 9997633 ；图 4(d) SPLINE ；相关系数 0.9997172。 应指出的是，所描述的实施例仅旨在便于对本发明的理解，而对其不起任何限定 作用。但本发明的保护范围并不局限于此，任何熟悉该技术的人在本发明所揭露的技术范 围内，可理解想到的变换或替换，都应属于本发明的范围。
一种对地球表面地形地貌离散表达的数字高程模型的构建方法，其特征在于，该方法的实现步骤如下步骤S1利用数据采集设备对地球表面原始数据进行高程采样，获得地球表面多个采样点的高程值；对多个采样点的高程值进行空间离散化处理，得到多个离散化的地形曲面网格点；将每个采样点的高程值赋给离其最近的离散化的地形曲面网格点，生成多个已知采样点高程值的地形曲面网格点；步骤S2利用插值方法，根据已知采样点高程值的地形曲面网格点对未采样点的地形曲面网格点进行空间预插值，求算出未采样点的地形曲面网格点高程值；步骤S3根据所有地形曲面网格点高程值计算每个地形曲面网格点的第一类基本量E，F，G和第二类基本量L，M，N，其中第一类基本量表示地形曲面上曲线的弧长、地形曲面上两个方向的夹角和地形曲面域的面积，第二类基本量刻画了地形曲面空间中的弯曲性；用高斯方程作为基本方程和高程采样点作为约束方程建立地形曲面方程组SUn+1＝Tn，其中，S为地形曲面方程组的系数矩阵，U为待求的数字高程模型数据，T为计算出的地形曲面常数项，n为地形曲面方程组迭代次数；步骤S4利用高斯 赛德尔迭代算法解算地形曲面方程组SUn+1＝Tn，整个解算过程为迭代计算，每次迭代结束得出数字高程模型的模拟结果；步骤S5判断数字高程模型模拟结果是否满足精度要求，如果满足要求，执行步骤6；如果不满足要求，执行步骤3；步骤S6停止迭代计算，输出一组有序数值阵列形式表示地面高程实体的数字高程模型，该数字高程模型用于对地球表面地形地貌进行离散表达。
2.根据权利要求1所述的数字高程模型的构建方法，其特征在于，设所述地形曲面ζ的 表达式为z = u(χ, y)，其中x，y分别表示地形曲面点的横、纵坐标；函数u表示地形曲面 点横、纵坐标(X，y)和地形曲面高程点ζ的映射关系。
3.根据权利要求1所述的数字高程模型的构建方法，其特征在于，所述第一类基本量 E，F，G和第二类基本量L，M，N的表达式为
其中 、分别表示对地形曲面点的横、纵坐标X，y求一阶偏导数；Uxy，Uxx, Uyy分别表示函数U对地形曲面点的横、纵 坐标x，y求二阶偏导数。
4.根据权利要求1所述的数字高程模型的构建方法，其特征在于，所述地形曲面方程 组SUn+i = τη为大型稀疏线性方程组，且该方程组系数矩阵S每行最多有九个非零元素。
5.根据权利要求1所述的数字高程模型的构建方法，其特征在于，所述解算地形曲面 方程组所需计算时间与地形曲面网格点总数呈线性关系。
6.根据权利要求1所述的对地球表面地形地貌离散表达的数字高程模型用作水文模 型、地貌模型的数据基础。
本发明为对地球表面地形地貌离散表达的数字高程模型的构建方法，首先对地球表面原始数据提取采样点的高程值及空间离散化处理，生成多个已知采样点高程值的地形曲面网格点；求算出未采样点的地形曲面网格点高程值；根据所有地形曲面网格点高程值计算每个地形曲面网格点的第一和第二类基本量，用高斯方程和高程采样点建立地形曲面方程组；利用高斯-赛德尔迭代算法解算地形曲面方程组，每次迭代结束得出数字高程模型的模拟结果；判断数字高程模型模拟结果是否满足精度要求，如果满足要求，停止迭代计算，输出一组有序数值阵列形式表示地面高程实体的数字高程模型用于对地球表面地形地貌进行离散表达；如果不满足要求，重复迭代计算。
文档编号G06T17/05GKSQ
公开日日 申请日期日 优先权日日
发明者岳天祥, 杜正平, 陈传法 申请人:中国科学院地理科学与资源研究所下载作业帮安装包
扫二维码下载作业帮
1.75亿学生的选择
请问Ln(x^2+y^2) 的Uxx和Uyy分别是啥
Uxx是U对x求偏导后再对x求偏导,也就是二阶偏导Uyy是U对y求二阶偏导
我是想知道导出来的是什么。。。
Ux=2x/(x^2+y^2)
Uxx=[2(x^2+y^2)-2x*2x]/(x^2+y^2)^2
Uy=2y/(x^2+y^2)
Uyy=[2(x^2+y^2)-2y*2y]/(x^2+y^2)^2
为您推荐：
其他类似问题
U = ln(x²+y²) ,
Uxx和Uyy 分别是 U 对x和y的二阶偏导数Ux = (2x) /(x²+y²)Uxx = 2[ (x²+y²) - x * 2x ] / (x²+y²)² = 2(y²- x²) / (x²+y²)²...
Uxx = 1/2xuyy = 1/2y
扫描下载二维码}