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 。



目录
相关文章
|
1月前
|
算法 C语言
C语言求方程解的隐藏技能,你不会就亏大了!
C语言求方程解的隐藏技能,你不会就亏大了!
|
6月前
|
C语言
C 语言实例 - 一元二次方程
求一元二次方程:ax2+bx+c=0 的根。输入三个实数a,b,c的值,且a不等于0。
23 1
|
7月前
|
算法 C语言
C语言实现牛顿迭代法解方程
C语言实现牛顿迭代法解方程
180 0
|
C语言
造轮子之-C语言实现ArrayList
造轮子之-C语言实现ArrayList
06【C语言 & 趣味算法】牛顿迭代法求方程根(可回看)
06【C语言 & 趣味算法】牛顿迭代法求方程根(可回看)
06【C语言 & 趣味算法】牛顿迭代法求方程根(可回看)
|
C语言
C语言实现学生成绩管理系统
本文提供一例C语言实现的命令行学生信息管理系统,供初学者参考。
290 0
|
存储 Linux C语言
生产者消费者模式保姆级教程 (阻塞队列解除耦合性) 一文帮你从C语言版本到C++ 版本, 从理论到实现 (一文足以)
生产者消费者模式保姆级教程 (阻塞队列解除耦合性) 一文帮你从C语言版本到C++ 版本, 从理论到实现 (一文足以)
生产者消费者模式保姆级教程 (阻塞队列解除耦合性) 一文帮你从C语言版本到C++ 版本, 从理论到实现 (一文足以)
|
存储 C语言
c语言实现扫雷(含循环递归展开)
本笔记通过c语言实现扫雷小游戏(包含递归展开) 游戏实现逻辑位于test.c文件,整个游戏头文件位于game.h,游戏进程的具体操作于game.c中实现。
149 0
c语言实现扫雷(含循环递归展开)
|
C语言
c语言实现三子棋(内含阅读思路,简单易实现)
本文如果按顺序来阅读可能不太好接受,建议阅读顺序为,由test.c的逻辑顺序读下去,遇见具体函数的实现跳转到game.c中来理解
114 0
c语言实现三子棋(内含阅读思路,简单易实现)
|
SQL 关系型数据库 MySQL
C语言连接并实现对MySQL的增删改查
C语言访问数据库并不如Java、Python那般容易。本文介绍C语言连接并实现对MySQL的增删改查的方法。
436 0