[發(fā)明專利]一種海域流動(dòng)點(diǎn)外部擾動(dòng)引力水平分量的無奇異計(jì)算方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201810054343.4
• 申請(qǐng)日: 2018-01-19
• 公開（公告）號(hào): CN108416082B
• 公開（公告）日: 2021-09-10
• 發(fā)明（設(shè)計(jì)）人: 鄧凱亮;黃謨濤;黃辰虎;吳太旗;李凱峰;王川;陳欣
• 申請(qǐng)（專利權(quán)）人: 中國人民解放軍92859部隊(duì)
• 主分類號(hào): G06F30/20
• 分類號(hào): G06F30/20
• 代理公司: 天津盛理知識(shí)產(chǎn)權(quán)代理有限公司 12209
• 代理人: 王利文

1.一種海域流動(dòng)點(diǎn)外部擾動(dòng)引力水平分量的無奇異計(jì)算方法，其特征在于包括以下步驟：

步驟1、利用局域泰勒級(jí)數(shù)展開實(shí)現(xiàn)海域格網(wǎng)點(diǎn)外部擾動(dòng)引力水平分量的無奇異計(jì)算；

步驟2、聯(lián)合局域泰勒級(jí)數(shù)展開和內(nèi)插方法得到計(jì)算點(diǎn)擾動(dòng)引力分量的對(duì)應(yīng)部分；

所述步驟1通過聯(lián)合以下兩組公式計(jì)算實(shí)現(xiàn)：

式中，ψ為計(jì)算點(diǎn)和流動(dòng)點(diǎn)之間的球面角距；l為計(jì)算點(diǎn)和流動(dòng)點(diǎn)之間的空間距離；μ稱為廣義面密度；R為地球橢球平均半徑；h為計(jì)算高度；dσ為單位球積分面積元；ψ₀₀為計(jì)算點(diǎn)網(wǎng)格數(shù)據(jù)塊半徑；和δg_λ(ψ₀₀)分別表示導(dǎo)致奇異問題的中央?yún)^(qū)塊對(duì)海域流動(dòng)點(diǎn)外部擾動(dòng)引力水平分量的貢獻(xiàn)；x軸表示指向正北坐標(biāo)軸，y軸表示正向東坐標(biāo)軸；μ_x′代表μ對(duì)x軸的偏導(dǎo)數(shù)，s₀代表數(shù)據(jù)網(wǎng)格大小的一半；

所述步驟2的實(shí)現(xiàn)方法為：在計(jì)算點(diǎn)為非網(wǎng)格點(diǎn)時(shí)，設(shè)計(jì)算點(diǎn)距4個(gè)網(wǎng)格點(diǎn)的距離分別為：s₁、s₂、s₃和s₄，由局域泰勒級(jí)數(shù)展開計(jì)算得到4個(gè)網(wǎng)格點(diǎn)外部擾動(dòng)引力水平分量影響量，其中，子午方向上的網(wǎng)格點(diǎn)外部擾動(dòng)引力水平分量影響量分別為：和卯酉方向的網(wǎng)格點(diǎn)外部擾動(dòng)引力水平分量影響量分別為δg_λ1、δg_λ2、δg_λ3和δg_λ4，則由下式內(nèi)插得到計(jì)算點(diǎn)擾動(dòng)引力分量的對(duì)應(yīng)部分：