Как сделать анимацию движения объекта по кривой (алгоритм Безье)

@Bronzor · Регистрация: 13.10.2018

Студворк — интернет-сервис помощи студентам

Смог написать программу для рисования кривой используя алгоритм Безье (anaconda), можете пожалуйста помочь, что нужно добавить, чтобы получилась анимация движения, например круга по этой кривой. Просят использовать ffmpeg, celluloid, matplotlib для анимации и Tkinter или что-то похожее для графического интерфейса, читал про библиотеки, но пока туго с пониманием того, как добиться нужного результата, надеюсь на вашу помощь.

Python
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
import numpy as np               
import matplotlib.pyplot as plt  
 
img = np.ones((600,600,3))
 
def DrawBezier(x0,y0,x1,y1,x2,y2,x3,y3,x4,y4,x5,y5):
    
    x = x0
    y = y0
    t = 0
    ts = 0.001
 
    while t <= 1.0:
        
        B0 = (1-t)*(1-t)*(1-t)*(1-t)*(1-t)
        B1 = 5*t*(1-t)*(1-t)*(1-t)*(1-t)
        B2 = 10*t*t*(1-t)*(1-t)*(1-t)
        B3 = 10*t*t*t*(1-t)*(1-t)
        B4 = 5*t*t*t*t*(1-t)
        B5 = t*t*t*t*t
        
        x = int(x0*B0 + x1*B1 + x2*B2 + x3*B3 + x4*B4 + x5*B5)
        y = int(y0*B0 + y1*B1 + y2*B2 + y3*B3 + y4*B4 + y5*B5)
        
        img[y,x] = (1,0,0)
        t=t+ts 
        
    return
   
DrawBezier(180,290,140,335,50,170,10,15,360,240,140,67)
    
%matplotlib notebook
plt.figure(figsize=(10,8), dpi = 100, facecolor='w')
plt.imshow(img)
plt.gca().invert_yaxis()
plt.show()

Название: jr.JPG
Просмотров: 237

Размер: 10.2 Кб

Хочется, чтобы было вот так:

Как сделать анимацию движения объекта по кривой (алгоритм Безье)

@Excalibur921 · 08.01.2021, 20:32

Сообщение от Bronzor

B0 = (1-t)*(1-t)*(1-t)*(1-t)*(1-t)
B1 = 5*t*(1-t)*(1-t)*(1-t)*(1-t)
B2 = 10*t*t*(1-t)*(1-t)*(1-t)
B3 = 10*t*t*t*(1-t)*(1-t)
B4 = 5*t*t*t*t*(1-t)
B5 = t*t*t*t*t

Никогда такого не видел. На википедии формула в одну строку, что это за монстр и зачем?
Это наверно высокой степени. Их вроде как никто не использует. Набирают кривую из кусков квадратичной кривой. Можно легко сделать сложной формы и замкнуть.
Как правильно задать систему урАвнений для кубических сплайнов?
Для анимации по таймеру на каждом кадре увеличиваете параметр t, подставляете в формулу кривой получаете координаты точки меняются по времени. Когда будет 2 кривых вместе то нужно будет завести два параметра, один общий t0 от 0 до 2 и второй t должен расти от 0 до 1 смотря какому диапазону принадлежит t0.Чтобы точка на кривой прошла сначала одну кривую затем вторую.

@Bronzor · 09.01.2021, 11:29 **[ТС]**

Это 6 контрольных точек, можно и с двумя, мне главное понять, желательно код с маленькими пояснениями, как на моей кривой нарисовать круг и заставить его двигаться по этой кривой.

@Excalibur921 · 09.01.2021, 13:56

Сообщение от Bronzor

желательно код с маленькими пояснениями

В гугле "point animation python"
Включил “картинки” смотрю где 1 точка а не 100500.
https://stackoverflow.com/ques... f-a-circle
Поменял строку

Python
1
 return np.array([r*np.cos(phi), r*np.sin(phi)])

