統計分布
仲上 分布 または 仲上 m 分布は、 ガンマ分布 に関連する 確率分布 です 。仲上分布のファミリーには、 形状パラメータ と スケールパラメータの 2 つのパラメータがあります。これは、医療用超音波画像診断、通信工学、気象学、水文学、マルチメディア、地震学などで見られる物理現象をモデル化するために使用されます。
m
≥
1
/
2
{\displaystyle m\geq 1/2}
Ω
>
0
{\displaystyle \Omega >0}
特徴づけ
その 確率密度関数 (pdf)は [1]
f
(
x
;
m
,
Ω
)
=
2
m
m
Γ
(
m
)
Ω
m
x
2
m
−
1
exp
(
−
m
Ω
x
2
)
for
x
≥
0.
{\displaystyle f(x;\,m,\Omega )={\frac {2m^{m}}{\Gamma (m)\Omega ^{m}}}x^{2m-1}\exp \left(-{\frac {m}{\Omega }}x^{2}\right){\text{ for }}x\geq 0.}
ここで 、および 。
m
≥
1
/
2
{\displaystyle m\geq 1/2}
Ω
>
0
{\displaystyle \Omega >0}
その 累積分布関数 (CDF)は [1]
F
(
x
;
m
,
Ω
)
=
γ
(
m
,
m
Ω
x
2
)
Γ
(
m
)
=
P
(
m
,
m
Ω
x
2
)
{\displaystyle F(x;\,m,\Omega )={\frac {\gamma \left(m,{\frac {m}{\Omega }}x^{2}\right)}{\Gamma (m)}}=P\left(m,{\frac {m}{\Omega }}x^{2}\right)}
ここで Pは 正規化された(下側の) 不完全ガンマ関数 です。
パラメータ化
パラメータ とは次 の通りである [2]
m
{\displaystyle m}
Ω
{\displaystyle \Omega }
m
=
(
E
[
X
2
]
)
2
Var
[
X
2
]
,
{\displaystyle m={\frac {\left(\operatorname {E} [X^{2}]\right)^{2}}{\operatorname {Var} [X^{2}]}},}
そして
Ω
=
E
[
X
2
]
.
{\displaystyle \Omega =\operatorname {E} [X^{2}].}
この分布の中央値 には閉じた形式の解は存在しませんが、 m = 1など の特殊なケースは存在します 。実用上、中央値は観測値の 50 パーセンタイルとして計算する必要があります。
Ω
ln
(
2
)
{\displaystyle {\sqrt {\Omega \ln(2)}}}
パラメータ推定
分布を近似する別の方法は、 σ = Ω/ m として再パラメータ化することです 。 [3]
Ω
{\displaystyle \Omega }
仲上分布からの 独立した 観測値 を考えると、 尤度関数 は
X
1
=
x
1
,
…
,
X
n
=
x
n
{\textstyle X_{1}=x_{1},\ldots ,X_{n}=x_{n}}
L
(
σ
,
m
)
=
(
2
Γ
(
m
)
σ
m
)
n
(
∏
i
=
1
n
x
i
)
2
m
−
1
exp
(
−
∑
i
=
1
n
x
i
2
σ
)
.
{\displaystyle L(\sigma ,m)=\left({\frac {2}{\Gamma (m)\sigma ^{m}}}\right)^{n}\left(\prod _{i=1}^{n}x_{i}\right)^{2m-1}\exp \left(-{\frac {\sum _{i=1}^{n}x_{i}^{2}}{\sigma }}\right).}
その対数は
ℓ
(
σ
,
m
)
=
log
L
(
σ
,
m
)
=
−
n
log
Γ
(
m
)
−
n
m
log
σ
+
(
2
m
−
1
)
∑
i
=
1
n
log
x
i
−
∑
i
=
1
n
x
i
2
σ
.
{\displaystyle \ell (\sigma ,m)=\log L(\sigma ,m)=-n\log \Gamma (m)-nm\log \sigma +(2m-1)\sum _{i=1}^{n}\log x_{i}-{\frac {\sum _{i=1}^{n}x_{i}^{2}}{\sigma }}.}
したがって
∂
ℓ
∂
σ
=
−
n
m
σ
+
∑
i
=
1
n
x
i
2
σ
2
and
∂
ℓ
∂
m
=
−
n
Γ
′
(
m
)
Γ
(
m
)
−
n
log
σ
+
2
∑
i
=
1
n
log
x
i
.
{\displaystyle {\begin{aligned}{\frac {\partial \ell }{\partial \sigma }}={\frac {-nm\sigma +\sum _{i=1}^{n}x_{i}^{2}}{\sigma ^{2}}}\quad {\text{and}}\quad {\frac {\partial \ell }{\partial m}}=-n{\frac {\Gamma '(m)}{\Gamma (m)}}-n\log \sigma +2\sum _{i=1}^{n}\log x_{i}.\end{aligned}}}
これらの派生語が消えるのは、
σ
=
∑
i
=
1
n
x
i
2
n
m
{\displaystyle \sigma ={\frac {\sum _{i=1}^{n}x_{i}^{2}}{nm}}}
そして、 m に関する導関数がゼロになる m の値は 、ニュートン・ラプソン法 などの数値計算によって求められます 。
臨界点では全体的な最大値が達成されることがわかるので、臨界点は ( m , σ ) の最大尤度推定値となります。最大尤度推定値の 等価性 により 、Ω の最大尤度推定値も得られます。
ランダム変量生成
仲上分布は ガンマ分布 と関連しています。特に、ランダム変数 が与えられた場合、 、 を設定し、 の平方根を取ることで 、 ランダム 変数 を得ることができます 。
Y
∼
Gamma
(
k
,
θ
)
{\displaystyle Y\,\sim {\textrm {Gamma}}(k,\theta )}
X
∼
Nakagami
(
m
,
Ω
)
{\displaystyle X\,\sim {\textrm {Nakagami}}(m,\Omega )}
k
=
m
{\displaystyle k=m}
θ
=
Ω
/
m
{\displaystyle \theta =\Omega /m}
Y
{\displaystyle Y}
X
=
Y
.
{\displaystyle X={\sqrt {Y}}.\,}
あるいは、パラメータを に 設定して カイ分布 から仲上分布 を生成し 、その後にランダム変数のスケーリング変換を行うこともできます。つまり、仲上ランダム変数は、以下のように カイ分布ランダム変数の単純なスケーリング変換によって生成されます 。
f
(
y
;
m
,
Ω
)
{\displaystyle f(y;\,m,\Omega )}
k
{\displaystyle k}
2
m
{\displaystyle 2m}
X
{\displaystyle X}
Y
∼
χ
(
2
m
)
{\displaystyle Y\sim \chi (2m)}
X
=
(
Ω
/
2
m
)
Y
.
{\displaystyle X={\sqrt {(\Omega /2m)Y}}.}
カイ分布の場合、自由度は 整数でなければなりませんが、仲上分布の場合、自由度は 1/2 より大きい任意の実数にすることができます。これが重要な違いであり、したがって、仲上 m 分布はカイ分布の一般化と見なされ、ガンマ分布がカイ 2 乗分布の一般化と見なされるのと同様です。
2
m
{\displaystyle 2m}
m
{\displaystyle m}
歴史と応用
仲上分布は比較的新しいもので、1960年に仲上実によって長距離高周波電波伝搬における小規模フェージングの数学的モデルとして初めて提案されました。 [4]これは 、複数の経路を通過する 無線 信号 の減衰をモデル化するために使用され [5] 、フェージングチャネルが無線通信に 与える 影響を研究するために使用されてきました 。 [6]
mを 単位区間( q = m ; 0 < q < 1)に 制限すると [ 疑わしい – 議論する ]仲上q 分布 が 定義 され 、これは ホイト 分布は 、1940年代にRSホイトによって初めて研究されました。 [7] [8] [9] 特に、 極座標 で書き直された 2変量正規 の真の平均の周りの 半径は 、ホイト分布に従います。同様に、 複素正規 の 係数 も同様です。
2 m = k の場合、仲上分布はスケールされた カイ2乗分布 を与えます。
の場合 、仲上分布はスケールされた 半正規分布 を与えます。
m
=
1
2
{\displaystyle m={\tfrac {1}{2}}}
仲上分布は、 p = 2 および d = 2 m である 一般化ガンマ分布 の特別な形式です 。
参照
数学ポータル
参考文献
^ ab Laurenson, Dave (1994). 「Nakagami Distribution」. レイトレーシング技術による屋内無線チャネル伝搬モデリング 。2007-08-04 に 取得 。
^ R. Kolar、R. Jirik、J. Jan (2004)「Nakagami-mパラメータの推定値の比較と心エコー検査への応用」、 Radioengineering 、13 (1)、8–12
^ Mitra, Rangeet; Mishra, Amit Kumar; Choubisa, Tarun (2012). 「Nakagami-m分布のパラメータの最大尤度推定」。 国際通信・デバイス・インテリジェントシステム会議 (CODIS)、2012年 :9~12ページ。
^ 中上正之 (1960)「m分布、急速なフェーディングの強度の一般的な公式」ウィリアム・C・ホフマン編『 電波伝搬の統計的手法:1958年6月18~20日に開催されたシンポジウムの議事録』 3~36ページ。ペルガモン・プレス、 doi :10.1016/B978-0-08-009306-2.50005-4
^ パーソンズ、JD(1992) モバイル無線伝搬チャネル 。ニューヨーク:ワイリー。
^ Ramon Sanchez-Iborra、Maria-Dolores Cano、Joan Garcia-Haro (2013)。「フェーディング チャネルでの VoIP トラフィックの QoE のパフォーマンス評価」。2013 World Congress on Computer and Information Technology (WCCIT) 。pp. 1–6。doi : 10.1109 /WCCIT.2013.6618721。ISBN 978-1-4799-0462-4 .S2CID 16810288 。
^ Paris, JF (2009). 「Nakagami-q (Hoyt) 分布関数とその応用」. Electronics Letters . 45 (4): 210–211. Bibcode :2009ElL....45..210P. doi :10.1049/el:20093427.
^ 「HoytDistribution」.
^ 「NakagamiDistribution」.