【数值分析】二分法求解-Python实现

最新推荐文章于 2025-04-12 18:56:39 发布

大菜市场

最新推荐文章于 2025-04-12 18:56:39 发布

阅读量2.5k

点赞数 2

分类专栏：数值分析 Python

本文链接：https://blog.youkuaiyun.com/qq_36440163/article/details/88761101

版权

本文档介绍了使用Python3.6及matplotlib、numpy库实现数值分析中的二分法，提供了代码示例，旨在帮助读者理解二分法的原理与应用。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

本篇为数值分析课程代码实现-二分法的实现

仅供参考

配置环境

Python3.6
matplotlib
numpy

# -*- coding: utf-8 -*-
"""
Created on Thu Nov 15 22:28:22 2018
二分法
@author: hhuaf
"""
import numpy as np
import matplotlib.pyplot as plt
# input
'''
a :下限
b :上限
theta:阈值
'''
#a=input('输入下限\n:')
#b=input('输入上限\n:')
#a=float(a)
#b=float(b)
a=1
b=5
theta=0.05

#可以显示中文
plt.rcParams["font.sans-serif"] = ["SimHei"]
plt.rcParams['axes.unicode_minus'] = False

# 设置风格
plt.style.use('ggplot')

# 定义函数，构造数值

fun = lambda x: x**3-3*x+1
#fun = lambda x: x**2-1
x = np.arange(a,b,0.05)
y = fun(x)

# 画出图像

fig = plt.figure(figsize = (15, 10))
plt.subplot(2,3,1)
plt.xlim(a, b)
plt.xlabel('X')
plt.ylabel('Y')
plt.title('$f(x)=x^3-3x+1$ 图像')