и добавил опорные точки.

Python
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
42
43
import numpy as np
import matplotlib.pyplot as plt
plt.rcParams["figure.figsize"] = 4,3
from matplotlib.animation import FuncAnimation
 
p0x=0.0
p0y=-0.5
p1x=0.0
p1y=0.5
p2x=0.5
p2y=0.0
 
r = 1 # radius of circle
def circle(phi):
    #return np.array([r*np.cos(phi), r*np.sin(phi)])
    return np.array([
        (1.0-phi)**2*p0x + 2*phi*(1.0-phi)*p1x + phi**2*p2x,
        (1.0-phi)**2*p0y + 2*phi*(1.0-phi)*p1y + phi**2*p2y,
        ])
 
# create a figure with an axes
fig, ax = plt.subplots()
# set the axes limits
ax.axis([-1.5,1.5,-1.5,1.5])
# set equal aspect such that the circle is not shown as ellipse
ax.set_aspect("equal")
# create a point in the axes
point, = ax.plot(0,1, marker="o")
 
# Updating function, to be repeatedly called by the animation
def update(phi):
    # obtain point coordinates 
    x,y = circle(phi)
    # set point's coordinates
    point.set_data([x],[y])
    return point,
 
# create animation with 10ms interval, which is repeated,
# provide the full circle (0,2pi) as parameters
ani = FuncAnimation(fig, update, interval=30, blit=True, repeat=True,
                    frames=np.linspace(0,2*np.pi,360, endpoint=False))
 
plt.show()

Тут другой с экспортом в mp4.
https://medium.com/@kirthi.tej... 0e075e5c6e
закоментил экспорт и добавил

Python
1
 plt.show()

Python
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
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.animation as animation
 
TWOPI = 2*np.pi
 
fig, ax = plt.subplots()
 
t = np.arange(0.0, TWOPI, 0.001)
s = np.sin(t)
l = plt.plot(t, s)
 
ax = plt.axis([0,TWOPI,-1,1])
 
redDot, = plt.plot([0], [np.sin(0)], 'ro')
 
def animate(i):
    redDot.set_data(i, np.sin(i))
    return redDot,
 
# create animation using the animate() function
myAnimation = animation.FuncAnimation(fig, animate, frames=np.arange(0.0, TWOPI, 0.1),interval=10, blit=True, repeat=True)
 
plt.show()
 
# myAnimation.save('sine.mp4', fps=30, extra_args=['-vcodec', 'libx264'])

@u235 · 09.01.2021, 15:35

Python
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
42
43
44
45
from matplotlib import pyplot as plt
from celluloid import Camera
import numpy as np
 
tau=np.linspace(0,1,1000)
P=np.array([[180, 290],
           [140, 335],
           [50,  170],
           [10,   15],
           [360, 240],
           [140,  67]
           ])
 
def Bezier(t, P):
    # Bezier Matrix:
    Mb=np.matrix([
        [ -1,   5, -10,  10, -5, 1],
        [  5, -20,  30, -20,  5, 0],
        [-10,  30, -30,  10,  0, 0],
        [ 10, -20,  10,   0,  0, 0],
        [ -5,   5,   0,   0,  0, 0],
        [  1,   0,   0,   0,  0, 0]
                 ])
    # T matrix
    xx, yy=np.meshgrid(np.arange(5, -1, -1), t)
    xx=np.matrix(xx)
    yy=np.matrix(yy)
    T_mat=np.power(yy, xx)
    B=T_mat*Mb*P
    return B
    
 
fig = plt.figure(figsize=(10,8), dpi = 100, facecolor='w')
plt.axis('equal')
plt.axis([0,600,0,600])
X,Y=Bezier(tau, P)[:,0],Bezier(tau, P)[:,1]
 
camera = Camera(fig)
for i in range(100):
    coord=Bezier(i/100, P)
    plt.plot(X,Y,color = '#ff0000')
    plt.scatter(coord[0,0], coord[0,1], color = '#00ff00' , marker='o')
    camera.snap()
