最美情侣中文字幕电影,在线麻豆精品传媒,在线网站高清黄,久久黄色视频

歡迎光臨散文網(wǎng) 會員登陸 & 注冊

衛(wèi)星數(shù)據(jù)處理(一)——HT分析

2021-06-17 21:18 作者:Berton9407  | 我要投稿

此模塊涉及相對專業(yè)的數(shù)理知識,借鑒《analysis Methods for Multi-Spacecraft Data》(1998)一書內(nèi)容,講述deHoffmann-Teller(HT)分析、最小方差分析(minimum variance analysis,MVA)、奇異值分解法(singular value decomposition,SVD)求解波動性質(zhì)三模塊內(nèi)容。

(一)HT分析

通常,對于背景介質(zhì)均一、磁場和速度呈現(xiàn)相干的準靜態(tài)模式(如:阿爾芬波等),就可以找到一個坐標系使得等離子體流v沿磁場B徑向,此時電場消失,稱作HT坐標系,這種分析方法叫做HT分析(deHoffmann-Teller,1950)。其中,最重要的參量是速度矢量V_%7BHT%7D,使得各點的“剩余電場(residual electric field)”E%5E%5Cprime基本為0,根據(jù)法拉第定律(Farady’s law)有:

%5Cbigtriangledown%5Ctimes%20E%5E%7B%5Cprime%7D%3D-%5Cfrac%7B%5Cpartial%20B%5E%7B%5Cprime%7D%7D%7B%5Cpartial%20t%7D%3D0.

這說明在此坐標系下,磁場不隨時間而改變,即其位形是相對穩(wěn)定的。關(guān)于V_%7BHT%7D需要滿足的關(guān)系式也就可以表示為:

E%5E%7B%5Cprime%7D%3DE%2BV_%7BHT%7D%5Ctimes%20B%3D-v%5Ctimes%20B%2BV_%7BHT%7D%5Ctimes%20B%3D-(v-V_%7BHT%7D)%5Ctimes%20B.

這樣,就可以用V_%7BHT%7D來“描述”一整段速度和磁場的觀測數(shù)據(jù)。Aggson et al.(1983)使用復(fù)分析的方法求解V_%7BHT%7D,而現(xiàn)在常用的是基于非復(fù)分析和最小二乘法得到(Sonnerup et al., 1987; 1990)。給出一段時間內(nèi)有N個觀測樣本數(shù)的E%5E%5Cprime的方差和D(E%5E%5Cprime),表達式為:

D(E%5E%7B%5Cprime%7D)%3D%5Cfrac%20%7B1%7D%7BN%7D%5Csum%5Cnolimits_%7B1%7D%5EN%7CE_i%5E%5Cprime-E(E_i%5E%5Cprime)%7C%20%5E2%3D%5Cfrac%20%7B1%7D%7BN%7D%5Csum%5Cnolimits_%7B1%7D%5EN%7C-(v-V_%7BHT%7D)%5Ctimes%20B_i%7C%20%5E2.

不難看出,上式也是關(guān)于V_%7BHT%7D的函數(shù),可記為D(E%5E%7B%5Cprime%7D)%5Cequiv%20f(V_%7BHT%7D)。因此,要使得D(E%5E%7B%5Cprime%7D)達到最小值,即通過對每個分量的偏導(dǎo)數(shù)為0得到最優(yōu)解V_%7BHT%7D,也就是%5Cbigtriangledown_%7BV_%7BHT%7D%7Df(V_%7BHT%7D)%3D0。當然,也可以表示成:

%5Cfrac%7B%5Cpartial%20f%7D%7B%5Cpartial%20V_%7BHT%7D(x)%7D%3D%5Cfrac%7B%5Cpartial%20f%7D%7B%5Cpartial%20V_%7BHT%7D(y)%7D%3D%5Cfrac%7B%5Cpartial%20f%7D%7B%5Cpartial%20V_%7BHT%7D(z)%7D%3D0.

