طريقة مونت كارلو

تقريب التوزيع الطبيعي بطريقة مونت كارلو

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

تُستخدم طرق مونت كارلو بشكل أساسي في ثلاث فئات مختلفة من المشكلات: التحسين، والتكامل العددي، وتوليد النتائج من توزيع الاحتمالات. ويمكن استخدامها أيضًا لنمذجة الظواهر ذات عدم اليقين الكبير في المدخلات، مثل حساب خطر فشل محطة الطاقة النووية. غالبًا ما يتم تنفيذ طرق مونت كارلو باستخدام المحاكاة الحاسوبية، ويمكنها تقديم حلول تقريبية للمشكلات التي يصعب حلها أو معقدة للغاية بحيث يصعب تحليلها رياضيًا.

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

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

ملخص

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

  1. تحديد مجال المدخلات الممكنة
  2. إنشاء مدخلات عشوائيًا من توزيع الاحتمالات على المجال
  3. إجراء حساب حتمي للمخرجات
  4. تجميع النتائج
طريقة مونت كارلو المطبقة لتقريب قيمة π

على سبيل المثال، ضع في اعتبارك ربع دائرة (قطاع دائري) محصور في مربع الوحدة . مع العلم أن نسبة مساحتيهما هيπ/4، يمكن تقريب قيمة π باستخدام طريقة مونت كارلو: [1]

  1. ارسم مربعًا، ثم ارسم ربعًا داخله
  2. توزيع عدد معين من النقاط بشكل موحد على المربع
  3. احسب عدد النقاط داخل الربع، أي التي تكون المسافة بينها وبين الأصل أقل من 1
  4. نسبة العدد الداخلي وعدد العينة الإجمالي هي تقدير لنسبة المنطقتين ،π/4. اضرب النتيجة في 4 لتقدير π .

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

هناك اعتباران مهمان:

  1. إذا لم يتم توزيع النقاط بشكل متساوٍ، فسيكون التقريب ضعيفًا.
  2. التقريب ضعيف بشكل عام إذا تم وضع عدد قليل من النقاط عشوائيًا في المربع بأكمله. في المتوسط، يتحسن التقريب مع وضع المزيد من النقاط.

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

طلب

تُستخدم طرق مونت كارلو غالبًا في المشكلات الفيزيائية والرياضية وتكون مفيدة للغاية عندما يكون من الصعب أو المستحيل استخدام طرق أخرى. تُستخدم طرق مونت كارلو بشكل أساسي في ثلاث فئات من المشكلات: [ 2] التحسين والتكامل العددي وتوليد السحوبات من توزيع الاحتمالات .

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

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

من حيث المبدأ، يمكن استخدام طرق مونت كارلو لحل أي مشكلة لها تفسير احتمالي. وفقًا لقانون الأعداد الكبيرة ، يمكن تقريب التكاملات الموصوفة بالقيمة المتوقعة لبعض المتغيرات العشوائية عن طريق أخذ المتوسط ​​التجريبي ( المعروف أيضًا باسم "متوسط ​​العينة") للعينات المستقلة للمتغير. عندما يتم تحديد معلمات توزيع احتمالات المتغير، يستخدم علماء الرياضيات غالبًا عينة سلسلة ماركوف مونت كارلو (MCMC). [4] [5] [6] الفكرة الأساسية هي تصميم نموذج سلسلة ماركوف حكيم بتوزيع احتمالي ثابت محدد . أي أنه في الحد الأقصى، ستكون العينات التي يتم إنشاؤها بواسطة طريقة MCMC عينات من التوزيع المطلوب (المستهدف). [7] [8] وفقًا لنظرية إرجوديك ، يتم تقريب التوزيع الثابت من خلال المقاييس التجريبية للحالات العشوائية لعينة MCMC.

في مشاكل أخرى، الهدف هو توليد سحوبات من سلسلة من توزيعات الاحتمالات التي تلبي معادلة التطور غير الخطية. يمكن دائمًا تفسير تدفقات توزيعات الاحتمالات هذه على أنها توزيعات للحالات العشوائية لعملية ماركوف التي تعتمد احتمالات انتقالها على توزيعات الحالات العشوائية الحالية (انظر عمليات ماكين-فلاسوف ، معادلة الترشيح غير الخطية ). [9] [10] في حالات أخرى، ينشأ تدفق لتوزيعات الاحتمالات بمستوى متزايد من تعقيد أخذ العينات (نماذج مساحات المسار ذات أفق زمني متزايد، ومقاييس بولتزمان-جيبس المرتبطة بمعلمات درجة الحرارة المتناقصة، وغيرها الكثير). يمكن أيضًا اعتبار هذه النماذج تطورًا لقانون الحالات العشوائية لسلسلة ماركوف غير الخطية. [10] [11] إحدى الطرق الطبيعية لمحاكاة عمليات ماركوف غير الخطية المعقدة هذه هي أخذ عينات من نسخ متعددة من العملية، واستبدال التوزيعات غير المعروفة للحالات العشوائية في معادلة التطور بالمقاييس التجريبية المأخوذة من العينات . وعلى النقيض من منهجيات مونت كارلو وMCMC التقليدية، تعتمد تقنيات الجسيمات ذات المجال المتوسط ​​على عينات متفاعلة متسلسلة. ويعكس مصطلح المجال المتوسط ​​حقيقة مفادها أن كل عينة ( المعروفة أيضًا باسم الجسيمات أو الأفراد أو السائرين أو الوكلاء أو المخلوقات أو النمط الظاهري) تتفاعل مع المقاييس التجريبية للعملية. وعندما يميل حجم النظام إلى اللانهاية، تتقارب هذه المقاييس التجريبية العشوائية مع التوزيع الحتمي للحالات العشوائية لسلسلة ماركوف غير الخطية، بحيث يختفي التفاعل الإحصائي بين الجسيمات.

مونت كارلو بسيطة

لنفترض أن شخصًا يريد معرفة القيمة المتوقعة μ لمجتمع ما (ويعلم أن μ موجودة)، ولكن ليس لديه صيغة متاحة لحسابها. تقدم طريقة مونت كارلو البسيطة تقديرًا لـ μ من خلال تشغيل n محاكاة ومتوسط ​​نتائج المحاكاة. لا توجد قيود على توزيع احتمالات المدخلات في المحاكاة، وتتطلب فقط أن يتم إنشاء المدخلات عشوائيًا وأن تكون مستقلة عن بعضها البعض وأن μ موجودة. ستنتج قيمة n كبيرة بما يكفي قيمة لـ m قريبة بشكل تعسفي من μ ؛ وبشكل أكثر رسمية، ستكون الحالة هي أنه لأي ε > 0، | μm | ≤ ε . [12]

عادةً، الخوارزمية للحصول على m هي

s = 0؛ 
 بالنسبة إلى  i = 1 إلى  قم 
   بتشغيل المحاكاة للمرة i ، مما يعطي النتيجة r i ؛
    s = s + r i ؛
 كرر 
m = s / n ؛

مثال

لنفترض أننا نريد أن نعرف عدد المرات التي ينبغي لنا أن نتوقع فيها رمي ثلاثة أحجار نرد ذات ثمانية أوجه حتى يكون مجموع رميات النرد T على الأقل . نحن نعلم أن القيمة المتوقعة موجودة. يتم توزيع رميات النرد بشكل عشوائي ومستقل عن بعضها البعض. لذا فإن طريقة مونت كارلو البسيطة قابلة للتطبيق:

س = 0؛ 
 بالنسبة إلى  i = 1 إلى  قم 
   برمي النرد الثلاثة حتى يتم الوصول إلى T أو تجاوزه أولاً ؛ ر i = عدد الرميات ؛
    س = س + ر i ؛ 
كرر 
م = س / ن ؛

إذا كان n كبيرًا بدرجة كافية، فسيكون m ضمن ε من μ لأي ε > 0.

تحديد حجم كبير بما فيه الكفايةن

الصيغة العامة

ليكن ε = | μm | > 0. اختر مستوى الثقة المطلوب - النسبة المئوية لاحتمالية أن يكون m ضمن ε من μ عند اكتمال خوارزمية مونت كارلو . ليكن z هو الدرجة z المقابلة لمستوى الثقة هذا.

لنفترض أن s 2 هو التباين المقدر، والذي يُطلق عليه أحيانًا تباين "العينة"؛ وهو تباين النتائج التي تم الحصول عليها من عدد صغير نسبيًا من عمليات محاكاة "العينة" k . اختر k ؛ لاحظ درييلز وشين أنه "حتى بالنسبة لأحجام العينات التي تقل بمقدار مرتبة واحدة عن العدد المطلوب، فإن حساب هذا العدد مستقر تمامًا". [13]

تحسب الخوارزمية التالية s 2 في تمريرة واحدة مع تقليل احتمالية أن يؤدي الخطأ العددي المتراكم إلى إنتاج نتائج خاطئة: [12]

س 1 = 0؛
قم بتشغيل المحاكاة للمرة الأولى، مما ينتج عنه النتيجة r 1 ؛
 m 1 = r 1 ؛ // m i هو متوسط ​​أول i محاكاة 
لـ i = 2 إلى k  قم
    بتشغيل المحاكاة للمرة i ، مما ينتج عنه النتيجة r i ؛ 
   δ i = r i - m i -1 ؛ 
   m i = m i-1  + (1/ i ) δ i ؛ 
   s i = s i-1  + (( i - 1)/ i )( δ i ) 2 ؛ 
كرر 
s 2 = s k /( k - 1)؛

لاحظ أنه عند اكتمال الخوارزمية، يكون m k هو متوسط ​​نتائج k .

n كبير بما فيه الكفاية عندما

[12] [13]

إذا كان nk ، فإن m k = m ؛ وقد تم إجراء محاكاة عينة كافية لضمان أن m k ضمن ε من μ . إذا كان n > k ، فيمكن إجراء n محاكاة "من الصفر"، أو نظرًا لأن k محاكاة قد تم إجراؤها بالفعل، فيمكن للمرء إجراء n - k محاكاة أخرى وإضافة نتائجها إلى تلك المحاكاة من العينات:

س = م ك * ك ؛
بالنسبة إلى i = k + 1 إلى  قم 
   بتشغيل المحاكاة للمرة i ، مما يعطي النتيجة r i ؛
    s = s + r i ؛
 m = s / n ؛

صيغة عندما تكون نتائج المحاكاة محدودة

