กลับไปหน้าบทความ

อ่าน 6 นาที

ฟังก์ชันความแปรปรวนร่วมของ Matérn

ความแปรปรวนร่วมและความสัมพันธ์/ธรณีสถิติ/การวิเคราะห์เชิงพื้นที่

ในทางสถิติค่าความแปรปรวนร่วมของ Matérnหรือที่เรียกว่าเคอร์เนล Matérn เป็นฟังก์ชันความแปรปรวนร่วมที่ใช้ในสถิติเชิงพื้นที่ ภูมิสถิติการเรียนรู้ของเครื่อง

ฟังก์ชันความแปรปรวนร่วมของ Matérn

ในทางสถิติค่าความแปรปรวนร่วมของ Matérnหรือที่เรียกว่าเคอร์เนล Matérn [ 1 ] เป็นฟังก์ชันความแปรปรวนร่วมที่ใช้ในสถิติเชิงพื้นที่ ภูมิสถิติการเรียนรู้ของเครื่อง การวิเคราะห์ภาพและการประยุกต์ใช้การวิเคราะห์ทางสถิติแบบหลายตัวแปรอื่นๆ บนปริภูมิเมตริก โดยตั้งชื่อตาม Bertil Matérnนักสถิติป่าไม้ชาวสวีเดน[ 2 ]ฟังก์ชันนี้ระบุความแปรปรวนร่วมระหว่างการวัดสองครั้งเป็นฟังก์ชันของระยะทาง{\displaystyle d}ระหว่างจุดที่ทำการเก็บข้อมูล เนื่องจากค่าความแปรปรวนร่วมขึ้นอยู่กับระยะห่างระหว่างจุดเท่านั้น จึงทำให้ค่าความแปรปรวนร่วมคงที่หากระยะห่างเป็นระยะทางแบบยุคลิด ค่าความแปรปรวนร่วมของ Matérn ก็จะเป็น แบบไอโซโทรปิกด้วย

คำนิยาม

ความแปรปรวนร่วมของ Matérn ระหว่างการวัดที่จุดสองจุดที่แยกจากกันด้วย ระยะทาง dหน่วย จะได้รับจาก[ 3 ]

ซีν()=σ221νΓ(ν)(2νρ)νเคν(2νρ),{\displaystyle C_{\nu }(d)=\sigma ^{2}{\frac {2^{1-\nu }}{\Gamma (\nu )}}{{\Bigg (}{\sqrt {2\nu }}{\frac {d}{\rho }}{\Bigg )}}^{\nu }K_{\nu }{\Bigg (}{\sqrt {2\nu }}{\frac {d}{\rho }}{\บิ๊ก )},}

ที่ไหนΓ{\displaystyle \Gamma }คือฟังก์ชันแกมมาเคν{\displaystyle K_{\nu }}คือฟังก์ชันเบสเซลแบบดัดแปลงชนิดที่สอง และρและν{\displaystyle \nu }เป็นพารามิเตอร์ บวก ของความแปรปรวนร่วม

กระบวนการเกาส์เซียนที่มีความแปรปรวนร่วมแบบ Matérn คือν1{\displaystyle \lceil \nu \rceil -1}หลายครั้งที่สามารถหาอนุพันธ์ได้ในแง่ของค่าเฉลี่ยกำลังสอง[ 3 ] [ 4 ]

ความหนาแน่นสเปกตรัม

สเปกตรัมกำลังของกระบวนการที่มีความแปรปรวนร่วมแบบ Matérn ที่กำหนดไว้บนอาร์n{\displaystyle \mathbb {R} ^{n}}คือการแปลงฟูริเยร์ ( nมิติ) ของฟังก์ชันความแปรปรวนร่วมของ Matérn (ดูทฤษฎีบท Wiener–Khinchin ) โดยเฉพาะอย่างยิ่ง กำหนดโดย

เอส(เอฟ)=σ22nπn/2Γ(ν+n2)(2ν)νΓ(ν)ρ2ν(2νρ2+4π2เอฟ2)(ν+n2).{\displaystyle S(f)=\sigma ^{2}{\frac {2^{n}\pi ^{n/2}\Gamma (\nu +{\frac {n}{2}})(2\nu )^{\nu }}{\Gamma (\nu )\rho ^{2\nu }}}\left({\frac {2\nu }{\rho ^{2}}}+4\pi ^{2}f^{2}\right)^{-\left(\nu +{\frac {n}{2}}\right)}.}[ 3 ]

การลดรูปสำหรับค่าν ที่เฉพาะเจาะจง

การลดรูปสำหรับνครึ่งจำนวนเต็ม

