没有合适的资源?快使用搜索试试~ 我知道了~
拟合算法-基于卫星高度计海面高度异常资料获取潮汐调和常数的方法及应用.pdf
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 86 浏览量
2022-05-02
06:49:32
上传
评论 1
收藏 6.05MB PDF 举报
温馨提示
试读
64页
拟合算法-基于卫星高度计海面高度异常资料获取潮汐调和常数的方法及应用.pdf
资源推荐
资源详情
资源评论
1
题目: 基于卫星高度计海面高度异常资料
获取潮汐调和常数的方法及应用
摘要:
本文以 TOPEX/POSEION 卫星高度计观测所得南海区域海面高度异常产品
为基础,利用调和分析方法得到了符合资料长度要求的 TOPEX/POSEION 星下
观测点处 M
2
、S
2
、K
1
和 O
1
四个主要分潮的调和常数;进而,在每个轨道上利用
多项式拟合方法分离了四个主要分潮的正压潮和内潮部分;最后,利用数据插值
法和车贝雪夫多项式拟合法两种方法对南海海域 M
2
、S
2
、K
1
和 O
1
四个主要分潮
的正压潮调和常数进行了二维插值或者拟合,得到了南海四个主要分潮的同潮图。
针对问题一,首先确定 TOPEX/POSEION 观测数据中 M
2
、S
2
、K
1
和 O
1
四个
主要分潮的混淆周期,然后根据 Rayleigh 准则,计算得到了分离四个主要分潮所
需的最短资料长度为 2.97 年,进而对题目中给出的所有 TOPEX/POSEION 星下
观测点进行了筛选,得到了符合要求的星下轨道点共 4309 个。利用调和分析方
法,对于筛选后的 4309 个星下轨道点处的观测数据进行了分析,得到了上述所
有星下轨道点处 M
2
、S
2
、K
1
和 O
1
四个主要分潮的调和常数(振幅 H 和迟角 G);
根据星下轨道点和验潮站之间的距离,选取了两者之间距离小于 15 千米的 40 个
验潮站对计算所得四个主要分潮的调和常数进行了检验:40 个验潮站处,M
2
、
S
2
、K
1
和 O
1
四个分潮计算值和观测值之间的矢量差的平均值分别为为 4.99 厘米,
3.43 厘米,4.42 厘米和 4.54 厘米,说明利用 TOPEX/POSEION 观测数据所提取
的调和常数结果较准。
针对第二问,首先基于 TOPEX/POSEION 卫星绕地球一圈的时间和轨道重
复周期,依据筛选后的 4309 个星下轨道点处数据的观测时间,对南海海域 20 条
轨道进行了分离和提取,确定了每条轨道上存在观测数据的星下轨道点的经纬度。
利用多项式拟合法对每一条轨道上 HcosG 和 HsinG 进行拟合,所得多项式拟合
部分为正压潮,而剩余部分为内潮。以第四问为启发,对于多项式拟合的次数进
行了遍历,以 4309 个星下点处调和常数拟合值和观测所得值之间矢量差的平均
值作为判别标准,确定了最优的多项式拟合次数为 8 次。最终,在所有 20 条轨
道上,利用 8 次多项式对四个主要分潮的 HcosG 和 HsinG 进行了拟合,然后反
算出四个分潮正压潮和内潮两个部分的调和常数,实现了对于正压潮和内潮的分
离。
针对第三问,以第二问所得 4309 个星下轨道点处正压潮的调和常数为基础,
以题目中所给水深数据的经纬度网格为目标区域,利用基于三角网格的 cubic 插
值法(方法 1)和二维车贝雪夫多项式方法(方法 2)分别对星下轨道点处的正
压潮调和常数进行了插值或拟合。在利用二维车贝雪夫多项式方法进行拟合时,
以第四问为启发,对于经度和纬度两个方向的一维车贝雪夫多项式的次数进行了
遍历,得到了最优的多项式次数:对于 M
2
和 S
2
分潮,两个方向的最优次数均分
别为 4 次和 4 次;对于 K
1
分潮,两个方向的最优次数分别为 4 次和 3 次;对于
O
1
分潮,两个方向的最优次数分别为 3 次和 3 次。基于三角网格的 cubic 插值法
(方法 1)所得 M
2
、S
2
、K
1
和 O
1
四个主要分潮的同潮图中,振幅插值结果较好,
2
但等迟角线不够平滑;二维车贝雪夫多项式方法(方法 2)所得 M
2
、S
2
、K
1
和 O
1
四个主要分潮的同潮图,四个分潮振幅和迟角的拟合结果均明显优于方法 1 所得
插值结果。需要指出的是,两种方法所得四个主要分潮的同潮图与实际同潮图还
有一定的差距,相关文献表明,利用伴随同化方法同化 TOPEX/POSEION 观测
数据所得调和常数,对南海四个主要分潮进行数值模拟与同化研究,将会进一步
提升所得同潮图的质量。
关键词: 调和分析、潮汐调和常数、卫星高度计、车贝雪夫多项式拟合
3
1. 问题重述
海洋潮汐是在天体引潮力作用下形成的长周期波动现象,在水平方向上表现
为潮流的涨落,在铅直方向上则表现为潮位的升降。潮汐潮流运动是海洋中的基
本运动之一,它是动力海洋学研究的重要组成部分,对它的研究直接影响着波浪、
风暴潮、环流、水团等其他海洋现象的研究,在大陆架浅海海洋中,对潮汐潮流
的研究更具重要性[1]。
区域海洋潮汐的数值模拟需要提供开边界的水位调和常数,而开边界的水位
调和常数,或者来源于观测、或者来源于全球海洋潮汐的数值模拟;而全球海洋
潮汐的数值模拟,相当耗费资源。虽然目前有国外学者或研究机构,能够提供区
域海洋潮汐的调和常数,但实质上的评价结果难以令人满意。
从各分潮的调和常数获取的发展史来说,通过对已有观测结果进行插值曾经
是首选,但发展过程中逐渐被数值模拟方法所取代。高度计资料的出现,引发部
分学者开展了插值方法的研究,本文以 TOPEX/POSEIDON 卫星高度计星与潮汐
相关的的资料为研究依据,通过已知各个经纬坐标点的南海海面高度异常数据,
计算潮汐调和常数,与潮汐验潮点的调和常数给予评价或检验,并给出评价结果
的分析或评价。
具体任务如下:
根据沿轨道的星下观测点的海面高度异常值,提取所有星下观测点各主要分
潮(M2、S2、K1、O1)的潮汐调和常数,注意能有效提取那些分潮的潮汐调和
常数取决于相应的资料长度;对提取的潮汐调和常数,应利用潮汐验潮点的调和
常数给予评价或检验,并给出评价结果的分析或评价。
得到所有星下观测点各主要分潮(M2、S2、K1、O1)的潮汐调和常数,沿
轨道作图后,可发现潮汐调和常数在沿轨道方向,在空间有细结构,而此细结构
是内潮对正压潮的调制;请设法对沿轨道的各分潮的潮汐调和常数进行正压潮和
内潮的分离。
设计数据插值或拟合方法给出南海的各主要分潮的同潮图,并利用潮汐验潮
点的调和常数给予评价或检验,并给出评价结果的分析或评价。
如果你们还有时间和兴趣,还可考虑下列:
如果在对沿轨道的潮汐调和常数分离、插值或拟合的过程中,利用了特定的
函数进行拟合,是否能够确定出需利用的特定函数的最佳(高)次数?上述结论
是否对第 3 问有启示或帮助。
更多数学建模资料请关注微店店铺“数学建模学习交流”
https://k.weidian.com/RHO6PSpA
4
2. 问题分析
根据本题主要任务可见,本题是按照分步设计的思路,首先完成对提供的
TOPEX/POSEIDON 卫星下高度计所测的南海海面高度异常数据的文件进行预处
理,在此基础上对所选点进行主要分潮(M2、S2、K1、O1)的潮汐调和常数分
析。对提取的潮汐调和常数,应利用潮汐验潮点的调和常数利用插值法给予评价
或检验,并探讨所求的潮汐调和常数与观测资料所得的实际值产生差异的原因;
然后根据后续方法探讨将正压潮与内压分离[2]。
2.1. 问题一的分析
问题一所给出的数据是 TOPEX/POSEION 卫星高度计所测的南海海面高度
异常数据文件,时间跨度为 1992 年到 2016 年。首先整理这 14 年的卫星高度计
资料得到南海海区各轨道观测点的海表面高度异常数据的时间序列,然后对此时
间序列进行调和分析得到各分潮调和常数。考虑到观测次数过少会降低调和分析
结果的可靠性,本文在计算调和分析的过程中,根据潮汐混淆周期和 Rayleigh 周
期进行分析、筛选、提取数据。根据对提取的各观测点数据,重新构建模型对潮
汐调和常数进行求解。最后利用潮汐验潮点的调和常数与所求的调和常数进行对
比分析,评价和检验,并探讨所求的潮汐调和常数与观测资料的差距及其产生差
距的原因。
2.2. 问题二的分析
通过构造该区域的轨道提取模型,将 TOPEX/POSEION 卫星高度计测量点
投影到该研究区域内,即可形成卫星运行的轨道图,并由此确定每条轨道上的观
测点的位置信息,根据所得的每条轨道的位置信息。得到所有星下观测点各主要
分潮(M2、S2、K1、O1)的潮汐调和常数,沿轨道作图后查找潮汐调和常数在
沿轨道方向的变化规律。继而利用多项式拟合的方法得到该轨道上任意点的潮汐
调和常数,最后利用同一位置处拟合所得的潮汐调和常数减去所测的潮汐调和常
数即可将正压潮与内压潮分离。
2.3. 问题三的分析
问题三要求得出南海的各主要分潮的同潮图,首先就需要求出南海区域内任
意点的各主要分潮的潮汐调和常数,我们初步分析有两种思路:一是根据问题一
所分析出来的 4309 个位置的潮汐调和常数进行差值处理得到整个平面研究区域
内的潮汐调和常数数值;二是根据问题二中多项式拟合所得到的所有轨道上任意
点处的潮汐调和常数进行差值处理得到整个平面研究区域内的潮汐调和常数数
值。
再依据验潮站观测数据所得的潮汐调和常数分别与两种插值方法所得的验
潮站处的潮汐调和常数最对比并分析其结果。
5
3. 问题一的求解
根据题意,问题一是让我们根据沿轨道的星下观测点的海面高度异常值,提
取所有星下观测点各主要分潮(
2
M
、
2
S
、
1
K
、
1
O
)的潮汐调和常数,同时要注
意能有效提取主要分潮(
2
M
、
2
S
、
1
K
、
1
O
)的分潮的潮汐调和常数与相应的资
料长度的关系,最后利用题目给出的潮汐验潮点的调和常数对我们所计算出来的
潮汐调和常数予评价或检验。
3.1. 调和分析的基本原理与编程实现
实际潮汐的分潮从其来源看可分为以下四种:天文分潮、气象分潮、天文
-
气
象分潮和浅水分潮。
从分潮的频率分布来看,分潮在频率上的分布是极不均匀的,而是分成族、
群和亚群。在
Doodson
展开中,按
Doodson
数
1
区分潮族,按
2
区分群,按
3
区分亚群。在潮族中一般分为长周期分潮族
)0(
1
=
、全日分潮族
)1(
1
=
、半日
分潮族
)2(
1
=
、三分日分潮族
)3(
1
=
直到十二分日分潮族
)12(
1
=
,共
13
个
潮族。在每一个潮族中,具有不同数量的群和亚群。
在亚群中的各个分潮的角速度是非常接近的,彼此之间只有微小的差异。因
此,在资料长度有限的情况下,亚群中的各个分潮是无法区分的。因此,在实际
的潮汐分析中,往往将一个亚群合成一个分潮,此时这一分潮的振幅和迟角不再
是常数,而是随着升交点的黄经十分缓慢地变化,一般在较短的时间内可近似看
作不变。这样的分潮实质上是准调和的,但习惯上仍叫做调和分潮。
实际水位可以看作是很多个调和分潮迭加的结果,但是在实际分析中只能选
取其中有限个较主要的分潮。假设我们选取了
J
个分潮,对于任一点的潮位表达
式为:
==
−+++=−++=
J
j
jjjjjj
J
j
jjjjj
guvthfSguvhfSh
1
00
1
0
)cos()cos(
,
其中,
0
S
为余水位,
j
f
为交点因子,
j
u
为交点订正角,
jj
gh ,
为分潮的调和
常数(振幅和迟角)。
3.1.1.
分潮角速度的计算
••••••
+
++
++= pNphs
654321
其中:
为分潮的角速度,
654321
,,,,,
为 Doodson 数,
剩余63页未读,继续阅读
资源评论
普通网友
- 粉丝: 12w+
- 资源: 9335
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功