توزيع مغلف

في نظرية الاحتمالات والإحصاء الاتجاهي ، يُعرف التوزيع الاحتمالي الملتف بأنه توزيع احتمالي مستمر يصف نقاط البيانات التي تقع على كرة الوحدة ذات البعد n . في بُعد واحد، يتكون التوزيع الملتف من نقاط على دائرة الوحدة .ϕ{\displaystyle \phi }هو متغير عشوائي في الفترة(-،){\displaystyle (-\infty ,\infty )}باستخدام دالة كثافة الاحتمال (PDF)ص(ϕ){\displaystyle p(\phi )}، ثمz=هـأناϕ{\displaystyle z=e^{i\phi }}هو متغير دائري موزع وفقًا للتوزيع الملتفصwz(θ){\displaystyle p_{wz}(\theta )}وθ=arg(z){\displaystyle \theta =\arg(z)}هو متغير زاوي في الفترة(-π،π]{\displaystyle (-\pi ,\pi ]}يتم التوزيع وفقًا للتوزيع المغلفصw(θ){\displaystyle p_{w}(\theta )}.

أي دالة كثافة احتماليةص(ϕ){\displaystyle p(\phi )}يمكن "لف" الخط حول محيط دائرة نصف قطرها وحدة واحدة. [ 1 ] أي أن دالة كثافة الاحتمال للمتغير الملفوف

θ=ϕمود2π{\displaystyle \theta =\phi \mod 2\pi }في فترة زمنية ما من الطول2π{\displaystyle 2\pi }

يكون

صw(θ)=ك=-ص(θ+2πك){\displaystyle p_{w}(\theta )=\sum _{k=-\infty }^{\infty }{p(\theta +2\pi k)}}

وهو مجموع دوري ذو فترة2π{\displaystyle 2\pi }الفترة الزمنية المفضلة عموماً هي(-π<θπ){\displaystyle (-\pi <\theta \leq \pi )}والتيln(هـأناθ)=arg(هـأناθ)=θ{\displaystyle \ln(e^{i\theta })=\arg(e^{i\theta })=\theta }.

نظرية