يمكن استخدام صيغة بديلة في الحالة الخاصة حيث تكون جميع نتائج المحاكاة محدودة من الأعلى والأسفل.

اختر قيمة لـ ε تساوي ضعف الحد الأقصى المسموح به للفرق بين μ و m. دع 0 < δ < 100 يكون مستوى الثقة المطلوب، معبرًا عنه كنسبة مئوية. دع كل نتيجة محاكاة r 1 , r 2 , … r i , … r n تكون بحيث ar ib لـ a و b المحدودتين . للحصول على ثقة لا تقل عن δ في أن | μm | < ε /2، استخدم قيمة لـ n بحيث

على سبيل المثال، إذا كانت δ = 99%، فإن n ≥ 2( ba ) 2 ln(2/0.01)/ ε 2 ≈ 10.6( ba ) 2 / ε 2 . [12]

التكاليف الحسابية

على الرغم من بساطتها المفاهيمية والخوارزمية، فإن التكلفة الحسابية المرتبطة بمحاكاة مونت كارلو يمكن أن تكون مرتفعة بشكل مذهل. بشكل عام، تتطلب الطريقة العديد من العينات للحصول على تقريب جيد، مما قد يتسبب في وقت تشغيل إجمالي كبير بشكل تعسفي إذا كان وقت معالجة عينة واحدة مرتفعًا. [14] على الرغم من أن هذا يمثل قيدًا شديدًا في المشكلات المعقدة للغاية، فإن الطبيعة المتوازية المحرجة للخوارزمية تسمح بتقليل هذه التكلفة الكبيرة (ربما إلى مستوى ممكن) من خلال استراتيجيات الحوسبة المتوازية في المعالجات المحلية والمجموعات والحوسبة السحابية ووحدة معالجة الرسومات ووحدة FPGA وما إلى ذلك. [15] [16] [17] [18]

تاريخ

قبل تطوير طريقة مونت كارلو، كانت عمليات المحاكاة تختبر مشكلة حتمية مفهومة مسبقًا، وتم استخدام أخذ العينات الإحصائية لتقدير أوجه عدم اليقين في عمليات المحاكاة. تعكس عمليات محاكاة مونت كارلو هذا النهج، حيث تحل المشكلات الحتمية باستخدام الخوارزميات الاحتمالية (انظر محاكاة التلدين ).

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

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

كانت الأفكار والمحاولات الأولى التي قمت بها لممارسة [طريقة مونت كارلو] قد اقترحها سؤال خطر ببالي في عام 1946 بينما كنت أتعافى من مرض وألعب لعبة سوليتير. كان السؤال هو ما هي احتمالات نجاح لعبة سوليتير كانفيلد الموضوعة باستخدام 52 بطاقة؟ بعد قضاء الكثير من الوقت في محاولة تقديرها من خلال الحسابات التوافقية البحتة، تساءلت عما إذا كانت الطريقة الأكثر عملية من "التفكير المجرد" قد لا تكون طرحها مائة مرة ومراقبة عدد اللعبات الناجحة وحسابها ببساطة. كان من الممكن بالفعل تصور ذلك مع بداية العصر الجديد من أجهزة الكمبيوتر السريعة، وفكرت على الفور في مشاكل انتشار النيوترونات وغيرها من أسئلة الفيزياء الرياضية، وبشكل عام كيفية تغيير العمليات الموصوفة بمعادلات تفاضلية معينة إلى شكل مكافئ يمكن تفسيره على أنه سلسلة من العمليات العشوائية. لاحقًا [في عام 1946]، وصفت الفكرة لجون فون نيومان ، وبدأنا في التخطيط للحسابات الفعلية. [20]

نظرًا لكونه سريًا، فقد تطلب عمل فون نيومان وأولام اسمًا رمزيًا. [21] اقترح أحد زملاء فون نيومان وأولام، نيكولاس متروبوليس ، استخدام اسم مونت كارلو ، الذي يشير إلى كازينو مونت كارلو في موناكو حيث كان عم أولام يقترض المال من الأقارب للمقامرة. [19] كانت طرق مونت كارلو أساسية في عمليات المحاكاة المطلوبة لمزيد من تطوير الأسلحة النووية بعد الحرب، بما في ذلك تصميم القنبلة الهيدروجينية، على الرغم من محدوديتها الشديدة بواسطة الأدوات الحسابية في ذلك الوقت. قام فون نيومان ونيكولاس متروبوليس وآخرون ببرمجة كمبيوتر إينياك لإجراء أول حسابات مونت كارلو آلية بالكامل، لنواة سلاح انشطار ، في ربيع عام 1948. [22] في الخمسينيات من القرن الماضي، استُخدمت طرق مونت كارلو في لوس ألاموس لتطوير القنبلة الهيدروجينية ، وأصبحت شائعة في مجالات الفيزياء والكيمياء الفيزيائية وبحوث العمليات . كانت مؤسسة راند والقوات الجوية الأمريكية اثنتين من المنظمات الرئيسية المسؤولة عن تمويل ونشر المعلومات حول أساليب مونت كارلو خلال هذا الوقت، وبدأت في العثور على تطبيق واسع في العديد من المجالات المختلفة.

بدأت نظرية طرق مونت كارلو للجسيمات من النوع الحقلي المتوسط ​​الأكثر تطورًا بالتأكيد بحلول منتصف الستينيات، مع عمل هنري ب. ماكين جونيور على تفسيرات ماركوف لفئة من المعادلات التفاضلية الجزئية المكافئة غير الخطية الناشئة في ميكانيكا الموائع. [23] [24] استخدمت مقالة رائدة سابقة لثيودور إي هاريس وهيرمان كان، نُشرت عام 1951، طرق مونت كارلو الجينية من النوع الحقلي المتوسط ​​لتقدير طاقات انتقال الجسيمات. [25] تُستخدم أيضًا منهجيات مونت كارلو الجينية من النوع الحقلي المتوسط ​​كخوارزميات بحث طبيعية استدلالية (المعروفة أيضًا باسم الاستدلالية ) في الحوسبة التطورية. يمكن تتبع أصول تقنيات الحوسبة الحقلية المتوسطة هذه إلى عامي 1950 و1954 مع عمل آلان تورينج على آلات التعلم لاختيار الطفرة الجينية [26] والمقالات التي كتبها نيلز آل باريسيلي في معهد الدراسات المتقدمة في برينستون، نيو جيرسي . [27] [28]

يمكن أيضًا تفسير معادلة مونت كارلو الكمومية ، وبشكل أكثر تحديدًا معادلة مونت كارلو الانتشارية، على أنها تقريب مونت كارلو لجسيمات المجال المتوسط ​​لتكاملات مسار فاينمان - كاك . [29] [30] [31] [32] [33] [34] [35] غالبًا ما تُنسب أصول طرق مونت كارلو الكمومية إلى إنريكو فيرمي وروبرت ريتشماير اللذين طورا في عام 1948 تفسيرًا متوسطًا لتفاعلات سلسلة النيوترون باستخدام جسيمات المجال المتوسط، [36] ولكن أول خوارزمية للجسيمات الشبيهة بالاستدلال والوراثية (المعروفة أيضًا باسم طرق مونت كارلو المعاد أخذ العينات منها أو إعادة تكوينها) لتقدير طاقات الحالة الأساسية للأنظمة الكمومية (في نماذج المصفوفة المخفضة) ترجع إلى جاك إتش هيثرينغتون في عام 1984. [35] في الكيمياء الجزيئية، يمكن إرجاع استخدام منهجيات الجسيمات الشبيهة بالاستدلال الوراثي (المعروفة أيضًا باسم استراتيجيات التقليم والإثراء) إلى عام 1955 مع العمل الرائد لمارشال إن. روزنبلوث وأريانا دبليو روزنبلوث . [37]

إن استخدام خوارزمية مونت كارلو المتسلسلة في معالجة الإشارات المتقدمة والاستدلال البايزي هو أكثر حداثة. ففي عام 1993، نشر جوردون وآخرون، في عملهم الرائد [38]، أول تطبيق لخوارزمية إعادة أخذ العينات مونت كارلو في الاستدلال الإحصائي البايزي. أطلق المؤلفون على خوارزميتهم اسم "مرشح التمهيد"، وأثبتوا أنه مقارنة بطرق الترشيح الأخرى، فإن خوارزمية التمهيد الخاصة بهم لا تتطلب أي افتراض حول مساحة الحالة أو ضوضاء النظام. كانت مقالة جينشيرو كيتاغاوا الرائدة الأخرى في هذا المجال، حول "مرشح مونت كارلو" ذي الصلة، [39] وتلك التي كتبها بيير ديل مورال [40] وهيميلكون كارفاليو وبيير ديل مورال وأندريه مونين وجيرارد سالوت [41] حول مرشحات الجسيمات المنشورة في منتصف التسعينيات. تم تطوير مرشحات الجسيمات أيضًا في معالجة الإشارات في الفترة من 1989 إلى 1992 بواسطة P. Del Moral و JC Noyer و G. Rigal و G. Salut في LAAS-CNRS في سلسلة من التقارير البحثية المقيدة والمصنفة مع STCAN (Service Technique des Constructions et Armes Navales) وشركة تكنولوجيا المعلومات DIGILOG و LAAS-CNRS (مختبر تحليل وهندسة الأنظمة) حول مشاكل معالجة إشارات الرادار / السونار ونظام تحديد المواقع العالمي (GPS). [42] [43] [44] [45] [46] [47] يمكن تفسير منهجيات مونت كارلو المتسلسلة هذه على أنها عينة قبول ورفض مزودة بآلية إعادة تدوير متفاعلة.

من عام 1950 إلى عام 1996، قدمت جميع المنشورات حول منهجيات مونت كارلو المتسلسلة، بما في ذلك طرق مونت كارلو للتقليم وإعادة أخذ العينات التي تم تقديمها في الفيزياء الحاسوبية والكيمياء الجزيئية، خوارزميات طبيعية وخوارزميات استدلالية مطبقة على مواقف مختلفة دون دليل واحد على اتساقها، ولا مناقشة حول تحيز التقديرات والخوارزميات القائمة على الشجرة الأنسابية والأسلاف. كتب بيير ديل مورال الأسس الرياضية والتحليل الدقيق الأول لهذه الخوارزميات الجسيمية في عام 1996. [40] [48]

