شبیه سازی پیاده روی تصادفی در پایتون — راهنمای گام به گام

۷۵۱ بازدید
آخرین به‌روزرسانی: ۰۳ اردیبهشت ۱۴۰۱
زمان مطالعه: ۴ دقیقه
شبیه سازی پیاده روی تصادفی در پایتون — راهنمای گام به گام

پیاده روی تصادفی (Random Walk) فرایندی است که در طول زمان رخ می‌دهد و هر قدم (Step) به‌صورت تصادفی ایجاد می‌شود و موقعیت بعدی را ایجاد می‌کند. در این آموزش، با روش شبیه سازی پیاده روی تصادفی در پایتون آشنا می‌شویم.

پیاده روی تصادفی چیست؟

در شکل زیر حرکت براونی (Brownian Motion) برای ذرات گاز نشان داده شده است که نقطه ابتدا و انتها نشان داده شده و نمونه‌ای از پیاده روی تصادفی است.

حرکت براونی

در پیاده‌روی تصادفی، از توزیع‌های مختلف مانند یکنواخت و نرمال می‌توان استفاده کرد. اندازه گام نیز می‌تواند پیوسته یا گسسته باشد.

ابتدا می‌خواهیم یک پیاده‌روی تصادفی در دو بعد ایجاد کنیم. بدین منظور، یک موقعیت برای ذره تعریف می‌کنیم:

$$ \begin{aligned}
&x_{0}=0 \\
&y_{0}=0
\end{aligned} $$

سپس به تعداد مراحلی مشخص (N) اعدادی تصادفی با توزیع یکنواخت از $$-1$$ تا $$+1$$ ایجاد می‌کنیم:

$$\begin{aligned}
&\Delta x_{t}=\text { uniform }(-1,+1) \\
&\Delta y_{t}=\text { uniform }(-1,+1)
\end{aligned}$$

سپس، می‌توانیم موقعیت ذره را اصلاح کنیم:

$$\begin{aligned}
&x_{t}=x_{t-1}+\Delta x_{t} \\
&y_{t}=y_{t-1}+\Delta y_{t}
\end{aligned}$$

به این ترتیب، موقعیت در هر لحظه، از روی موقعیت قبلی و با تغییراتی تصادفی ایجاد می‌شود.

پیاده روی تصادفی در پایتون

برای پیاده‌سازی، وارد محیط برنامه‌نویسی می‌شویم:

1import numpy as np
2import matplotlib.pyplot as plt

حال تنظیمات زیر را اعمال می‌کنیم:

1np.random.seed(0)
2plt.style.use('ggplot')

در ابتدا، موقعیت‌های اولیه را تعریف می‌کنیم:

1x0 = 0
2y0 = 0

سپس تعداد گام‌های شبیه‌سازی را تعریف می‌کنیم:

1N = 100

حال یک فهرست برای ذخیره Xها و یک فهرست دیگر برای ذخیره Yها ایجاد می‌کنیم:

1X = [x0]
2Y = [y0]

حال می‌توانیم یک حلقه ایجاد می‌کنیم و در هر مرحله، به‌صورت تصادفی اندازه گام محاسبه و به آخرین موقعیت اضافه کنیم:

1for _ in range(N):
2    dx = np.random.uniform(low=-1, high=+1)
3    dy = np.random.uniform(low=-1, high=+1)
4    newx = X[-1] + dx
5    newy = Y[-1] + dy
6    X.append(newx)
7    Y.append(newy)

به این ترتیب، مقادیر محاسبه و در دو فهرست ذخیره می‌شود. حال می‌توانیم یک نموداری براساس دو لیست ایجاد کنیم تا حرکات به‌خوبی مشاهده شود:

1plt.plot(X, Y, ls='--', lw=0.9, c='crimson', marker='o', ms=4)
2plt.title('Random Walk')
3plt.xlabel('X')
4plt.ylabel('Y')
5plt.show()

