newtons-method相关内容

牛顿-拉夫森在帕斯卡,结果不是很好

我用Pascal语言实现了牛顿-拉夫森算法。这很奇怪,因为同样的代码在C++中会产生很好的结果(9的结果是3),但在Pascal中9的结果是3.25,为什么呢? 帕斯卡: Program NewtonRaphsonIter(output); {$mode objFPC} function newton_raphson_iter(a: real; p: real; eps: real ..
发布时间:2022-09-02 10:34:37 C/C++开发

二分法(数值分析)

在找到每个根之前进行了多少次递归?还有,哪些是根? 这是我的代码: e=0.000001;f1=@(x) 14.*x.*exp(x-2)-12.*exp(x-2)-7.*x.^3+20.*x.^2-26.*x+12;a=0;c=3;而 abs(c-a)>eb=(c+a)/2;如果 f1(a)*f1(b) 解决方案 二分法的工作原理是取某个初始区间 [a,b] 的端点,并找出区间的哪一半 ..
发布时间:2021-11-27 12:03:37 其他开发

SciPy 优化:Newton-CG vs BFGS vs L-BFGS

我正在使用 Scipy 做一个优化问题,我正在使用一个由顶点和大小为 NNxNN 的键组成的平面网络,连接它的两侧(即使其具有周期性),并最小化能量函数,使其卷曲形成圆柱体.(请参阅下面的链接.) 由于我有函数energy(xyz-position)而且是梯度,我决定使用Scipy手册中推荐的三种方法——Newton-CG,BFGS、L-BFGS-B -- 并比较它们的表现. 我调用优 ..
发布时间:2021-06-12 20:48:03 Python

在Fortran中以牛顿方法传递其他参数

我在实现一种在Fortran程序中调用牛顿方法的方法时遇到了麻烦. 因此,我想使用牛顿方法来解决链接但是,我的程序与上面的示例略有不同.在我的情况下,方程式需要一些其他信息,这些信息是在运行时生成的. subroutine solve(f, fp, x0, x, iters, debug) 这意味着f不仅基于x进行计算,还基于其他一些变量(但x是未知数)计算. 我有一个解决方案,该 ..
发布时间:2020-11-10 06:50:37 其他开发

牛顿拉夫森迭代陷入无限循环

我是这个主题的初学者,无法找出原因:有时程序可以工作,有时不能(在问了问题之后,它根本不想接受我的答案,我可以尽可能多地写下来,但它没有响应,只列出我输入的数字) #include float abszolut (float szam) { float abszoluterteke; if (szam >=0) abszoluter ..
发布时间:2020-06-12 19:25:22 其他开发

Newton Raphson迭代-无法迭代

我不确定这个问题是在这里还是在其他地方(或者根本不在任何地方). 我继承了执行Newton Raphson插值的Fortran 90代码,其中将温度的对数与压力的对数进行插值. 插值类型 t = a ln(p) + b 其中a,b定义为 a = ln(tup/tdwn)/(alogpu - alogpd) 和 b = ln T - a * ln P 这是测 ..
发布时间:2020-06-12 19:25:02 其他开发

如何判断牛顿法是否失败

我正在为无约束的优化问题创建一个基本的牛顿法算法,但是我从算法中得到的结果与我的预期不同。它是一个简单的目标函数,因此很明显该算法应收敛于(1,1)。这是由我之前在这里创建的梯度下降算法所证实的: def grad_descent(x,t,count,幅值): xvalues.append(x) gradvalues.append(np.array([dfx1(x),dfx2(x)] ..
发布时间:2020-06-03 21:39:41 Python

J:牛顿法的默示副词

我在布伦特方法的'addons/math/misc/brent.ijs'实现中发现它是副词.我也想将牛顿的方法构建为副词,但比构建隐式动词要困难得多. 这是牛顿迭代的显式版本: newton_i =: 1 : '] - u % u d.1' 这种用法: 2&o. newton_i^:_ (1) NB. (-: 1p1) must be found 1.5708 ..
发布时间:2020-05-21 20:09:24 其他开发

使用MATLAB编写在二维中实现牛顿方法的函数

我正在尝试编写一个在二维上实现牛顿方法的函数,尽管这样做,我现在必须调整脚本,以便函数的输入参数在列向量中必须为f(x), f(x)的雅可比矩阵,初始猜测x0以及函数f(x)及其雅可比矩阵位于单独的.m文件中的容差. 作为我编写的实现牛顿方法的脚本的示例,我有: n=0; %initialize iteration counter eps=1; ..
发布时间:2020-05-06 13:32:06 其他开发

Matlab中的Newton-Raphson方法

我是Matlab的新手,我需要创建一个函数,该函数以x = a开始逼近,对Newton-Raphson方法进行n次迭代.此开始近似值不算作中间项,另一个要求是需要for循环.我看过其他类似的问题,但就我而言,我不想使用while循环. 这是我的输入应该是的: mynewton(f,a,n) which takes three inputs: f: A function handle ..
发布时间:2020-05-06 12:21:29 其他开发

Matlab中的Newton Raphsons方法?

Newtons-Raphsons方法在Mathematica中易于实现,但在Matlab中似乎有点困难.我不知道是否可以将函数传递给函数,以及如何将导数用作函数. newtonRaphson[f_, n_, guess_] := If[n == 0, guess, newtonRaphson[f, n - 1, guess - f[guess]/f'[guess]]] newtonRap ..
发布时间:2020-05-06 12:20:44 其他开发

使用牛顿法求平方根(错误!)

我正在使用Python中的牛顿的猜测和检查方法来完成一个数学问题,该问题近似于数字的平方根.用户应该输入一个数字,对该数字的初始猜测,以及他们想返回多少次才能检查答案.为了使事情变得更容易并了解Python(几个月前我才刚刚开始学习该语言),我将其分解为许多较小的函数.现在的问题是,我在调用每个函数并传递数字时遇到麻烦. 这是我的代码,提供注释以帮助您(每个功能按使用顺序排列): # ..
发布时间:2020-05-06 10:35:22 Python

RGB和RYB颜色空间之间的转换

我目前正在尝试在RGB(红色,绿色,蓝色)颜色空间和RYB(红色,黄色,蓝色)颜色空间之间转换颜色,然后再次转换. 基于下面的详细信息,我能够使用三线性插值将RYB转换为RGB-其中参数权重(s,t,u)是RYB颜色,而立方体的顶点是3d点在RGB空间中. 绘画启发的颜色混合和合成以实现可视化-Gossett和Chen -第2.1节-实现细节 我的困难在于逆转转换过程. 第 ..
发布时间:2020-05-06 10:32:00 其他开发