-
Notifications
You must be signed in to change notification settings - Fork 82
Expand file tree
/
Copy pathgrad1_2.py
More file actions
41 lines (30 loc) · 1.08 KB
/
grad1_2.py
File metadata and controls
41 lines (30 loc) · 1.08 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
import time
import numpy as np
import matplotlib.pyplot as plt
def f(x):
return np.sin(x)+0.5*x
def df(x):
return np.cos(x) + 0.5
N = 20 # число итераций
xx = 2.5 # начальное значение
lmd = 0.9 # шаг сходимости
x_plt = np.arange(-5.0, 5.0, 0.1)
f_plt = [f(x) for x in x_plt]
plt.ion() # включение интерактивного режима отображения графиков
fig, ax = plt.subplots() # Создание окна и осей для графика
ax.grid(True) # отображение сетки на графике
ax.plot(x_plt, f_plt) # отображение параболы
point = ax.scatter(xx, f(xx), c='red') # отображение точки красным цветом
mn = 100
for i in range(N):
lmd = 1/min(i+1, mn)
xx = xx - lmd*np.sign(df(xx)) # изменение аргумента на текущей итерации
point.set_offsets([xx, f(xx)]) # отображение нового положения точки
# перерисовка графика и задержка на 20 мс
fig.canvas.draw()
fig.canvas.flush_events()
time.sleep(0.02)
plt.ioff() # выключение интерактивного режима отображения графиков
print(xx)
ax.scatter(xx, f(xx), c='blue')
plt.show()