1

У меня есть метод золотого сечения, мне нужно показать как меняется интервал на каждой итерации метода. Нужно сделать анимацию того как меняется интервал на оси X, но пока только вот такой статический вариант есть. Нужна помощь

import math
import matplotlib as mpl
import matplotlib.pyplot as plt
import numpy as np

x_graf1 = []
y_graf1 = []

x_points1 = []
y_points1 = []

x_graf2 = []
y_graf2 = []

x_points2 = []
y_points2 = []


def draw():
    plt.plot(x_graf1, y_graf1)
    plt.plot(x_graf2, y_graf2)
    plt.scatter(x_points1, y_points1)
    plt.scatter(x_points2, y_points2)
    plt.show()


def f(x):
    return 2 * (x - 2) ** 2


def gold_section(a, b, l):
    max_iteration = 100000
    iteration = 0
    left = a + ((3 - math.sqrt(5)) / 2) * (b - a)
    right = a + b - left
    print(f"Текущий интервал: [{a} ; {b}]")
    x_graf1.append(a)
    y_graf1.append(f(a))
    x_graf2.append(b)
    y_graf2.append(f(b))
    x_points1.append(a)
    y_points1.append(0)
    x_points2.append(b)
    y_points2.append(0)
    while iteration < max_iteration:
        Fy = f(left)
        Fz = f(right)
        if Fy <= Fz:
            b = right
            right = left
            left = a + b - left
        else:
            a = left
            left = right
            right = a + b - right


        x_graf1.append(a)
        y_graf1.append(f(a))

        x_graf2.append(b)
        y_graf2.append(f(b))

        x_points1.append(a)
        y_points1.append(0)

        x_points2.append(b)
        y_points2.append(0)

        print(f"Текущий интервал: [{a} ; {b}], Итерация: {iteration + 1}")

        if math.fabs(a - b) <= l:
            print(f"Текущий интервал: [{a} ; {b}], кол-во итераций {iteration + 1}")
            return (a + b) / 2
        iteration += 1
    print("Привышено кол-во итераций")
    return (a + b) / 2


# a, b = map(float, input("Введите интервал(через пробел): ").split())
# l = float(input("Точность: "))
a = -10
b = 10
l = 0.000001
# print(f"Вы выбрали интервал от {a} до {b}, точность: {l}")
print(gold_section(a, b, l))
draw()


введите сюда описание изображения

1
  • Пожалуйста, уточните вашу конкретную проблему или приведите более подробную информацию о том, что именно вам нужно. В текущем виде сложно понять, что именно вы спрашиваете. 4 дек 2022 в 15:52

1 ответ 1

1

не знаю, насколько это то, что вы хотите, но можно накопить концы интервала в массивы A и B, тогда

        df1 = pd.DataFrame({'A':A, 'B':B})
        plt.bar(range(len(df1)), df1.B-df1.A, bottom=A)
        plt.show()

даст примерно такую картинку

введите сюда описание изображения

возможно, надо заморочаться с логарифмической шкалой по ординате

4
  • Я думал по другому, чтобы по точкам двигалось на шкале X. Но это интересный вариант тоже, где ось икс будет номер итераций а Y как раз отражение (a и b) 4 дек 2022 в 15:43
  • именно так. только вот про нелинейный масштаб надо подумать, а то быстро картинка вырождается
    – splash58
    4 дек 2022 в 15:45
  • у нас алгоритм вычисляет конечную точку к которой сходится интервал, можно ли ее указать на графике? 4 дек 2022 в 16:46
  • да, конечно, нарисуйте точку на графике, какие проблемы примерно так plt.plot([x], [y], 'ro')
    – splash58
    4 дек 2022 в 17:31

Ваш ответ

By clicking “Отправить ответ”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.