计算机组成原理:牛顿迭代法(求根的估计)(python实现

今天无意间发现了nympy这个库
于是就想着能够使用它来实现以下牛顿迭代法:
但是发现利用nympy如果写到函数里边,总是出现cannot unpack non-iterable Float object的错误
调试无果后 一气之下 直接都不要函数了 直接写
(生气.jpg
牛顿迭代法流程图如下:
在这里插入图片描述
代码如下:

#coding=gbk;
from sympy import *;
from future import *;
import math;
import sys;
x=symbols("x");
#f=x*exp(x)-1;  #这个地方的表达式可以根据需要进行更改
#注意到这里,就是先将f对x进行不定积分,然后再求导。
#因为一开始input时候f也不过就是个表达式
f=input("请输入函数的表达式:");
f=integrate(f, x);
f=diff(f);

list_in=input("请输入初始值x,极小值sigma以及想要计算的次数: ").split(" ");
x0=float(list_in[0]);
sigma=float(list_in[1]);
N=int(list_in[2]);

if diff(f).subs(x,x0)==0:
    print(x0)
    sys.exit();
account=1;
x_new=x0;
while account<=N:
    x_new=x0-f.subs(x,x0)/diff(f).subs(x,x0);
        #x_new=float(x_new);
    if math.fabs(x_new-x0)<sigma:
        print(x_new)
        sys.exit();
    x0=x_new;
    account+=1;
print("无法解出,不满足条件")

代码运行结果:
注意这里所谓的输入函数表达式,一定是我们用python写代码的那种表达式,而不能将exp(x)写成e^x。
在这里插入图片描述
书上结果:
在这里插入图片描述

遇事不决,可问春风。

  • 2
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值