به این ترتیب، نمودار زیر مشاهده می‌شود.

شبیه سازی پیاده روی تصادفی در پایتون

به این ترتیب، حرکات به‌خوبی نشان داده می‌شود.

برای یادگیری برنامه‌نویسی با زبان پایتون، پیشنهاد می‌کنیم به مجموعه آموزش‌های مقدماتی تا پیشرفته پایتون فرادرس مراجعه کنید که لینک آن در ادامه آورده شده است.

  • برای مشاهده مجموعه آموزش‌های برنامه نویسی پایتون (Python) — مقدماتی تا پیشرفته + اینجا کلیک کنید.

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

1P = np.zeros((N+1, 2))
2
3for i in range(N):
4    P[i+1] = P[i] + np.random.uniform(low=-1, high=+1, size=(2, ))
5
6plt.plot(P[:, 0], P[:, 1], ls='--', lw=0.9, c='crimson', marker='o', ms=4)
7plt.title('Random Walk')
8plt.xlabel('X')
9plt.ylabel('Y')
10plt.show()

در این حالت نیز نتایج قبلی ایجاد خواهد شد. می‌توان گام‌ها را به‌صورت تصادفی با توزیع نرمال ایجاد کرد:

1for i in range(N):
2    P[i+1] = P[i] + np.random.normal(loc=0, scale=1, size=(2, ))

که در این صورت، نموداری به شکل زیر خواهد بود.

پیاده روی تصادفی

توجه داشته باشید که با توجه به رفتار هر سیستم، باید از توزیع متناسبی استفاده کرد.

به این ترتیب، شبیه‌سازی پیاده‌روی تصادفی دوبعدی پیاده‌سازی شد.

حال می‌خواهیم در بخش دوم از مطلب، به شبیه‌سازی پیاده‌روی تصادفی قیمت یک نماد بپردازیم. توجه داشته باشید که در این حالت، پیشفرضمان را به این صورت در نظر می‌گیریم که رفتار سری زمانی کاملاً تصادفی بوده و توزیع حرکات را می‌دانیم.

به این منظور، کتابخانه‌های مورد نیاز را فراخوانی می‌کنیم و تنظمات را اعمال می‌کنیم:

1import numpy as np
2import yfinance as yf
3import matplotlib.pyplot as plt
4
5np.random.seed(0)
6plt.style.use('ggplot')

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

1DF = yf.download('BTC-USD', start='2019-01-01', end='2022-01-01', interval='1d')
2C = DF['Close'].to_numpy()
3R = C[1:] / C[:-1] - 1

توجه داشته باشید که می‌توان به‌شکل زیر نیز عمل کرد:

1DF = yf.download('BTC-USD', start='2019-01-01', end='2022-01-01', interval='1d')
2C = DF['Close'].to_numpy()
3R = DF['Close'].pct_change().dropna().to_numpy()

حال می‌توانیم یک نمودار هیستوگرام (Histogram Plot) برای تغییرات نسبی قیمت رسم کنیم:

1plt.hist(R, bins=51, color='crimson', alpha=0.8)
2plt.title('Relative Change Histogram')
3plt.xlabel('Relative Change')
4plt.ylabel('Frequency')
5plt.show()

که شکل زیر را خواهیم داشت.

نمودار هیستوگرام

مشاهده می‌کنیم که توزیع داده‌ها به توزیع نرمال نزدیک است. به همین دلیل، می‌توانیم از توزیع نرمال برای تولید گام‌ها استفاده کنیم. برای ایجاد این توزیع، میانگین (Mean) و انحراف معیار (Standard Deviation) داده‌ها را محاسبه می‌کنیم:

1m = np.mean(R)
2s = np.std(R)

حال می‌توانیم موقعیت اولیه را تعریف کنیم:

1p = C[-1]

حال تعداد گام را تعیین و یک آرایه برای ذخیره مقادیر جدید ایجاد می‌کنیم:

