RError.com

RError.com Logo RError.com Logo

RError.com Navigation

  • 主页

Mobile menu

Close
  • 主页
  • 系统&网络
    • 热门问题
    • 最新问题
    • 标签
  • Ubuntu
    • 热门问题
    • 最新问题
    • 标签
  • 帮助
主页 / 问题 / 1516631
Accepted
Alex_Kazantsev
Alex_Kazantsev
Asked:2023-04-29 23:05:09 +0000 UTC2023-04-29 23:05:09 +0000 UTC 2023-04-29 23:05:09 +0000 UTC

为什么在 Python 脚本中计算时只得到一次迭代?

  • 772

在用二分法求解一个变量的方程时,根据问题的条件,我需要计算求根时的迭代次数。对于一次迭代,正如我所读,您可以将段的划分数减半。对于每个根至少一次,段被划分。即迭代次数显然不会少于根数。但是,我的计数器只显示一个分区,即一次迭代。怎么了?

代码:

import numpy as np
from decimal import *
import pandas as pd

def print_dichotomy(funcs):
    def wrapper(a,b, eps_func, eps_arg):     
        *res, func, eps_list = list(funcs(a,b, eps_func, eps_arg))
        result_table = pd.DataFrame()
        result_table['Корни'] = pd.Series(res)
        print('Количество итераций: ', iter_count)
        print(result_table)
    return wrapper 

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

@print_dichotomy
def dichotomy(a, b, eps_func, eps_arg):      
    # 1. определим глобальную переменную (расчет количества итераций)
    global iter_count
    
    # 2. определим ограничения на границы отрезка 
    assert a != 0, 'a равно 0'
    assert b != 0, 'b равно 0'

    # 3. запросим количество частей n, на которые будет дробиться отрезок 
    n = int( input('Введите параметр n (целое число): ') )
    
    # 5. отделим корни
    # для этого разделим заданный отрезок на n частей 
    grid = np.linspace(a, b, n)
    
    #  6. создадим счетчик для количества итераций 
    iter_count = 0
    
    # 7. уточним наличие корней на каждом отрезке
    for x, y in zip(grid, grid[1:]):
        # 7.1. если на отрезке нет корня, смотрим следующий отрезок
        if f(x) * f(y) > 0:  
            continue
        root = None
        # 7.2. уточним корни
        while abs(f((x+y)/2)) > eps_func and ((y - x)/2) > eps_arg:
            if y != x:
                # 7.2.1.1. получаем середину отрезка
                mid = (y + x) / 2
                # 7.2.1.2. пополним счетчик итераций на один (было одно деление пополам)
                iter_count =+ 1
                # 7.2.1.3. корень примем равным серединному значению:
                root = mid  
                break
                # 7.2.1.5. иначе если произведение функции в середине отрезка на функцию в т. а <= 0:
                if Decimal(f(mid) * f(x)) <= 0:
                    # 7.2.1.5.1. серединой отрезка признаем точку b:
                    y = mid  
                else:
                    # 7.2.1.5.3. в другом случае — точку а:
                    x = mid 
        # 8. возвращаем рассчитанный корень 
        if root:
            yield root

dichotomy(0.0001, 50, 0.000001, 0.000001)
python
  • 1 1 个回答
  • 36 Views

1 个回答

  • Voted
  1. Best Answer
    sLoNcE
    2023-04-30T02:24:50Z2023-04-30T02:24:50Z

    在您的代码中,您使用 iter_count 全局变量来计算迭代次数。iter_count =+ 1但是,在您使用 =+ 运算符的行中,它将值 1 分配给 iter_count 而不是将其递增 1。更改为iter_count += 1

    • 2

相关问题

  • 是否可以以某种方式自定义 QTabWidget?

  • telebot.anihelper.ApiException 错误

  • Python。检查一个数字是否是 3 的幂。输出 无

  • 解析多个响应

  • 交换两个数组的元素,以便它们的新内容也反转

Sidebar

Stats

  • 问题 10021
  • Answers 30001
  • 最佳答案 8000
  • 用户 6900
  • 常问
  • 回答
  • Marko Smith

    我看不懂措辞

    • 1 个回答
  • Marko Smith

    请求的模块“del”不提供名为“default”的导出

    • 3 个回答
  • Marko Smith

    "!+tab" 在 HTML 的 vs 代码中不起作用

    • 5 个回答
  • Marko Smith

    我正在尝试解决“猜词”的问题。Python

    • 2 个回答
  • Marko Smith

    可以使用哪些命令将当前指针移动到指定的提交而不更改工作目录中的文件?

    • 1 个回答
  • Marko Smith

    Python解析野莓

    • 1 个回答
  • Marko Smith

    问题:“警告:检查最新版本的 pip 时出错。”

    • 2 个回答
  • Marko Smith

    帮助编写一个用值填充变量的循环。解决这个问题

    • 2 个回答
  • Marko Smith

    尽管依赖数组为空,但在渲染上调用了 2 次 useEffect

    • 2 个回答
  • Marko Smith

    数据不通过 Telegram.WebApp.sendData 发送

    • 1 个回答
  • Martin Hope
    Alexandr_TT 2020年新年大赛! 2020-12-20 18:20:21 +0000 UTC
  • Martin Hope
    Alexandr_TT 圣诞树动画 2020-12-23 00:38:08 +0000 UTC
  • Martin Hope
    Air 究竟是什么标识了网站访问者? 2020-11-03 15:49:20 +0000 UTC
  • Martin Hope
    Qwertiy 号码显示 9223372036854775807 2020-07-11 18:16:49 +0000 UTC
  • Martin Hope
    user216109 如何为黑客设下陷阱,或充分击退攻击? 2020-05-10 02:22:52 +0000 UTC
  • Martin Hope
    Qwertiy 并变成3个无穷大 2020-11-06 07:15:57 +0000 UTC
  • Martin Hope
    koks_rs 什么是样板代码? 2020-10-27 15:43:19 +0000 UTC
  • Martin Hope
    Sirop4ik 向 git 提交发布的正确方法是什么? 2020-10-05 00:02:00 +0000 UTC
  • Martin Hope
    faoxis 为什么在这么多示例中函数都称为 foo? 2020-08-15 04:42:49 +0000 UTC
  • Martin Hope
    Pavel Mayorov 如何从事件或回调函数中返回值?或者至少等他们完成。 2020-08-11 16:49:28 +0000 UTC

热门标签

javascript python java php c# c++ html android jquery mysql

Explore

  • 主页
  • 问题
    • 热门问题
    • 最新问题
  • 标签
  • 帮助

Footer

RError.com

关于我们

  • 关于我们
  • 联系我们

Legal Stuff

  • Privacy Policy

帮助

© 2023 RError.com All Rights Reserve   沪ICP备12040472号-5