في معظم الحالات، تنتج عملية تتضمن إحصاءات دائرية زوايا (ϕ{\displaystyle \phi }) التي تقع في الفترة(-،){\displaystyle (-\infty ,\infty )}ويتم وصفها بواسطة دالة كثافة احتمالية "غير ملفوفة".ص(ϕ){\displaystyle p(\phi )}ومع ذلك، فإن القياس سيعطي زاويةθ{\displaystyle \theta }والتي تقع في فترة طول معينة2π{\displaystyle 2\pi }(على سبيل المثال، من 0 إلى2π{\displaystyle 2\pi }بمعنى آخر، لا يمكن للقياس أن يحدد ما إذا كانت الزاوية الحقيقيةϕ{\displaystyle \phi }أو زاوية ملفوفةθ=ϕ+2πأ{\displaystyle \theta =\phi +2\pi a}، أينأ{\displaystyle a}هو عدد صحيح غير معروف، وقد تم قياسه.

إذا أردنا حساب القيمة المتوقعة لدالة ما تعتمد على الزاوية المقاسة، فستكون كالتالي:

و(θ)=-ص(ϕ)و(ϕ+2πأ)دϕ{\displaystyle \langle f(\theta )\rangle =\int _{-\infty }^{\infty }p(\phi )f(\phi +2\pi a)d\phi }.

يمكننا التعبير عن التكامل كمجموع تكاملات على فترات زمنية محددة.2π{\displaystyle 2\pi }:

و(θ)=ك=-2πك2π(ك+1)ص(ϕ)و(ϕ+2πأ)دϕ{\displaystyle \langle f(\theta )\rangle =\sum _{k=-\infty }^{\infty }\int _{2\pi k}^{2\pi (k+1)}p(\phi )f(\phi +2\pi a)d\phi }.

تغيير متغير التكامل إلىθ=ϕ-2πك{\displaystyle \theta '=\phi -2\pi k}وبتبادل ترتيب التكامل والجمع، نحصل على

و(θ)=02πصw(θ)و(θ+2πأ)دθ{\displaystyle \langle f(\theta )\rangle =\int _{0}^{2\pi }p_{w}(\theta ')f(\theta '+2\pi a')d\theta '}

أينصw(θ){\displaystyle p_{w}(\theta ')}ملف PDF الخاص بالتوزيع المغلف وأ{\displaystyle a'}هو عدد صحيح مجهول آخر(أ=أ+ك){\displaystyle (a'=a+k)}العدد الصحيح المجهولأ{\displaystyle a'}يُدخل ذلك غموضًا في القيمة المتوقعة لـو(θ){\displaystyle f(\theta )}يشبه هذا مشكلة حساب المتوسط ​​الزاوي . ويمكن حل هذه المشكلة عن طريق إدخال المعاملz=هـأناθ{\displaystyle z=e^{i\theta }}، منذz{\displaystyle z}لها علاقة لا لبس فيها بالزاوية الحقيقيةϕ{\displaystyle \phi }:

z=هـأناθ=هـأناϕ{\displaystyle z=e^{i\theta }=e^{i\phi }}.

حساب القيمة المتوقعة لدالة منz{\displaystyle z}سيؤدي ذلك إلى إجابات لا لبس فيها:

و(z)=02πصw(θ)و(هـأناθ)دθ{\displaystyle \langle f(z)\rangle =\int _{0}^{2\pi }p_{w}(\theta ')f(e^{i\theta '})d\theta '}.

ولهذا السبب، فإنz{\displaystyle z}يُفضل استخدام المعلمة على الزوايا المقاسةθ{\displaystyle \theta }في التحليل الإحصائي الدائري. يشير هذا إلى أن دالة التوزيع الملتف يمكن التعبير عنها كدالة لـz{\displaystyle z}بحيث:

و(z)=صwz(z)و(z)دz{\displaystyle \langle f(z)\rangle =\oint p_{wz}(z)f(z)\,dz}

أينصw(z){\displaystyle p_{w}(z)}يتم تعريفها بحيثصw(θ)|دθ|=صwz(z)|دz|{\displaystyle p_{w}(\theta )\,|d\theta |=p_{wz}(z)\,|dz|}يمكن توسيع هذا المفهوم ليشمل السياق متعدد المتغيرات من خلال توسيع المجموع البسيط إلى عدد منF{\displaystyle F}المجاميع التي تغطي جميع الأبعاد في مساحة الميزات:

صw(θ)=ك1،...،كF=-ص(θ+2πك1هـ1++2πكFهـF){\displaystyle p_{w}({\vec {\theta }})=\sum _{k_{1},...,k_{F}=-\infty }^{\infty }{p({\vec {\theta }}+2\pi k_{1}\mathbf {e} _{1}+\dots +2\pi k_{F}\mathbf {e} _{F})}}

أينهـك=(0،...،0،1،0،...،0)تي{\displaystyle \mathbf {e} _{k}=(0,\dots ,0,1,0,\dots ,0)^{\mathsf {T}}}هوك{\displaystyle k}متجه الأساس الإقليدي رقم 1.

التعبير بدلالة الدوال المميزة

التوزيع الملتف الأساسي هو مشط ديراك ، وهو دالة دلتا ديراك الملتفة :

Δ2π(θ)=ك=-دلتا(θ+2πك){\displaystyle \Delta _{2\pi }(\theta )=\sum _{k=-\infty }^{\infty }{\delta (\theta +2\pi k)}}.

باستخدام دالة دلتا، يمكن كتابة توزيع ملتف عام

صw(θ)=ك=--ص(θ)دلتا(θ-θ+2πك)دθ{\displaystyle p_{w}(\theta )=\sum _{k=-\infty }^{\infty }\int _{-\infty }^{\infty }p(\theta ')\delta (\theta -\theta '+2\pi k)\,d\theta '}.

بتغيير ترتيب الجمع والتكامل، يمكن كتابة أي توزيع ملفوف على أنه التفاف التوزيع غير الملفوف مع مشط ديراك:

صw(θ)=-ص(θ)Δ2π(θ-θ)دθ{\displaystyle p_{w}(\theta )=\int _{-\infty }^{\infty }p(\theta ')\Delta _{2\pi }(\theta -\theta ')\,d\theta '}.

يمكن أيضًا التعبير عن مشط ديراك كمجموع دوال أسية، لذلك يمكننا كتابة ما يلي:

صw(θ)=12π-ص(θ)ن=-هـأنان(θ-θ)دθ{\displaystyle p_{w}(\theta )={\frac {1}{2\pi }}\,\int _{-\infty }^{\infty }p(\theta ')\sum _{n=-\infty }^{\infty }e^{in(\theta -\theta ')}\,d\theta '}.

مع تبديل ترتيب الجمع والتكامل مرة أخرى:

صw(θ)=12πن=--ص(θ)هـأنان(θ-θ)دθ{\displaystyle p_{w}(\theta )={\frac {1}{2\pi }}\,\sum _{n=-\infty }^{\infty }\int _{-\infty }^{\infty }p(\theta ')e^{in(\theta -\theta ')}\,d\theta '}.

باستخدام تعريفϕ(s){\displaystyle \phi (s)}، الوظيفة المميزة لـص(θ){\displaystyle p(\theta )}ينتج عنه متسلسلة لوران حول الصفر للتوزيع الملفوف بدلالة الدالة المميزة للتوزيع غير الملفوف:

صw(θ)=12πن=-ϕ(ن)هـ-أنانθ{\displaystyle p_{w}(\theta )={\frac {1}{2\pi }}\,\sum _{n=-\infty }^{\infty }\phi (n)\,e^{-in\theta }}

أو

صwz(z)=12πن=-ϕ(ن)z-ن{\displaystyle p_{wz}(z)={\frac {1}{2\pi }}\,\sum _{n=-\infty }^{\infty }\phi (n)\,z^{-n}}

على غرار التوزيعات الخطية،ϕ(م){\displaystyle \phi (m)}يُشار إليها باسم الدالة المميزة للتوزيع الملتف (أو بدقة أكبر، المتتالية المميزة ). [ 2 ] هذه حالة من حالات صيغة جمع بواسون ، ويمكن ملاحظة أن معاملات متسلسلة فورييه للتوزيع الملتف هي ببساطة معاملات تحويل فورييه للتوزيع غير الملتف عند القيم الصحيحة.

لحظات

لحظات التوزيع المغلفصw(z){\displaystyle p_{w}(z)}تُعرَّف على النحو التالي:

zم=صwz(z)zمدz{\displaystyle \langle z^{m}\rangle =\oint p_{wz}(z)z^{m}\,dz}.

التعبيرصw(z){\displaystyle p_{w}(z)}من حيث الدالة المميزة وتبديل ترتيب التكامل والجمع ينتج ما يلي:

zم=12πن=-ϕ(ن)zم-ندz{\displaystyle \langle z^{m}\rangle ={\frac {1}{2\pi }}\sum _{n=-\infty }^{\infty }\phi (n)\oint z^{m-n}\,dz}.

من نظرية الباقي لدينا

zم-ندz=2πدلتام-ن{\displaystyle \oint z^{m-n}\,dz=2\pi \delta _{m-n}}

أيندلتاك{\displaystyle \delta _{k}}هي دالة دلتا كرونكر . ويترتب على ذلك أن العزوم تساوي ببساطة الدالة المميزة للتوزيع غير المغلف بالنسبة للوسائط الصحيحة:

zم=ϕ(م){\displaystyle \langle z^{m}\rangle =\phi (m)}.

توليد المتغيرات العشوائية

لوX{\displaystyle X}هو متغير عشوائي يتم سحبه من توزيع احتمالي خطيP{\displaystyle P}، ثمZ=هـأناX{\displaystyle Z=e^{iX}}هو متغير دائري موزع وفقًا للتغليفP{\displaystyle P}التوزيع، وθ=arg(Z){\displaystyle \theta =\arg(Z)}هل يتم توزيع المتغير الزاوي وفقًا للالتفاف؟P{\displaystyle P}التوزيع، مع-π<θπ{\displaystyle -\pi <\theta \leq \pi }.

إنتروبيا

إنتروبيا المعلومات لتوزيع دائري ذي كثافة احتماليةصw(θ){\displaystyle p_{w}(\theta )}يُعرَّف على النحو التالي:

ح=-Γصw(θ)ln(صw(θ))دθ{\displaystyle H=-\int _{\Gamma }p_{w}(\theta )\,\ln(p_{w}(\theta ))\,d\theta }

أينΓ{\displaystyle \Gamma }أي فاصل زمني بطول2π{\displaystyle 2\pi }[ 1 ] إذا كان من الممكن التعبير عن كل من كثافة الاحتمال ولوغاريتمها كمتسلسلة فورييه (أو بشكل عام، أي تحويل تكاملي على الدائرة)، فيمكن استخدام الأساس المتعامد للمتسلسلة للحصول على تعبير مغلق للإنتروبيا.

لحظات التوزيعϕ(ن){\displaystyle \phi (n)}هي معاملات فورييه لتوسيع متسلسلة فورييه لكثافة الاحتمال:

صw(θ)=12πن=-ϕنهـ-أنانθ{\displaystyle p_{w}(\theta )={\frac {1}{2\pi }}\sum _{n=-\infty }^{\infty }\phi _{n}e^{-in\theta }}.

إذا كان من الممكن أيضًا التعبير عن لوغاريتم كثافة الاحتمال كمتسلسلة فورييه:

ln(صw(θ))=م=-جمهـأنامθ{\displaystyle \ln(p_{w}(\theta ))=\sum _{m=-\infty }^{\infty }c_{m}e^{im\theta }}

أين

جم=12πΓln(صw(θ))هـ-أنامθدθ{\displaystyle c_{m}={\frac {1}{2\pi }}\int _{\Gamma }\ln(p_{w}(\theta ))e^{-im\theta }\,d\theta }.

ثم، بتغيير ترتيب التكامل والجمع، يمكن كتابة الإنتروبيا على النحو التالي:

ح=-12πم=-ن=-جمϕنΓهـأنا(م-ن)θدθ{\displaystyle H=-{\frac {1}{2\pi }}\sum _{m=-\infty }^{\infty }\sum _{n=-\infty }^{\infty }c_{m}\phi _{n}\int _{\Gamma }e^{i(m-n)\theta }\,d\theta }.

باستخدام خاصية التعامد في أساس فورييه، يمكن اختزال التكامل إلى:

ح=-ن=-جنϕن{\displaystyle H=-\sum _{n=-\infty }^{\infty }c_{n}\phi _{n}}.

في الحالة الخاصة التي يكون فيها توزيع الاحتمال متناظرًا حول المتوسط،ج-م=جم{\displaystyle c_{-m}=c_{m}}ويمكن كتابة اللوغاريتم على النحو التالي:

ln(صw(θ))=ج0+2م=1جمكوس(مθ){\displaystyle \ln(p_{w}(\theta ))=c_{0}+2\sum _{m=1}^{\infty }c_{m}\cos(m\theta )}

و

جم=12πΓln(صw(θ))كوس(مθ)دθ{\displaystyle c_{m}={\frac {1}{2\pi }}\int _{\Gamma }\ln(p_{w}(\theta ))\cos(m\theta )\,d\theta }

وبما أن التطبيع يتطلب ذلكϕ0=1{\displaystyle \phi _{0}=1}ويمكن كتابة الإنتروبيا على النحو التالي:

ح=-ج0-2ن=1جنϕن{\displaystyle H=-c_{0}-2\sum _{n=1}^{\infty }c_{n}\phi _{n}}.

انظر أيضاً

مراجع

  1. 1 2 مارديا, كانتيلال ; جوب، بيتر إي. (1999). إحصائيات الاتجاه . وايلي. رقم ISBN 978-0-471-95333-3.
  2. مارديا، ك. (1972). إحصاءات البيانات الاتجاهية . نيويورك: دار النشر الأكاديمية. ISBN 978-1-4832-1866-3.