第十一周—牛顿迭代法求根


/* 

Copyright(c)2016,烟台大学计算机学院 

 All rights reserved. 

文件名称:tset.cpp 

作者:尤文敏 

完成日期:2016年11月12日 

版本号:v1.0 

题目描述 :

 牛顿迭代法求根

问题及代码:

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
int main()

{
    double solut(double ,double ,double ,double );
    double a,b,c,d;
    scanf("%lf %lf %lf %lf",&a,&b,&c,&d);
    printf("%.2f",solut(a,b,c,d));
    return 0;
}
double solut(double a,double b,double c,double d )
{
    double x=1,x0;
   do
    {
        x0=x;
        x=x0-(((a*x+b)*x+c)*x+d)/((3*a*x+2*b)*x+c);
    }
     while(fabs(x-x0)>=1e-5);
    return x;
}
运行结果:
知识点总结:
  运用牛顿迭代法求根时,要判断两者之间的绝对值是否满足精确度,所以运用到fabs这个函数。
学习心得:
  我在运用fabs函数的时候一开始忘记了子啊前面要加上#include <math.h>,因此大家千万不要忘记喽!


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
### 回答1: 牛顿迭代法是一种求解非线性方程组的方法,可以使用Matlab进行实现。具体步骤如下: 1. 定义非线性方程组的函数表达式,例如: function F = myfun(x) F(1) = x(1)^2 + x(2)^2 - 1; F(2) = x(1) - x(2)^2; 2. 定义牛顿迭代法的迭代公式,例如: function [x, k] = newton(fun, x0, tol, maxiter) k = 0; x = x0; while k < maxiter F = fun(x); J = jacobian(fun, x); dx = -J\F'; x = x + dx'; if norm(F) < tol break; end k = k + 1; end 3. 调用函数进行求解,例如: [x, k] = newton(@myfun, [1, 1], 1e-6, 100); 其中,@myfun表示使用myfun函数进行求解,[1, 1]表示初始值,1e-6表示误差容限,100表示最大迭代次数。 4. 输出结果,例如: disp(['Solution: x = [', num2str(x(1)), ', ', num2str(x(2)), ']']); disp(['Iterations: ', num2str(k)]); 这样就可以使用Matlab实现牛顿迭代法求解非线性方程组了。 ### 回答2: 牛顿迭代法是求解非线性方程组的一种有效方法,它通过一系列迭代公式逼近方程组的根。在matlab中,我们可以使用该方法求解非线性方程组。 首先,我们需要定义一个函数句柄来表示非线性方程组,比如: ``` f = @(x) [x(1)^2 + x(2)^2 - 4; x(1)*x(2) - 1]; ``` 这里定义的函数句柄f表示一个含有两个未知变量的非线性方程组,其中第一个方程表示一个以原点为圆心,半径为2的圆,第二个方程表示一个过点(1,1)的直线与x轴的交点。 接下来,我们需要设定初始值x0和迭代终止条件tol,比如: ``` x0 = [1;1]; tol = 1e-6; ``` x0表示迭代的起点,tol表示迭代的终止条件,通常设置为一个较小的正数,如1e-6,表示当两个相邻迭代结果的差值小于等于1e-6时停止迭代。 然后,我们可以使用牛顿迭代公式对方程组进行迭代求解,具体公式如下: ``` x = x - J\f(x); ``` 其中,x表示当前迭代点的值,J表示方程组f在当前迭代点的雅可比矩阵,f(x)表示当前迭代点对应的方程组的函数值,\表示矩阵的左除,即求解如下线性方程组: ``` J*dx = -f(x) ``` 其中,dx表示当前迭代点相对于上一个迭代点的增量,即: ``` dx = x - x_prev; ``` 我们可以使用一个循环来实现牛顿迭代的过程,如下: ``` x = x0; x_prev = x0; while norm(x - x_prev) > tol J = [2*x(1) 2*x(2); x(2) x(1)]; dx = J\-f(x); x_prev = x; x = x + dx; end ``` 其中,norm函数用来计算向量的2-范数,表示向量的长度。迭代过程中,我们先计算当前点的雅可比矩阵J和函数值f(x),然后求解线性方程组得到增量dx,最后更新迭代点的值。 最后,我们可以使用disp函数输出最终的迭代结果,如下: ``` disp(['x = (' num2str(x(1)) ', ' num2str(x(2)) ')']); ``` 通过以上步骤,我们就可以成功地使用牛顿迭代法求解非线性方程组。 ### 回答3: 牛顿迭代法是一种求解非线性方程组的常用方法,它是基于牛顿-拉夫逊迭代法的思想,通过不断迭代逼近非线性方程组的解。在matlab中,可以使用牛顿迭代法求解非线性方程组,其步骤如下: 1. 首先定义非线性方程组的函数表达式,如:f = @(x) [x(1)^2+x(2)-11;x(1)+x(2)^2-7]; 2. 然后定义非线性方程组的雅可比矩阵,即f的偏导数矩阵,如:df = @(x) [2*x(1),1;1,2*x(2)]; 3. 初始化解向量,如:x = [1;1]; 4. 设置收敛条件,如:tol = 1e-6; 5. 开始迭代,如:for i=1:100 f_val = f(x); df_val = df(x); dx = -df_val\f_val; x = x + dx; if(norm(dx)<tol) break; end end 以上就是用牛顿迭代法求解非线性方程组的基本步骤,通过不断迭代可以逼近方程组的解。需要注意的是,初始解向量的设置、收敛条件的确定以及迭代次数的控制都会影响迭代结果的精度和速度,需要根据具体需要进行调整。此外,在matlab中还可以使用fsolve函数来实现牛顿迭代法求解非线性方程组,其使用方法更加方便快捷。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值