كما تم تطوير منهجيات الجسيمات المتفرعة ذات أحجام السكان المتنوعة في نهاية التسعينيات بواسطة دان كريسان وجيسيكا جاينز وتيري ليونز، [49] [50] [51] ودان كريسان وبيير ديل مورال وتيري ليونز. [52] وتم وصف المزيد من التطورات في هذا المجال في الفترة من 1999 إلى 2001 بواسطة ب. ديل مورال وأ. جوينيت ول. مايكرو. [30] [53] [54]

التعاريف

لا يوجد إجماع حول كيفية تعريف مونت كارلو . على سبيل المثال، يعرّف ريبلي [55] معظم النمذجة الاحتمالية بأنها محاكاة عشوائية ، مع حجز مونت كارلو لتكامل مونت كارلو واختبارات مونت كارلو الإحصائية. يميز ساويلوفسكي [56] بين المحاكاة وطريقة مونت كارلو ومحاكاة مونت كارلو: المحاكاة هي تمثيل خيالي للواقع، وطريقة مونت كارلو هي تقنية يمكن استخدامها لحل مشكلة رياضية أو إحصائية، وتستخدم محاكاة مونت كارلو أخذ العينات المتكررة للحصول على الخصائص الإحصائية لبعض الظواهر (أو السلوك).

وفيما يلي بعض الأمثلة:

  • المحاكاة: يمكن استخدام سحب متغير عشوائي موحد من الفاصل الزمني [0,1] لمحاكاة رمي عملة معدنية: إذا كانت القيمة أقل من أو تساوي 0.50، فحدد النتيجة على أنها وجه العملة، ولكن إذا كانت القيمة أكبر من 0.50، فحدد النتيجة على أنها ذيل العملة. هذه محاكاة، ولكنها ليست محاكاة مونت كارلو.
  • طريقة مونت كارلو: سكب صندوق من العملات المعدنية على طاولة، ثم حساب نسبة العملات المعدنية التي تسقط على الوجه مقابل الظهر هي طريقة مونت كارلو لتحديد سلوك رمي العملة المعدنية المتكرر، ولكنها ليست محاكاة.
  • محاكاة مونت كارلو: رسم عدد كبير من المتغيرات العشوائية المنتظمة من الفاصل [0,1] في وقت واحد، أو مرة واحدة في أوقات مختلفة عديدة، وتعيين قيم أقل من أو تساوي 0.50 كرؤوس وأكبر من 0.50 كذيول، هي محاكاة مونت كارلو لسلوك رمي العملة المعدنية بشكل متكرر.

يشير كالوس وويتلوك [57] إلى أن مثل هذه التمييزات ليس من السهل الحفاظ عليها دائمًا. على سبيل المثال، فإن انبعاث الإشعاع من الذرات هو عملية عشوائية طبيعية. يمكن محاكاتها بشكل مباشر، أو يمكن وصف سلوكها المتوسط ​​من خلال معادلات عشوائية يمكن حلها بنفسها باستخدام طرق مونت كارلو. "في الواقع، يمكن النظر إلى نفس الكود الحاسوبي في نفس الوقت باعتباره "محاكاة طبيعية" أو كحل للمعادلات عن طريق أخذ العينات الطبيعية."

يمكن التحقق من تقارب محاكاة مونت كارلو باستخدام إحصائية جلمان-روبين .

مونت كارلو والأرقام العشوائية

الفكرة الرئيسية وراء هذه الطريقة هي أن النتائج تُحسب بناءً على أخذ عينات عشوائية متكررة وتحليل إحصائي. محاكاة مونت كارلو هي في الواقع تجارب عشوائية، في حالة عدم معرفة نتائج هذه التجارب جيدًا. تتميز محاكاة مونت كارلو عادةً بالعديد من المعلمات غير المعروفة، والتي يصعب الحصول على العديد منها تجريبيًا. [58] لا تتطلب طرق محاكاة مونت كارلو دائمًا أرقامًا عشوائية حقيقية لتكون مفيدة (على الرغم من أن عدم القدرة على التنبؤ أمر حيوي لبعض التطبيقات مثل اختبار البدائية ). [59] تستخدم العديد من التقنيات الأكثر فائدة تسلسلات شبه عشوائية حتمية ، مما يجعل من السهل اختبار وإعادة تشغيل المحاكاة. الجودة الوحيدة اللازمة عادةً لإجراء محاكاة جيدة هي أن يظهر التسلسل شبه العشوائي "عشوائيًا بدرجة كافية" بمعنى معين.

يعتمد ما يعنيه هذا على التطبيق، ولكن عادةً ما يجب أن تجتاز سلسلة من الاختبارات الإحصائية. يعد اختبار توزيع الأرقام بشكل موحد أو اتباع توزيع مرغوب آخر عند النظر في عدد كبير بما يكفي من عناصر التسلسل أحد أبسط الاختبارات وأكثرها شيوعًا. غالبًا ما تكون الارتباطات الضعيفة بين العينات المتعاقبة مرغوبة/ضرورية أيضًا.

يسرد ساويلوفسكي خصائص محاكاة مونت كارلو عالية الجودة: [56]

  • يحتوي مولد الأرقام (شبه العشوائي) على خصائص معينة (على سبيل المثال "فترة" طويلة قبل تكرار التسلسل)
  • يقوم مولد الأرقام (شبه العشوائي) بإنتاج قيم تجتاز اختبارات العشوائية
  • هناك عينات كافية لضمان نتائج دقيقة
  • يتم استخدام تقنية أخذ العينات المناسبة
  • الخوارزمية المستخدمة صالحة لما يتم نمذجته
  • فهو يحاكي الظاهرة المعنية.

تُستخدم خوارزميات أخذ العينات من الأرقام العشوائية الزائفة لتحويل الأرقام العشوائية الزائفة الموزعة بشكل موحد إلى أرقام موزعة وفقًا لتوزيع احتمالي معين .

غالبًا ما تُستخدم تسلسلات التباين المنخفض بدلاً من أخذ العينات العشوائية من مساحة ما لأنها تضمن تغطية متساوية وعادةً ما يكون ترتيب التقارب أسرع من محاكاة مونت كارلو باستخدام تسلسلات عشوائية أو شبه عشوائية. تُسمى الطرق القائمة على استخدامها بأساليب مونت كارلو شبه العشوائية .

في محاولة لتقييم تأثير جودة الأرقام العشوائية على نتائج محاكاة مونت كارلو، اختبر الباحثون في الفيزياء الفلكية أرقامًا شبه عشوائية آمنة تشفيريًا تم إنشاؤها عبر مجموعة تعليمات RDRAND من Intel ، مقارنةً بتلك المشتقة من الخوارزميات، مثل Mersenne Twister ، في محاكاة مونت كارلو للتوهجات الراديوية من الأقزام البنية . لم يتم العثور على فرق ذي دلالة إحصائية بين النماذج التي تم إنشاؤها باستخدام مولدات الأرقام شبه العشوائية النموذجية وRDRAND للتجارب المكونة من توليد 10 7 أرقام عشوائية. [60]

محاكاة مونت كارلو مقابل سيناريوهات "ماذا لو"

هناك طرق لاستخدام الاحتمالات التي لا تعتبر بالتأكيد محاكاة مونت كارلو - على سبيل المثال، النمذجة الحتمية باستخدام تقديرات النقطة الواحدة. يتم تعيين تقدير "أفضل تخمين" لكل متغير غير مؤكد داخل النموذج. يتم اختيار السيناريوهات (مثل أفضل أو أسوأ أو الحالة الأكثر احتمالاً) لكل متغير إدخال وتسجيل النتائج. [61]

على النقيض من ذلك، تقوم محاكاة مونت كارلو بأخذ عينات من توزيع الاحتمالات لكل متغير لإنتاج مئات أو آلاف النتائج المحتملة. يتم تحليل النتائج للحصول على احتمالات حدوث نتائج مختلفة. [62] على سبيل المثال، تُظهر مقارنة نموذج بناء تكلفة جدول بيانات يتم تشغيله باستخدام سيناريوهات "ماذا لو" التقليدية، ثم تشغيل المقارنة مرة أخرى باستخدام محاكاة مونت كارلو وتوزيعات الاحتمالات المثلثية أن تحليل مونت كارلو له نطاق أضيق من تحليل "ماذا لو". [ مطلوب مثال ] وذلك لأن تحليل "ماذا لو" يعطي وزنًا متساويًا لجميع السيناريوهات (انظر تحديد عدم اليقين في التمويل المؤسسي )، بينما نادرًا ما تأخذ طريقة مونت كارلو عينات في مناطق الاحتمالية المنخفضة للغاية. تسمى العينات في مثل هذه المناطق "أحداث نادرة".

التطبيقات

تعتبر طرق مونت كارلو مفيدة بشكل خاص لمحاكاة الظواهر ذات عدم اليقين الكبير في المدخلات والأنظمة ذات درجات الحرية المتعددة . تشمل مجالات التطبيق ما يلي:

العلوم الفيزيائية

تعتبر طرق مونت كارلو مهمة جدًا في الفيزياء الحاسوبية والكيمياء الفيزيائية والمجالات التطبيقية ذات الصلة، ولها تطبيقات متنوعة من حسابات الديناميكا اللونية الكمومية المعقدة إلى تصميم الدروع الحرارية والأشكال الديناميكية الهوائية وكذلك في نمذجة نقل الإشعاع لحسابات قياس جرعات الإشعاع. [63] [64] [65] في الفيزياء الإحصائية ، تعد النمذجة الجزيئية لمونت كارلو بديلاً للديناميكيات الجزيئية الحاسوبية ، وتُستخدم طرق مونت كارلو لحساب نظريات المجال الإحصائي لأنظمة الجسيمات والبوليمر البسيطة. [37] [66] تحل طرق مونت كارلو الكمومية مشكلة الأجسام المتعددة للأنظمة الكمومية. [9] [10] [29] في علم المواد الإشعاعية ، يعتمد تقريب الاصطدام الثنائي لمحاكاة زرع الأيونات عادةً على نهج مونت كارلو لاختيار الذرة المتصادمة التالية. [67] في فيزياء الجسيمات التجريبية ، تُستخدم طرق مونت كارلو لتصميم الكواشف وفهم سلوكها ومقارنة البيانات التجريبية بالنظرية. في الفيزياء الفلكية ، يتم استخدامها بطرق متنوعة لنمذجة تطور المجرات [68] وانتقال الإشعاع الميكروي عبر سطح كوكبي خشن. [69] تُستخدم طرق مونت كارلو أيضًا في نماذج المجموعة التي تشكل أساس التنبؤ بالطقس الحديث .

هندسة

