توزيع كوشي المغلف

في نظرية الاحتمالات والإحصاء الاتجاهي ، يُعرف توزيع كوشي الملتف بأنه توزيع احتمالي ملتف ينتج عن التفاف توزيع كوشي حول دائرة الوحدة . يُعرف توزيع كوشي أحيانًا بتوزيع لورنتزي، وقد يُشار إلى توزيع كوشي الملتف أحيانًا بتوزيع لورنتزي ملتف.

غالبًا ما يتم العثور على توزيع كوشي الملفوف في مجال التحليل الطيفي حيث يتم استخدامه لتحليل أنماط الحيود (على سبيل المثال، انظر مقياس التداخل فابري-بيرو ).

وصف

دالة كثافة الاحتمال لتوزيع كوشي الملفوف هي: [ 1 ]

ومرحاض(θ;μ،γ)=ن=-γπ(γ2+(θ-μ+2πن)2)-π<θ<π{\displaystyle f_{\text{WC}}(\theta ;\mu ,\gamma )=\sum _{n=-\infty}^{\infty }{\frac {\gamma }{\pi (\gamma ^{2}+(\theta -\mu +2\pi n)^{2})}}\qquad -\pi <\theta <\pi }

أينγ{\displaystyle \gamma }هو عامل المقياس وμ{\displaystyle \mu }يمثل موضع الذروة للتوزيع "غير الملفوف". ويؤدي التعبير عن دالة كثافة الاحتمال المذكورة أعلاه بدلالة الدالة المميزة لتوزيع كوشي إلى ما يلي:

ومرحاض(θ;μ،γ)=12πن=-هـأنان(θ-μ)-|ن|γ=12πسينهγضرب بالعصاγ-كوس(θ-μ){\displaystyle f_{\text{WC}}(\theta ;\mu ,\gamma )={\frac {1}{2\pi }}\sum _{n=-\infty}^{\infty }e^{in(\theta -\mu ) -|n|\gamma }={\frac {1}{2\pi }}\,\,{\frac {\sinh \gamma }{\cosh \gamma -\cos(\theta -\mu )}}}

يمكن أيضًا التعبير عن دالة كثافة الاحتمال بدلالة المتغير الدائري z = e i θ والمعامل المركب ζ = e i ( μ + i γ )

ومرحاض(z;ζ)=12π1-|ζ|2|z-ζ|2{\displaystyle f_{\text{WC}}(z;\zeta )={\frac {1}{2\pi }}\,\,{\frac {1-|\zeta |^{2}}{|z-\zeta |^{2}}}}

حيث، كما هو موضح أدناه، ζ = ⟨ z ⟩.

فيما يتعلق بالمتغير الدائريz=هـأناθ{\displaystyle z=e^{i\theta }}اللحظات الدائرية لتوزيع كوشي الملفوف هي الدالة المميزة لتوزيع كوشي عند قيم عددية صحيحة:

zن=Γهـأنانθومرحاض(θ;μ،γ)دθ=هـأنانμ-|ن|γ.{\displaystyle \langle z^{n}\rangle =\int _{\Gamma}e^{in\theta }\,f_{\text{WC}}(\theta ;\mu ,\gamma )\,d\theta =e^{in\mu -|n|\gamma }.}

أينΓ{\displaystyle \Gamma \,}هي فترة زمنية ما طولها2π{\displaystyle 2\pi }. اللحظة الأولى هي القيمة المتوسطة لـ z ، والمعروفة أيضًا باسم المحصلة المتوسطة، أو متجه المحصلة المتوسطة:

z=هـأناμ-γ{\displaystyle \langle z\rangle =e^{i\mu -\gamma }}

الزاوية المتوسطة هي

θ=أرزz=μ{\displaystyle \langle \theta \rangle =\mathrm {Arg} \langle z\rangle =\mu }

وطول المحصلة المتوسطة هو

R=|z|=هـ-γ{\displaystyle R=|\langle z\rangle |=e^{-\gamma }}

مما ينتج عنه تباين دائري مقداره 1 − R.

تقدير المعلمات

سلسلة من N قياساتzن=هـأناθن{\displaystyle z_{n}=e^{i\theta _{n}}}يمكن استخدام البيانات المستمدة من توزيع كوشي الملتف لتقدير بعض معلمات التوزيع. متوسط ​​السلسلةz¯{\displaystyle {\overline {z}}}يُعرَّف بأنه