animation = camera.animate()
animation.save('animation.mp4')

https://youtu.be/WSB_Smh6aIw

@Bronzor · 09.01.2021, 23:50 **[ТС]**

Excalibur921, В anaconda ни один код не заработал, нету анимации просто статичная картинка. Как я понимаю, первый код это движение точки по кругу, а второй по синусоиде, а как передвигать по точкам кривой Безье?
u235, Всё работает, большое спасибо за код, но можно побольше объяснений, откуда взялась матрица 6 на 6, почему именно такие значения, зачем нужна T матрица, meshgrid? какой параметр отвечает за ускорение/замедление передвижения маркера?

Добавлено через 17 минут
u235, Анимация почему-то работает, только если добавить строчку %matplotlib qt, а затем удалить её, если просто скопировать код и запустить, отображается просто координатная система без линии и анимации.
Запускаю этот код

Кликните здесь для просмотра всего текста

Python
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
42
43
from matplotlib import pyplot as plt
from celluloid import Camera
import numpy as np
import ffmpeg
tau=np.linspace(0,1,1000)
P=np.array([[180, 290],
           [140, 335],
           [50,  170],
           [10,   15],
           [360, 240],
           [140,  67]
           ])
 
def Bezier(t, P):
    # Bezier Matrix:
    Mb=np.matrix([
        [ -1,   5, -10,  10, -5, 1],
        [  5, -20,  30, -20,  5, 0],
        [-10,  30, -30,  10,  0, 0],
        [ 10, -20,  10,   0,  0, 0],
        [ -5,   5,   0,   0,  0, 0],
        [  1,   0,   0,   0,  0, 0]
                 ])
    # T matrix
    xx, yy=np.meshgrid(np.arange(5, -1, -1), t)
    xx=np.matrix(xx)
    yy=np.matrix(yy)
    T_mat=np.power(yy, xx)
    B=T_mat*Mb*P
    return B
 
fig = plt.figure(figsize=(10,8), dpi = 100, facecolor='w')
plt.axis('equal')
plt.axis([0,600,0,600])
X,Y=Bezier(tau, P)[:,0],Bezier(tau, P)[:,1]
camera = Camera(fig)
for i in range(100):
    coord=Bezier(i/100, P)
    plt.plot(X,Y,color = '#ff0000')
    plt.scatter(coord[0,0], coord[0,1], color = '#00ff00' , marker='o')
    camera.snap()
animation = camera.animate(interval = 200, repeat = True,
                           repeat_delay = 500)

Всё работает, перезагружаю kernel, всё работает, анимация показывается, вставляю ваш код, анимация работает, перезагружаю kernel снова просто координатная плоскость, иначе говоря, если в новый файл вставить ваш код анимация не работает, почему?

@Excalibur921 · 09.01.2021, 23:51

Сообщение от Bronzor

В anaconda ни один код не заработал, нету анимации просто статичная картинка. Как я понимаю, первый код это движение точки по кругу,

Не знаю что такое anaconda.
Я настроил Visual Studio Code для Python 3.8 и коды работают. Первый код движение по кривой Безье. Второй по синусоиде с показом кривой. Плохо что одну и туже задачу можно рисовать 100500 способами .

@Bronzor · 09.01.2021, 23:58 **[ТС]**

Да, как оказалось, та же проблема, если вставить этот код и затем вставить ваш, анимация работает, но если вставить код в новый файл, будет статичная картинка, пока не понятно как это исправить.

Кликните здесь для просмотра всего текста

Python
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
42
43
from matplotlib import pyplot as plt
from celluloid import Camera
import numpy as np
import ffmpeg
tau=np.linspace(0,1,1000)
P=np.array([[180, 290],
           [140, 335],
           [50,  170],
           [10,   15],
           [360, 240],
           [140,  67]
           ])
 