تُستخدم طرق مونت كارلو على نطاق واسع في الهندسة لتحليل الحساسية والتحليل الاحتمالي الكمي في تصميم العمليات . تنشأ الحاجة من السلوك التفاعلي والخطي وغير الخطي لمحاكاة العمليات النموذجية. على سبيل المثال،

تغير المناخ والقوة الإشعاعية

تعتمد الهيئة الحكومية الدولية المعنية بتغير المناخ على أساليب مونت كارلو في تحليل دالة كثافة الاحتمال للقوة الإشعاعية . [73]

علم الأحياء الحاسوبي

تُستخدم طرق مونت كارلو في مجالات مختلفة من علم الأحياء الحسابي ، على سبيل المثال للاستدلال البايزي في علم النشوء والتطور ، أو لدراسة الأنظمة البيولوجية مثل الجينومات أو البروتينات، [74] أو الأغشية. [75] يمكن دراسة الأنظمة في الأطر ذات الحبيبات الخشنة أو الأطر الأولية اعتمادًا على الدقة المطلوبة. تسمح المحاكاة الحاسوبية بمراقبة البيئة المحلية لجزيء معين لمعرفة ما إذا كان يحدث تفاعل كيميائي على سبيل المثال. في الحالات التي لا يكون فيها من الممكن إجراء تجربة فيزيائية، يمكن إجراء تجارب فكرية (على سبيل المثال: كسر الروابط، أو إدخال شوائب في مواقع محددة، أو تغيير البنية المحلية/العالمية، أو إدخال مجالات خارجية).

رسومات الحاسوب

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

الاحصاء التطبيقي

تم وضع معايير تجارب مونت كارلو في الإحصاء بواسطة ساويلوفسكي. [76] في الإحصاء التطبيقي، يمكن استخدام طرق مونت كارلو لأربعة أغراض على الأقل:

  1. لمقارنة الإحصاءات المتنافسة للعينات الصغيرة في ظل ظروف بيانات واقعية. على الرغم من أنه يمكن حساب خطأ النوع الأول وخصائص القوة للإحصاءات للبيانات المأخوذة من التوزيعات النظرية الكلاسيكية ( على سبيل المثال ، المنحنى الطبيعي ، وتوزيع كوشي ) للظروف المقاربة ( أي حجم العينة اللانهائي وتأثير المعالجة الصغير للغاية)، فإن البيانات الحقيقية غالبًا لا تحتوي على مثل هذه التوزيعات. [77]
  2. توفير تنفيذات لاختبارات الفرضيات التي تكون أكثر كفاءة من الاختبارات الدقيقة مثل اختبارات التباديل (والتي غالبًا ما يكون من المستحيل حسابها) بينما تكون أكثر دقة من القيم الحرجة للتوزيعات المقاربة .
  3. لتوفير عينة عشوائية من التوزيع الخلفي في الاستدلال البايزي . تقوم هذه العينة بعد ذلك بتقريب وتلخيص جميع السمات الأساسية للتوزيع الخلفي.
  4. لتوفير تقديرات عشوائية فعالة لمصفوفة هيسيان لدالة الاحتمال اللوغاريتمي السالبة التي يمكن حساب متوسطها لتشكيل تقدير لمصفوفة معلومات فيشر . [78] [79]

إن أساليب مونت كارلو تشكل أيضاً حلاً وسطاً بين الاختبارات العشوائية التقريبية واختبارات التباديل. ويستند اختبار العشوائية التقريبية إلى مجموعة فرعية محددة من كل التباديل (وهو ما يستلزم قدراً هائلاً من العناية والتنظيم فيما يتصل بالتباديل التي تم أخذها في الاعتبار). أما أسلوب مونت كارلو فيستند إلى عدد محدد من التباديل المرسومة عشوائياً (مع استبدال الخسارة الطفيفة في الدقة إذا تم رسم التباديل مرتين ـ أو أكثر من مرة ـ بكفاءة عدم الاضطرار إلى تتبع التباديل التي تم اختيارها بالفعل).

الذكاء الاصطناعي للألعاب

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

تتضمن طريقة البحث الشجري لمونت كارلو (MCTS) أربع خطوات: [81]

  1. ابدأ من عقدة الجذر للشجرة، ثم حدد العقد الفرعية المثالية حتى الوصول إلى عقدة الورقة.
  2. قم بتوسيع عقدة الورقة واختر أحد أبنائها.
  3. قم بلعب لعبة محاكاة بدءًا من تلك العقدة.
  4. استخدم نتائج تلك اللعبة المحاكاة لتحديث العقدة وأسلافها.

إن التأثير الصافي، على مدار العديد من الألعاب المحاكاة، هو أن قيمة العقدة التي تمثل حركة سترتفع أو تنخفض، ومن المأمول أن يتوافق ذلك مع ما إذا كانت تلك العقدة تمثل حركة جيدة أم لا.

تم استخدام Monte Carlo Tree Search بنجاح لتشغيل ألعاب مثل Go ، [82] Tantrix ، [83] Battleship ، [84] Havannah ، [85] و Arimaa . [86]

التصميم والمرئيات

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

البحث والإنقاذ

يستخدم خفر السواحل الأمريكي طرق مونت كارلو ضمن برنامج النمذجة الحاسوبية SAROPS من أجل حساب المواقع المحتملة للسفن أثناء عمليات البحث والإنقاذ . يمكن لكل محاكاة توليد ما يصل إلى عشرة آلاف نقطة بيانات موزعة عشوائيًا بناءً على المتغيرات المقدمة. [88] ثم يتم إنشاء أنماط البحث بناءً على استقراء هذه البيانات من أجل تحسين احتمال الاحتواء (POC) واحتمال الكشف (POD)، والتي تساوي معًا احتمال النجاح الإجمالي (POS). في النهاية، يعمل هذا كتطبيق عملي لتوزيع الاحتمالات من أجل توفير أسرع وأكثر الطرق ملاءمة للإنقاذ، مما ينقذ الأرواح والموارد. [89]

المالية والأعمال

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

غالبًا ما تُستخدم طرق مونت كارلو في التمويل لتقييم الاستثمارات في المشاريع على مستوى وحدة الأعمال أو الشركة، أو التقييمات المالية الأخرى. يمكن استخدامها لنمذجة جداول المشروع ، حيث تجمع المحاكاة التقديرات لأسوأ الحالات وأفضل الحالات والأكثر احتمالية لفترات كل مهمة لتحديد النتائج للمشروع الإجمالي. [90] تُستخدم طرق مونت كارلو أيضًا في تسعير الخيارات وتحليل مخاطر التخلف عن السداد. [91] [92] بالإضافة إلى ذلك، يمكن استخدامها لتقدير التأثير المالي للتدخلات الطبية. [93]

قانون

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

علم المكتبات

كما تم استخدام نهج مونت كارلو لمحاكاة عدد منشورات الكتب بناءً على نوع الكتاب في ماليزيا. استخدمت محاكاة مونت كارلو بيانات منشورات الكتب الوطنية المنشورة سابقًا وسعر الكتاب وفقًا لنوع الكتاب في السوق المحلية. تم استخدام نتائج مونت كارلو لتحديد نوع الكتاب الذي يحبه الماليزيون وتم استخدامها لمقارنة منشورات الكتب بين ماليزيا واليابان . [ 95]

آخر

يكتب نسيم نيكولاس طالب عن مولدات مونت كارلو في كتابه الصادر عام 2001 بعنوان "مخدوعون بالعشوائية" كمثال حقيقي لاختبار تورينج العكسي : يمكن اعتبار الإنسان غير ذكي إذا لم يكن من الممكن التمييز بين كتابته وبين الكتابة المولدة.

الاستخدام في الرياضيات

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

اندماج

تعمل عملية تكامل مونت كارلو عن طريق مقارنة النقاط العشوائية بقيمة الدالة.
يتم تقليل الأخطاء بعامل .

تعمل خوارزميات التكامل العددي الحتمي بشكل جيد في عدد صغير من الأبعاد، ولكنها تواجه مشكلتين عندما تحتوي الدوال على العديد من المتغيرات. أولاً، يزداد عدد تقييمات الدوال المطلوبة بسرعة مع عدد الأبعاد. على سبيل المثال، إذا وفرت 10 تقييمات دقة كافية في بُعد واحد، فستكون هناك حاجة إلى 10 100 نقطة لـ 100 بُعد - وهو عدد كبير جدًا بحيث لا يمكن حسابه. وهذا ما يسمى بلعنة الأبعاد . ثانيًا، قد تكون حدود المنطقة متعددة الأبعاد معقدة للغاية، لذلك قد لا يكون من الممكن تقليل المشكلة إلى تكامل متكرر . [96] 100 بُعد ليس أمرًا غير عادي بأي حال من الأحوال، لأنه في العديد من المشكلات الفيزيائية، يكون "البعد" معادلًا لدرجة الحرية .

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

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

يستخدم نهج مماثل، وهو طريقة مونت كارلو شبه ، تسلسلات ذات تباين منخفض . تملأ هذه التسلسلات المنطقة بشكل أفضل وتأخذ عينات من النقاط الأكثر أهمية بشكل أكثر تكرارًا، وبالتالي يمكن لطرق مونت كارلو شبه أن تتقارب غالبًا مع التكامل بشكل أسرع.

فئة أخرى من الطرق لأخذ عينات من النقاط في حجم ما هي محاكاة المشي العشوائي عليها ( سلسلة ماركوف مونت كارلو ). تتضمن هذه الطرق خوارزمية متروبوليس-هاستينجز ، وعينات جيبس ، وخوارزمية وانج ولانداو ، ومنهجيات MCMC من النوع المتفاعل مثل أجهزة أخذ العينات المتسلسلة مونت كارلو . [99]

المحاكاة والتحسين

هناك تطبيق آخر قوي وشائع جدًا للأرقام العشوائية في المحاكاة العددية وهو التحسين العددي . تكمن المشكلة في تقليل (أو تعظيم) وظائف بعض المتجهات التي غالبًا ما يكون لها أبعاد عديدة. يمكن صياغة العديد من المشكلات بهذه الطريقة: على سبيل المثال، يمكن اعتبار برنامج شطرنج الكمبيوتر وكأنه يحاول العثور على مجموعة من 10 حركات على سبيل المثال والتي تنتج أفضل دالة تقييم في النهاية. في مشكلة بائع السفر ، الهدف هو تقليل المسافة المقطوعة. هناك أيضًا تطبيقات لتصميم الهندسة، مثل تحسين التصميم متعدد التخصصات . تم تطبيقه مع نماذج شبه أحادية البعد لحل مشاكل ديناميكيات الجسيمات من خلال استكشاف مساحة التكوين الكبيرة بكفاءة. المرجع [100] هو مراجعة شاملة للعديد من القضايا المتعلقة بالمحاكاة والتحسين.

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

