بهبود کارایی شبیه‌سازی مونت‌کارلو برای تولید داده‌های همبسته: کاربردی در تغییر اقلیم

نوع مقاله : مقاله پژوهشی

نویسندگان

1 استاد، گروه، مهندسی آب، دانشکده کشاورزی، دانشگاه فردوسی مشهد

2 دانشجوی مهندسی برق و کامپیوتر، دانشکده مهندسی، دانشگاه مانیتوبا، وینیپگ، کانادا

چکیده

از شبیه‌سازی مونت‌کارلو به‌خوبی برای تولید گروه داده‌هایی که از توزیعی مشخص پیروی می‌کنند استفاده می‌شود. با این‌حال، تولید داده‌های همبسته تنها برای توزیعهای خاص (بهنجار دو- متغیره) امکان‌پذیر است، در حالی‌که در پاره‌ای از مواقع، همچون تحلیل فراوانی منطقه‌ای، پرداختن به توزیعهای غیربهنجار مورد نیاز است. برای حل مساله­ی مزبور، در این مقاله استفاده از الگوریتم ژنتیک به‌عنوان روش ابتکاری پیشنهاد گردیده است. برای تطابق الگوریتم ژنتیک با مساله‌ی در دست بررسی، به‌جای تلاقی دو کروموزوم از دست‌کاری در هر کروموزوم استفاده شد، زیرا قرارگیری مقادیر در یک گروه زمانی به‌گونه‌ای است که تکرار اعداد مجاز نبوده، و تمامی اعداد باید در یک گروه مورد استفاده قرار گیرند. این شیوه برای یک متغیر تصادفی انجام گرفت و نشان داده شد که الگوریتم ژنتیک منجر به پاسخهای چندگانه می­گردد، در حالی‌که در روش مونت‌کارلو، تنها پاسخی یکتا به‌دست می‌آید. با ارائه‌ی تابع هدف مناسب، روش پیشنهادی برای توزیعهای دو- متغیره نیز به‌کار رفت و شمول گسترده‌تر پاسخها در مقایسه با روش مونت‌کارلو (تنها برای توزیعهای بهنجار) ارایه شد. دلیل گستردگی پاسخ در روش پیشنهادی به این حقیقت باز می­گردد که یک گروه زمانی را می‌توان به چندین شیوه‌ی گوناگون مرتب کرد که بتواند با یک گروه زمانی ثابت ضریب همبستگی یکسان را تولید کند. نشان داده شد که هرچه ضریب همبستگی کمتر باشد، شمار حالتهای گوناگون افزایش می‌یابد.

کلیدواژه‌ها


عنوان مقاله [English]

Enhancing the Performance of the Mont Carlo Simulation for Generating Correlated Data: An Application to Climate Change

نویسندگان [English]

  • B. Ghahraman 1
  • A. Ghahraman 2
چکیده [English]

The Mont Carlo simulation is widely used for data generation based on specified distribution. However, generation of correlated data is only possible for specific distributions (bivariate normal), while it is required to deal with non-normal distributions, e.g. in regional frequency analysis. To tackle the problem, an innovative approach was proposed by using the genetic algorithm. Since it is possible only to change the order of numbers in an array, our problem is not completely in agreement with the classic genetic algorithm. Therefore, we reordered the algorithm; the internal of each chromosome was changed instead of crossing two different ones. The proposed method was applied for one random variable, and it was shown that the genetic algorithm produces multi-solutions, while unique solution is obtained using the Mont Carlo simulation procedure. We presented an appropriate objective function to include any bivariate distribution and showed its hither versatility as compared with the Mont Carlo simulation (only for the bivariate normal distribution). The multi-solutions in our proposed methodology is due to the fact that any array can be reordered in different possible combinations, each of which having a unique correlation coefficient with another fixed array. It was shown that these different combinations increase as the correlation coefficient decreases.

کلیدواژه‌ها [English]

  • Optimization
  • non-normal distribution
  • normal distribution
  • heuristic methods
  • regional frequency analysis
  1. Ang, A.H-S., and W.H. Tang. 1984. Probability concepts in engineering and design. Volume II- Decision, Risk, and Reliability. John Wiley and Sons, Inc. 562 p.
  2. Ang, A.H-S., and W.H. Tang. 2007. Probability concepts in engineering. Empasis on applications to civil and environmental engineering. John Wiley and Sons, Inc. 2nd edition. 406 p.
  3. Casterllarin, A., D.H., Burn, and A. Brath, 2008. Homogeniety testing: how homogeneous do heterogeneous cross-correlated regions seem? J. Hydrol.360,67-76.
  4. Hosking, J.R.M. 1994. The four- parameter kappa distribution. IBM J. Res. Develop.38, 251-258.
  5. Hosking, J.R.M., and J.R. Wallis. 1997. Regional frequency alalysis. An approach based on L-moments. Cambridge University Press. 224 p.
  6. Kleidon, A. 2010. Non-equilibrium thermodynamics, maximum entropy production and earth-system evolution. Philos. Tran. Royal Soc. A. 368(1910), 181-196. doi: 10.1098/ rsta. 2009. 01188.
  7. Rodriguez-Iturbe, I., A. Rinaldo, R. Rigon, R.L. Bras, A. Marani, and E. Ijjasz-Vasquez. 1992. Fractal structure as least energy patterns: The case of river networks. Geophys. Res. Lett. 19:889-892.