เมื่อไร ν=พี+1/2, พีเอ็น+{\displaystyle \nu =p+1/2,\ p\in \mathbb {N} ^{+}}ค่าความแปรปรวนร่วมของ Matérnสามารถเขียนได้ในรูปผลคูณของเลขชี้กำลังและพหุนามดีกรีพี{\displaystyle p}[ 5 ] [ 6 ] ฟังก์ชันเบสเซลที่ดัดแปลง ของอันดับเศษส่วนกำหนดโดยสมการ 10.1.9 และ 10.2.15 [ 7 ]ดังนี้

π2zเคพี+1/2(z)=π2zอีzเค=0n(n+เค)!เค!Γ(nเค+1)(2z)เค{\displaystyle {\sqrt {\frac {\pi }{2z}}}K_{p+1/2}(z)={\frac {\pi }{2z}}e^{-z}\sum _{k=0}^{n}{\frac {(n+k)!}{k!\Gamma (n-k+1)}}\left(2z\right)^{-k}}.

สิ่งนี้ทำให้สามารถใช้ความแปรปรวนร่วมของ Matérn สำหรับค่าครึ่งจำนวนเต็ม ได้ν{\displaystyle \nu }เพื่อแสดงออกมาเป็น

ซีพี+1/2()=σ2เอ็กซ์(2พี+1ρ)พี!(2พี)!ฉัน=0พี(พี+ฉัน)!ฉัน!(พีฉัน)!(22พี+1ρ)พีฉัน,{\displaystyle C_{p+1/2}(d)=\sigma ^{2}\exp \left(-{\frac {{\sqrt {2p+1}}d}{\rho }}\right){\frac {p!}{(2p)!}}\sum _{i=0}^{p}{\frac {(p+i)!}{i!(p-i)!}}\left({\frac {2{\sqrt {2p+1}}d}{\rho }}\right)^{p-i},}

ซึ่งให้ผลลัพธ์ดังนี้:

  • สำหรับν=1/2 (พี=0){\displaystyle \nu =1/2\ (p=0)}:ซี1/2()=σ2เอ็กซ์(ρ),{\displaystyle C_{1/2}(d)=\sigma ^{2}\exp \left(-{\frac {d}{\rho }}\right),}
  • สำหรับν=3/2 (พี=1){\displaystyle \nu =3/2\ (p=1)}:ซี3/2()=σ2(1+3ρ)เอ็กซ์(3ρ),{\displaystyle C_{3/2}(d)=\sigma ^{2}\left(1+{\frac {{\sqrt {3}}d}{\rho }}\right)\exp \left(-{\frac {{\sqrt {3}}d}{\rho }}\right),}
  • สำหรับν=5/2 (พี=2){\displaystyle \nu =5/2\ (p=2)}:ซี5/2()=σ2(1+5ρ+523ρ2)เอ็กซ์(5ρ).{\displaystyle C_{5/2}(d)=\sigma ^{2}\left(1+{\frac {{\sqrt {5}}d}{\rho }}+{\frac {5d^{2}}{3\rho ^{2}}}\right)\exp \left(-{\frac {{\sqrt {5}}d}{\rho }}\right).}

กรณีเกาส์เซียนในขีดจำกัดของν ที่อนันต์

เช่นν{\displaystyle \nu \rightarrow \infty }ฟังก์ชันความแปรปรวนร่วมของ Matérnจะลู่เข้าสู่ฟังก์ชันความแปรปรวน ร่วมแบบเลขชี้กำลังกำลังสอง

ลิมνซีν()=σ2เอ็กซ์(22ρ2).{\displaystyle \lim _{\nu \rightarrow \infty }C_{\nu }(d)=\sigma ^{2}\exp \left(-{\frac {d^{2}}{2\rho ^{2}}}\right).}

อนุกรมเทย์เลอร์ที่ศูนย์และโมเมนต์สเปกตรัม

จากความสัมพันธ์พื้นฐานที่ฟังก์ชันแกมมาเป็นไปตามนั้นΓ(z)Γ(1z)=πบาป(πz){\displaystyle \Gamma (z)\Gamma (1-z)={\frac {\pi }{\sin(\pi z)}}} และความสัมพันธ์พื้นฐานที่ฟังก์ชันเบสเซลดัดแปลงของตัวที่สองเป็นไปตามนั้น

เคν(x)=π2ฉันν(x)ฉันν(x)บาป(πν){\displaystyle K_{\nu }(x)={\frac {\pi }{2}}{\frac {I_{-\nu }(x)-I_{\nu }(x)}{\sin(\pi \nu )}}}

และนิยามของฟังก์ชันเบสเซลที่ดัดแปลงของข้อแรก ฉันν(x)==01!Γ(+ν+1)(x2)2+ν,{\displaystyle I_{\nu }(x)=\sum _{m=0}^{\infty }{\frac {1}{m!\,\Gamma (m+\nu +1)}}\left({\frac {x}{2}}\right)^{2m+\nu },}

