云题海 - 专业文章范例文档资料分享平台

当前位置:首页 > 计算方法上机作业

计算方法上机作业

  • 62 次阅读
  • 3 次下载
  • 2025/7/10 0:24:48

方程求根

一、目的和意义

非线性方程在科学研究与工程实践中广泛出现,例如,优化问题、特征值问题、微分方程问题等。但是,除少量方程外,大多数非线性方程求根相当困难,常见的几个简单、有效的数值求根方法,包括二分法、迭代法、牛顿法、割线法等,本实验旨在比较各种算法的计算性能和使用范围。

二、计算公式 1.二分法 2.不动点迭代法 三、结构程序设计代码 1.二分法

1).定义所求解函数

function [ y ] = f( x ) y = x^3 + 4*x^2 - 10; end

2).执行算法

%初始化,设置区间端点a、b,误差限tol a = 1; b = 2; tol = 0.5*10^(-6); k = 0; fa = f(a); %设置最大二分次数为30 for k = 1:50

p = (a + b)/2; fp = f(p); if(fp == 0 || (b - a)/2 < tol) break end

if(fa * fp < 0) b = p; else a = p;

end

disp('近似解p = ');disp(vpa(p,10)); disp('迭代次数k = ');disp(k); end

2.不动点迭代法

1).定义不动点方程g(x)

function [ y ] = g( x )

y = x^3 + 4*x^2 + x - 10; end

2).执行算法

%初始化,设置误差限,设置初值p0 tol = 0.5*10^(-6); k = 1; p0 = 1.5; %迭代次数为10次 while k <= 10 p = g(p0);

if abs(p - p0) < tol break end

disp('近似解p = ');disp(vpa(p,10)); disp('迭代次数k = ');disp(k); k = k + 1; p0 = p; end

四、结果及其讨论 1.二分法结果

由于结果较长,只取了一部分,从图中可以看出,迭代20次可得到误差限范围内的近似解p=1.36522。

2.不动点迭代法结果

由结果可以看出,函数不收敛,当迭代次数为6时,已经近似趋向于无穷大。

3.讨论

上面实验通过两种方法近似计算了函数的解,二分法可在迭代20次的情况下求出近似解,而由于所设的不动点方程发散,所以使用不动点迭代法无法求出近似解。因此,在实际计算过程中,使用不动点迭代法时,在建立不动点方程后,要分析其收敛性。

搜索更多关于: 计算方法上机作业 的文档
  • 收藏
  • 违规举报
  • 版权认领
下载文档10.00 元 加入VIP免费下载
推荐下载
本文作者:...

共分享92篇相关文档

文档简介:

方程求根 一、目的和意义 非线性方程在科学研究与工程实践中广泛出现,例如,优化问题、特征值问题、微分方程问题等。但是,除少量方程外,大多数非线性方程求根相当困难,常见的几个简单、有效的数值求根方法,包括二分法、迭代法、牛顿法、割线法等,本实验旨在比较各种算法的计算性能和使用范围。 二、计算公式 1.二分法 2.不动点迭代法 三、结构程序设计代码 1.二分法 1).定义所求解函数 function [ y ] = f( x ) y = x^3 + 4*x^2 - 10; end 2).执行算法 %初始化,设置区间端点a、b,误差限tol a = 1; b = 2; tol = 0.5*10^(-6); k = 0; fa = f(a); %设置最大二分次数为30 for k = 1

× 游客快捷下载通道(下载后可以自由复制和排版)
单篇付费下载
限时特价:10 元/份 原价:20元
VIP包月下载
特价:29 元/月 原价:99元
低至 0.3 元/份 每月下载150
全站内容免费自由复制
VIP包月下载
特价:29 元/月 原价:99元
低至 0.3 元/份 每月下载150
全站内容免费自由复制
注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信:fanwen365 QQ:370150219
Copyright © 云题海 All Rights Reserved. 苏ICP备16052595号-3 网站地图 客服QQ:370150219 邮箱:370150219@qq.com