C语言实现牛顿迭代法解方程详解
C语言实现牛顿迭代法解方程详解
利用迭代算法解决问题，需要做好以下三个方面的工作：
一、确定迭代变量
在可以用迭代算法解决的问题中，我们可以确定至少存在一个可直接或间接地不断由旧值递推出新值的变量，这个变量就是迭代变量。
二、建立迭代关系式
所谓迭代关系式，指如何从变量的前一个值推出其下一个值的公式（或关系）。迭代关系式的建立是解决迭代问题的关键，通常可以使用递推或倒推的方法来完成。
三、对迭代过程进行控制
在什么时候结束迭代过程？这是编写迭代程序必须考虑的问题。不能让迭代过程无休止地执行下去。迭代过程的控制通常可分为两种情况：一种是所需的迭代次数是个确定的值，可以计算出来；另一种是所需的迭代次数无法确定。对于前一种情况，可以构建一个固定次数的循环来实现对迭代过程的控制；对于后一种情况，需要进一步分析得出可用来结束迭代过程的条件。
接下来，我介绍一种迭代算法的典型案例----牛顿-拉夫逊（拉弗森）方法
牛顿-拉夫逊（拉弗森）方法，又称牛顿迭代法，也称牛顿切线法：先任意设定一个与真实的根接近的值x0作为第一次近似根，由x0求出f(x0),过(x0，f(x0))点做f(x)的切线，交x轴于x1，把它作为第二次近似根，再由x1求出f(x1),过(x1，f(x1))点做f(x)的切线，交x轴于x2，……如此继续下去，直到足够接近（比如|x- x0|<1e-6时）真正的根x*为止。
而f '(x0)=f(x0)/( x1- x0)
所以 x1= x0- f(x0)/ f ' (x0)。
我们来看一副从网上找到的图:
接下来，我们来看一个例子:
我们还是直接上代码:
例子:用牛顿迭代法求下列方程在值等于2.0附近的根：2x3-4x2+3x-6=0。
#include <stdio.h>
#include <math.h>
int main(void)
{
float x,x0,f,f1;
x = 2.0;
do{
x0=x;
f=2*x0*x0*x0-4*x0*x0+3*x0-6;
f1=6*x0*x0-8*x0+3;
x=x0-f/f1;
//函数fabs：求浮点数x的绝对值
//说明：计算|x|, 当x不为负时返回 x，否则返回 -x
}while(fabs(x-x0)>=1e-5);
printf ("%f\n",x);
return 0 ;
}
执行结果:
当x=1.5时，方程2x3-4x2+3x-6=0。附近的根为2.000000 。
感谢阅读，希望能帮助到大家，谢谢大家对本站的支持！
相关文章
这篇文章主要介绍了C++11的右值引用的具体使用，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2020-02-02
这篇文章主要为大家详细介绍了C语言实现弹跳小球，设置了小球分数及过关难度，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2020-03-03
哈夫曼编码是一种编码方式，又称“霍夫曼编码”，其是可变字长的编码（VCL）的一种，这篇文章主要介绍了基于C语言利用哈夫曼树实现文件压缩,需要的朋友可以参考下2021-08-08
这篇文章主要为大家详细介绍了C语言实现个人财务管理，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2021-11-11
pthread是C++98接口且只支持Linux,使用时需要包含头文件#include <pthread.h>,编译时需要链接pthread库,其中p是POSIX的缩写,而POSIX是Portable Operating System Interface的缩写,这篇文章主要介绍了C++ pthread简介,需要的朋友可以参考下2024-05-05
C++矩阵运算矩阵运算包括矩阵相加、相减、相乘、转置、求逆矩阵等等,用计算机程序实现矩阵运算的方法算法很多,这篇文章主要给大家介绍了关于C/C++如何实现两矩阵相乘之模拟法的相关资料,需要的朋友可以参考下2023-02-02
这篇文章主要为大家详细介绍了如何利用C语言自制华容道游戏，文中的示例代码讲解详细，对我们学习C语言有一定的帮助，感兴趣的小伙伴可以跟随小编一起了解一下2023-03-03
这篇文章主要介绍了数据结构 双机调度问题的实例详解的相关资料,双机调度问题,又称独立任务最优调度：用两台处理机A和B处理n个作业的实例，需要的朋友可以参考下2017-08-08
大家好，本篇文章主要讲的是C++实现堆排序实例介绍，感兴趣的同学赶快来看一看吧，对你有帮助的话记得收藏一下，方便下次浏览2021-12-12
这篇文章主要为大家详细介绍了C++如何利用类实现矩阵的数乘,乘法以及点乘，文中的示例代码讲解详细，对我们学习C++有一定帮助，需要的可以参考一下2022-11-11
最新评论