def Bezier(t, P):
    # Bezier Matrix:
    Mb=np.matrix([
        [ -1,   5, -10,  10, -5, 1],
        [  5, -20,  30, -20,  5, 0],
        [-10,  30, -30,  10,  0, 0],
        [ 10, -20,  10,   0,  0, 0],
        [ -5,   5,   0,   0,  0, 0],
        [  1,   0,   0,   0,  0, 0]
                 ])
    # T matrix
    xx, yy=np.meshgrid(np.arange(5, -1, -1), t)
    xx=np.matrix(xx)
    yy=np.matrix(yy)
    T_mat=np.power(yy, xx)
    B=T_mat*Mb*P
    return B
 
fig = plt.figure(figsize=(10,8), dpi = 100, facecolor='w')
plt.axis('equal')
plt.axis([0,600,0,600])
X,Y=Bezier(tau, P)[:,0],Bezier(tau, P)[:,1]
camera = Camera(fig)
for i in range(100):
    coord=Bezier(i/100, P)
    plt.plot(X,Y,color = '#ff0000')
    plt.scatter(coord[0,0], coord[0,1], color = '#00ff00' , marker='o')
    camera.snap()
animation = camera.animate(interval = 200, repeat = True,
                           repeat_delay = 500)

@u235 · 10.01.2021, 00:41

Bronzor, Матрица 6х6 это матрица Безье, если из ваших B0-B5 сделать полиномы (раскрыть скобки ), то матрица Безье это коэффициенты при t^5.. t^0 степени.
Матрица T - матрица степеней вектора времени. Если вектор t=0.5 (одна точка) T=[0.5**5, 0.5**4, 0.5**3,..0.5**0]
Если вектор t=[0, 0.5, 1] (три точки), то матрица T=[[0**5, 0**4,..0**0], [0.5**5, 0.5**4,..0.5**0],[1**5, 1**4,..1**0]].
См. https://ru.wikipedia.org/wiki/... 1%8C%D0%B5
Т.о. мы ускоряем код, используя быстрое умножение матриц numpy, а не циклы.

Добавлено через 20 минут
Попробуйте добавить в конец скрипта

Python
1
 plt.show()

@Bronzor · 10.01.2021, 12:20 **[ТС]**

u235, После добавления ничего не изменилось. При этом анимация сохраняется и с ней всё хорошо.

@u235 · 10.01.2021, 12:33

Так вы через Jupyter запускаете.. Там работать не будет. Запустите просто из комм.строки анаконды.

@Bronzor · 10.01.2021, 12:34 **[ТС]**

u235, Notebook, а почему другая анимация работает, а эта нет?

@u235 · 10.01.2021, 12:48

Bronzor, откуда же я знаю.

У меня нет кода другой анамации чтобы сравнить, а у вас есть. Вот и найдите различия.

@Bronzor · 10.01.2021, 13:59 **[ТС]**

Оказывается два раза сверху не то скопировал, тоже самое всё, но эта анимация работает...

Python
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
#Importing required libraries
%matplotlib notebook
import matplotlib.pyplot as plt
import numpy as np
from celluloid import Camera
import ffmpeg
 
 
#Creating Data
x = np.linspace(0, 10, 100)
 
#defining a function to return sine of input values.
def fun(i):
    y = np.sin(i)
    return y
 
x_y = fun(x)
 
#Creating matplotlib figure and camera object
fig = plt.figure()
plt.xlim(0,10)
plt.ylim(-2,2)
camera = Camera(fig)
 
#Looping the data and capturing frame at each iteration
for i in x:
    plt.plot(x,x_y , color = 'green' , lw = 0.8)
    f = plt.scatter(i, fun(i) , color = 'red' , s = 200)
    plt.title('tracing a sin function')
    camera.snap()
 
#Creating the animation from captured frames
animation = camera.animate(interval = 200, repeat = True,
                           repeat_delay = 500)