المشاكل العكسية

تؤدي الصياغة الاحتمالية للمسائل العكسية إلى تعريف توزيع الاحتمالات في فضاء النموذج. يجمع توزيع الاحتمالات هذا بين المعلومات السابقة والمعلومات الجديدة التي تم الحصول عليها من خلال قياس بعض المعلمات القابلة للملاحظة (البيانات). وكما هو الحال في الحالة العامة، فإن النظرية التي تربط البيانات بمعلمات النموذج غير خطية، وقد لا يكون من السهل وصف الاحتمال اللاحق في فضاء النموذج (قد يكون متعدد الوسائط، وقد لا يتم تعريف بعض اللحظات، وما إلى ذلك).

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

يمكن تعميم طريقة أخذ العينات ذات الأهمية الأكثر شهرة، وهي خوارزمية متروبوليس، وهذا يوفر طريقة تسمح بتحليل المشكلات العكسية (التي قد تكون غير خطية بدرجة كبيرة) مع معلومات وبيانات مسبقة معقدة مع توزيع ضوضاء عشوائي. [101] [102]

فلسفة

أجرى ماكراكين عرضًا شعبيًا لطريقة مونت كارلو. [103] وناقش إليشاكوف [104] وجرون-يانوف وويريتش الفلسفة العامة للطريقة. [105]

انظر أيضا

مراجع

