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.0001m=(a+b)/2;if f(m)==0break;endif f(a)*f(m)>0a=m;elseb=m;end
end
m

单点弦截法

单点弦截法适用于导数值计算较为困难，且f(x)在区间(a,b)内凹向不变的函数。

%p228task3_2
%单点弦截法求x^2+x-3=0近似解
clc,clear
f=@(x)x^2+x-3
ezplot(f,[-3,2])
grid%画出网格线
%根在[-2.5,2],[1,1.5]之间
a=-2.5;
b=-2;digits 8
while abs(f(a))>0.00001a=b-f(b)*(b-a)/(f(b)-f(a));
end
a

牛顿切线迭代法

牛顿切点迭代法的条件是：f(x)在闭区间[a,b]上有二阶导数，f(a)*f(b)<0，且f’(x)与f”(x)在[a,b]上不变号。

%p228task4
%牛顿切线迭代法求x^2-2*x*exp(x)+exp(-x)=0的根
clc,clear
syms x;
f=@(x)x^2-2*x*exp(x)+exp(-x)
d1f=diff(f(x))
d2f=diff(d1f)
% ezplot(@(x)0)
figure(1)
ezplot(f(x),[0,0.5])
grid
figure(2)
ezplot(d1f,[0,0.5])
grid
figure(3)
ezplot(d2f,[0,0.5])
grid
%作图可知，一阶导数、二阶导数都小于0；符合牛顿切线弦截法的要求
%根在0.3到0.35之间
a=0.3;
b=0.35
if subs(f,a)*subs(d2f,a)>0x0=a;
elsex0=b;
end
dlt=1.0e-5;%设置精度
m=min(abs(subs(d1f,a)),abs(subs(d1f,b)));%m的作用就是使达到精度要求后停止循环
k=1;
while abs(subs(f,x0))>m*dltx1=x0-subs(f,x0)/subs(d1f,x0);x0=x1;fprintf('k=%d,x=%.7f\n',k,x0);k=k+1;
end
solve(f)

此处用到一个新知识，用subs求符号函数在某一点的函数值，
如:

>> f=@(x)x^2-2*x*exp(x)+exp(-x)f =包含以下值的 function_handle:@(x)x^2-2*x*exp(x)+exp(-x)>> subs(f,2)ans =exp(-2) - 4*exp(2) + 4>> vpa(subs(f,2))ans =-25.420889112485988217027710347328

可见subs的用法：subs(f,x0)可求出在x0处的f的值，但对于sin,exp等函数，不会直接显示近似数值，调用vpa函数即可令其显示数值。

matlab二分法，单点弦截法，牛顿切线迭代法相关推荐

非线性方程求解——单点弦截法
单点弦截法是双电弦截法的特殊情况,迭代的情况可能在有根区间的两个端点之间变换.如果迭代的结果始终在一侧的端点,这个端点会收敛于精确根,而有根区间另一个端点固定不变,那么双点弦截法变为单点弦截法. 单点 ...
用弦截法求方程的根matlab,matlab 语言用弦截法任意实数方程求实根
满意答案 dwgg2n0das8 2013.04.26 采纳率:45% 等级:12 已帮助:8202人弦截法这个方法一般用作学习,实际用的很少.这里我提供一个较完整的弦截法求根的函数.func ...
单点、双点弦截法求解方程根
弦截法: (a) 用牛顿法解方程f(x)=0,虽然在单根附近具有较快的收敛速度,但它有个明显的缺点,就是需要计算导数f'(x),当f(x)比较复杂时,计算f'(x)可能有困难. (b)弦截法和牛顿迭代 ...
牛顿迭代法是一种速度很快的迭代方法，但是它需要预先求得导函数。若用差商代替导数，可得下列弦截法
牛顿迭代法是一种速度很快的迭代方法,但是它需要预先求得导函数.若用差商代替导数,可得下列弦截法这一迭代法需要两个初值,迭代过程和牛顿法类似,当时停止迭代.编写程序实现弦截法,要求输出方程的根.函数在 ...
MATLAB弦截法求解非线性方程
MATLAB弦截法求解非线性方程用Newton法解非线性方程时,当f比较复杂时就比较难以实现了,弦解法的好处可以用f(x)在两点上的值构造一次插值函数来回避微商的计算.其迭代格式: 1.弦截法的MA ...
数值计算大作业：非线性方程求根（二分法、牛顿法、弦截法在Matlab实现）
作为研究生的入门课,数值计算的大作业算是所有研究生开学的重要编程作业. 我把二分法.牛顿法.弦截法求解非线性方程求根的数值计算作业在MATLAB中编程实现.具体的程序详细标注后放在文章附录了,算法数学 ...
数值分析-基本算法流程图：二分法、简单迭代法、复化梯形公式、弦截法、辛普森
1.二分法 2.复化Simpson公式 3.复化梯形公式 4.简单迭代法 5.牛顿迭代法 6.弦截法
弦截法求方程根例题c语言,弦截法求方程根.ppt
弦截法求方程根弦截法求方程弦截法求方程根弦截法求方程根 12 §3 迭代收敛的加速法取g(x)=x3-1 , 则: 程序设计 function [f,k]= Steffensen(eps,x0) % ...
弦截法c语言程序,高数介质定理——弦截法求根代码实践（C语言）
在高等数学中,我们一开始接触概念时就接受了ε-δ(epsilon-delta)语言的洗礼,但即使到课程的结束,许多人依然会对各种抽象的数学符号.定理证明感到无所适从,我也不例外,尽管在写这篇博客以前已 ...

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

二分法

单点弦截法

牛顿切线迭代法

matlab二分法，单点弦截法，牛顿切线迭代法相关推荐

最新文章

热门文章