数值分析实验
数值分析综合实验报告
一、实验目的通过本次综合实验，掌握数值分析中常用的插值方法、方程求根方法以及数值积分方法，了解这些方法在实际问题中的应用，提高数值计算能力。
二、实验内容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）矩形法：计算简单，但精度较低。
土木数值分析实验报告(3篇)
第1篇一、实验背景随着我国土木工程领域的不断发展，数值分析方法在工程设计和施工中的应用越来越广泛。
为了更好地理解和掌握数值分析方法，本实验报告以某典型土木工程问题为背景，通过数值分析软件对问题进行模拟，分析结果并得出结论。
二、实验目的1. 熟悉数值分析软件的基本操作和功能。
2. 建立合理的数值模型，对土木工程问题进行模拟分析。
3. 分析模拟结果，验证理论计算的准确性，为实际工程提供参考。
三、实验内容1. 问题背景：某桥梁工程中，需要进行桥梁结构的稳定性分析。
2. 数值模型建立：- 选择合适的数值分析软件（如ANSYS、ABAQUS等）。
- 建立桥梁结构的几何模型，包括桥梁的梁、板、柱等构件。
- 确定材料属性，如弹性模量、泊松比等。
- 设置边界条件和加载方式。
3. 数值模拟：- 进行网格划分，确保网格质量满足分析要求。
- 运行模拟，获取桥梁结构的应力、应变等数据。
4. 结果分析：- 分析桥梁结构的应力分布情况，确定结构的安全性。
- 分析桥梁结构的变形情况，评估结构的舒适性。
- 将模拟结果与理论计算结果进行对比，验证数值方法的准确性。
四、实验步骤1. 模型建立：- 使用CAD软件绘制桥梁结构的几何模型。
- 将几何模型导入数值分析软件。
2. 材料属性设置：- 根据设计规范和实际材料性能，设置材料的弹性模量、泊松比等参数。
3. 边界条件和加载方式设置：- 根据实际工程情况，设置边界条件和加载方式。
4. 网格划分：- 选择合适的网格划分方法，确保网格质量满足分析要求。
5. 模拟运行：- 运行模拟，获取桥梁结构的应力、应变等数据。
6. 结果分析：- 分析桥梁结构的应力分布情况，确定结构的安全性。
- 分析桥梁结构的变形情况，评估结构的舒适性。
- 将模拟结果与理论计算结果进行对比，验证数值方法的准确性。
五、实验结果与分析1. 应力分布情况：- 模拟结果显示，桥梁结构的最大应力出现在梁的支座处，符合理论计算结果。
数值分析积分实验报告(3篇)
第1篇一、实验目的本次实验旨在通过数值分析的方法，研究几种常见的数值积分方法，包括梯形法、辛普森法、复化梯形法和龙贝格法，并比较它们在计算精度和效率上的差异。
通过实验，加深对数值积分理论和方法的理解，提高编程能力和实际问题解决能力。
二、实验内容1. 梯形法梯形法是一种基本的数值积分方法，通过将积分区间分割成若干个梯形，计算梯形面积之和来近似积分值。
实验中，我们选取了几个不同的函数，对积分区间进行划分，计算积分近似值，并与实际积分值进行比较。
2. 辛普森法辛普森法是另一种常见的数值积分方法，它通过将积分区间分割成若干个等距的区间，在每个区间上使用二次多项式进行插值，然后计算多项式与x轴围成的面积之和来近似积分值。
实验中，我们对比了辛普森法和梯形法的计算结果，分析了它们的精度差异。
3. 复化梯形法复化梯形法是对梯形法的一种改进，通过将积分区间分割成多个小区间，在每个小区间上使用梯形法进行积分，然后计算所有小区间积分值的和来近似积分值。
实验中，我们对比了复化梯形法和辛普森法的计算结果，分析了它们的精度和效率。
4. 龙贝格法龙贝格法是一种通过外推加速提高计算精度的数值积分方法。
它通过比较使用不同点数（n和2n）的积分结果，得到更高精度的积分结果。
实验中，我们使用龙贝格法对几个函数进行积分，并与其他方法进行了比较。
三、实验步骤1. 编写程序实现梯形法、辛普森法、复化梯形法和龙贝格法。
2. 选取几个不同的函数，对积分区间进行划分。
3. 使用不同方法计算积分近似值，并与实际积分值进行比较。
4. 分析不同方法的精度和效率。
四、实验结果与分析1. 梯形法梯形法在计算精度上相对较低，但当积分区间划分足够细时，其计算结果可以接近实际积分值。
2. 辛普森法辛普森法在计算精度上优于梯形法，但当积分区间划分较细时，计算量较大。
3. 复化梯形法复化梯形法在计算精度上与辛普森法相当，但计算量较小。
4. 龙贝格法龙贝格法在计算精度上优于复化梯形法，且计算量相对较小。
数值分析方法实验报告
一、实验目的通过本次实验，掌握数值分析方法的基本原理和应用，熟悉MATLAB编程环境，学会使用MATLAB进行数值计算，并分析不同数值方法的优缺点。
二、实验内容1. 二分法求方程根（1）原理：二分法是一种迭代方法，通过不断缩小根所在的区间，直到满足精度要求为止。
（2）步骤：①给定初始区间[a, b]，使得f(a) f(b) < 0；②计算区间中点c = (a + b) / 2；③判断f(c)的符号：a. 若f(c) = 0，则c为方程的根；b. 若f(c) f(a) < 0，则新的区间为[a, c]；c. 若f(c) f(b) < 0，则新的区间为[c, b]；④重复步骤②和③，直到满足精度要求。
（3）代码实现：```MATLABfunction root = bisection_method(f, a, b, tol)while (b - a) / 2 > tolc = (a + b) / 2;if f(c) == 0break;elseif f(a) f(c) < 0b = c;elsea = c;endendroot = (a + b) / 2;end```2. Newton法求方程根（1）原理：Newton法是一种基于切线逼近的迭代方法，通过不断逼近函数的零点。
（2）步骤：①给定初始值x0；②计算导数f'(x)和f(x)在x0处的值；③计算新的近似值x1 = x0 - f(x0) / f'(x0)；④重复步骤②和③，直到满足精度要求。
（3）代码实现：```MATLABfunction root = newton_method(f, df, x0, tol)while abs(f(x0)) > tolx1 = x0 - f(x0) / df(x0);x0 = x1;endroot = x0;end```3.不动点迭代法求方程根（1）原理：不动点迭代法是一种迭代方法，通过不断逼近不动点，即方程的根。
数值分析追赶法实验报告(3篇)
第1篇一、实验目的通过本次实验，掌握追赶法的基本原理和计算步骤，了解追赶法在解三对角线性方程组中的应用，并学会利用C++编程实现追赶法，提高编程能力。
二、实验原理追赶法是一种解三对角线性方程组的迭代方法，其基本原理是利用递推公式逐步求解未知数。
对于形如Ax=b的三对角线性方程组，其中系数矩阵A具有如下形式：A = [a_00, a_01, 0, ..., 0;a_10, a_11, a_12, ..., 0;0, a_21, a_22, ..., a_2n-1;...;0, ..., 0, a_n2]追赶法将系数矩阵A分解为两个因子L和U，其中L为下三角矩阵，U为上三角矩阵，即：A = LU其中：L = [1, 0, ..., 0;a_10/a_11, 1, 0, ..., 0;..., ..., ..., ...;a_n1/a_n2, ..., ..., ..., 1]U = [a_11, a_12, ..., a_1n;0, a_22, ..., a_2n-1;..., ..., ..., ...;0, ..., ..., a_n2]通过递推公式求解L和U中的元素，进而得到解向量x：x_1 = b_1 / a_11x_i = (b_i - ∑(j=1 to i-1) l_ij x_j) / u_ij, i = 2, ..., n三、实验步骤1. 编写C++程序，实现追赶法的基本算法。
2. 生成三对角线性方程组的系数矩阵A和解向量b。
3. 调用C++程序，计算追赶法的结果，并输出解向量x。
4. 分析追赶法的计算过程，验证结果是否正确。
四、实验数据及结果1. 生成三对角线性方程组的系数矩阵A和解向量b。
假设A为：A = [2, -1, 0, 0, 0;-1, 2, -1, 0, 0;0, -1, 2, -1, 0;0, 0, -1, 2, -1;0, 0, 0, -1, 2]b = [1, 1, 1, 1, 1]2. 追赶法计算结果。
数值分析实验报告心得(3篇)
第1篇在数值分析这门课程的学习过程中，我深刻体会到了理论知识与实践操作相结合的重要性。
通过一系列的实验，我对数值分析的基本概念、方法和应用有了更加深入的理解。
以下是我对数值分析实验的心得体会。
一、实验目的与意义1. 巩固数值分析理论知识：通过实验，将课堂上学到的理论知识应用到实际问题中，加深对数值分析概念和方法的理解。
2. 培养实际操作能力：实验过程中，我学会了使用Matlab等软件进行数值计算，提高了编程能力。
3. 增强解决实际问题的能力：实验项目涉及多个领域，通过解决实际问题，提高了我的问题分析和解决能力。
4. 培养团队协作精神：实验过程中，我与同学们分工合作，共同完成任务，培养了团队协作精神。
二、实验内容及方法1. 实验一：拉格朗日插值法与牛顿插值法（1）实验目的：掌握拉格朗日插值法和牛顿插值法的原理，能够运用这两种方法进行函数逼近。
（2）实验方法：首先，我们选择一组数据点，然后利用拉格朗日插值法和牛顿插值法构造插值多项式。
最后，我们将插值多项式与原始函数进行比较，分析误差。
2. 实验二：方程求根（1）实验目的：掌握二分法、Newton法、不动点迭代法、弦截法等方程求根方法，能够运用这些方法求解非线性方程的根。
（2）实验方法：首先，我们选择一个非线性方程，然后运用二分法、Newton法、不动点迭代法、弦截法等方法求解方程的根。
最后，比较不同方法的收敛速度和精度。
3. 实验三：线性方程组求解（1）实验目的：掌握高斯消元法、矩阵分解法等线性方程组求解方法，能够运用这些方法求解线性方程组。
（2）实验方法：首先，我们构造一个线性方程组，然后运用高斯消元法、矩阵分解法等方法求解方程组。
最后，比较不同方法的计算量和精度。
4. 实验四：多元统计分析（1）实验目的：掌握多元统计分析的基本方法，能够运用这些方法对数据进行分析。
（2）实验方法：首先，我们收集一组多元数据，然后运用主成分分析、因子分析等方法对数据进行降维。
数值分析实验报告
一、实验目的1. 理解数值分析的基本概念和常用算法；2. 掌握数值方法在求解实际问题中的应用；3. 培养编程能力，提高对数值分析软件的使用熟练度。
二、实验内容本次实验主要涉及以下内容：1. 拉格朗日插值法；2. 牛顿插值法；3. 线性方程组的求解方法；4. 方程求根的数值方法；5. 最小二乘法曲线拟合。
三、实验步骤1. 拉格朗日插值法（1）输入数据：给定一组数据点（x1, y1）、（x2, y2）、...、（xn, yn）。
（2）计算拉格朗日插值多项式L(x)。
（3）利用L(x)计算待求点x0的函数值y0。
2. 牛顿插值法（1）输入数据：给定一组数据点（x1, y1）、（x2, y2）、...、（xn, yn）。
（2）计算牛顿插值多项式N(x)。
（3）利用N(x)计算待求点x0的函数值y0。
3. 线性方程组的求解方法（1）输入数据：给定线性方程组的系数矩阵A和常数向量b。
（2）采用高斯消元法求解线性方程组Ax=b。
4. 方程求根的数值方法（1）输入数据：给定函数f(x)和初始值x0。
（2）采用二分法求解方程f(x)=0的根。
5. 最小二乘法曲线拟合（1）输入数据：给定一组数据点（x1, y1）、（x2, y2）、...、（xn, yn）。
（2）建立线性最小二乘模型y=F(x)。
（3）利用最小二乘法求解模型参数。
四、实验结果与分析1. 拉格朗日插值法与牛顿插值法的比较通过实验，我们发现牛顿插值法的精度高于拉格朗日插值法。
这是因为牛顿插值法在计算过程中考虑了前一项的导数信息，从而提高了插值多项式的平滑性。
2. 线性方程组的求解方法高斯消元法在求解线性方程组时，计算过程较为繁琐，但稳定性较好。
在实际应用中，可根据具体问题选择合适的方法。
3. 方程求根的数值方法二分法在求解方程时，收敛速度较慢，但具有较好的稳定性。
对于初始值的选择，应尽量接近真实根。
4. 最小二乘法曲线拟合最小二乘法在拟合曲线时，误差较小，适用于数据点较多的情况。
数值分析原理实验报告
一、实验目的通过本次实验，掌握数值分析的基本原理和方法，了解数值分析在科学和工程领域的应用，培养动手能力和分析问题的能力。
二、实验内容1. 二分法求方程根（1）原理：二分法是一种在实数域上寻找函数零点的算法。
对于连续函数f(x)，如果在区间[a, b]上f(a)f(b)<0，则存在一个根在区间(a, b)内。
二分法的基本思想是将区间[a, b]不断二分，缩小根所在的区间，直到满足精度要求。
（2）实验步骤：① 输入函数f(x)和精度要求；② 初始化区间[a, b]和中间点c=a+(b-a)/2；③ 判断f(c)与f(a)的符号，若符号相同，则将区间缩小为[a, c]，否则缩小为[c,b]；④ 重复步骤②和③，直到满足精度要求；⑤ 输出根的近似值。
2. 牛顿法求方程根（1）原理：牛顿法是一种在实数域上寻找函数零点的算法。
对于可导函数f(x)，如果在点x0附近，f(x0)f'(x0)≠0，则存在一个根在点x0附近。
牛顿法的基本思想是通过泰勒展开近似函数，然后求解近似方程的根。
（2）实验步骤：① 输入函数f(x)和精度要求；② 初始化迭代次数n=0，近似根x0；③ 计算导数f'(x0)；④ 求解近似方程x1=x0-f(x0)/f'(x0)；⑤ 判断|x1-x0|是否满足精度要求，若满足，则停止迭代；否则，将x0更新为x1，n=n+1，返回步骤③。
3. 雅可比迭代法解线性方程组（1）原理：雅可比迭代法是一种解线性方程组的迭代算法。
对于线性方程组Ax=b，雅可比迭代法的基本思想是利用矩阵A的对角线元素将方程组分解为多个一元线性方程，然后逐个求解。
（2）实验步骤：① 输入系数矩阵A和常数向量b；② 初始化迭代次数n=0，近似解向量x0；③ 计算对角线元素d1, d2, ..., dn；④ 更新近似解向量x1=x0-A/d1, x2=x0-A/d2, ..., xn=x0-A/dn；⑤ 判断|x1-x0|是否满足精度要求，若满足，则停止迭代；否则，将x0更新为x1, x2, ..., xn，n=n+1，返回步骤③。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
时间(t) 0 1 2 3 4 5 6 7 8 9 10 11 12 温度 15 14 14 14 14 15 16 18 20 22 23 25 28 (x(t))
时间(t) 13 14 15 16 17 18 19 20 21 22 23 24
温度 (x(t))
31 32 31 29 27 25 24 22 20 18 17 16
（1）
1x 0 4 x2 dx,
将区间8等分；
（2）
9
xdx,
将区间4等分；
1
（3） 6 4sin2 xdx, 将区间6等分； 0
实验作业5
1. 考虑迭代公式 x k 1(x k)0 .9 9x kx k 2,
取初始点 x0 0.5 。分别考察一般迭代法与Steffensen加速 迭代法的收敛情况，以及各自的迭代次数。
的值, 输出函数值.
(4)对(3)中的函数，画出其在区域 2x3 , 上1的等y 值7 线，等值 线的值分别为0,1,2,3,4,5,10,15,20,30,40,50, 共12条，并在等 值线上标出等值线的值。
实验作业2
对以下问题,编写相应的M文件:
1.研究高次插值的龙格现象。考虑函数
在 [ 1 , 1上] 取 式。
Sa 2 0
1ac2sin2d,
式中，a 是椭圆的半长轴，c是地球中心与轨道中心（椭圆 中心）的距离。令h为近地点距离，H为远地点距离， R=6371（km)为地球半径，则
a (2RHh)/ 2 c (Hh)/ 2
我国第一颗人造地球卫星近地点距离h=439(km),远地点距离 H=3484(km),试求卫星轨道长度。
试用最小二乘法确定这一天的气温变化规律，考虑用下列 类型的函数，计算误差平方和，并作图比较效果。
（1）二次函数 （2）三次函数 （3）四次函数 （4）函数 x(t)aeb(tc)2, 式中 a , b , c 为常数（拟合参数）。
实验作业4
对以下问题,编写相应的M文件:
1. 分别用复化梯形公式与复化Simpson公式计算下列积分， 2. 并与精确积分值相比较，探讨两类积分公式的精度。
2. （机翼加工问题）
已知机翼轮廓上的数据如下表所示：
x / m 0 3 5 7 9 11 12 13 14 15
0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6
y/m
加工时需要x每改变0.1m时的y值，画出相应的轮廓曲线。
实验作业3
对以下问题,编写相应的M文件: 1. 数据的最小二乘拟合 假定某天的气温变化记录如下表所示
2. 考虑一般迭代公式 x k 1 x k ( 1 x k ) ,k 1 ,2 L .
取 [0.2,4] 中不同的值，x0 (0,1) 进行迭代，画出不同
情况下的 x k 的图形，并分析 取值与 x k 图形的关系。
2. 定积分的应用（选用适当的求积函数计算定积分）。 地球卫星飞行轨道是一个椭圆，椭圆周长的计算公式是
பைடு நூலகம்
f(x)1215x2,x[1,1]
个n 等1距节点
xk1n 2 构k造,k拉0 格,1朗,2.日..,n 插值多项
Ln ( x )
（1）分别画出 n 2 ,4 ,6 ,8 ,1 0 ,1 的2 ,1 拉6 ,格2 0 朗日插值函数和函数
y 的f (图x) 形。 （2）画出 n 2 ,4 ,6 ,8 ,1 0 ,1 2 ,1 6 ,2 0情况下的分段线性插值函数图形。
数值分析 实验
实验作业1
对以下问题,编写相应的M文件:
(1)编写一个求任意 m矩n阵的最大值及其所处的位置的通用程序, 并计算任意随机 产生的 矩阵4 的5最大值及其所处位置。
20
(2)编程求 n 的! 值。 n 1
(3)有一函数 f(x ,y ) x 4 2 x 2 y x 2 2 x y , 写2 一y 2 程 序9 x , 输4 入y 自4 变量 2