พฤติกรรมสำหรับ0{\displaystyle d\rightarrow 0}สามารถหาได้จากอนุกรมเทย์เลอร์ ต่อไปนี้ (เมื่อ ν{\displaystyle \nu }ไม่ใช่จำนวนเต็มและมากกว่า 2):

ซีν()=σ2(1+ν2(1ν)(ρ)2+ν28(23ν+ν2)(ρ)4+โอ(6(2ν))),ν>2.{\displaystyle C_{\nu }(d)=\sigma ^{2}\left(1+{\frac {\nu }{2(1-\nu )}}\left({\frac {d}{\rho }}\right)^{2}+{\frac {\nu ^{2}}{8(2-3\nu +\nu ^{2})}}\left({\frac {d}{\rho }}\right)^{4}+{\mathcal {O}}\left(d^{6\wedge (2\nu )}\right)\right),\,\,\nu >2.}[ 8 ]

เมื่อกำหนดนิยามแล้ว โมเมนต์สเปกตรัมต่อไปนี้สามารถหาได้จากอนุกรมเทย์เลอร์:

λ0=ซีν(0)=σ2,λ2=2ซีν()2|=0=σ2νρ2(ν1).{\displaystyle {\begin{aligned}\lambda _{0}&=C_{\nu }(0)=\sigma ^{2},\\[8pt]\lambda _{2}&=-\left.{\frac {\partial ^{2}C_{\nu }(d)}{\partial d^{2}}}\right|_{d=0}={\frac {\sigma ^{2}\nu }{\rho ^{2}(\nu -1)}}.\end{aligned}}}

สำหรับกรณีของν(0,1)(1,2){\displaystyle \nu \in (0,1)\cup (1,2)}สามารถสร้าง อนุกรมเทย์เลอร์ที่คล้ายกัน ได้:ซีν()=σ2(1+ν2(1ν)(ρ)2Γ(1ν)Γ(1+ν)(ν2)ν(ρ)2ν+โอ(4(2ν+2))),ν(0,1)(1,2).{\displaystyle C_{\nu }(d)=\sigma ^{2}\left(1+{\frac {\nu }{2(1-\nu )}}\left({\frac {d}{\rho }}\right)^{2}-{\frac {\Gamma (1-\nu )}{\Gamma (1+\nu )}}\left({\frac {\nu }{2}}\right)^{\nu }\left({\frac {d}{\rho }}\right)^{2\nu }+{\mathcal {O}}\left(d^{4\wedge (2\nu +2)}\right)\right),\,\,\nu \in (0,1)\cup (1,2).} เมื่อไรν{\displaystyle \nu }หากเป็นจำนวนเต็ม ควรใช้ค่าจำกัด (ดู[ 8 ] )

ดูเพิ่มเติม

ดึงข้อมูลมาจาก " https://en.wikipedia.org/w/index.php?title=Matérn_covariance_function&oldid=1314480248 "

สรุปเนื้อหา

ข้อมูลสำคัญจากบทความ

ข้อมูลสำคัญเกี่ยวกับ ฟังก์ชันความแปรปรวนร่วมของ Matérn

ในทางสถิติค่าความแปรปรวนร่วมของ Matérnหรือที่เรียกว่าเคอร์เนล Matérn เป็นฟังก์ชันความแปรปรวนร่วมที่ใช้ในสถิติเชิงพื้นที่ ภูมิสถิติการเรียนรู้ของเครื่อง

คำนิยาม

ความแปรปรวนร่วมของ Matérn ระหว่างการวัดที่จุดสองจุดที่แยกจากกันด้วย ระยะทาง d หน่วย จะได้รับจาก [ 3 ]

ความหนาแน่นสเปกตรัม

สเปกตรัมกำลังของกระบวนการที่มีความแปรปรวนร่วมแบบ Matérn ที่กำหนดไว้บน อาร์ n {\displaystyle \mathbb {R} ^{n}} คือ การแปลงฟูริเยร์ ( n มิติ) ของฟังก์ชันความแปรปรวนร่วมของ Matérn (ดู ทฤษฎีบท Wiener–Khinchin ) โดยเฉพาะอย่างยิ่ง กำหนดโดย

การลดรูปสำหรับ ν ครึ่งจำนวนเต็ม

เมื่อไร ν = พี + 1 / 2 , พี ∈ เอ็น + {\displaystyle \nu =p+1/2,\ p\in \mathbb {N} ^{+}} ค่า ความแปรปรวนร่วมของ Matérn สามารถเขียนได้ในรูปผลคูณของเลขชี้กำลังและพหุนามดีกรี พี {\displaystyle p} [ 5 ] [ 6 ] ฟังก์ชันเบสเซลที่ดัดแปลง ของ อันดับเศษส่วนกำหนดโดยสมการ...