z¯=1شمالن=1شمالzن{\displaystyle {\overline {z}}={\frac {1}{N}}\sum _{n=1}^{N}z_{n}}

وستكون قيمتها المتوقعة هي اللحظة الأولى فقط:

z¯=هـأناμ-γ{\displaystyle \langle {\overline {z}}\rangle =e^{i\mu -\gamma }}

بعبارة أخرى،z¯{\displaystyle {\overline {z}}}هو مُقدِّر غير متحيز للعزم الأول. إذا افترضنا أن موضع الذروةμ{\displaystyle \mu }يقع في الفترة[-π،π){\displaystyle [-\pi ,\pi )}ثم أرج(z¯){\displaystyle ({\overline {z}})}سيكون تقديرًا (متحيزًا) لموقع الذروةμ{\displaystyle \mu }.

مشاهدةzن{\displaystyle z_{n}}باعتبارها مجموعة من المتجهات في المستوى المركب،R¯2{\displaystyle {\overline {R}}^{2}}الإحصائية هي طول المتجه المتوسط:

R¯2=z¯z*¯=(1شمالن=1شمالكوسθن)2+(1شمالن=1شمالالخطيئةθن)2{\displaystyle {\overline {R}}^{2}={\overline {z}}\,{\overline {z^{*}}}=\left({\frac {1}{N}}\sum _{n=1}^{N}\cos \theta _{n}\right)^{2}+\left({\frac {1}{N}}\sum _{n=1}^{N}\sin \theta _{n}\right)^{2}}

وقيمتها المتوقعة هي

R¯2=1شمال+شمال-1شمالهـ-2γ.{\displaystyle \langle {\overline {R}}^{2}\rangle ={\frac {1}{N}}+{\frac {N-1}{N}}e^{-2\gamma }.}

بمعنى آخر، الإحصائية

Rهـ2=شمالشمال-1(R¯2-1شمال){\displaystyle R_{e}^{2}={\frac {N}{N-1}}\left({\overline {R}}^{2}-{\frac {1}{N}}\right)}

سيكون مقدراً غير متحيز لـهـ-2γ{\displaystyle e^{-2\gamma }}، وln(1/Rهـ2)/2{\displaystyle \ln(1/R_{e}^{2})/2}سيكون مقدراً (متحيزاً) لـγ{\displaystyle \gamma }.

إنتروبيا

يتم تعريف إنتروبيا المعلومات لتوزيع كوشي الملتف على النحو التالي: [ 1 ]

ح=-Γومرحاض(θ;μ،γ)ln(ومرحاض(θ;μ،γ))دθ{\displaystyle H=-\int _{\Gamma }f_{\text{WC}}(\theta ;\mu ,\gamma )\,\ln(f_{\text{WC}}(\theta  ;\mu ,\gamma ))\,d\theta }

أينΓ{\displaystyle \Gamma }أي فاصل زمني بطول2π{\displaystyle 2\pi }يمكن كتابة لوغاريتم كثافة توزيع كوشي الملتف على شكل متسلسلة فورييه فيθ{\displaystyle \theta \,}:

ln(ومرحاض(θ;μ،γ))=ج0+2م=1جمكوس(مθ){\displaystyle \ln(f_{\text{WC}}(\theta ;\mu ,\gamma ))=c_{0}+2\sum _{m=1}^{\infty }c_{m}\cos(m\theta )}

أين

جم=12πΓln(سينهγ2π(ضرب بالعصاγ-كوسθ))كوس(مθ)دθ{\displaystyle c_{m}={\frac {1}{2\pi }}\int _{\Gamma }\ln \left({\frac {\sinh \gamma }{2\pi (\cosh \gamma -\cos \theta )}}\right)\cos(m\theta )\,d\theta }

مما ينتج عنه:

ج0=ln(1-هـ-2γ2π){\displaystyle c_{0}=\ln \left({\frac {1-e^{-2\gamma }}{2\pi }}\right)}

(راجع Gradshteyn و Ryzhik [ 2 ] 4.224.15) و

جم=هـ-مγموoرم>0{\displaystyle c_{m}={\frac {e^{-m\gamma }}{m}}\qquad \mathrm {for} \,m>0}

(انظر Gradshteyn و Ryzhik [ 2 ] 4.397.6). تمثيل الدالة المميزة لتوزيع كوشي الملتف في الجانب الأيسر من التكامل هو:

