C语言实现牛顿迭代法解方程详解

首页 > 软件编程 > C 语言 更新日期: 2017-09-09
这篇文章主要介绍了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)。

我们来看一副从网上找到的图:

C语言实现牛顿迭代法解方程详解

接下来,我们来看一个例子:

我们还是直接上代码:

例子:用牛顿迭代法求下列方程在值等于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语言实现牛顿迭代法解方程详解

 感谢阅读,希望能帮助到大家,谢谢大家对本站的支持!


> 本站内容系网友提交或本网编辑转载,其目的在于传递更多信息,并不代表本网赞同其观点和对其真实性负责。如涉及作品内容、版权和其它问题,请及时与本网联系,我们将在第一时间删除内容!

相关文章
  • JavaScript语言对Unicode字符集的支持详解
    JavaScript语言对Unicode字符集的支持详解
    这篇文章主要介绍了JavaScript语言对Unicode字符集的支持详解,需要的朋友可以参考下上个月,我做了一次分享,详细介绍了Unicode字符集,以及JavaScript语言对它的支持.下面就是这次分享的讲稿. 一.Unicode是什么? Unicode源于一个很简单的想法:将全世界所有的字符包含在一个集合里,计算机只要支持这一个字符集,就能显示所有的 ...
  • c语言中基于随机函数的使用详解
    本篇文章对c语言的随机函数进行了详细的分析介绍.需要的朋友参考下在C语言中,rand()函数可以用来产生随机数,但是这不是真真意义上的随机数,是一个伪随机数,是根据一个数,我们可以称它为种子,为基准以某个递推公式推算出来的一系数,当这系列数很大的时候,就符合正态公布,从而相当于产生了随机数,但这不是真正的随机数,当计算机正常开机后,这个种子的值是定了的,除非 ...
  • 详解C语言中条件判断语句if和switch的用法
    这篇文章主要介绍了详解C语言中条件判断语句if和switch的用法,是C语言入门学习中的基础知识,需要的朋友可以参考下if 语句 用 if 语句可以构成分支结构,它根据给的条件进行判定,以决定执行哪个分支程序段. C 语言的 if 语句有三种基本形式 第一种形式: if(条件表达式) { 语句1: } if(条件表达式) { 语句1: } 这种形式运行顺序为 ...
  • c语言stack(栈)和heap(堆)的使用详解
    c语言stack(栈)和heap(堆)的使用详解
    这篇文章主要介绍了c语言stack栈和heap堆的使用详解,需要的朋友可以参考下一.预备知识-程序的内存分配 一个由C/C++编译的程序占用的内存分为以下几个部分 1.栈区(stack)-由编译器自动分配释放,存放函数的参数值,局部变量的值等.其操作方式类似于数据结构中的栈.2.堆区(heap)-一般由程序员分配释放,若程序员不释放,程序结束时可能由OS回收 ...
  • 详解C语言结构体中的函数指针
    这篇文章主要介绍了详解C语言结构体中的函数指针,文中对函数指针的基本概念也有讲解,需要的朋友可以参考下结构体是由一系列具有相同类型或不同类型的数据构成的数据集合.所以,标准C中的结构体是不允许包含成员函数的,当然C++中的结构体对此进行了扩展.那么,我们在C语言的结构体中,只能通过定义函数指针的方式,用函数指针指向相应函数,以此达到调用函数的目的. 函数指针 ...
  • 详解C语言中的内存四区模型及结构体对内存的使用
    这篇文章主要介绍了C语言中的内存四区模型及结构体对内存的使用,包括结构体中内存泄漏情况的注意点提醒,需要的朋友可以参考下内存四区1.代码区代码区code,程序被操作系统加载到内存的时候,所有的可执行代码都加载到代码区,也叫代码段,这块内存是不可以在运行期间修改的.2.静态区所有的全局变量以及程序中的静态变量都存储到静态区.3.栈区栈stack是一种先进后出的 ...
猜你喜欢