1N = 200
2
3P = np.zeros(N)

حال می‌توانیم حلقه اصلی را ایجاد کنیم و مقدایر گام را تعیین کنیم:

1for i in range(N):
2    r = np.random.normal(loc=m, scale=s)

عدد $$r$$ نشان‌دهنده تغییرات نسبی خواهد بود:

$$ r=\frac{P_{t+1}}{P_{t}}-1 \Rightarrow r+1=\frac{P_{t+1}}{P_{t}} \Rightarrow P_{t+1}=P_{t} \times(r+1) $$

اکنون مقدار بعدی را می‌توانیم محاسبه کنیم:

1for i in range(N):
2    r = np.random.normal(loc=m, scale=s)
3    p *= (r + 1)
4    P[i] = p

به این ترتیب، مقادیر محاسبه شد. برای رسم نمودار، به شکل زیر عمل می‌کنیم:

1plt.semilogy(T1, C, lw=0.8, c='crimson', label='Real Data')
2plt.semilogy(T2, P, lw=0.8, c='teal', label='Random Walk')
3plt.title('Real Data + Random Walk')
4plt.xlabel('Time (Day)')
5plt.ylabel('Price ($)')
6plt.legend()
7plt.show()

که در این صورت، نمودار به شکل زیر خواهد بود.

نمودار پیاده روی تصادفی

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

اگر بخواهیم چندین سناریو را به‌صورت همزمان در نمودار نشان دهیم، می‌توانیم به شکل زیر کد را تغییر دهیم:

1N = 200
2nScenario = 10
3
4T1 = np.arange(start=0, stop=C.size)
5T2 = np.arange(start=C.size, stop=C.size + N)
6
7plt.semilogy(T1, C, lw=0.8, c='crimson', label='Real Data')
8for i in range(nScenario):
9    p = C[-1]
10    P = np.zeros(N)
11    for i in range(N):
12        r = np.random.normal(loc=m, scale=s)
13        p *= (r + 1)
14        P[i] = p
15    plt.semilogy(T2, P, lw=0.6)
16plt.title('Real Data + Random Walk')
17plt.xlabel('Time (Day)')
18plt.ylabel('Price ($)')
19plt.legend()
20plt.show()

که در این صورت، به تعداد ۱۰ بار، پیاده‌روی تصادفی انجام می‌شود و نتایج زیر حاصل می‌شود.

پیاده‌روی تصادفی

به این ترتیب، مشاهده می‌کنیم که هر سناریو از یک نقطه شروع شده ولی رفته‌رفته فاصله هرکدام از سایرین افزایش می‌یابد. با انجام یک مدل‌سازی مناسب و دقیق‌تر، می‌توان پیش‌بینی‌هایی انجام داد که به واقعیت نزدیک هستند و در نتیجه، با بررسی سناریو‌های مختلف، احتمال مشاهده قیمتی مشخص در تاریخی مشخص را محاسبه کرد.

جمع‌بندی

در این مطلب، با پیاده‌روی تصادفی آشنا شدیم و دو شبیه‌سازی با استفاده از آن انجام دادیم. برای مطالعه بیشتر، می‌توان به موارد زیر مراجعه کرد:

  1. روند (Trend) در داده‌های سری زمانی
  2. اجزای تشکیل‌دهنده یک سری زمانی
  3. توزیع‌های احتمال پرکاربرد
  4. خودهمبستگی در سری‌های زمانی
بر اساس رای ۱۶ نفر
آیا این مطلب برای شما مفید بود؟
اگر بازخوردی درباره این مطلب دارید یا پرسشی دارید که بدون پاسخ مانده است، آن را از طریق بخش نظرات مطرح کنید.
منابع:
مجله فرادرس
نظر شما چیست؟

نشانی ایمیل شما منتشر نخواهد شد. بخش‌های موردنیاز علامت‌گذاری شده‌اند *