ومرحاض(θ;μ،γ)=12π(1+2ن=1ϕنكوس(نθ)){\displaystyle f_{\text{WC}}(\theta ;\mu ,\gamma )={\frac {1}{2\pi }}\left(1+2\sum _{n=1}^{\infty }\phi _{n}\cos(n\theta )\right)}

أينϕن=هـ-|ن|γ{\displaystyle \phi _{n}=e^{-|n|\gamma }}بإدخال هذه التعبيرات في تكامل الإنتروبيا، وتبديل ترتيب التكامل والجمع، واستخدام تعامد جيوب التمام، يمكن كتابة الإنتروبيا على النحو التالي:

ح=-ج0-2م=1ϕمجم=-ln(1-هـ-2γ2π)-2م=1هـ-2نγن{\displaystyle H=-c_{0}-2\sum _{m=1}^{\infty }\phi _{m}c_{m}=-\ln \left({\frac {1-e^{-2\gamma }}{2\pi }}\right)-2\sum _{m=1}^{\infty }{\frac {e^{-2n\gamma }}{n}}}

المتسلسلة هي ببساطة توسيع تايلور للوغاريتم(1-هـ-2γ){\displaystyle (1-e^{-2\gamma })}لذا يمكن كتابة الإنتروبيا بصيغة مغلقة على النحو التالي:

ح=ln(2π(1-هـ-2γ)){\displaystyle H=\ln(2\pi (1-e^{-2\gamma }))\,}

توزيع كوشي الدائري

إذا كان المتغير X يتبع توزيع كوشي بمتوسط ​​وسيط μ ومعامل مقياس γ ، فإن المتغير المركب

Z=X-أناX+أنا{\displaystyle Z={\frac {X-i}{X+i}}}

له معامل وحدة ويتم توزيعه على دائرة الوحدة بكثافة: [ 3 ]

ونسخة(θ،μ،γ)=12π1-|ζ|2|هـأناθ-ζ|2{\displaystyle f_{\text{CC}}(\theta ,\mu ,\gamma )={\frac {1}{2\pi }}{\frac {1-|\zeta |^{2}}{|e^{i\theta }-\zeta |^{2}}}}

أين

ζ=ψ-أناψ+أنا{\displaystyle \zeta ={\frac {\psi -i}{\psi +i}}}

و ψ تعبر عن المعلمتين لتوزيع كوشي الخطي المرتبط بـ x كعدد مركب :

ψ=μ+أناγ{\displaystyle \psi =\mu +i\gamma \,}

يتضح أن توزيع كوشي الدائري له نفس الشكل الوظيفي لتوزيع كوشي الملتف في z و ζ (أي f WC ( z , ζ )). توزيع كوشي الدائري هو توزيع كوشي ملتف معاد تحديد معالمه.

ونسخة(θ،م،γ)=ومرحاض(هـأناθ،م+أناγ-أنام+أناγ+أنا){\displaystyle f_{\text{CC}}(\theta ,m,\gamma )=f_{\text{WC}}\left(e^{i\theta },\,{\frac {m+i\gamma -i}{m+i\gamma +i}}\right)}

التوزيعونسخة(θ;μ،γ){\displaystyle f_{\text{CC}}(\theta يُطلق على التوزيع ;\mu ,\gamma )} اسم توزيع كوشي الدائري [ 3 ] [ 4 ] (وأيضًا توزيع كوشي المركب [ 3 ] ) بمعاملات μ و γ . (انظر أيضًا إلى معلمات ماكولاغ لتوزيعات كوشي ونواة بواسون للمفاهيم ذات الصلة.)

يتميز توزيع كوشي الدائري، المعبر عنه بالصيغة المركبة، بعزوم محدودة من جميع الرتب.

هـ[Zن]=ζن،هـ[Z¯ن]=ζ¯ن{\displaystyle \operatorname {E} [Z^{n}]=\zeta ^{n},\quad \operatorname {E} [{\bar {Z}}^{n}]={\bar {\zeta }}^{n}}

بالنسبة للأعداد الصحيحة n ≥ 1. أما بالنسبة لـ | φ | < 1، فإن التحويل

يو(z،ϕ)=z-ϕ1-ϕ¯z{\displaystyle U(z,\phi )={\frac {z-\phi }{1-{\bar {\phi }}z}}}

يكون هولومورفي على قرص الوحدة، ويتم توزيع المتغير المحول U ( Z , φ ) كدالة كوشي معقدة مع المعلمة U ( ζ , φ ).