Добавлено через 8 минут
u235, Вставил %matplotlib qt5 перед созданием координатной плоскости (верно? fig = plt.figure()) и теперь анимация работает в отдельном окне.

Python
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
42
43
44
45
from matplotlib import pyplot as plt
from celluloid import Camera
import numpy as np
 
tau=np.linspace(0,1,1000)
P=np.array([[180, 290],
           [140, 335],
           [50,  170],
           [10,   15],
           [360, 240],
           [140,  67]
           ])
 
def Bezier(t, P):
    # Bezier Matrix:
    Mb=np.matrix([
        [ -1,   5, -10,  10, -5, 1],
        [  5, -20,  30, -20,  5, 0],
        [-10,  30, -30,  10,  0, 0],
        [ 10, -20,  10,   0,  0, 0],
        [ -5,   5,   0,   0,  0, 0],
        [  1,   0,   0,   0,  0, 0]
                 ])
    # T matrix
    xx, yy=np.meshgrid(np.arange(5, -1, -1), t)
    xx=np.matrix(xx)
    yy=np.matrix(yy)
    T_mat=np.power(yy, xx)
    B=T_mat*Mb*P
    return B
%matplotlib qt5
fig = plt.figure()
plt.xlim(0,600)
plt.ylim(0,600)
X,Y=Bezier(tau, P)[:,0],Bezier(tau, P)[:,1]
camera = Camera(fig)
 
for i in range(100):
    coord=Bezier(i/100, P)
    plt.plot(X,Y,color = '#ff0000')
    plt.scatter(coord[0,0], coord[0,1], color = '#00ff00' , marker='o')
    camera.snap()
    
animation = camera.animate(interval = 200, repeat = True,
                           repeat_delay = 500)

Добавлено через 53 минуты
Осталось понять как работает код.
tau (что это за сокращение, кстати?) это одномерный массив состоящий из 1000 элементов от 0 до 1?
P это двумерный массив, где хранятся координаты x,y моих контрольных точек.
Mb это двумерный массив, где хранятся коэффициенты при t,
Затем создаём прямоугольную сетку используя одномерные массивы xx и yy (от -5 до -1 с шагом -1), и зачем нужен параметр t здесь? не понял 25 строку.
xx = матрица размерностью хх
yy = матрица размерностью yy
28 снова не понял, но это возведение элемента массива yy в степень элемента массива xx?
Затем создаём объект figure без осей
Затем устанавливаем пределы x и y осей
35 не понял что происходит, но что-то присваивается переменным x и y
Затем в цикле получаем координаты точек, отображаем их на графике красным цветом
Добавляем маркер зелёного цвета (что значит coord[0,0], coord[0,1] ? и какие параметры отвечают за размер кружка, можно ли вместо маркера создать круг и его двигать?)
Затем снимок камеры
После того, как все кадры были получены, создаём анимацию?

@u235 · 10.01.2021, 15:15

tau - параметр кривой Безье - меняется от 0 до 1 , 1000 - число точек из которых она состоит. Если сделать не 1000, а, например, 20, то на кривой будут изломы.
25 строка - (от 5 до 0 включительно) это степени, в которые будет возводится t. Посмотрите примеры np.meshgrid на оф. сайте numpy. На входе - два вектора, на выходе две матрицы-сетки.
28 - да, так и есть.

35 - Получаем два вектора координат X и Y из вектор tau - для статичного графика
coord[0,0], coord[0,1] - координаты точки кривой, соответствующие параметру i/100 - центр кружка.

Python
1
plt.scatter(coord[0,0], coord[0,1], color = '#00ff00' , marker='o', s=10)

s - размер кружка

@Bronzor · 10.01.2021, 15:31 **[ТС]**

u235, вроде стало понятнее, а как переделать код в 3 контрольные точки? попробовал удалить координаты 3 точек, и последние три строки матрицы, но не заработало, видимо не всё так просто, также а как отобразить закрасить контрольные точки?
Нашёл вот такой идеальный код, мне бы анимировать точку и добавить возможность ввода конечной точки всё.