通過推導(dǎo),引入每個時刻點的矩陣K_i,其內(nèi)元素K_i%5E%7B%5Cmu%20%5Cnu%20%7D(%5Cmu%2C%5Cnu%3Dx%2Cy%2Cz%20)的值由下式得到(Khrabrov & Sonnerup, 1998):

K_i%5E%7B%5Cmu%20%5Cnu%20%7D%3D%7CB_i%7C%5E2(%5Cdelta%5E%7B%5Cmu%20%5Cnu%20%7D%20-%5Cfrac%7BB_i%5E%5Cmu%20B_i%5E%5Cnu%7D%7B%7CB_i%7C%5E2%7D)%5Cequiv%20%7CB_i%7C%5E2P_i%5E%7B%5Cmu%5Cnu%7D.

那么,V_%7BHT%7D就滿足下式:

%3CK_i%3EV_%7BHT%7D%5Cequiv%20%5Csum%5Cnolimits_%7B1%7D%5EN%20K_iV_%7BHT%7D%3D%5Csum%5Cnolimits_%7B1%7D%5EN%20K_iv_i%5Cequiv%20%3CK_iv_i%3E.

所以,容易得到:

V_%7BHT%7D%3D%3CK_i%3E%5E%7B-1%7D%3CK_iv_i%3E%3D%3CK_i%3E%5E%7B-1%7D%5Cfrac%7B1%7D%7BN%7D%5Csum%5Cnolimits_%7B1%7D%5ENB_i%5E2(v_%5Cperp)_i.%20

同時,給出對應(yīng)V_%7BHT%7D的誤差%5CDelta%20V_%7BHT%7D表達式為:

%5CDelta%20V_%7BHT%7D%3D%3CK_i%3E%5E%7B-1%7D%3CE_i%5E%5Cprime%5Ctimes%20B_i%3E.

從而,也得到了總體評估矩陣S%5E%7B%5Cmu%20%5Cnu%7D%3D%3C%3C%5CDelta%20V_%7BHT%7D%5E%5Cmu%5CDelta%20V_%7BHT%7D%5E%5Cnu%3E%3E,經(jīng)過推導(dǎo),其滿足下式:

S%5E%7B%5Cmu%5Cnu%7D%3D%5Cfrac%7BD(E%5E%5Cprime)%7D%7B2N-1%7D(%3CK_i%3E%5E%7B-1%7D)%5E%7B%5Cmu%5Cnu%7D.

則關(guān)于S%3CK_i%3E的特征矢量%5CLambda_m%3D%5Cfrac%7BD(E%5E%5Cprime)%7D%7B2N-1%7D%5Cfrac%7B1%7D%7B%5Clambda_m%7D.相同,特征值分別為%5CLambda_m(m%3D1%2C2%2C3)%5Clambda_m%5CLambda_m表示V_%7BHT%7De_m上的方差,有:

%5CLambda_m%3D%5Cfrac%7BD(E%5E%5Cprime)%7D%7B2N-1%7D%5Cfrac%7B1%7D%7B%5Clambda_m%7D.

同時,V_%7BHT%7D在單位矢量%5Chat%7Bn%7D%20的投影方差為%5Csigma_n%5E2%3D%5Chat%7Bn%7D%5E%5Cmu%20S%5E%7B%5Cmu%5Cnu%7D%5Chat%7Bn%7D%5E%5Cnu。如此,假如存在速度差%5Cdelta%20V%3DV-V_%7BHT%7D,則關(guān)于V的方差f(V)可由下式表示:

f(V)%3Df(V_%7BHT%7D)%2B%5Cdelta%20V%3CK_i%3E%5Cdelta%20V%5ET%3DD(E%5E%5Cprime)%2B%5Cdelta%20V%3CK_i%3E%5Cdelta%20V%5ET.

因此,在特征值所約束的特征向量空間中,可以得到:

