احتمال ضئيل

في الإحصاء ، يُعدّ احتمال ويتل تقريبًا لدالة الاحتمال لسلسلة زمنية غاوسية ثابتة . سُمّي هذا الاحتمال نسبةً إلى عالم الرياضيات والإحصاء بيتر ويتل ، الذي قدّمه في أطروحته للدكتوراه عام 1951. [ 1 ] ويُستخدم هذا الاحتمال بشكل شائع في تحليل السلاسل الزمنية ومعالجة الإشارات لتقدير المعلمات والكشف عن الإشارات.

سياق

في نموذج السلاسل الزمنية الغاوسية الثابتة ، تكون دالة الاحتمال (كما هو معتاد في النماذج الغاوسية) دالة لمعاملات المتوسط ​​والتباين المرتبطة بها. مع عدد كبير (شمال{\displaystyle N}) من الملاحظات، الـ (شمال×شمال{\displaystyle N\times N}قد تصبح مصفوفة التغاير كبيرة جدًا، مما يجعل الحسابات مكلفة للغاية عمليًا. ومع ذلك، نظرًا لثبات البيانات، فإن مصفوفة التغاير لها بنية بسيطة نسبيًا، وباستخدام تقريب، يمكن تبسيط الحسابات بشكل كبير (منيا(شمال2){\displaystyle O(N^{2})}ليا(شمالسجل(شمال)){\displaystyle O(N\log(N))}[ 2 ] تتلخص الفكرة في افتراض نموذج غاوسي غير متجانس التباين ذي متوسط ​​صفري في مجال فورييه ؛ ويستند صياغة النموذج إلى تحويل فورييه المنفصل للسلسلة الزمنية وكثافة طيف القدرة الخاصة بها . [ 3 ] [ 4 ] [ 5 ]

تعريف

يتركX1،...،Xشمال{\displaystyle X_{1},\ldots ,X_{N}}لتكن سلسلة زمنية غاوسية ثابتة ذات كثافة طيفية للقدرة ( من جانب واحد )S1(و){\displaystyle S_{1}(f)}، أينشمال{\displaystyle N}يكون التوزيع منتظماً، ويتم أخذ العينات على فترات زمنية ثابتة.Δت{\displaystyle \Delta _{t}}. يتركX~1،...،X~شمال/2+1{\displaystyle {\tilde {X}}_{1},\ldots ,{\tilde {X}}_{N/2+1}}ليكن تحويل فورييه المنفصل (DFT) ذو القيم المركبة للسلسلة الزمنية. عندئذٍ، بالنسبة لدالة احتمال ويتل، يُفترض فعليًا وجود توزيعات غاوسية مستقلة ذات متوسط ​​صفري لجميعX~ج{\displaystyle {\tilde {X}}_{j}}مع اختلافات الأجزاء الحقيقية والخيالية المعطاة بواسطة

متغير(يكرر(X~ج))=متغير(أنا(X~ج))=S1(وج){\displaystyle \operatorname {Var} \left(\operatorname {Re} ({\tilde {X}}_{j})\right)=\operatorname {Var} \left(\operatorname {Im} ({\tilde {X}}_{j})\right)=S_{1}(f_{j})}

أينوج=جشمالΔت{\displaystyle f_{j}={\frac {j}{N\,\Delta _{t}}}}هوج{\displaystyle j}التردد التقريبي لفورييه. يؤدي هذا النموذج التقريبي مباشرةً إلى دالة الاحتمال (اللوغاريتمية).

سجل(P(x1،...،xشمال))-ج(سجل(S1(وج))+|x~ج|2شمال2ΔتS1(وج)){\displaystyle \log \left(P(x_{1},\ldots ,x_{N})\right)\propto -\sum _{j}\left(\log \left(S_{1}(f_{j})\right)+{\frac {|{\tilde {x}}_{j}|^{2}}{{\frac {N}{2\,\Delta _{t}}}S_{1}(f_{j})}}\right)}

أين||{\displaystyle |\cdot |}يشير إلى القيمة المطلقة باستخدام مخطط الفترة|x~ج|2=(يكرر(x~ج))2+(أنا(x~ج))2{\displaystyle |{\tilde {x}}_{j}|^{2}=\left(\operatorname {Re} ({\tilde {x}}_{j})\right)^{2}+\left(\operatorname {Im} ({\tilde {x}}_{j})\right)^{2}}[ 3 ] [ 4 ] [ 6 ]

حالة خاصة من طيف الضوضاء المعروف

في حال افتراض معرفة طيف الضوضاء مسبقًا ، وعدم استنتاج خصائص الضوضاء من البيانات، يمكن تبسيط دالة الاحتمال بشكل أكبر عن طريق تجاهل الحدود الثابتة، مما يؤدي إلى تعبير مجموع المربعات

سجل(P(x1،...،xشمال))-ج|x~ج|2شمال2ΔتS1(وج){\displaystyle \log \left(P(x_{1},\ldots ,x_{N})\right)\;\propto \;-\sum _{j}{\frac {|{\tilde {x}}_{j}|^{2}}{{\frac {N}{2\,\Delta _{t}}}S_{1}(f_{j})}}}

هذا التعبير هو أيضاً أساس المرشح المطابق الشائع .

دقة التقريب

إن دالة ويتل الاحتمالية، بشكل عام، هي مجرد تقريب، ولا تكون دقيقة إلا إذا كان الطيف ثابتًا، أي في الحالة البسيطة للضوضاء البيضاء . وتعتمد كفاءة تقريب ويتل دائمًا على الظروف الخاصة. [ 7 ] [ 8 ]

لاحظ أنه نظرًا لخطية تحويل فورييه، فإن التوزيع الغاوسي في مجال فورييه يستلزم التوزيع الغاوسي في المجال الزمني والعكس صحيح. ما يجعل دالة ويتل دقيقة تقريبًا فقط يرتبط بنظرية أخذ العينات - تأثير تحويل فورييه لعدد محدود فقط من نقاط البيانات، والذي يتجلى أيضًا في صورة تسرب طيفي في مسائل ذات صلة (والذي يمكن تحسينه باستخدام نفس الأساليب، أي التقطيع ). في هذه الحالة، يفترض الفرض الضمني للدورية وجود ارتباط بين العينة الأولى والأخيرة.x1{\displaystyle x_{1}}وxشمال{\displaystyle x_{N}}والتي تُعامل فعلياً على أنها عينات "مجاورة" (مثلx1{\displaystyle x_{1}}وx2{\displaystyle x_{2}}).

التطبيقات

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

يُستخدم احتمال ويتل عادةً لتقدير معلمات الإشارة للإشارات المدفونة في ضوضاء غير بيضاء. ويمكن افتراض أن طيف الضوضاء معروف، [ 9 ] أو يمكن استنتاجه مع معلمات الإشارة. [ 4 ] [ 6 ]

كشف الإشارة

يُجرى الكشف عن الإشارة عادةً باستخدام المرشح المطابق ، الذي يعتمد على دالة ويتل الاحتمالية في حالة معرفة كثافة طيف قدرة الضوضاء. [ 10 ] [ 11 ] يقوم المرشح المطابق فعليًا بملاءمة الإشارة مع البيانات المشوشة باستخدام أقصى احتمالية ، ويستخدم نسبة الاحتمالية الناتجة كإحصائية للكشف. [ 12 ]

يمكن تعميم المرشح المطابق إلى إجراء مماثل يعتمد على توزيع Student-t، وذلك من خلال مراعاة عدم اليقين (مثل عدم اليقين في التقدير ) في طيف الضوضاء. من الناحية التقنية، يمكن استخدام خوارزمية EM هنا، مما يؤدي فعليًا إلى ترشيح مطابق متكرر أو تكراري. [ 12 ]

تقدير الطيف

يُمكن تطبيق دالة الاحتمال ويتل أيضًا لتقدير طيف الضوضاء ، سواء بمفردها أو بالاقتران مع معلمات الإشارة. [ 13 ] [ 14 ]

انظر أيضاً

مراجع

  1. ^ ويتل، ب. (1951). اختبار الفرضيات في تحليل السلاسل الزمنية . أوبسالا: المقفيست وويكسلز بوكتريكيري AB.
  2. هورفيتش، سي. (2002). "تقريب ويتل لدالة الاحتمال" (ملف PDF) . جامعة نيويورك ستيرن .
  3. 1 2 كالدير، م.؛ ديفيس، ر.أ. (1997)، "مقدمة إلى ويتل (1953) "تحليل السلاسل الزمنية الثابتة المتعددة"", في كوتز، س.؛ جونسون، ن. ل. (محرران)، اختراقات في الإحصاء ، سلسلة سبرينغر في الإحصاء، نيويورك: سبرينغر-فيرلاغ، ص 141-169 ، doi : 10.1007/978-1-4612-0667-5_7 ، ISBN  978-0-387-94989-5 انظر أيضًا: كالدير، م.؛ ديفيس، ر.أ. (1996)، "مقدمة إلى ويتل (1953) "تحليل السلاسل الزمنية الثابتة المتعددة"« ، التقرير الفني 1996/41 ، قسم الإحصاء، جامعة ولاية كولورادو
  4. 1 2 3 هانان، إي جيه (1994)، "تقدير احتمالية وتكرار ويتل"، في كيلي، إف بي (محرر)، الاحتمالات والإحصاء والتحسين؛ تكريم لبيتر ويتل ، تشيتشستر: وايلي
  5. باويتان، ي. (1998)، "احتمالية ويتل"، في كوتز، س.؛ ريد، س.ب.؛ بانكس، د.ل. (محررون)، موسوعة العلوم الإحصائية ، المجلد 2، نيويورك: وايلي وأولاده، ص 708-710 ، doi : 10.1002/0471667196.ess0753 ، ISBN   978-0-471-66719-3
  6. 1 2 روفر، سي.؛ ماير، ر.؛ كريستنسن، ن. (2011). "نمذجة الضوضاء المتبقية الملونة في معالجة إشارات موجات الجاذبية". الجاذبية الكلاسيكية والكمية . 28 (1): 025010. arXiv : 0804.3853 . Bibcode : 2011CQGra..28a5010R . doi : 10.1088/0264-9381/28/1/015010 . S2CID 46673503 . 
  7. تشودري، ن.؛ غوسال، س.؛ روي، أ. (2004). "تجاور مقياس ويتل لسلسلة زمنية غاوسية" . بيومتريكا . 91 (4): 211-218 . doi : 10.1093/biomet/91.1.211 .
  8. كونتريراس-كريستان، أ.؛ غوتيريز-بينيا، إ.؛ ووكر، إس. جي. (2006). "ملاحظة حول احتمالية ويتل". الاتصالات في الإحصاء - المحاكاة والحساب . 35 (4): 857-875 . doi : 10.1080/03610910600880203 . S2CID 119395974 . 
  9. فين، إل إس (1992). "الكشف والقياس والإشعاع الجاذبي". مجلة Physical Review D. 46 ( 12): 5236–5249 . arXiv : gr-qc/9209010 . Bibcode : 1992PhRvD..46.5236F . doi : 10.1103/ PhysRevD.46.5236 . PMID 10014913. S2CID 19004097 .  
  10. تورين، جي إل (1960). "مقدمة في المرشحات المتطابقة" . معاملات معهد مهندسي الراديو في نظرية المعلومات . 6 (3): 311-329 . doi : 10.1109/TIT.1960.1057571 . S2CID 5128742 . 
  11. واينشتاين، إل إيه؛ زوباكوف، في دي (1962). استخلاص الإشارات من الضوضاء . إنجلوود كليفس، نيوجيرسي: برنتيس هول.
  12. 1 2 روفر، سي. (2011). "مرشح قائم على توزيع t للطالب للكشف القوي عن الإشارات". مجلة Physical Review D. 84 ( 12) 122004. arXiv : 1109.0442 . Bibcode : 2011PhRvD..84l2004R . doi : 10.1103/PhysRevD.84.122004 .
  13. شودري، ن.؛ غوسال، س.؛ روي، أ. (2004). "التقدير البايزي للكثافة الطيفية لسلسلة زمنية" (ملف PDF) . مجلة الجمعية الإحصائية الأمريكية . 99 (468): 1050-1059 . CiteSeerX 10.1.1.212.2814 . doi : 10.1198/016214504000000557 . S2CID 17906077 .  
  14. إدواردز، إم سي؛ ماير، آر؛ كريستنسن، إن. (2015). "تقدير كثافة القدرة الطيفية شبه البارامترية البايزية في تحليل بيانات الموجات الثقالية". مجلة Physical Review D. 92 ( 6) 064011. arXiv : 1506.00185 . Bibcode : 2015PhRvD..92f4011E . doi : 10.1103/PhysRevD.92.064011 . S2CID 11508218 .