Python
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
import matplotlib.pyplot as plt
import numpy as np
from celluloid import Camera
 
tau=np.linspace(0,1,1000)
P0, P1, P2=np.array([[0, 0],
           [2, 4],
           [5, 3],
           ])
 
P = lambda t: (1 - t)**2 * P0 + 2 * t * (1 - t) * P1 + t**2 * P2
points = np.array([P(t) for t in np.linspace(0, 1, 50)])
 
# get x and y coordinates of points separately
x, y = points[:,0], points[:,1]
 
%matplotlib qt5
fig = plt.figure()
plt.xlim(-1,6)
plt.ylim(-1,6)
camera = Camera(fig)
 
for i in range(100):
    plt.plot(x,y,color = '#ff0000')
    plt.plot(*P0, 'b.')
    plt.plot(*P1, 'b.')
    plt.plot(*P2, 'b.')
    camera.snap()
    
animation = camera.animate(interval = 20, repeat = True,
                           repeat_delay = 0)

@u235 · 10.01.2021, 15:49

Bronzor, Для 3х контрольных точек нужна другая матрица Безье, 3x3, а не 5x5. Как ее считать я уже говорил.
Контрольные точки лучше хранить единой матрицей P, , а не в виде P0, P1... тогда
достаточно вызвать что-то типа:

Python
1
plt.scatter(P[:,0], P[:,1], color = '#ffff00' , marker='o')

чтобы отобразить сразу все контрольные точки за одну команду без циклов.
Работайте с матрицами и векторами как с единым целым, а не поэлементно.
Почитайте хорошую книжку Роджерс Д., Адамс Дж. Математические основы машинной графики.

@Bronzor · 10.01.2021, 16:20 **[ТС]**

u235, Спасибо за разъяснения, пока что такой результат, осталось добавить возможность, чтобы пользователь вводил координаты последней контрольной точки P3.

Python
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
from matplotlib import pyplot as plt
from celluloid import Camera
import numpy as np
 
tau=np.linspace(0,1,1000)
P=np.array([[0, 0],
           [2, 4],
           [5, 3],
           ])
 
def Bezier(t, P):
    # Bezier Matrix:
    Mb=np.matrix([
        [1, -2, 1],
        [-2, 2, 0],
        [1, 0, 0],
                 ])
    # T matrix
    xx, yy=np.meshgrid(np.arange(2, -1, -1), t)
    xx=np.matrix(xx)
    yy=np.matrix(yy)
    T_mat=np.power(yy, xx)
    B=T_mat*Mb*P
    return B
%matplotlib qt5
fig = plt.figure()
plt.xlim(-2,6)
plt.ylim(-2,6)
X,Y=Bezier(tau, P)[:,0],Bezier(tau, P)[:,1]
camera = Camera(fig)
 
for i in range(100):
    coord=Bezier(i/100, P)
    plt.plot(X,Y,color = '#ff0000')
    plt.scatter(P[:,0], P[:,1], color = '#0000FF' , marker='o')
    plt.scatter(coord[0,0], coord[0,1], color = '#00ff00' , marker='o', s=100)
    camera.snap()
    
animation = camera.animate(interval = 20, repeat = True,
                           repeat_delay = 0)
animation.save('1.mp4')

@Bronzor · 10.01.2021, 16:58 **[ТС]**

u235, Добавил ввод, может перемудрил и как-то иначе можно? Правильно ли написал комментарии?

Python
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
42
43
44
45
46
47
48
49
## Импортирование библиотек
from matplotlib import pyplot as plt
from celluloid import Camera
import numpy as np
 
## Создание одномерного массива, значения от 0 до 1, количество точек - 1000 (параметр кривой Безье)
tau=np.linspace(0,1,1000)
print('Введите х координату:')
x_end = input()
print('Введите y координату:')
y_end = input()
P=np.array([[0, 0],  # Координаты контрольны точек записанные в виде двумерного массива
           [2, 4],
           [x_end, y_end],
           ])