بفرض وجود عينة z 1 ، ... ، z n بحجم n > 2، فإن معادلة الاحتمال الأقصى

ن-1يو(z،ζ^)=ن-1يو(zج،ζ^)=0{\displaystyle n^{-1}U\left(z,{\hat {\zeta }}\right)=n^{-1}\sum U\left(z_{j},{\hat {\zeta }}\right)=0}

يمكن حلها من خلال تكرار بسيط للنقطة الثابتة:

ζ(ر+1)=يو(ن-1يو(z،ζ(ر))،-ζ(ر)){\displaystyle \zeta ^{(r+1)}=U\left(n^{-1}U(z,\zeta ^{(r)}),\,-\zeta ^{(r)}\right)\,}

بدءًا من ζ (0) = 0. سلسلة قيم الاحتمالية غير متناقصة، والحل فريد للعينات التي تحتوي على ثلاث قيم مميزة على الأقل. [ 5 ]

تقدير الاحتمال الأقصى للوسيط (μ^{\displaystyle {\hat {\mu }}}) ومعامل المقياس (γ^{\displaystyle {\hat {\gamma }}}يتم الحصول على ) لعينة كوشي حقيقية عن طريق التحويل العكسي:

μ^±أناγ^=أنا1+ζ^1-ζ^.{\displaystyle {\hat {\mu }}\pm i{\hat {\gamma }}=i{\frac {1+{\hat {\zeta }}}{1-{\hat {\zeta }}}}.}

بالنسبة لـ n ≤ 4، توجد صيغ مغلقة معروفة لـζ^{\displaystyle {\hat {\zeta }}}[ 6 ] كثافة مقدر الاحتمال الأقصى عند t في القرص الواحد تكون بالضرورة على الشكل التالي:

14πصن(χ(ت،ζ))(1-|ت|2)2،{\displaystyle {\frac {1}{4\pi }}{\frac {p_{n}(\chi (t,\zeta ))}{(1-|t|^{2})^{2}}},}

أين

χ(ت،ζ)=|ت-ζ|24(1-|ت|2)(1-|ζ|2){\displaystyle \chi (t,\zeta )={\frac {|t-\zeta |^{2}}{4(1-|t|^{2})(1-|\zeta |^{2})}}}.

تتوفر صيغ p 3 و p 4. [ 7 ]

انظر أيضاً

مراجع

  1. 1 2 مارديا, كانتيلال ; جوب، بيتر إي. (1999). إحصائيات الاتجاه . وايلي. رقم ISBN 978-0-471-95333-3.
  2. 1 2 Gradshteyn, إزرائيل سولومونوفيتش ; ريجيك، يوسف مويسيفيتش ; جيرونيموس، يوري فينيامينوفيتش ؛ تسيتلين ، ميخائيل يوليفيتش (فبراير 2007). جيفري، آلان؛ زويلينجر، دانيال (محرران). جدول التكاملات والمتسلسلات والمنتجات . تمت الترجمة بواسطة Scripta Technica، Inc. (7 ed.). الصحافة الأكاديمية، وشركة ISBN  0-12-373637-4. إل سي سي إن 2010481177 . 
  3. 1 2 3 ماكولاغ، بيتر (يونيو 1992). "الاستدلال الشرطي ونماذج كوشي" (ملف PDF) . Biometrika . 79 (2): 247–259 . doi : 10.1093/biomet/79.2.247 . تاريخ الاسترجاع: 26 يناير 2016 .
  4. كي في مارديا (1972). إحصاءات البيانات الاتجاهية . دار النشر الأكاديمية .
  5. ج. كوباس (1975). "حول أحادية نمط دالة الاحتمال لتوزيع كوشي". Biometrika . 62 (3): 701–704 . doi : 10.1093/biomet/62.3.701 .
  6. فيرغسون، توماس س. (1978). "تقديرات الاحتمال الأقصى لمعلمات توزيع كوشي لعينات بحجم 3 و4". مجلة الجمعية الإحصائية الأمريكية . 73 (361): 211-213 . doi : 10.1080/01621459.1978.10480031 . JSTOR 2286549 . 
  7. ^ ب. ماكولا (1996). “تحويل موبيوس وتقدير معلمة كوشي”. حوليات الإحصاء . 24 (2): 786-808 . دوى : 10.1214/آوس/1032894465 . جستور 2242674 .