华科《科学计算引论》编程作业

发布 2023-05-16 23:46:28 阅读 7000

科学计算引论第五章编程作业。

一。 本章总结。

第五章,非线性方程的数值解法主要介绍用几何方法(包括二分法、弦截法、方法)、迭代法、迭代法来解得非线性方程的数值解。并比较了各种数值解法的收敛速度,提出了收敛阶的概念。

二.实际问题的提出。

如图,小明一天将一个密度,半径的木球放置在水中,已知水的密度,重力加速度,试求系统稳定后球体浸入水中的深度。

三.问题分析。

由阿基米德定理:,而即半径为,高度为的球冠体积,所以。 再由而力平衡条件:。

综上所述得方程:,化简后代入数值得:

故只须求解上述非线性方程即可。

四.问题求解。

令,,,单调递减,而,所以在中一定有唯一根。先用函数图像找初值。**为:

fplot('12500*x^3-3750*x^2+3',[0,0.2]);grid on; %画出函数f(x)的图像。

图像经放大后为:

由图易知,零点在附近,故取为迭代初值,用迭代法进行求解,迭代格式为。要求迭代至满足。

**为:e=1e-8;format long; %定义e作为判断迭代停止的误差。

x1=0.03将x1定义为0.03

x0=x1+1让while语句成立。

while(abs(x0-x1)>e)

x0=x1,x1=x0-(12500*x0^3-3750*x0^2+3)/(37500*x0^2-7500*x0) %迭代过程。

end;format

运行程序得。x1 =

x0 =x1 =

x0 =x1 =

x0 =x1 =

综上,经3次迭代,得到满足误差要求的数值解。

五.结论。本题经牛顿迭代法得数值解。

计算机引论大作业

学生信息管理系统。班级 电子1202 学号 20123685 姓名 张欣。学生信息管理系统。一 概述。该程序是一个关于学生信息管理的程序,它包括了新建学生信息 增加学生信息 修改学生信息 删除学生信息 查找学生信息 学生信息排序 查看学生信息等七大功能。如下图所示 每个学生的信息又包括学号 姓名 性...

华师计算机作业

1 第1题。在windows中,打开一个窗口后,通常在其顶部是一个 a.标题栏。b.任务栏。c.状态栏。d.工具栏。您的答案 c 题目分数 1.0 此题得分 0.0 2 第2题。计算机中的所有信息都是以二进制方式表示的,主要理由是 a.运算速度快。b.节约元件。c.所需的物理元件最简单。d.信息处理...

华科电气MATLAB大作业

华中科技大学。电气与电子工程学院。matlab课程作业 班级。学号。姓名。时间 2014年12月25日。目录。一 概述 2 二 设计要求 2 三 设计分析 2 1.系统的稳态误差理论分析 3 2.系统稳态误差 分析 3 3.阶跃响应 分析 4 四 根轨迹法设计相位滞后环节 9 1.相位滞后环节设计 ...