P=P.astype(float)
# Алгоритм Безье
def Bezier(t, P):
    # Базисная матрица Безье (коэффициенты при t^2,t^1,t^0)
    Mb=np.matrix([
        [1, -2, 1],
        [-2, 2, 0],
        [1, 0, 0],
                 ])
    # Т матрица (матрица степеней вектора времени)
    xx, yy=np.meshgrid(np.arange(2, -1, -1), t) # Степени, в которые будет возводится t (На входе - два вектора, на выходе две матрицы-сетки)
    xx=np.matrix(xx)
    yy=np.matrix(yy)
    T_mat=np.power(yy, xx) # Возведение элемента массива yy в степень элемента массива xx
    B=T_mat*Mb*P # Получаем координаты точки
    return B # Возвращаем координаты точки
 
# Открытие всплывающего окна (GUI)
%matplotlib qt 
fig = plt.figure() # Создание объекта-фигуры без осей
plt.xlim(-2,6) #Создание и задание пределов оси Х
plt.ylim(-2,6) #Создание и задание пределов оси Y
X,Y=Bezier(tau, P)[:,0],Bezier(tau, P)[:,1] # Получаем два вектора координат X и Y из вектор tau - для статичного графика
camera = Camera(fig) 
 
for i in range(100):
    coord=Bezier(i/100, P)
    plt.plot(X,Y,color = '#ff0000') # Отображение координатной плоскости
    plt.scatter(P[:,0], P[:,1], color = '#0000FF') # Отображение контрольных точек
    plt.scatter(coord[0,0], coord[0,1], color = '#00ff00' , marker='o', s=100) # Координаты точки кривой, соответствующие параметру i/100 - центр кружка.
    camera.snap() # Снимок камеры
    
animation = camera.animate(interval = 20, repeat = True, repeat_delay = 0) # Создание анимации
animation.save('1.mp4') # Сохранение анимации

Введите х координату:
5
Введите y координату:
-1

@u235 · 10.01.2021, 22:43

Все нормально, только код можно еще сократить. В цикле for не считать coord заново, а просто брать скажем каждый 10ый элемент из X , Y

Добавлено через 1 минуту
типа

Python
1
plt.scatter( X[i*10, 0], Y[i*10, 0], color = '#00ff00' , marker='o', s=200)

а 42 строку убрать.

Новые блоги и статьи Все статьи Все блоги /
Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .
SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .

@Bronzor 2 / 2 / 0 Регистрация: 13.10.2018 Сообщений: 245
	09.01.2021, 11:29 [ТС]
	Это 6 контрольных точек, можно и с двумя, мне главное понять, желательно код с маленькими пояснениями, как на моей кривой нарисовать круг и заставить его двигаться по этой кривой. 0

@Bronzor 2 / 2 / 0 Регистрация: 13.10.2018 Сообщений: 245
	10.01.2021, 12:20 [ТС]
	u235, После добавления ничего не изменилось. При этом анимация сохраняется и с ней всё хорошо. 0

@u235 5516 / 2869 / 571 Регистрация: 07.11.2019 Сообщений: 4,759
	10.01.2021, 12:33
	Так вы через Jupyter запускаете.. Там работать не будет. Запустите просто из комм.строки анаконды. 1

@Bronzor 2 / 2 / 0 Регистрация: 13.10.2018 Сообщений: 245
	10.01.2021, 12:34 [ТС]
	u235, Notebook, а почему другая анимация работает, а эта нет? 0

@u235 5516 / 2869 / 571 Регистрация: 07.11.2019 Сообщений: 4,759
	10.01.2021, 12:48
	Bronzor, откуда же я знаю. У меня нет кода другой анамации чтобы сравнить, а у вас есть. Вот и найдите различия. 1