f(V)%3DD(E%5E%5Cprime)%2B%5Csum%5Cnolimits_%7B1%7D%5EN%20%5Clambda_m%20(%5Cdelta%20V_m)%5E2%3B

%5Csum%5Cnolimits_%7B1%7D%5EN%20%5Cfrac%7B(%5Cdelta%20V_m)%5E2%7D%7B%5CLambda_m%7D%3D1%3B

%5Csum%5Cnolimits_%7B1%7D%5EN%20%5Clambda_m(%5Cdelta%20V_m)%5E2%3D%5Cfrac%7BD(E%5E%5Cprime)%7D%7B2N-1%7D.%20

D(E%5E%5Cprime)最小值類似的,引入等離子體越場速度v_%5Cperp%5E%5Cprime%3DE_%5Cperp%5E%5Cprime%5Ctimes%20B,使得其均方差Q(v_%5Cperp%5E%5Cprime)最小,表達式為:

Q%3D%5Cfrac%7B1%7D%7BN%7D%5Csum%5Cnolimits_%7B1%7D%5EN(v_%5Cperp%5E%5Cprime)%5E2%3D%5Cfrac%7B1%7D%7BN%7D%5Csum%5Cnolimits_%7B1%7D%5EN%5Cfrac%7B%7C-(v_i-V_%7BHT%7D)%5Ctimes%20B_i%7C%5E2%7D%7BB_i%5E2%7D.

不難看出,此時的V_%7BHT%7D就滿足(即用P矩陣代替K矩陣):%3CP_i%3EV_%7BHT%7D%3D%3CP_iv_i%3E,也就是說V_%7BHT%7D%3D%3CP_i%3E%5E%7B-1%7D%3CP_iv_i%3E

當然,對于MHD和混雜(hybrid)模擬結(jié)果都顯示在HT坐標系下顯現(xiàn)完美的Walén關(guān)系(Lin et al., 2009)。在實際數(shù)據(jù)處理中,為了避免引入HT坐標系和背景磁場帶來的不確定性,Li et al.(2016)展示了一種用多次帶通濾波的方式來檢驗Walén關(guān)系。近來,源自希爾伯特-黃氏變換(Hilbert-Hung transform,HHT)(Huang et al., 1998)方法,Liu et al.(2020)基于總體平均經(jīng)驗?zāi)J椒纸猓╡nsemble empirical mode decomposition,EEMD)的方法,認為V_%7BHT%7D隨時間而發(fā)生相應(yīng)改變。另外,Yang et al.(2019; 2020)基于觀測樣本,以均值來替代得到V_%7BHT%7D,即表達式為:

V_%7BHT%7D%3D%5Cfrac%7B1%7D%7BN%7D%5Csum%5Cnolimits_%7B1%7D%5EN(v_i%2B%7C%5Cfrac%7B%5CDelta%20v%5E%5Cmu%7D%7B%5CDelta%20b%5E%5Cmu%7D%7Cb_i).%20

但要注意此種情況下,不一定是最小二乘的結(jié)果,此式被稱為定速度下的波場(wave frame with constant velocity,WFCV)。如果考慮速度不是定常值,則可加入滑動平均(也可用濾波代替)而引出變速度下的波場(wave frame with varying velocity,WFFV)。

衛(wèi)星數(shù)據(jù)處理(一)——HT分析的評論 (共 條)

分享到微博請遵守國家法律
乐至县| 泸州市| 乌海市| 平度市| 德江县| 富川| 枣阳市| 酉阳| 宾阳县| 吉隆县| 林周县| 晋中市| 赤水市| 定远县| 中西区| 石渠县| 姜堰市| 内乡县| 定襄县| 晋城| 北京市| 保康县| 宁阳县| 伊金霍洛旗| 青冈县| 壤塘县| 渭源县| 调兵山市| 略阳县| 乐清市| 攀枝花市| 德惠市| 泰州市| 桦川县| 盘锦市| 潢川县| 施甸县| 赣榆县| 和静县| 沧源| 灵宝市|