a华南理工大学数值分析报告A
_____________.
10． 设{}k k 0q x ()∞
=为区间[0,1]上带权x ρ=且首项系数为1的k 次正
交多项式序列, 其中0q x 1()=, 则1q x ()=_________. 二.(10分) 用直接三角分解方法解下列线形方程组
123215x 114112x 27245x 12⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪--⎝⎭⎝⎭⎝⎭
三. (12分) 对于线性方程组
123142x 202310x 3521x 12-⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭
写出其Jacobi 迭代法及其Guass-Seidel 迭代法的分量形式, 并判断它们的收敛性.
四. (12分)
, 若将其视为22x 30()-=的根, (1). 写出相应的Newton 迭代公式.
(2). 指出其收敛阶(需说明依据).
五. (12分) 依据如下函数值表
(1). 构造插值多项式满足以上插值条件
(2). 推导出插值余项.
六.(10分) 已知离散数据表
若用形如2
y ax bx =+进行曲线拟合, 求出该拟合曲线.
七. (12分) 构造带权x
()ρ=
的Guass 型求积公式.
1
00110
1x dx A f x A f x ()()()≈+⎰
八. (12分) 对于常微分方程的初值问题
dy
2y dx y 02
()⎧=-⎪⎨⎪=⎩ (1). 若用改进的欧拉方法求解, 证明该方法的收敛性. (2). 讨论改进欧拉方法的稳定条件.。
数值分析综合实验报告
一、实验目的通过本次综合实验，掌握数值分析中常用的插值方法、方程求根方法以及数值积分方法，了解这些方法在实际问题中的应用，提高数值计算能力。
二、实验内容1. 插值方法（1）拉格朗日插值法：利用已知数据点构造多项式，以逼近未知函数。
（2）牛顿插值法：在拉格朗日插值法的基础上，通过增加基函数，提高逼近精度。
2. 方程求根方法（1）二分法：适用于函数在区间内有正负值的情况，通过不断缩小区间来逼近根。
（2）Newton法：利用函数的导数信息，通过迭代逼近根。
（3）不动点迭代法：将方程转化为不动点问题，通过迭代逼近根。
3. 数值积分方法（1）矩形法：将积分区间等分，近似计算函数值的和。
（2）梯形法：将积分区间分成若干等分，用梯形面积近似计算积分。
（3）辛普森法：在梯形法的基础上，将每个小区间再等分，提高逼近精度。
三、实验步骤1. 拉格朗日插值法（1）输入已知数据点，构造拉格朗日插值多项式。
（2）计算插值多项式在未知点的函数值。
2. 牛顿插值法（1）输入已知数据点，构造牛顿插值多项式。
（2）计算插值多项式在未知点的函数值。
3. 方程求根方法（1）输入方程和初始值。
（2）选择求解方法（二分法、Newton法、不动点迭代法）。
（3）迭代计算，直到满足精度要求。
4. 数值积分方法（1）输入被积函数和积分区间。
（2）选择积分方法（矩形法、梯形法、辛普森法）。
（3）计算积分值。
四、实验结果与分析1. 插值方法（1）拉格朗日插值法：通过构造多项式，可以较好地逼近已知数据点。
（2）牛顿插值法：在拉格朗日插值法的基础上，增加了基函数，提高了逼近精度。
2. 方程求根方法（1）二分法：适用于函数在区间内有正负值的情况，计算简单，但收敛速度较慢。
（2）Newton法：利用函数的导数信息，收敛速度较快，但可能存在数值不稳定问题。
（3）不动点迭代法：将方程转化为不动点问题，收敛速度较快，但可能存在初始值选择不当的问题。
3. 数值积分方法（1）矩形法：计算简单，但精度较低。
华工计算机计算方法(数值分析)考试试卷
考完试了，顺便把记得地题目背下来，应该都齐全了.我印象中也就只有这些题，题目中地数字应该是对地，我也验证过，不过也不一定保证是对地，也有可能我也算错了.还有就是试卷上面地题目可能没有我说地这么短，但是我也不能全把文字背下来，大概意思就是这样吧.每个部分地题目地顺序可能不是这样，但总体就是这四大块.至于每道题目地分值，我记得地就写出来了，有些题目没注意.我题目后面写地结果都是我考试时算出来地，考完了也懒得验证了，可能不一定对，自己把握吧，仅供参考.华南理工大学计算机计算方法（数值分析）考试试卷一填空题（分）1.（分）* ，准确值，求绝对误差(*) ，相对误差(*) ，有效数位是.（分）当插值函数地越大时，会出现龙格现象，为解决这个问题，分段函数不一个不错地办法，请写出分段线性插值、分段三次插值和三次样条插值各自地特点.3.（分）已知和相近，将–变换成可以使其计算结果更准确.4.（分）已知–，求牛顿迭代法地迭代式子.解题思路：. 这里地绝对误差和相对误差是没有加绝对值地，而且要注意是用哪个数减去哪个数得到地值，正负号会不一样；. 可以从它们函数地连续性方面来说明；. 只要满足课本所说地那几个要求就可以；这个记得迭代公式就可以直接写，记不住可以自己推导，就是用泰勒展开式来近似求值得到地迭代公式.我最终地结果是：1.2.分段线性插值保证了插值函数地连续性，但是插值函数地一次导数不一定连续；分段三次既保证了插值函数地连续性，也保证了其一次导数地连续性；三次样条插值保证了插值函数及其一次导数和二次导数地连续性3.()4.– ( –)( )二计算题（分）已知() –，用对分法求其在[ , ]区间内地根，误差要满小于，需要对分多少次？请写出最后地根结果.解题思路：每次求区间地中值并计算其对应地函数值，然后再计算下一个区间中值及函数值，一直到两次区间中值地绝对值小于为止.我最终算得地对分次数是，根地结果为.2.根据以下数据回答相应问题：(1)请根据以上数据构造三次插值函数；(2)请列出差商表并写出三次插值函数.解题思路：() 直接按照书本地定义把公式列出来就可以了，这个要把公式记住了才行，不然也写不了；()差商表就是计算三次插值函数过程中计算到地中间值及结果值，可以先在草稿上按照公式地计算过程把公式写出来，然后把中间用到地值整理成一个表格，这个表格就是差商表了，最后再把公式和表格都写到试卷上就行了.当然也可以先把表格写出来，再用表格地数据写出公式都可以.因为我考试地时候也是先写表格，但是我感觉算地时候容易错，特别是除数地位置，很容易搞错相减地两个地值.所以我想如果直接按照公式用到地值来算，可能没那么容易混乱，因为需要哪个就算哪个，地值比较明确，最后再把中间算出来地值填到表格里就可以了.当然这要看个人喜好了.这里地结果有点长，不好写在这里，自己搞定吧，不难，只是直接套公式就可以了.3. 请用分解法求解以下方程组地解⎪⎩⎪⎨⎧3- = x3 - 9x2 + 6x17 = 3x3+ x2 - 4x11- = x3 - x2 + 2x1解题思路：这个直接套公式算就好了，只要数没有算错，基本都是对地.有时候要注意看是列主元还是直接法，我当时好像没注意，这里应该没有要求用列主元.我最终算得地结果是, , ，其中算出来地矩阵分别是： ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-123121 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--12531124. （分）已知下列矩阵方程，根据以下要求回答问题： ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡210131012⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321x x x ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-111 (1) 求该矩阵方程地高斯赛达尔()迭代法地收敛性；(2) 求该矩阵方程地高斯赛达尔()迭代法地迭代公式；(3) 已知() ()，求()?解题思路：() 这个证明可以有两种方法，第一种用课本地定义来算，就是将系数矩阵地下三角系数全都乘上一个λ值，然后计算行列式，把所有地λ求出来，只要所有地λ都小于，那么就收敛；第二种方法就是用课本地定理证明，如果系数矩阵是强对角占优地，那么简单迭代法（）和迭代法都收敛，这道题刚好满足条件；() 这个迭代公式只要把矩阵和矩阵求出来就可以写出迭代公式了；() 把()代入()中地迭代公式就可以求出来.我地最终结果是：我直接用强对角占优证明，只写了两句话，不知道老师是不是要求我们用算地...至于强对角占优地判定，书上有，大概意思就是每一行中在主对角线上地那个数地绝对值比旁边所有数地绝对值加起来都要大就是强对角占优了.弱对角就是可以等于.详细定义翻书吧.(2) 我算出来地和矩阵如下：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--02/1003/10，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--03/1002/10，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-2/13/12/1迭代公式就是() () ()(3) () (, , )5. 已知以下方程，请利用最小二乘法求解：⎪⎪⎩⎪⎪⎨⎧0 = 7x2 + 2x1-13= 6x2 + 3x12 = 5x2 + x1-5 = 2x2+ x1解题思路：首先构造一个多变量拟合函数() ，可以把，看成是系数来求解，按照多变量拟合函数求解方法就可以得到结果.我最终算得地结果是：方程组为：⎪⎩⎪⎨⎧⨯=⨯+⨯⨯=⨯+⨯∑∑∑∑∑∑y t t t x t t x yt t t x t t x 22222111212111计算值并代入：⎩⎨⎧=+=+9821141422115x x x x计算地结果为：,请用复化梯形求积公式求出积分dx ⎰10x -e （注：里面地函数是）地近似值，要求误差限满足，请问需要将区间[]分成多少份？解题思路：首先是先把复化梯形求积公式地误差公式写出来，这个要记得，利用误差公式计算出满足精度要求地即可.我最终算得地结果是：误差公式为’’(ŋ)ŋŋ≤≤，≥√≈，也就是满足条件.三证明题（分）已知函数()，其在区间[]内地三个插值点为,（）. 请证明函数()在[]区间内满足下列关系： 6/)]()2/)((4)()[()(b f b a f a f a b dx x f b a +++-≈⎰解题思路：利用这三个插值点写出插值函数，原函数约等于插值函数，所以原函数地积分也约等于插值函数地积分，然后算出插值函数地积分结果就是证明地公式，其实这个就是课本地公式地证明.这个证明过程看课本吧.四程序题（分）前面有一段介绍列主元高斯消元法地步骤地说明（没背下来，都是文字，参考课本吧） 请按照列主元高斯消元法地思路将代码中地空格填写完整：1. 输入系数矩阵，右端项及ε；2. 选主元及消元：选主元： ≤≤若 <ε，则打印“求解失败”，停机；否则若≠，则交换地第行和行，交换行和行；消元：––3. 回代若≤ε，则打印“求解失败”，停机，否则(∑+=nijaijxj1)4.打印(…)解题思路：这个直接按照列主元高斯消去法地计算过程去写就好了.结果我写在代码里面了，是按照课本写地，我考试地时候写地应该也是这样.。
数值分析实验报告4
function f=gg(x) f=sqrt(2.5-(x^3)/4); end
则三种方法运行出的结果分别为： >> BDD(1.5) k= 20 ans = 1.365229578333959 >> WA(1.5) k= 6 ans = 1.365230351032824 >> ATJ(1.5) k= 4 ans = 1.365230013413594 （5）在 MATLAB 的 Editor 中建立一个 M-文件，输入程序代码，实现 matlab 自带函数求根 的程序代码如下：
>>format long BDD(1.5) k= 7 ans = 1.365230575673434
（2）在 MATLAB 的 Editor 中建立一个 M-文件，输入程序代码，实现加权加速求根的程序 代码如下：
function [y,n]=WA(a,eps) if nargin==1 eps=1e-16; end syms x L=subs(diff(sqrt(10/(4+x))),a); xl=gg(a,L); n=1; E=abs(xl-a); while (E>=eps)&(n<=10000) x=xl; xl=gg(x,L); n=n+1; E=abs(xl-x); end
在 command Windows 中输入命令：XJF(1.5)，得出的结果为：
>> XJF(1.5) k= 5 ans = 1.365230020178121 （8）在 MATLAB 的 Editor 中建立一个 M-文件，输入程序代码，实现抛物线法求根的程序
代码如下：
function [y,n] = parabola(fx,a,b,c,eps) if(nargin == 4) eps=1e-6; syms x fa = subs(fx,a); fb = subs(fx,b); fc = subs(fx,c); fb_a = (fb-fa)/(b-a); fc_b = (fc-fb)/(c-b); fc_b_a = (fc_b - fb_a)/(c-a); w = fc_b + fc_b_a*(c-b); x0 = c - 2*fc/(w+(w^2 - 4*fc*fc_b_a)^.5); n=1; if(abs(x0-c) > eps && n <= 10000) a=b; b=c; c=x0; fa = subs(fx,a); fb = subs(fx,b); fc = subs(fx,c); fb_a = (fb-fa)/(b-a); fc_b = (fc-fb)/(c-b); fc_b_a = (fc_b - fb_a)/(c-a); w = fc_b + fc_b_a*(c-b); x0=c - 2*fc/(w+(w^2 - 4*fc*fc_b_a)^.5); n=n+1; end y=x0; n; end end
华南理工大学数值分析
二.(10分)用直接三角分解方法解下列线形方程组
三. (12分)对于线性方程组
写出其Jacobi迭代法及其Guass-Seidel迭代法的分量形式,并判断它们的收敛性.
四. (12分)对于求 的近似值,若将其视为 的根,
(1).写出相应的Newton迭代公式.
(2).指出其收敛阶(需说明依据).
五.(12分)依据如下函数值表
0
1
1
2
0
(1).构造插值多项式满足以上插值条件
(2).推导出插值余项.
六.(10分)已知离散数据表
x
1
2
3
4
y=f(x)
0.8
1.5
1.8
2.0
若用形如 进行曲线拟合,求出该拟合曲线.
七. (12分)构造带权 的Guass型求积公式.
诚信应考,考试作弊将带来严重后果！华南理工ຫໍສະໝຸດ 学期末考试《数值分析》试卷A卷
注意事项：1.考前请将密封线内各项信息填写清楚；
2.可使用计算器，解答就答在试卷上；
3．考试形式：闭卷；
4.本试卷共八大题，满分100分。考试时间120分钟。
题号
一
二
三
四
五
六
七
八
总分
得分
评卷人
一．填空题（每小题2分,共20分）
八. (12分)对于常微分方程的初值问题
(1).若用改进的欧拉方法求解,证明该方法的收敛性.
(2).讨论改进欧拉方法的稳定条件.
1．已知自然数e＝2.718281828459045…，取e≈2.71828，那么e具有的有效数字是____________.
数值分析实验报告心得(3篇)
第1篇在数值分析这门课程的学习过程中，我深刻体会到了理论知识与实践操作相结合的重要性。
通过一系列的实验，我对数值分析的基本概念、方法和应用有了更加深入的理解。
以下是我对数值分析实验的心得体会。
一、实验目的与意义1. 巩固数值分析理论知识：通过实验，将课堂上学到的理论知识应用到实际问题中，加深对数值分析概念和方法的理解。
2. 培养实际操作能力：实验过程中，我学会了使用Matlab等软件进行数值计算，提高了编程能力。
3. 增强解决实际问题的能力：实验项目涉及多个领域，通过解决实际问题，提高了我的问题分析和解决能力。
4. 培养团队协作精神：实验过程中，我与同学们分工合作，共同完成任务，培养了团队协作精神。
二、实验内容及方法1. 实验一：拉格朗日插值法与牛顿插值法（1）实验目的：掌握拉格朗日插值法和牛顿插值法的原理，能够运用这两种方法进行函数逼近。
（2）实验方法：首先，我们选择一组数据点，然后利用拉格朗日插值法和牛顿插值法构造插值多项式。
最后，我们将插值多项式与原始函数进行比较，分析误差。
2. 实验二：方程求根（1）实验目的：掌握二分法、Newton法、不动点迭代法、弦截法等方程求根方法，能够运用这些方法求解非线性方程的根。
（2）实验方法：首先，我们选择一个非线性方程，然后运用二分法、Newton法、不动点迭代法、弦截法等方法求解方程的根。
最后，比较不同方法的收敛速度和精度。
3. 实验三：线性方程组求解（1）实验目的：掌握高斯消元法、矩阵分解法等线性方程组求解方法，能够运用这些方法求解线性方程组。
（2）实验方法：首先，我们构造一个线性方程组，然后运用高斯消元法、矩阵分解法等方法求解方程组。
最后，比较不同方法的计算量和精度。
4. 实验四：多元统计分析（1）实验目的：掌握多元统计分析的基本方法，能够运用这些方法对数据进行分析。
（2）实验方法：首先，我们收集一组多元数据，然后运用主成分分析、因子分析等方法对数据进行降维。
数值分析实验报告——非线性方程求根
数值分析实验报告——非线性方程求根一、实验目的：1.掌握求解非线性方程的常用方法；2.了解非线性方程求根问题的数值解法；3.熟悉使用数值分析软件进行非线性方程求根的实现。
二、实验原理：非线性方程指的是形如f(x)=0的方程，其中f(x)是一个非线性函数。
非线性方程求根的常用方法包括二分法、割线法和牛顿法等。
其中，二分法是通过不断缩小区间范围来逼近方程的解；割线法是通过使用割线来逼近方程的解；牛顿法则是通过使用切线来逼近方程的解。
对于给定的非线性方程，可以根据实际情况选择合适的方法进行求根。
三、实验内容：1.编写求解非线性方程的函数，包括二分法、割线法和牛顿法；2.使用编写的函数求解给定的非线性方程，比较各个方法的收敛速度和精确程度；3.根据实际情况分析和选择合适的方法进行求根。
四、实验步骤：1.针对给定的非线性方程，编写二分法的函数实现：（1）首先确定方程的解存在的区间；（2）根据方程的解存在的区间，使用二分法逐步缩小区间范围；（3）根据设定的精度要求，不断循环迭代，直至满足要求或达到迭代次数限制；2.针对给定的非线性方程，编写割线法的函数实现：（1）首先需要确定方程的解存在的初始点；（2）根据方程的解存在的初始点，根据割线的定义进行迭代；（3）设定迭代的精度要求和限制次数，结束迭代；3.针对给定的非线性方程，编写牛顿法的函数实现：（1）首先需要确定方程的解存在的初始点；（2）根据方程的解存在的初始点，根据牛顿法的定义进行迭代；（3）设定迭代的精度要求和限制次数，结束迭代；4.根据给定的非线性方程，分别使用二分法、割线法和牛顿法进行求解，并比较各个方法的收敛速度和精确程度；5.分析实际情况，选择合适的方法进行求解。
五、实验结果：4.通过比较，发现割线法和牛顿法的收敛速度较快，精确程度较高，因此选择割线法进行求解。
六、实验总结：通过本次实验，我掌握了求解非线性方程的常用方法，并使用数值分析软件实现了二分法、割线法和牛顿法。
数值分析
3. 数值计算中的误差
• 来源及种类 --- 模型误差、参数误差、 截断误差、舍入误差。
1. 模型误差（也称描述误差）
模型误差是在建立数学模型时，由于忽略了一些次要因素
而产生的误差，它是数学建模阶段要考虑的误差，不是计算 方法可以解决的。
2. 参数误差（也称观测误差）
测量已知参数时，数据带来的误差 ，它也不是计算方法 能解决的问题。
插值：已知[a, b]上的函数y= f (x)在n+1个互异点处的函数值:
x
x0
x1
x2
xn
f(x)
f0
f1
f2
fn
求简单函数 P (x)，使得
P(xi ) fi
i 0,1,L ,n L L (*)
计算 f (x)可通过计算 P (x)来近似代替。如下图所示。
y
f(x)
▪ Matlab 简介
第一章 绪 论
主要内容： 一些常用概念； 数值算法的复杂度与稳定性。 数值计算中的误差； 数值算法设计的若干原则；
1.数值分析中常用的一些概念
数值问题 数值解 算法 计算量 病态问题 算法数值稳定性
数值问题、数值解 、算法
➢ 由一组已知数据（输入数据），求出一组结果 数据（输出数据），使得这两组数据之间满足 预先制定的某种关系的问题，称为数值问题。
舍入误差是由于计算机只能表示有限位数字，因而只能取 有限位数进行计算所得的误差，它也是计算方法关注的内容。
（举例： 3.141592653L ）
数值计算中的误差
• 误差的基本概念
➢ 绝对误差 --- 近似数 x * 关于准确数 x 的绝对误差：
数值分析实验报告5篇
1.69376699767424 0.92310666706964 0.08471614569741 0.40804026409411
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
讨论：
利用这种方法进行这类实验，可以很精确的扰动敏感性的一般规律。即 当对扰动项的系数越来越小时，对其多项式扰动的结果也就越来越小， 即扰动敏感性与扰动项的系数成正比，扰动项的系数越大，对其根的扰 动敏感性就越明显，当扰动的系数一定时，扰动敏感性与扰动的项的幂 数成正比，扰动的项的幂数越高，对其根的扰动敏感性就越明显。
解线性方程组的直接方法
实验 （主元的选取与算法的稳定性） 问题提出：Gauss消去法是我们在线性代数中已经熟悉的。但由于计算 机的数值运算是在一个有限的浮点数集合上进行的，如何才能确保 Gauss消去法作为数值算法的稳定性呢？Gauss消去法从理论算法到数值 算法，其关键是主元的选择。主元的选择从数学理论上看起来平凡，它 却是数值分析中十分典型的问题。 实验内容：考虑线性方程组 编制一个能自动选取主元，又能手动选取主元的求解线性方程组的 Gauss消去过程。 实验要求： （1）取矩阵，则方程有解。取n=10计算矩阵的条件数。让程序自动选 取主元，结果如何？ （2）现选择程序中手动选取主元的功能。每步消去过程总选取按模最 小或按模尽可能小的元素作为主元，观察并记录计算结果。若每步消去 过程总选取按模最大的元素作为主元，结果又如何？分析实验的结果。 （3）取矩阵阶数n=20或者更大，重复上述实验过程，观察记录并分析 不同的问题及消去过程中选择不同的主元时计算结果的差异，说明主元
数值分析实验报告总结
一、实验背景数值分析是研究数值计算方法及其理论的学科，是计算机科学、数学、物理学等领域的重要基础。
为了提高自身对数值分析理论和方法的理解，我们进行了数值分析实验，通过实验加深对理论知识的掌握，提高实际操作能力。
二、实验目的1. 理解数值分析的基本理论和方法；2. 掌握数值分析实验的基本步骤和技巧；3. 培养实验设计和数据分析能力；4. 提高编程和计算能力。
三、实验内容本次实验主要分为以下几个部分：1. 线性方程组求解实验：通过高斯消元法、LU分解法等求解线性方程组，并分析算法的稳定性和误差；2. 矩阵特征值问题计算实验：利用幂法、逆幂法等计算矩阵的特征值和特征向量，分析算法的收敛性和精度；3. 非线性方程求根实验：运用二分法、牛顿法、不动点迭代法等求解非线性方程的根，比较不同算法的优缺点；4. 函数插值实验：运用拉格朗日插值、牛顿插值等方法对给定的函数进行插值，分析插值误差；5. 常微分方程初值问题数值解法实验：运用欧拉法、改进的欧拉法、龙格-库塔法等求解常微分方程初值问题，比较不同算法的稳定性和精度。
四、实验过程1. 线性方程组求解实验：首先，编写程序实现高斯消元法、LU分解法等算法；然后，对给定的线性方程组进行求解，记录计算结果；最后，分析算法的稳定性和误差。
2. 矩阵特征值问题计算实验：编写程序实现幂法、逆幂法等算法；然后，对给定的矩阵进行特征值和特征向量的计算，记录计算结果；最后，分析算法的收敛性和精度。
3. 非线性方程求根实验：编写程序实现二分法、牛顿法、不动点迭代法等算法；然后，对给定的非线性方程进行求根，记录计算结果；最后，比较不同算法的优缺点。
4. 函数插值实验：编写程序实现拉格朗日插值、牛顿插值等方法；然后，对给定的函数进行插值，记录计算结果；最后，分析插值误差。
5. 常微分方程初值问题数值解法实验：编写程序实现欧拉法、改进的欧拉法、龙格-库塔法等算法；然后，对给定的常微分方程初值问题进行求解，记录计算结果；最后，比较不同算法的稳定性和精度。
数值分析实验报告模板
数值分析实验报告模板篇一：数值分析实验报告(一)(完整)数值分析实验报告12345篇二：数值分析实验报告实验报告一题目：非线性方程求解摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。
本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。
利用二分法求解给定非线性方程的根，在给定的范围内，假设f(x,y)在[a,b]上连续，f(a)xf(b) 直接影响迭代的次数甚至迭代的收敛与发散。
即若x0 偏离所求根较远，Newton法可能发散的结论。
并且本实验中还利用利用改进的Newton法求解同样的方程，且将结果与Newton法的结果比较分析。
前言：（目的和意义）掌握二分法与Newton法的基本原理和应用。
掌握二分法的原理，验证二分法，在选对有根区间的前提下，必是收敛，但精度不够。
熟悉Matlab语言编程，学习编程要点。
体会Newton使用时的优点，和局部收敛性，而在初值选取不当时，会发散。
数学原理：对于一个非线性方程的数值解法很多。
在此介绍两种最常见的方法：二分法和Newton法。
对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b) Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式xk?1?xk?f(xk) f'(xk)产生逼近解x*的迭代数列{xk}，这就是Newton法的思想。
当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。
另外，若将该迭代公式改进为xk?1?xk?rf(xk) 'f(xk)其中r为要求的方程的根的重数，这就是改进的Newton 法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。
程序设计：本实验采用Matlab的M文件编写。
其中待求解的方程写成function的方式，如下function y=f(x);y=-x*x-sin(x);写成如上形式即可，下面给出主程序。