不动点迭代法的matlab实现
不动点迭代法(Fixed point iteration method)
原理网上已经很多了,这里不再赘述。
需要注意的是不动点迭代法的使用条件。
举个简单的例子,利用不动点法求解方程f(x) = 1 + 0.5*sin(x) − x = 0
在(1,2)区间的根。
令g(x)=1+0.5*sin(x)
g = @(x) fun2(x);
% Initialization
x0 = 0;
tol = 1e-5;
maxIter = 40;% Test biSection function
[xStar,xRoot] = fixPoint(g,x0,tol,maxIter);
fprintf('The fixed point is: %d\n',xStar);
fprintf('The root of the equation is: %d\n',xRoot);function [xStar,xRoot] = fixPoint(fun2,x0,tol,maxIter)
% Inputs:
% fun2: a function handle, standing for the function written above
% x0: the initial guess of the fixed point
% tol: the tolerance within which the program can stop
% maxIter: the maximum number of iterations the program is allowed to run
% Outputs:
% xStar: the numerical value of the fixed point
% xRoot: the numerical value of the rootx = zeros(maxIter,1);x(1) = fun2(x0);i = 1;while abs(fun2(x(i))-x(i))>tol && i<maxIterx(i+1) = fun2(x(i));xStar = x(i);i = i+1;endxRoot = x(i);endfunction gx = fun2(x)gx = 1+0.5*sin(x);
end
不动点迭代法的matlab实现相关推荐
- 不动点迭代法matlab视频,不动点迭代法matlab
不动点迭代法的 MATLAB 程序代码如下: Function [root,n]=StablePoint(f,x0,eps) %用不动点迭代法求函数的一个零点 %初始迭代向量:x0 %根的精度:eps ...
- matlab牛顿法求区间根程序,MATLAB用二分法、不动点迭代法及Newton迭代(切线)法求非线性方程的根...
一.实验原理 二.实验步骤 三.实验过程 1.(程序) (1)二分法:求 在区间(1,2)之间的根,取 (a)bipart.m: function [x,m]=bipart(fun,a0,b0,to ...
- 重根的二阶迭代法matlab,MATLAB用二分法、不动点迭代法及Newton迭代(切线)法求非线性方程的根...
作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 一.实验原理 二.实验步骤 三.实验过程 1.(程序) (1)二分法:求 在区间(1,2)之间的根,取 ...
- 【MATLAB实验】数学实验实验求根问题三种方法、积分画图及古典密码设计(二分法、牛顿法、不动点迭代法)
目录 特征值与特征向量 求最简型 编程题 求积分问题(quad.trapz.int) 求根问题(二分法.不动点迭代法.牛顿法) 1.二分法: 2.不动点迭代法 3.牛顿切线法 画图题 画螺旋曲线图 画 ...
- 数值分析2-解非线性方程的不动点迭代法,斯蒂芬森(steffensen)迭代法,牛顿法以及简化牛顿法的matlab程序
对于形如f(x)=0的单变量非线性方程,可以构造不同的迭代函数进行迭代求根,以f(x)=x3-x-1为例,我们可以简单的通过等式变形构造出x=x3-1和x=(x+1)(1/3)两种等价形式:也可以自己 ...
- 计算方法实验:方程求根二分法、不动点迭代法、牛顿法
计算方法实验一,方程求根 分别用Matlab和C写了一下,初学Matlab,如有不足还请指正. 实验内容: 1)在区间[0,1]内用二分法求方程 e x + 10 ∗ x − 2 e^x+10*x-2 ...
- 二分法、试位法、不动点迭代法、牛顿法、割线法
二分法.试位法.不动点迭代法.牛顿法.割线法 问题回顾 问题分析 1.二分法 2.试位法 3.不动点迭代 4.Newton-Raphson法 5.割线法 小结 问题回顾 一段质量均匀分布的电缆线悬挂在 ...
- jacobi matlab程序,jacobi迭代法实验MATLAB程序数值分析
jacobi迭代法实验MATLAB程序数值分析 例1. 求线性方程组 得近似解.精确解为x*=[3,2,1]'. 解:对方程进行移项就得 记为Ax=b,或写为x=B0 x+f,其中 取初始值,代入原方 ...
- FPI(Fixed-point Iteration)不动点迭代法——迭代求方程的方法
一 不动点是什么? 不动点,其实定义比较简单,对于一些方程,例如f(x)=cosx,那么令cosx=x的点就是函数的不动点,说白了,就是y=x这条直线与函数曲线的交点.这个不动点有什么用呢?请继续往下 ...
最新文章
- 服务应该去版本化,不管是微服务还是SOA
- 模板模式(TemplateMethod)和策略(StrategyMethod)模式
- vxworks 学习和windows azure 学习
- javascript模式——Mixin
- Java 多线程(三) 线程的生命周期及优先级
- python分句_Python循环中的分句,继续和其他子句
- Python基础知识有哪些?你都知道吗
- php !$rs-gt;bof and !$rs-gt;eof,关于rs.bof rs.eof的一些解释,请大家帮我完善一下,谢谢...
- 合并多个js,css文件的方法 - 在服务端合并和压缩JavaScript和CSS文件
- java 解压文件夹_java新人,在学习解压缩文件时,解压文件后,里面的文件变文件夹...
- C# ComBox 垂直滚动条
- RTSP-传送ACC音频文件
- cad添加自己线性_如何在CAD中添加自行创建的线型.doc
- 99定时器设计c语言,51单片机写的曝光定时器C语言程序
- 字符集编码(一):Unicode 之前
- Linux环境变量PATH
- java拼图游戏(未补全)
- IDEA显示树状目录结构
- OnWorks-免费Windows/Linux/MAC等系统云服务器
- 推流式搅拌器选型功率计算方法_潜水搅拌机如何计算功率