الاستشهادات

  1. ^ كالوس وويتلوك 2008.
  2. ^ كروس، دي بي؛ بريريتون، تي؛ تايمري، تي؛ بوتيف، زي (2014). "لماذا تعتبر طريقة مونت كارلو مهمة للغاية اليوم". مجلة WIREs Comput Stat . 6 (6): 386–392. doi :10.1002/wics.1314. S2CID  18521840.
  3. ^ هوبارد، دوغلاس؛ صامويلسون، دوغلاس أ. (أكتوبر 2009). "النمذجة بدون قياسات". OR/MS Today : 28–33.
  4. ^ Metropolis, Nicholas; Rosenbluth, Arianna W.; Rosenbluth, Marshall N.; Teller, Augusta H.; Teller, Edward (June 1, 1953). "معادلة حسابات الحالة بواسطة آلات الحوسبة السريعة". مجلة الفيزياء الكيميائية . 21 (6): 1087–1092. Bibcode :1953JChPh..21.1087M. doi :10.1063/1.1699114. ISSN  0021-9606. OSTI  4390578. S2CID  1046577.
  5. ^ هاستينجز، دبليو كيه (1 أبريل 1970). "طرق أخذ العينات مونت كارلو باستخدام سلاسل ماركوف وتطبيقاتها". بيومتريكا . 57 (1): 97-109. رمز Bibcode :1970Bimka..57...97H. doi :10.1093/biomet/57.1.97. ISSN  0006-3444. S2CID  21204149.
  6. ^ ليو، جون إس.؛ ليانج، فامينج؛ وونج، وينج هونج (1 مارس 2000). "طريقة المحاولة المتعددة والتحسين المحلي في أخذ العينات الحضرية". مجلة الجمعية الإحصائية الأمريكية . 95 (449): 121-134. doi :10.1080/01621459.2000.10473908. ISSN  0162-1459. S2CID  123468109.
  7. ^ Spall, JC (2003). "Estimation via Markov Chain Monte Carlo". مجلة أنظمة التحكم IEEE . 23 (2): 34–45. doi :10.1109/MCS.2003.1188770.
  8. ^ هيل، ستايسي د.؛ سبال، جيمس س. (2019). "الثبات والتقارب لخوارزمية متروبوليس-هاستينجز: رؤى حول الجوانب النظرية". مجلة أنظمة التحكم IEEE . 39 : 56-67. doi :10.1109/MCS.2018.2876959. S2CID  58672766.
  9. ^ ab Kolokoltsov, Vassili (2010). Nonlinear Markov processes . Cambridge University Press . p. 375.
  10. ^ abc Del Moral, Pierre (2013). Mean field simulator for Monte Carlo integration. Chapman & Hall/ CRC Press . p. 626. Monographs on Statistics & Applied Probability
  11. ^ ديل مورال ، ص. دوسيه، أ.؛ جسرا، أ. (2006). “أخذ عينات مونت كارلو المتسلسلة”. مجلة الجمعية الإحصائية الملكية، السلسلة ب . 68 (3): 411-436. أرخايف : cond-mat/0212648 . دوى :10.1111/j.1467-9868.2006.00553.x. S2CID  12074789.
  12. ^ abcd Owen, Art B. (2013). نظرية مونت كارلو، الأساليب والأمثلة (PDF) . العمل قيد التنفيذ. ص 15-36.
  13. ^ ab Driels, Morris R.; Shin, Young S. (April 2004). "تحديد عدد التكرارات لمحاكاة مونت كارلو لفعالية الأسلحة". تقرير فني لمدرسة الدراسات العليا البحرية (مارس 2003 - مارس 2004): 10-11.
  14. ^ Shonkwiler, RW; Mendivil, F. (2009). Explorations in Monte Carlo Methods . Springer.
  15. ^ Atanassova, E.; Gurov, T.; Karaivanova, A.; Ivanovska, S.; Durchova, M.; Dimitrov, D. (2016). "حول أساليب التوازي لهندسة Intel MIC". وقائع مؤتمر AIP . 1773 (1): 070001. Bibcode :2016AIPC.1773g0001A. doi :10.1063/1.4964983.
  16. ^ Cunha Jr, A.; Nazir, R.; Sampaio, R.; Lopes, H.; Breitman, K. (2014). "قياس عدم اليقين من خلال طريقة مونت كارلو في بيئة الحوسبة السحابية". Computer Physics Communications . 185 (5): 1355–1363. arXiv : 2105.09512 . Bibcode :2014CoPhC.185.1355C. doi :10.1016/j.cpc.2014.01.006. S2CID  32376269.
  17. ^ Wei, J.; Kruis, FE (2013). "طريقة مونت كارلو المتوازية المعتمدة على وحدة معالجة الرسوميات لتخثر الجسيمات باستخدام استراتيجية القبول والرفض". علم الهندسة الكيميائية . 104 : 451–459. Bibcode :2013ChEnS.104..451W. doi :10.1016/j.ces.2013.08.008.
  18. ^ لين، ي.؛ وانج، ف.؛ ليو، ب. (2018). "مولدات الأرقام العشوائية لمحاكاة مونت كارلو المتوازية واسعة النطاق على FPGA". مجلة الفيزياء الحاسوبية . 360 : 93-103. رمز Bibcode : 2018JCoPh.360...93L. doi : 10.1016/j.jcp.2018.01.029.
  19. ^ abc Metropolis 1987.
  20. ^ إيكهارت 1987.
  21. ^ أب مازدراكوف، بينوف وفالكانوف 2018، ص. 250.
  22. ^ هايغ، توماس؛ بريستلي، مارك؛ روب، كريسبين (2014). "رهانات لوس ألاموس على إينياك: محاكاة مونت كارلو النووية، 1947-1948". حوليات معهد مهندسي الكهرباء والإلكترونيات لتاريخ الحوسبة . 36 (3): 42-63. doi :10.1109/MAHC.2014.40. S2CID  17470931.
  23. ^ ماكين، هنري ب. (1967). "انتشار الفوضى لفئة من المعادلات المكافئة غير الخطية". سلسلة محاضرات في المعادلات التفاضلية، الجامعة الكاثوليكية . 7 : 41-57.
  24. ^ McKean, Henry P. (1966). "فئة من عمليات ماركوف المرتبطة بالمعادلات المكافئة غير الخطية". Proc. Natl. Acad. Sci. USA . 56 (6): 1907–1911. Bibcode : 1966PNAS ...56.1907M. doi : 10.1073/pnas.56.6.1907 . PMC 220210. PMID  16591437. 
  25. ^ هيرمان، كان؛ ثيودور، هاريس إي. (1951). "تقدير انتقال الجسيمات عن طريق أخذ العينات العشوائية" (PDF) . مجلة المكتبات الوطنية. سلسلة الرياضيات التطبيقية . 12 : 27-30.
  26. ^ تورينج، آلان م. (1950). "آلات الحوسبة والذكاء". العقل . LIX (238): 433-460. doi :10.1093/mind/LIX.236.433.
  27. ^ باريتشيلي ، نيلز آل (1954). "نماذج رقمية من عمليات التطور". الطرق : 45-68.
  28. ^ باريسيلي، نيلز آل (1957). "عمليات التطور التكافلي التي تحققت بالطرق الاصطناعية". Methodos : 143–182.
  29. ^ ab Del Moral, Pierre (2004). Feynman–Kac formulae. Genealogical and interacting particle approximations. Probability and Its Applications. Springer. p. 575. ISBN 9780387202686. السلسلة: الاحتمالات والتطبيقات
  30. ^ ab Del Moral, P.; Miclo, L. (2000). "تقريبات أنظمة الجسيمات المتفرعة والمتفاعلة لصيغ فاينمان-كاتش مع التطبيقات على الترشيح غير الخطي". Séminaire de Probabilités XXXIV. Lecture Notes in Mathematics. المجلد 1729. برلين: سبرينغر. ص. 1-145. doi :10.1007/BFb0103798. ISBN 978-3-540-67314-9. السيد  1768060.
  31. ^ ديل مورال، بيير؛ ميكلو، لوران (2000). "تقريب نظام جسيمات موران لصيغ فاينمان-كاتش". العمليات العشوائية وتطبيقاتها . 86 (2): 193-216. doi : 10.1016/S0304-4149(99)00094-0 .
  32. ^ ديل مورال، بيير (2003). "تقريبات الجسيمات لأسس ليابونوف المتصلة بمشغلات شرودنجر وشبه مجموعات فاينمان-كاتش". مجلة الاحتمالات والإحصاءات الأوروبية . 7 : 171-208. doi : 10.1051/ps:2003001 .
  33. ^ Assaraf, Roland; Caffarel, Michel; Khelif, Anatole (2000). "Diffusion Monte Carlo Methods with a fixed number of walkers" (PDF) . Phys. Rev. E. 61 ( 4): 4566–4575. Bibcode :2000PhRvE..61.4566A. doi :10.1103/physreve.61.4566. PMID  11088257. مؤرشف من الأصل (PDF) في 7 نوفمبر 2014.
  34. ^ كافاريل، ميشيل؛ سيبرلي، ديفيد؛ كالوس، مالفين (1993). "تعليق على حساب تكامل مسار فاينمان-كاتش لطاقات الحالة الأرضية للذرات". Phys. Rev. Lett . 71 (13): 2159. Bibcode :1993PhRvL..71.2159C. doi :10.1103/physrevlett.71.2159. PMID  10054598.
  35. ^ ab Hetherington, Jack H. (1984). "ملاحظات حول التكرار الإحصائي للمصفوفات". Phys. Rev. A. 30 ( 2713): 2713–2719. Bibcode :1984PhRvA..30.2713H. doi :10.1103/PhysRevA.30.2713.
  36. ^ فيرمي، إنريكي؛ ريتشماير، روبرت د. (1948). "ملاحظة حول إجراء التعداد في حسابات مونت كارلو" (PDF) . لام . 805 (أ). تقرير غير سري أرشيف لوس ألاموس
  37. ^ ab Rosenbluth, Marshall N.; Rosenbluth, Arianna W. (1955). "حسابات مونت كارلو لمتوسط ​​امتداد السلاسل الجزيئية الكبيرة". J. Chem. Phys . 23 (2): 356–359. Bibcode :1955JChPh..23..356R. doi : 10.1063/1.1741967 . S2CID  89611599.
  38. ^ جوردون، نيوجيرسي؛ سالموند، دي جي؛ سميث، إيه إف إم (أبريل 1993). "نهج جديد لتقدير الحالة البايزية غير الخطية/غير الغاوسية". وقائع معهد الهندسة الكهربائية والإلكترونية ف - الرادار ومعالجة الإشارات . 140 (2): 107-113. doi :10.1049/ip-f-2.1993.0015. ISSN  0956-375X. S2CID  12644877.
  39. ^ كيتاغاوا، ج. (1996). "مرشح مونت كارلو وأداة التنعيم لنماذج الفضاء غير الخطية غير الغاوسية". مجلة الإحصاءات الحسابية والرسومية . 5 (1): 1-25. doi :10.2307/1390750. JSTOR  1390750.
  40. ^ ab Del Moral, Pierre (1996). "Non Linear Filtering: Interacting Particle Solution" (PDF) . Markov Processes and Related Fields . 2 (4): 555–580. مؤرشف من الأصل (PDF) في 4 مارس 2016. تم الاسترجاع في 11 يونيو 2015 .
  41. ^ كارفاليو، هيميلكون؛ ديل مورال، بيير؛ مونين، أندريه؛ سالوت، جيرارد (يوليو 1997). "الترشيح غير الخطي الأمثل في تكامل GPS/INS" (PDF) . معاملات معهد مهندسي الكهرباء والإلكترونيات في مجال الفضاء والأنظمة الإلكترونية . 33 (3): 835-850. رمز Bibcode :1997ITAES..33..835C. doi :10.1109/7.599254. S2CID  27966240. مؤرشف من الأصل (PDF) في 10 نوفمبر 2022. تم الاسترجاع في 11 يونيو 2015 .
  42. ^ ب. ديل مورال، ج. ريجال، و ج. سالوت. “التقدير والتحكم الأمثل غير الخطي: إطار موحد لحلول الجسيمات”. LAAS-CNRS، تولوز، تقرير البحث رقم. 91137، عقد DRET-DIGILOG- LAAS/CNRS، أبريل (1991).
  43. ^ P. Del Moral, G. Rigal, and G. Salut. "Nonlinear and non Gaussian particle candidates applied to inertial platform repositioning." LAAS-CNRS, Toulouse, Research Report no. 92207, STCAN/DIGILOG-LAAS/CNRS Convention STCAN no. A.91.77.013, (94p.) September (1991).
  44. ^ P. Del Moral, G. Rigal, and G. Salut. "Estimation and nonlinear optimal control: Particle resolution in filtering and estimation: Experimental results". Convention DRET no. 89.34.553.00.470.75.01, Research report no.2 (54p.), January (1992).
  45. ^ P. Del Moral, G. Rigal, and G. Salut. "Estimation and nonlinear optimal control: Particle resolution in filtering and estimation: Theoretical results". Convention DRET no. 89.34.553.00.470.75.01, Research report no.3 (123p.), October (1992).
  46. ^ P. Del Moral, J.-Ch. Noyer, G. Rigal, and G. Salut. "مرشحات الجسيمات في معالجة إشارات الرادار: الكشف والتقدير والتعرف على الأهداف الجوية". LAAS-CNRS، تولوز، تقرير بحثي رقم 92495، ديسمبر (1992).
  47. ^ P. Del Moral, G. Rigal, and G. Salut. "Estimation and nonlinear optimal control: Particle resolution in filtering and estimation". دراسات حول: الترشيح والتحكم الأمثل وتقدير أقصى احتمال. اتفاقية DRET رقم 89.34.553.00.470.75.01. تقرير بحثي رقم 4 (210 صفحة)، يناير (1993).
  48. ^ ديل مورال ، بيير (1998). “قياس العمليات القيمة وأنظمة الجسيمات المتفاعلة. التطبيق على مشاكل التصفية غير الخطية”. حوليات الاحتمالية التطبيقية . 8 (2) (منشورات مختبر الإحصاء والاحتمالات، 96-15 (1996) طبعة): 438-495. سيتيسيركس 10.1.1.55.5257 . دوى :10.1214/aoap/1028903535. 
  49. ^ كريسان، دان؛ جينز، جيسيكا؛ ليونز، تيري (1998). "تقارب طريقة الجسيمات المتفرعة لحل زاكاي". مجلة سيام للرياضيات التطبيقية . 58 (5): 1568-1590. doi :10.1137/s0036139996307371. S2CID  39982562.
  50. ^ كريسان، دان؛ ليونز، تيري (1997). "الترشيح غير الخطي والعمليات ذات القيمة المقاسة". نظرية الاحتمالات والمجالات ذات الصلة . 109 (2): 217-244. doi : 10.1007/s004400050131 . S2CID  119809371.
  51. ^ كريسان، دان؛ ليونز، تيري (1999). "تقريب جسيم لحل معادلة كوشنر-ستراتونوفيتش". نظرية الاحتمالات والمجالات ذات الصلة . 115 (4): 549-578. doi : 10.1007/s004400050249 . S2CID  117725141.
  52. ^ كريسان، دان؛ ديل مورال، بيير؛ ليونز، تيري (1999). "الترشيح المنفصل باستخدام أنظمة الجسيمات المتفرعة والمتفاعلة" (PDF) . عمليات ماركوف والحقول ذات الصلة . 5 (3): 293-318.
  53. ^ ديل مورال، بيير؛ جوينيت، أليس (1999). "حول استقرار العمليات ذات القيمة المقاسة مع التطبيقات على الترشيح". CR Acad. Sci. Paris . 39 (1): 429–434.
  54. ^ ديل مورال، بيير؛ جوينيت، أليس (2001). "حول استقرار العمليات المتفاعلة مع التطبيقات على الترشيح والخوارزميات الجينية". حوليات معهد هنري بوانكاريه . 37 (2): 155-194. رمز Bibcode :2001AIHPB..37..155D. doi :10.1016/s0246-0203(00)01064-5.
  55. ^ ريبلي 1987
  56. ^ بواسطة ساويلوفسكي 2003
  57. ^ كالوس وويتلوك 2008
  58. ^ Shojaeefard, MH; Khalkhali, A.; Yarmohammadisatri, Sadegh (2017). "طريقة تحليل حساسية فعّالة للهندسة المعدلة لنظام تعليق ماكفرسون استنادًا إلى معامل ارتباط بيرسون". ديناميكيات نظام المركبات . 55 (6): 827–852. رمز Bibcode :2017VSD....55..827S. doi :10.1080/00423114.2017.1283046. S2CID  114260173.
  59. ^ دافنبورت 1992
  60. ^ Route, Matthew (10 أغسطس 2017). "Radio-flaring Ultracool Dwarf Population Synthesis". مجلة الفيزياء الفلكية . 845 (1): 66. arXiv : 1707.02212 . Bibcode :2017ApJ...845...66R. doi : 10.3847/1538-4357/aa7ede . S2CID  118895524.
  61. ^ فوزي 2008، ص 13.
  62. ^ فوزي 2008، ص 16.
  63. ^ جيا، شون؛ زيجينهاين، بيتر؛ جيانج، ستيف ب (2014). "الحوسبة عالية الأداء القائمة على وحدة معالجة الرسوميات للعلاج الإشعاعي". الفيزياء في الطب والأحياء . 59 (4): R151–R182. رمز Bibcode :2014PMB....59R.151J. doi :10.1088/0031-9155/59/4/R151. PMC 4003902. PMID  24486639 . 
  64. ^ هيل، ر.؛ هيلي، ب.؛ هولواي، ل.؛ كونسيك، ز.؛ ثويتس، د.؛ بالدوك، س. (مارس 2014). "التطورات في قياس جرعات شعاع الأشعة السينية بالكيلوفولت". الفيزياء في الطب والأحياء . 59 (6): R183–R231. رمز Bibcode :2014PMB....59R.183H. doi :10.1088/0031-9155/59/6/R183. PMID  24584183. S2CID  18082594.
  65. ^ Rogers, DWO (2006). "خمسون عامًا من محاكاة مونت كارلو للفيزياء الطبية". الفيزياء في الطب والأحياء . 51 (13): R287–R301. Bibcode :2006PMB....51R.287R. doi :10.1088/0031-9155/51/13/R17. PMID  16790908. S2CID  12066026.
  66. ^ بيورل 2009
  67. ^ مولر، دبليو؛ إيكشتاين، دبليو (1 مارس 1984). "تريديين - كود محاكاة تريم يتضمن تغييرات في التركيب الديناميكي". الأجهزة والطرق النووية في أبحاث الفيزياء القسم ب: تفاعلات الحزمة مع المواد والذرات . 2 (1): 814-818. رمز Bibcode :1984NIMPB...2..814M. doi :10.1016/0168-583X(84)90321-5.
  68. ^ ماكجيليفراي ودود 1982
  69. ^ جولدن 1979
  70. ^ GA Bird، ديناميكيات الغاز الجزيئية، كلارندون، أكسفورد (1976)
  71. ^ Dietrich, S.; Boyd, I. (1996). "A Scalar optimized parallel implementation of the DSMC technique". مجلة الفيزياء الحاسوبية . 126 (2): 328–42. Bibcode :1996JCoPh.126..328D. doi : 10.1006/jcph.1996.0141 .
  72. ^ تشن، شانغ ينغ؛ هسو، كيو تشين؛ فان، تشيا مينغ (15 مارس 2021). "تحسين طريقة الفرق المحدودة المعممة لنمذجة التدفق تحت السطحي العشوائي". مجلة الفيزياء الحاسوبية . 429 : 110002. رمز Bibcode : 2021JCoPh.42910002C. doi : 10.1016/J.JCP.2020.110002. S2CID  228828681.
  73. ^ تغير المناخ 2013 الأساس العلمي الفيزيائي (PDF) . مطبعة جامعة كامبريدج . 2013. ص. 697. ISBN 978-1-107-66182-0تم الاسترجاع بتاريخ 6 يوليو 2023 .
  74. ^ أوجيدا وآخرون. 2009.
  75. ^ ميليك وسكولينك 1993.
  76. ^ كاسي؛ سميث (2014). "محاكاة الثقة لمؤشر إليسون-جلاسر". مجلة الاقتصاد الحضري . 81 : 93. doi :10.1016/j.jue.2014.02.005.
  77. ^ ساويلوفسكي وفاهوم 2003
  78. ^ Spall, James C. (2005). "حساب مونت كارلو لمصفوفة معلومات فيشر في الإعدادات غير القياسية". مجلة الإحصاءات الحسابية والرسومية . 14 (4): 889–909. CiteSeerX 10.1.1.142.738 . doi :10.1198/106186005X78800. S2CID  16090098. 
  79. ^ Das, Sonjoy; Spall, James C.; Ghanem, Roger (2010). "حساب مونت كارلو الفعّال لمصفوفة معلومات فيشر باستخدام المعلومات السابقة". إحصاءات حسابية وتحليل البيانات . 54 (2): 272–289. doi :10.1016/j.csda.2009.09.018.
  80. ^ تشاسلوت ، غيوم. باكس، ساندر؛ سزيتا، إستفان؛ سبرونك، بيتر. "بحث شجرة مونت كارلو: إطار عمل جديد للعبة الذكاء الاصطناعي" (PDF) . Sander.landofsand.com . تم الاسترجاع 28 أكتوبر، 2017 .
  81. ^ "Monte Carlo Tree Search - About". مؤرشف من الأصل في 29 نوفمبر 2015. تم الاسترجاع في 15 مايو 2013 .
  82. ^ تشاسلوت، غيوم إم جي -ب؛ ويناندز، مارك جلالة؛ فان دن هيريك، هـ. جاب (2008). “بحث موازي لشجرة مونت كارلو”. أجهزة الكمبيوتر والألعاب . ملاحظات محاضرة في علوم الكمبيوتر. المجلد. 5131. ص 60-71. سيتيسيركس 10.1.1.159.4373 . دوى :10.1007/978-3-540-87608-3_6. رقم ISBN  978-3-540-87607-6.
  83. ^ برونز، بيت. البحث عن شجرة مونت كارلو في لعبة تانتركس: التقرير النهائي لـ Cosc490 (PDF) (تقرير).
  84. ^ سيلفر، ديفيد؛ فينيس، جويل. "تخطيط مونت كارلو في مشاريع POMDP الكبيرة" (PDF) . 0.cs.ucl.ac.uk. مؤرشف من الأصل (PDF) في 18 يوليو 2016. تم الاسترجاع في 28 أكتوبر 2017 .
  85. ^ Lorentz, Richard J. (2011). "تحسين البحث الشجري باستخدام طريقة مونت كارلو في هافانا". Computers and Games . Lecture Notes in Computer Science. المجلد 6515. ص 105-115. Bibcode :2011LNCS.6515..105L. doi :10.1007/978-3-642-17928-0_10. ISBN 978-3-642-17927-3.
  86. ^ جاكل، توماس. "تحدي أريما - دراسة مقارنة لطرق MCTS مقابل طرق ألفا-بيتا" (PDF) . Arimaa.com . تم الاسترجاع في 28 أكتوبر 2017 .
  87. ^ Szirmay-Kalos 2008.
  88. ^ "كيف يستخدم خفر السواحل التحليلات للبحث عن المفقودين في البحر". Dice Insights . 3 يناير 2014.
  89. ^ Stone, Lawrence D.; Kratzke, Thomas M.; Frost, John R. "Search Modeling and Optimization in USCG's Search and Rescue Optimal Planning System (SAROPS)" (PDF) . Ifremer.fr . تم الاسترجاع في 28 أكتوبر 2017 .
  90. ^ "محاكاة مخاطر المشروع (BETA)". risk.octigo.pl . تم الاسترجاع في 21 مايو 2024 .
  91. ^ كارمونا، رينيه؛ ديل مورال، بيير؛ هو، بينج؛ أوجان، ناديا (2012). "مقدمة إلى أساليب الجسيمات مع التطبيقات المالية". في كارمونا، رينيه أ.؛ مورال، بيير ديل؛ هو، بينج؛ وآخرون (المحررون). الأساليب العددية في التمويل . وقائع سبرينغر في الرياضيات. المجلد 12. سبرينغر برلين هايدلبرغ. ص. 3-49. CiteSeerX 10.1.1.359.7957 . doi :10.1007/978-3-642-25746-9_1. ISBN  978-3-642-25745-2.
  92. ^ كروس، دي بي؛ تايمري، تي؛ بوتيف، زي (2011). دليل طرق مونت كارلو . جون وايلي وأولاده.
  93. ^ أريناس، دانيال جيه؛ ليت، لانير أ؛ كلوزاريتز، هيذر؛ تيتلمان، آن م. (2017). "نهج محاكاة مونت كارلو لتقدير التأثير الصحي والاقتصادي للتدخلات المقدمة في عيادة يديرها الطلاب". PLOS ONE . 12 (12): e0189718. رمز Bibcode : 2017PLoSO..1289718A. doi : 10.1371/journal.pone.0189718 . PMC 5746244. PMID  29284026 . 
  94. ^ إلوارت، ليز؛ إيمرسون، نينا؛ إندرز، كريستينا؛ فوميا، داني؛ مورفي، كيفن (ديسمبر 2006). "زيادة فرص الحصول على أوامر تقييدية لضحايا العنف المنزلي من ذوي الدخل المنخفض: تحليل التكلفة والفائدة لبرنامج المنح المقترحة للعنف المنزلي" (PDF) . نقابة المحامين في ولاية ويسكونسن . مؤرشف من الأصل (PDF) في 6 نوفمبر 2018. تم الاسترجاع في 12 ديسمبر 2016 .
  95. ^ دحلان، هادي أكبر (29 أكتوبر 2021). "Perbandingan Penerbitan dan Harga Buku Mengikut Genre di Malaysia dan Jepun Menggunakan Data Akses Terbuka dan Simulasi Monte Carlo" (PDF) . كاجيان ماليزيا . 39 (2): 179-202. دوى :10.21315/km2021.39.2.8. S2CID  240435973.
  96. ^ ab Press et al. 1996
  97. ^ MEZEI, M (31 ديسمبر 1986). "العينة المظلة التكيفية: التحديد الذاتي المتسق للتحيز غير بولتزمان". مجلة الفيزياء الحاسوبية . 68 (1): 237-248. رمز Bibcode :1987JCoPh..68..237M. doi :10.1016/0021-9991(87)90054-4.
  98. ^ بارتلز، كريستيان؛ كاربلس، مارتن (31 ديسمبر 1997). "توزيعات الاحتمالات للأنظمة المعقدة: أخذ العينات المظلة التكيفية للطاقة الكامنة". مجلة الكيمياء الفيزيائية ب . 102 (5): 865-880. doi :10.1021/jp972280j.
  99. ^ ديل مورال ، بيير. دوسيه، أرنو؛ جسرا، أجاي (2006). “أخذ عينات مونت كارلو المتسلسلة”. مجلة الجمعية الإحصائية الملكية، السلسلة ب . 68 (3): 411-436. أرخايف : cond-mat/0212648 . دوى :10.1111/j.1467-9868.2006.00553.x. S2CID  12074789.
  100. ^ Spall, JC (2003)، مقدمة إلى البحث العشوائي والتحسين: التقدير والمحاكاة والتحكم ، وايلي، هوبوكين، نيوجيرسي. http://www.jhuapl.edu/ISSO
  101. ^ موسيجارد وتارانتولا 1995
  102. ^ تارانتولا 2005
  103. ^ McCracken, DD, (1955) طريقة مونت كارلو، ساينتفك أمريكان، 192(5)، ص 90-97
  104. ^ إليشاكوف، آي. (2003) ملاحظات حول فلسفة طريقة مونت كارلو، الميكانيكا التطبيقية الدولية، 39(7)، ص.753-762
  105. ^ Grüne-Yanoff, T., & Weirich, P. (2010). فلسفة ونظرية المحاكاة: مراجعة، المحاكاة والألعاب، 41(1)، ص 20-50

