不动点迭代方程收敛判据及MATLAB实现
%函数 fx=x.^3-3.*x+1;
%收敛性判断
xl=0;
xr=1;
%不动点所在区间
x=xl:0.001:xr;
y=(x.^3+1)./3;%需要修改判断的方程
if max(y)>xr||min(y)<xldisp('迭代方程不收敛,定理1不符合');return
end
for x0=xl:0.001:xrfor y0=xl:0.001:xrfor L=0.001:0.001:0.999if abs((x0.^3+1)./3-(y0.^3+1)./3)<=L*abs(x0-y0)%需要修改判断的方程disp('迭代方程收敛');breakelsedisp('迭代方程不收敛,定理2不符合')returnendendbreakendbreak
end
%设置不动点方程 x=(x.^3+1)./3
%设置初始值x,最大迭代次数N,终止迭代条件p
k=0;
while k<Nk=k+1;x0=(x.^3+1)./3;%迭代方程式if abs(x0-x)<=pbreakendx=x0;
end
if k>Ndisp(['迭代次数= ',num2str(k),',算法超过最大迭代次数。']);
elsedisp(['迭代次数= ',num2str(k)]);
end
不动点迭代方程收敛判据及MATLAB实现相关推荐
- 非线性方程(组):一维非线性方程(一)二分法、不动点迭代、牛顿法 [MATLAB]...
1. 二分法(Bisection) 1) 原理 [介值定理] 对于连续的一元非线性函数,若其在两个点的取值异号,则在两点间必定存在零点. [迭代流程] 若左右两端取值不同,则取其中点,求其函数值,取中 ...
- 不动点迭代求解方程数值解
求解方程 2x−x3=0 解:可以采用不动点迭代的方式求出数值解. 其不动点是 x=3log2x ,故可以采用下式进行不动点迭代: xn+1=3log2xn 初始值为 x0=2 ,迭代终止条件是 |x ...
- 不动点迭代(Fixed Point Iteration)
题目:不动点迭代(Fixed Point Iteration) 本篇介绍不动点迭代(Fixed Point Iteration).之所以学习不动点迭代是由于近来看到了FPC算法,即Fixed Poin ...
- 最优化方法 26:不动点迭代
前面讲了很多具体的算法,比如梯度.次梯度.近似点梯度.加速近似点梯度.PPA.DR方法.ADMM.ALM等,对这些方法的迭代过程有了一些了解.这一节则主要是针对算法的收敛性进行分析,试图从一个更加抽象 ...
- 不动点迭代以及其收敛性
不动点迭代以及其收敛性 对于迭代的理解 不动点迭代 迭代的收敛性 区间收敛 局部收敛 对于迭代的理解 所谓迭代就是反复使用执行某一个过程,并且用本次执行该过程的结果作为下一次执行的起点,不断推进, ...
- 数值分析——求方程解的不动点迭代法和斯特芬森法(Python实现)
一.不动点迭代法求方程的解 import sympy #迭代的方程为 f(x)=pow(x,3)-x-1 def psi(x):return pow(x+1,1/3)def dif(x0): #求导函 ...
- matlab不动点迭代代码
%不动点迭代function [xc,k] = fpi(g,x0,tol,N) x1=g(x0); k=1; while abs(x1-x0) > tolx0 =x1;x1=g(x0);if k ...
- 用matlab求不动点迭代,matlab实现不动点迭代、牛顿法、割线法
不动点迭代 function xc = fpi( g, x0, tol ) x(1) = x0; i = 1; while 1 x(i + 1) = g(x(i)); if(abs(x(i+1) - ...
- 不动点迭代 开平方 Excel演示
不动点迭代的精髓在于可以从任意值收敛. 难点在于配这个方程.详细技巧请参阅 数值分析 内容. X 代表初始值
最新文章
- Javascript鼠标滚轮事件兼容写法
- 在Intellij Idea中使用jstl标签库
- NumberUtils源码分析
- 链表LinkedList的Java实现
- java就是我的幸福,嫁给幸福
- 南方科技大学计算机交换生,国际合作 – 合作交流分类 – 南方科技大学生物医学工程系...
- 美观又实用,10 款强大的开源 Javascript 图表库
- SugarNMS分布式网络运维部署攻略
- 关于get请求的乱码问题
- 算法合集之《信息学中守恒法的应用》
- TP-Link TD-W89841N 增强型无线路由器快速设置指南
- 计算机绘图培训心得,学习cad制图心得体会
- 【虹科分享】什么是 RFC 2544?网络设备的性能基准测试方法
- webfreer去广告
- Android开发自定义相机,自定义拍照界面
- 基于tensorflow的手写数字识别
- 周报8.22-8.28
- MySQL修改数据库名字
- pycharm报错warning: iCCP: known incorrect sRGB profile
- 安装tensorflow