2

Реализую функцию нахождения корней методом дихотомии. Не стану здесь объяснять его суть - код реализован, причем достаточно корректно. Однако вот в чем загвоздка.

Функция dichotomy() содержит логику расчета корней уравнения, и мне надо получить их вывод в удобочитаемом виде. Сейчас когда ставлю print(root), я получаю вывод в столбик. Мне же требуется списком в строку, причем с нужным количеством цифр после запятой.

По отдельности это реализовано:

import numpy as np
import timeit

start_time = time.time()
res = list(dichotomy(0.0001,184,10000, 0.000001)) # написанная нами функция с заданными пользователем аргументами
print('Корни по методу дихотомии находятся в точках:')
print(', '.join(map(lambda x: f'{x:.4f}', res)))
print(f'Время счета по модулю timeit: ', timeit.timeit('dichotomy(0.0001,184,10000, 0.000001)', 'from __main__ import dichotomy'), 'seconds')

Однако занести код вывода в тело функции не получается. Это первая трудность.

Вторая: если объявляю время счета timeit.timeit() вне функции, то метод работает, внутри - уже нет. В общем сейчас этот код дает сбои в финальных строках:

import numpy as np
import timeit

def f(x):
    return 1.2-np.log(x)-4*np.cos(2*x)

def dichotomy (a,b,n, eps): # отрезок от a до b делим на n частей, погрешность eps
    """
    Функция отделения и уточнения корня
    """
    assert a!=0,  'a равно 0'
    assert b!=0, 'b равно 0'
    
    # сначала отделим корни
    grid=np.linspace(a, b, n)
    
    # далее уточним корни
    for x,y in zip(grid, grid[1:]):
        if f(x) * f(y) > 0: 
            continue
        root = None
        while ( abs(f(y)-f(x)) )>eps:    
            mid = (y+x)/2                   
            if f(mid) == 0 or f(mid)<eps:   
                root = mid                 
                break
            elif (f(mid) * f(x)) < 0:       
                y = mid                     
            else:
                x = mid                     
        if root:
            yield root
         #   print(root) # дает вывод столбиком - неудобно, потому пробуем ниже:
            res = list(root)   # ОШИБКА здесь!
            print(', '.join(map(lambda x: f'{x:.4f}', res)))

    print(f'Время счета по модулю timeit: ', timeit.timeit('dichotomy(0.0001,184,10000, 0.000001)', 'from __main__ import dichotomy'), 'seconds')

print(list(dichotomy(0.0001,184,10000, 0.000001))) # здесь юзер пытается подставить в функцию значения своих аргументов и получить распечатку корней

На выходе сейчас ошибка с указанием на res = list(root):

TypeError: 'numpy.float64' object is not iterable

1 ответ 1

2
+50

Ну тут возможны разные варианты. Но если вы делаете именно генератор, то лучше выводить результат уже после работы генератора:

# получаем результат
res = list(dichotomy(0.0001,184,10000, 0.000001))

# выводим результат
print(', '.join(map(lambda x: f'{x:.4f}', res)))

А из генератора две строки после yield root просто уберите.

Если же вы хотите делать печать именно в генераторе, то будет не совсем красиво - придётся следить за запятыми, либо в конце будет лишняя запятая. Примерно можно сделать как-то так в конце функции:

        if root:
            yield root
            print(f'{root:.4f}', end=', ')
    print()
    print(f'Время счета по модулю timeit: ', timeit.timeit('dichotomy(0.0001,184,10000, 0.000001)', 'from __main__ import dichotomy'), 'seconds')
5
  • сделал по второму варианту (первый-то был изначально) - и корни посчитались, только заменил внутри функции одну из компонент принта на timeit.timeit(dichotomy(a,b,n, eps). Но... потом поймалась ошибка NameError: name 'a' is not defined, и время счета не вышло 16 мар 2022 в 17:18
  • @АлексейКазанцев А что за a, откуда она эта переменная?
    – CrazyElf
    16 мар 2022 в 17:34
  • фактически это аргумент функции dichotomy(), задающий начало отрезка для поиска корня 16 мар 2022 в 18:03
  • @АлексейКазанцев Ну так его надо как-то задать, это должна быть конкретная переменная при вызове функции. Она у вас есть?
    – CrazyElf
    16 мар 2022 в 18:34
  • вот при вызове и задаю: print(list(dichotomy(0.0001,184,10000, 0.000001))). А в теле функции прописываю print(f'Время счета по модулю timeit: ', timeit.timeit('dichotomy(a,b,n, eps)', 'from __main__ import dichotomy'), 'seconds'). Однако почему-то ошибка 16 мар 2022 в 18:39

Ваш ответ

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

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