matlab二分法,单点弦截法,牛顿切线迭代法

二分法

%p222task2_3
%二分法求f=@(x)1-x-sin(x)零点
clc,clear;
f=@(x)1-x-sin(x)
b=1;a=0;
f(0)
f(1)
ezplot(f,[0,1])
hold on
ezplot(@(x)0,[0,1])

while b-a>0.0001
    m=(a+b)/2;
    if f(m)==0
        
  • 2
    点赞
  • 29
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
1. 目的: (1)通过采用牛顿迭代法弦截法二分法根的程序设计,使学生更加系统地理和掌握C语言函数间参数传递方法、数组和指针的应用等编程技巧。培养学生综合利用C语言进行科学计算,使学生将所学知识转化为分析和设计数学中的实际问题的能力,学会查资料和工具书。 (2)提高学生建立程序文档、归纳总结的能力。 (3)进一步巩固和灵活运用先修课程《计算机文化基础》有关文字处理、图表分析、数据归整、应用软件之间图表、数据共享等信息技术处理的综合能力。 2. 基本要: (1)要用模块化设计和C语言的思想来完成程序的设计; (2)要分别编写牛顿迭代法弦截法二分法根的函数,分别存到不同的.CPP文件中; (3)在VC++6.0环境中,学会调试程序的方法,及时查究错误,独立调试完成。 (4)程序调试通过后,完成程序文档的整理,加必要的注释。 一般一元方程,常用采用的方法有:牛顿迭代法弦截法二分法等。 牛顿迭代法根 〖〖f(x)=a〗_0 x〗^n 〖〖 + a〗_1 x〗^(n-1) +⋯+〖 a〗_(n-2) x^2 +〖 a〗_(n-1) x +〖 a〗_n=0 f(x)在〖 x〗_0附近的根。 计算公式:〖 x〗_(n+1)=〖 x〗_n- f(〖 x〗_n )/(f(〖 x〗_n)) ́ 精度:ε=|〖 x〗_(n+1)-〖 x〗_n|<1.0e-m ,m=6。 牛顿迭代法的根:满足精度的〖 x〗_n 二分法 任取两点〖 x〗_1和〖 x〗_2,判断(〖 x〗_1, 〖 x〗_2)有无实根。如下图所示,如果f(〖 x〗_1 )和f(〖 x〗_2 )符号相反,说明(〖 x〗_1, 〖 x〗_2)之间有一实根。取(〖 x〗_1, 〖 x〗_2)的中点x,检查f(x)和f(〖 x〗_1 )是否同符号,如果不同号,说明实根在(〖 x〗_1,x)区间,x作为新的〖 x〗_2,舍弃(x, 〖 x〗_2)区间;若同号,则实根在(x, 〖 x〗_2)区间,x作为新的〖 x〗_1, 舍弃(〖 x〗_1,x)区间。再根据新的〖 x〗_1 、 〖 x〗_2,找中点,重复上述步骤。直到|〖 x〗_1-〖 x〗_2|〖<10〗^(-6)时,x =(〖 x〗_1+〖 x〗_2)/2为所。 (3)弦截法 取f(〖 x〗_1 )与f(〖 x〗_2 )连线与x轴的交点x,从(〖 x〗_1, x)和(x, 〖 x〗_2)两个区间中取舍的方法与二分法相同。 计算公式为: 判断f(〖 x〗_1 )与f(〖 x〗_2 )是否同符号的方法与二分法采用的方法相同。直到先后两次出的x的值之差小于〖10〗^(-6)为止。 分别用牛顿迭代法弦截法二分法下列方程的根,分析比较各种方法的迭代次数及精度。 〖f(x)=x〗^3 〖- 2x〗^2 +7x +4=0 牛顿迭代法的初值:x=0.5; 弦截法〖 x〗_1,〖 x〗_2的初值:-1,1 二分法〖 x〗_1,〖 x〗_2的初值:-1,0 精度要:|〖 x〗_1-〖 x〗_2| 〖<10〗^(-6)
二元牛顿迭代法 Matlab 实现: 假设要方程组 f(x,y)=0,初始值为 (x0,y0)。二元牛顿迭代法的迭代公式为: [x(k+1);y(k+1)] = [x(k);y(k)] - [J(x(k),y(k))]^-1 * [f(x(k),y(k))] 其中,J(x,y) 为 f(x,y) 的雅克比矩阵,即: J(x,y) = [df/dx(x,y), df/dy(x,y)] 按照上述公式,可以用 Matlab 编写二元牛顿迭代法的代码: function [x,y] = newton2(f,J,x0,y0,tol,maxiter) % f: 二元函数句柄 % J: 二元函数的雅克比矩阵句柄 % x0, y0: 初始值 % tol: 容忍误差 % maxiter: 最大迭代次数 % 初始化 k = 0; x = x0; y = y0; while k < maxiter % 计算 J 和 f Jk = J(x,y); fk = f(x,y); % 计算迭代步长 s = - Jk \ fk; % 更新 x 和 y x = x + s(1); y = y + s(2); % 判断是否达到容忍误差 if norm(fk) < tol break; end % 更新迭代次数 k = k + 1; end end Matlab 编写二分法牛顿迭代法二分法 Matlab 实现: 假设要方程 f(x)=0,在区间 [a,b] 内有唯一二分法的思路是不断缩小区间 [a,b] 的范围,直到满足容忍误差要。具体实现如下: function [x] = bisection(f,a,b,tol) % f: 函数句柄 % a,b: 区间端点 % tol: 容忍误差 % 初始化 fa = f(a); fb = f(b); % 判断是否满足 f(a)*f(b)<0 if fa*fb >= 0 error('Error: f(a)*f(b)>=0'); end while (b-a)/2 > tol % 计算中点 c = (a + b) / 2; fc = f(c); % 判断在哪个区间 if fc == 0 x = c; return; elseif fa*fc < 0 b = c; fb = fc; else a = c; fa = fc; end end % 返回 x = (a+b) / 2; end 牛顿迭代法 Matlab 实现: 假设要方程 f(x)=0,初始值为 x0。牛顿迭代法的迭代公式为: x(k+1) = x(k) - f(x(k))/f'(x(k)) 按照上述公式,可以用 Matlab 编写牛顿迭代法的代码: function [x] = newton(f,df,x0,tol,maxiter) % f: 函数句柄 % df: 导数函数句柄 % x0: 初始值 % tol: 容忍误差 % maxiter: 最大迭代次数 % 初始化 k = 0; x = x0; while k < maxiter % 计算 f 和 df fk = f(x); dfk = df(x); % 计算迭代步长 s = - fk / dfk; % 更新 x x = x + s; % 判断是否达到容忍误差 if abs(fk) < tol break; end % 更新迭代次数 k = k + 1; end end

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值