مصادر

  • أندرسون، هربرت ل. (1986). "ميتروبوليس، مونت كارلو والمجنون" (PDF) . لوس ألاموس ساينس . 14 : 96-108.
  • بينوف، دوبريان م. (2016). "مشروع مانهاتن، أول حاسوب إلكتروني وطريقة مونت كارلو". طرق وتطبيقات مونت كارلو . 22 (1): 73-79. doi :10.1515/mcma-2016-0102. S2CID  30198383.
  • باورل، ستيفان أ. (2009). "النمذجة متعددة المقاييس للمواد البوليمرية باستخدام منهجيات نظرية المجال: دراسة استقصائية حول التطورات الأخيرة". مجلة الكيمياء الرياضية . 46 (2): 363-426. doi :10.1007/s10910-008-9467-3. S2CID  117867762.
  • بيرج، بيرند أ. (2004). محاكاة سلسلة ماركوف لمونت كارلو وتحليلها الإحصائي (باستخدام كود فورتران المستند إلى الويب) . هاكنساك، نيوجيرسي: وورلد ساينتيفيك. رقم ISBN 978-981-238-935-0.
  • بايندر، كورت (1995). طريقة مونت كارلو في فيزياء المادة المكثفة . نيويورك: سبرينغر. ISBN 978-0-387-54369-7.
  • كافليش، ري (1998). طرق مونت كارلو وشبه مونت كارلو . اكتا نوميريكا. المجلد. 7. مطبعة جامعة كامبريدج. ص 1-49.
  • دافنبورت، جيه إتش (1992). "إعادة النظر في اختبار البدائية". أوراق من الندوة الدولية حول الحوسبة الرمزية والجبرية - ISSAC '92 . ص 123-129. CiteSeerX  10.1.1.43.9296 . doi :10.1145/143242.143290. ISBN 978-0-89791-489-5. S2CID  17322272.
  • دوسيت، أرنو؛ فريتاس، ناندو دي؛ جوردون، نيل (2001). أساليب مونت كارلو المتسلسلة في الممارسة العملية . نيويورك: سبرينغر. ISBN 978-0-387-95146-1.
  • إيكهاردت، روجر (1987). “ستان أولام، جون فون نيومان، وطريقة مونت كارلو” (PDF) . علوم لوس ألاموس (15): 131-137.
  • فيشمان، ج. س. (1995). مونت كارلو: المفاهيم والخوارزميات والتطبيقات . نيويورك: سبرينغر. رقم ISBN 978-0-387-94527-9.
  • C. Forastero و L. Zamora و D. Guirado و A. Lallena (2010). "أداة مونت كارلو لمحاكاة برامج فحص سرطان الثدي". Phys. Med. Biol . 55 (17): 5213–5229. Bibcode :2010PMB....55.5213F. doi :10.1088/0031-9155/55/17/021. PMID  20714045. S2CID  30021759.
  • جولدن، ليزلي م. (1979). "تأثير خشونة السطح على انتقال إشعاع الموجات الدقيقة عبر سطح كوكبي". إيكاروس . 38 (3): 451-455. رمز Bibcode :1979Icar...38..451G. doi :10.1016/0019-1035(79)90199-4.
  • جولد، هارفي؛ توبوشنيك، جان (1988). مقدمة إلى أساليب المحاكاة الحاسوبية، الجزء 2، تطبيقات على الأنظمة الفيزيائية . القراءة: أديسون ويسلي. رقم ISBN 978-0-201-16504-3.
  • جرينستيد، تشارلز؛ سنيل، ج. لوري (1997). مقدمة في الاحتمالات. الجمعية الرياضية الأمريكية . ص 10-11.
  • هامرسلي، جيه إم؛ هاندسكومب، دي سي (1975). طرق مونت كارلو . لندن: ميثيون. رقم ISBN 978-0-416-52340-9.
  • هارتمان، أيه كيه (2009). دليل عملي لمحاكاة الكمبيوتر. مجلة وورلد ساينتفيك. رقم ISBN 978-981-283-415-7. تم أرشفة النسخة الأصلية في 11 فبراير 2009.
  • هوبارد، دوغلاس (2007). كيفية قياس أي شيء: إيجاد قيمة الأشياء غير الملموسة في الأعمال . جون وايلي وأولاده . ص. 46. ISBN 9780470110126.
  • هوبارد، دوغلاس (2009). فشل إدارة المخاطر: لماذا هو معطل وكيفية إصلاحه . جون وايلي وأولاده .
  • كانيمان، د.؛ تفيرسكي، أ. (1982). الحكم في ظل عدم اليقين: الاستدلالات والتحيزات . مطبعة جامعة كامبريدج.
  • كالوس، مالفين H .؛ ويتلوك، باولا أ. (2008). طرق مونت كارلو . وايلي-VCH . رقم ISBN 978-3-527-40760-6.
  • كروسي، دي بي؛ تيمري، ت.؛ بوتيف، زي (2011). دليل أساليب مونت كارلو. نيويورك: جون وايلي وأولاده . ص. 772. ردمك 978-0-470-17793-8.
  • MacGillivray, HT; Dodd, RJ (1982). "محاكاة مونت كارلو لأنظمة المجرات". Astrophysics and Space Science . 86 (2): 419–435. doi :10.1007/BF00683346. S2CID  189849365.
  • ماكيون، ب. كيفين (1997). المحاكاة العشوائية في الفيزياء . نيويورك: سبرينغر. ISBN 978-981-3083-26-4.
  • متروبوليس، ن. (1987). “بداية طريقة مونت كارلو” (PDF) . علوم لوس ألاموس (عدد خاص لعام 1987 مخصص لستانيسلاف أولام): 125-130.
  • Metropolis, N .؛ Rosenbluth, Arianna W.؛ Rosenbluth, Marshall N.؛ Teller, Augusta H.؛ Teller, Edward (1953). "معادلة حسابات الحالة بواسطة آلات الحوسبة السريعة" . مجلة الفيزياء الكيميائية . 21 (6): 1087. رمز Bibcode :1953JChPh..21.1087M. doi :10.1063/1.1699114. OSTI  4390578. S2CID  1046577.
  • متروبوليس، نأولام، س. (1949). "طريقة مونت كارلو". مجلة الجمعية الإحصائية الأمريكية . 44 (247): 335-341. doi :10.1080/01621459.1949.10483310. JSTOR  2280232. PMID  18139350.
  • ميليك، م.؛ سكولنيك، ج. (يناير 1993). "إدخال سلاسل الببتيد في الأغشية الدهنية: نموذج ديناميكيات مونت كارلو خارج الشبكة". البروتينات . 15 (1): 10-25. doi :10.1002/prot.340150104. PMID  8451235. S2CID  7450512.
  • موسيجارد، كلاوس؛ تارانتولا، ألبرت (1995). "أخذ العينات من مونت كارلو للحلول للمشكلات العكسية" (PDF) . J. Geophys. Res . 100 (B7): 12431–12447. Bibcode :1995JGR...10012431M. doi :10.1029/94JB03097. مؤرشف من الأصل (PDF) في 10 مارس 2021. تم الاسترجاع في 1 نوفمبر 2017 .
  • أوجيدا، ب.؛ جارسيا، م.؛ لندنو، أ.؛ تشين، نيويورك (فبراير 2009). "محاكاة مونت كارلو للبروتينات في الأقفاص: تأثير الحبس على استقرار الحالات الوسيطة". مجلة بيوفيس . 96 (3): 1076-1082. رمز Bibcode : 2009BpJ....96.1076O. doi : 10.1529/biophysj.107.125369. PMC  2716574. PMID  18849410 .
  • إنت بانيس، ل.؛ دي نوكر، ل.؛ دي فليجر، ل.؛ تورفس، ر. (2001). "الاتجاهات وعدم اليقين في تأثيرات تلوث الهواء والتكاليف الخارجية لحركة مرور سيارات الركاب البلجيكية". المجلة الدولية لتصميم المركبات . 27 (1-4): 183-194. doi :10.1504/IJVD.2001.001963.
  • إنت بانيس، إل؛ ربل، أ. دي نوكر، L.؛ تورفس، ر. (2002). شتورم، P. (محرر). “الديزل أم البنزين؟ مقارنة بيئية يعوقها عدم اليقين “. معهد Mitteilungen لتصنيع الآلات والديناميكية الحرارية . Heft 81 Vol 1. الجامعة التقنية في غراتس النمسا: 48-54.
  • Press, William H.; Teukolsky, Saul A.; Vetterling, William T.; Flannery, Brian P. (1996) [1986]. وصفات رقمية في فورتران 77: فن الحوسبة العلمية . وصفات رقمية في فورتران. المجلد 1 (الطبعة الثانية). مطبعة جامعة كامبريدج . ISBN 978-0-521-43064-7.
  • ريبلي، بي دي (1987). المحاكاة العشوائية . وايلي وأولاده .
  • روبرت، C .؛ كاسيلا، ج. (2004). الأساليب الإحصائية لمونت كارلو (الطبعة الثانية). نيويورك: سبرينغر. رقم ISBN 978-0-387-21239-5.
  • روبنشتاين، آر واي؛ كروس، دي بي (2007). المحاكاة وطريقة مونت كارلو (الطبعة الثانية). نيويورك: جون وايلي وأولاده. رقم ISBN 978-0-470-17793-8.
  • سافيدس، سافاكاس سي. (1994). "تحليل المخاطر في تقييم الاستثمار" (PDF) . مجلة تقييم المشروع . 9 (1). doi :10.2139/ssrn.265905. S2CID  2809643.
  • Sawilowsky, Shlomo S.; Fahoome, Gail C. (2003). إحصائيات عبر محاكاة مونت كارلو باستخدام فورتران . روتشستر هيلز، ميشيغان: JMASM. ISBN 978-0-9740236-0-1.
  • Sawilowsky, Shlomo S. (2003). "هل تعتقد أن لديك أشياء تافهة؟". مجلة الأساليب الإحصائية التطبيقية الحديثة . 2 (1): 218-225. doi : 10.22237/jmasm/1051748460 .
  • Silver, David; Veness, Joel (2010). "Monte-Carlo Planning in Large POMDPs" (PDF) . في Lafferty, J.; Williams, CKI; Shawe-Taylor, J.; Zemel, RS; Culotta, A. (المحررون). Advances in Neural Information Processing Systems 23. Neural Information Processing Systems 2010. Neural Information Processing Systems Foundation. مؤرشف من الأصل (PDF) في 25 مايو 2012. تم الاسترجاع في 15 مارس 2011 .
  • زيرماي كالوس، لازلو (2008). طرق مونت كارلو في الإضاءة العالمية – تقديم صور واقعية مع التوزيع العشوائي . VDM Verlag Dr. Mueller eK ISBN 978-3-8364-7919-6.
  • تارانتولا، ألبرت (2005). نظرية المشكلة العكسية. فيلادلفيا: جمعية الرياضيات الصناعية والتطبيقية. رقم ISBN 978-0-89871-572-9.
  • فوزي، ديفيد (2008). تحليل المخاطر، دليل كمي (الطبعة الثالثة). جون وايلي وأولاده . رقم ISBN 9780470512845.
  • مازدراكوف، ميتودي؛ بينوف، دوبريان؛ فالكانوف، نيكولاي (2018). طريقة مونت كارلو. التطبيقات الهندسية. مطبعة ACMO الأكاديمية. رقم ISBN 978-619-90684-3-4.
  • الوسائط المتعلقة بطريقة مونت كارلو في ويكيميديا ​​كومنز
تم الاسترجاع من "https://en.wikipedia.org/w/index.php?title=طريقة_مونت_كارلو&oldid=1262423747"
Original text
Rate this translation
Your feedback will be used to help improve Google Translate