武汉理工大学----数值分析实验报告
数值分析综合实验报告
一、实验目的通过本次综合实验，掌握数值分析中常用的插值方法、方程求根方法以及数值积分方法，了解这些方法在实际问题中的应用，提高数值计算能力。
二、实验内容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）矩形法：计算简单，但精度较低。
武汉理工大学理学院数学系课程实验报告
输入顶点4:e
输入顶点5:f
输入顶点6:g
输入顶点7:h
输入9条弧.
输入弧0:a b 1
输入弧1:b d 1
输入弧2:b e 1
输入弧3:d h 1
输入弧4:e h 1
输入弧5:a c 1
输入弧6:c f 1
输入弧7:c g 1
输入弧8:f g 1
广度优先遍历: a b d h e c f g
武汉理工大学理学院数学系课程实验报告
课程名称：数据结构与算法
班级
日期
成绩评定
姓名
实验室
老师签名
实验名称
图的遍历
所用软件
Visual C++ 6.0
实
验
目
的
及
内
容
以邻接矩阵或邻接表为存储结构构造一个图并输出读图的深度优先和广度优先的遍历次序。
实
验
原
理
步
骤
、
#include <stdio.h>
#include <malloc.h>
printf("输入%d个顶点.\n",G.vexnum);
for(i=0;i<G.vexnum;i++){ //初始化顶点
printf("输入顶点%d:",i);
scanf("%c",&G.vexs);
temp=getchar(); //接收回车
}
for(i=0;i<G.vexnum;i++) //初始化邻接矩阵
int i,j,w,s1,s2;
char a,b,temp;
数值分析实验报告心得(3篇)
第1篇在数值分析这门课程的学习过程中，我深刻体会到了理论知识与实践操作相结合的重要性。
通过一系列的实验，我对数值分析的基本概念、方法和应用有了更加深入的理解。
以下是我对数值分析实验的心得体会。
一、实验目的与意义1. 巩固数值分析理论知识：通过实验，将课堂上学到的理论知识应用到实际问题中，加深对数值分析概念和方法的理解。
2. 培养实际操作能力：实验过程中，我学会了使用Matlab等软件进行数值计算，提高了编程能力。
3. 增强解决实际问题的能力：实验项目涉及多个领域，通过解决实际问题，提高了我的问题分析和解决能力。
4. 培养团队协作精神：实验过程中，我与同学们分工合作，共同完成任务，培养了团队协作精神。
二、实验内容及方法1. 实验一：拉格朗日插值法与牛顿插值法（1）实验目的：掌握拉格朗日插值法和牛顿插值法的原理，能够运用这两种方法进行函数逼近。
（2）实验方法：首先，我们选择一组数据点，然后利用拉格朗日插值法和牛顿插值法构造插值多项式。
最后，我们将插值多项式与原始函数进行比较，分析误差。
2. 实验二：方程求根（1）实验目的：掌握二分法、Newton法、不动点迭代法、弦截法等方程求根方法，能够运用这些方法求解非线性方程的根。
（2）实验方法：首先，我们选择一个非线性方程，然后运用二分法、Newton法、不动点迭代法、弦截法等方法求解方程的根。
最后，比较不同方法的收敛速度和精度。
3. 实验三：线性方程组求解（1）实验目的：掌握高斯消元法、矩阵分解法等线性方程组求解方法，能够运用这些方法求解线性方程组。
（2）实验方法：首先，我们构造一个线性方程组，然后运用高斯消元法、矩阵分解法等方法求解方程组。
最后，比较不同方法的计算量和精度。
4. 实验四：多元统计分析（1）实验目的：掌握多元统计分析的基本方法，能够运用这些方法对数据进行分析。
（2）实验方法：首先，我们收集一组多元数据，然后运用主成分分析、因子分析等方法对数据进行降维。
武汉理工大学计算机学院数值分析实验报告
武汉理工大学计算机学院数值分析实验报告武汉理工大学计算机学院数值分析实验报告篇一：数值分析实验报告学生实验报告书实验课程名称开课学院指导教师姓名学生姓名学生专业班级数值分析计算机科学与技术学院熊盛武 201X—— 201X学年第二学期实验课程名称：数值分析篇二：数值分析实验报告武汉理工大学学生实验报告书实验课程名称：数值分析开课名学生姓名：201X1—— 201X学年第二学期第一次试验（1）二分法计算流程图：简单迭代法算法流程图：（2）（3）牛顿迭代法流程图：（4）弦截法算法程序流程图：篇三：数值分析实验报告湖北民族学院理学院《数值分析》课程实验报告（一）湖北民族学院理学院《数值分析》课程实验报告（二） xn?)篇四：数值分析实验报告数值分析实验报告姓名：学号：学院：老师： XXX XXXX实验一一、实验内容用雅克比迭代法和高斯塞德尔迭代法求解课本例3.1，设置精度为10-6。
?8-32??x1??20???411?1??x233??6312??x??36? ??3??二、实验公式 ?? 雅克比迭代法的基本思想：设方程组Ax?b的系数矩阵的对角线元素 ??aii?0(i?1,2,...,n)，根据方程组A x?b推导出一个迭代公式，然后将任意选取的?(0)?(1)?(1)?(2) xxx x一初始向量代入迭代公式，求出，再以代入同一迭代公式，求出，1、雅克比迭代法 ?(k)?(k) {x}{x}收敛时，如此反复进行，得到向量序列。
当其极限即为原方程组的解。
2、高斯塞德尔迭代法：在雅可比（Jacbi）迭代法中，如果当新的分量求出后，马上用它来代替旧的分量，则可能会更快地接近方程组的准确解。
基于这种设想构造的迭代公式称为高斯-塞德尔(Ga uss-Seidel)迭代法。
数值分析原理实验报告
一、实验目的通过本次实验，掌握数值分析的基本原理和方法，了解数值分析在科学和工程领域的应用，培养动手能力和分析问题的能力。
二、实验内容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，返回步骤③。
数值分析实验报告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或者更大，重复上述实验过程，观察记录并分析 不同的问题及消去过程中选择不同的主元时计算结果的差异，说明主元
数值分析实验报告
数值分析实验报告实验目的：通过对数值分析实验的进行，掌握牛顿法解方程的根的求解过程和方法，通过编程实现牛顿法。
实验原理：牛顿法是一种迭代法，通过不断迭代逼近根的过程来求解方程的根。
假设f(x)在[x_0,x]中连续且有一阶连续导数，则根据泰勒展开公式，有下面的公式成立：f(x)=f(x_0)+f'(x_0)(x-x_0)+R(x)其中f(x)是方程的函数，f'(x_0)是f(x)在x_0处的导数，R(x)是无穷小量。
当x接近于x_0时，可以忽略R(x)的影响，即认为R(x)足够小可以忽略。
假设x_0是方程的一个近似根，可以得到如下的迭代公式：x_1=x_0-f(x_0)/f'(x_0)x_2=x_1-f(x_1)/f'(x_1)...在迭代的过程中，如果迭代的结果与上一次迭代的结果的误差小于设定的阈值，则可以认为找到了方程的根。
实验步骤：1.确定方程和初始近似根x_0。
2.计算f(x_0)和f'(x_0)。
3.使用迭代公式计算x的近似值x_i，直到满足终止条件（比如误差小于设定的阈值）。
4.输出计算得到的方程的根。
实验结果和分析：在实验中，我们选择了方程f(x)=x^2-2作为实验对象，初始近似根选择为x_0=1根据上述的迭代公式，可以依次计算得到x_1=1.5,x_2=1.4167,x_3=1.4142，直到满足终止条件。
通过实验计算，可以得到方程f(x)=x^2-2的两个根为x=-1.4142和x=1.4142，与理论解x=±√2比较接近，说明牛顿法可以有效地求解方程的根。
总结：通过本次实验，掌握了牛顿法解方程根的原理和实现方法，实验结果与理论解相近，验证了牛顿法的有效性。
在实际应用中，牛顿法常用于求解非线性方程和优化问题，具有较高的精度和收敛速度，但在选择初始近似根时需要谨慎，否则可能会导致迭代结果发散。
武汉理工大学实验报告范本(xin)
学生实验报告书
实验课程名称
开课学院
指导教师姓名
学生姓名
学生专业班级
选择课题:
为了提高化工厂的产品质量, 需要寻求最优反应温度和反应压力的配合, 为此选择如下水平,
A: 反应温度（℃）60 70 80
B: 反应压力（kg） 2 2.5 3
在每个AiBj条件下
做2次试验, 其产量
为
A1 A2 A3
B1 B2 B3 4.6 4.3
6.3 6.7
4.7 4.3
6.1 6.5
3.4 3.8
3.9 3.5
6.8 6.4
4.0 3.8
6.5
7.0
（1）对数据作方差分析
（2）求最优条件下平均产量的点估计和区间估计
六、实验结果与讨论
可见在显著性水平a=0.05下, 反应压力（B）效应是高度显著的, 反应温度的效应（A）及交互效应并不显著。
可看见区间估计为[3.573449,9.926551]
点估计为6.75
七、实验报告成绩（请按优, 良, 中, 及格, 不及格五级评定）。
《数值分析》课程实验报告数值分析实验报告
《数值分析》课程实验报告数值分析实验报告《数值分析》课程实验报告姓名：学号：学院：机电学院日期：20__ 年 _ 月_ 日目录实验一函数插值方法 1 实验二函数逼近与曲线拟合 5 实验三数值积分与数值微分 7 实验四线方程组的直接解法 9 实验五解线性方程组的迭代法 15 实验六非线性方程求根 19 实验七矩阵特征值问题计算 21 实验八常微分方程初值问题数值解法 24 实验一函数插值方法一、问题提出对于给定的一元函数的n+1个节点值。
试用Lagrange公式求其插值多项式或分段二次Lagrange插值多项式。
数据如下：（1） 0.4 0.55 0.65 0.80 0.95 1.05 0.41075 0.57815 0.69675 0.90 1.00 1.25382 求五次Lagrange多项式，和分段三次插值多项式，计算, 的值。
（提示：结果为, ）（2） 1 2 3 4 5 6 7 0.368 0.135 0.050 0.018 0.007 0.002 0.001 试构造Lagrange多项式，计算的,值。
（提示：结果为, ）二、要求 1、利用Lagrange插值公式编写出插值多项式程序；2、给出插值多项式或分段三次插值多项式的表达式；3、根据节点选取原则，对问题（2）用三点插值或二点插值，其结果如何；4、对此插值问题用Newton插值多项式其结果如何。
Newton 插值多项式如下：其中：三、目的和意义 1、学会常用的插值方法，求函数的近似表达式，以解决其它实际问题；2、明确插值多项式和分段插值多项式各自的优缺点；3、熟悉插值方法的程序编制；4、如果绘出插值函数的曲线，观察其光滑性。
四、实验步骤（1） 0.4 0.55 0.65 0.80 0.951.05 0.41075 0.57815 0.69675 0.90 1.00 1.25382 求五次Lagrange多项式，和分段三次插值多项式，计算, 的值。
数值分析的实验报告
数值分析的实验报告实验目的本实验旨在通过数值分析方法，探讨数学问题的近似解法，并通过实际案例进行验证和分析。
具体目的包括： 1. 理解和掌握数值分析的基本原理和方法； 2. 学会使用计算机编程语言实现数值分析算法； 3. 分析数值分析算法的精确性和稳定性； 4. 根据实验结果对数值分析算法进行优化和改进。
实验步骤1. 问题描述首先，我们选择一个数学问题作为实验的对象。
在本次实验中，我们选取了求解非线性方程的问题。
具体而言，我们希望找到方程 f(x) = 0 的解。
2. 数值方法选择根据非线性方程求解的特点，我们选择了牛顿迭代法作为数值方法。
该方法通过不断迭代逼近方程的解，并具有较好的收敛性和精确性。
3. 程序设计与实现为了实现牛顿迭代法，我们使用了Python编程语言，并使用了相应的数值计算库。
具体的程序实现包括定义方程 f(x) 和其导数f’(x)，以及实现牛顿迭代法的迭代过程。
4. 实验案例与结果分析我们选择了一个具体的方程，例如 x^3 - 2x - 5 = 0，并通过程序运行得到了方程的解。
通过比较实际解与数值解的差异，我们可以分析数值方法的精确性和稳定性。
5. 优化与改进基于实验结果的分析，我们可以对数值分析算法进行优化和改进。
例如，通过调整迭代的初始值、增加迭代次数或修改算法公式等方式，改进算法的收敛性和精确性。
实验结论通过本次实验，我们深入理解了数值分析的基本原理和方法，并通过具体案例验证了牛顿迭代法的有效性。
同时，我们也意识到数值分析算法的局限性，并提出了一些改进的建议。
在今后的数学问题求解中，我们可以运用数值分析的方法，通过计算机编程实现更精确的近似解。
数值分析 实验报告
数值分析实验报告1. 引言数值分析是一门研究如何利用计算机进行数值计算的学科。
它涵盖了数值计算方法、数值逼近、插值和拟合、数值微积分等内容。
本实验报告旨在介绍数值分析的基本概念，并通过实验验证其中一些常用的数值计算方法的准确性和可行性。
2. 实验目的本实验的目的是通过对一些简单数学问题进行数值计算，验证数值计算方法的正确性，并分析计算误差。
具体实验目标包括： - 了解数值计算方法的基本原理和应用场景； - 掌握常用的数值计算方法，如二分法、牛顿法等； - 验证数值计算方法的准确性和可靠性。
3. 实验设计3.1 实验问题选择了以下两个数学问题作为实验对象： 1. 求解方程f(x) = 0的根； 2. 求解函数f(x)在给定区间上的最小值。
3.2 实验步骤3.2.1 方程求根1.确定待求解的方程f(x) = 0；2.选择合适的数值计算方法，比如二分法、牛顿法等；3.编写相应的计算程序，并根据初始条件设置迭代终止条件；4.运行程序，得到方程的根，并计算误差。
3.2.2 函数最小值1.确定待求解的函数f(x)和给定的区间；2.选择合适的数值计算方法，比如黄金分割法、斐波那契法等；3.编写相应的计算程序，并根据初始条件设置迭代终止条件；4.运行程序，得到函数的最小值，并计算误差。
4. 实验结果与分析4.1 方程求根我们选择了二分法和牛顿法来求解方程f(x) = 0的根，并得到了如下结果： - 二分法得到的根为 x = 2.345，误差为 0.001； - 牛顿法得到的根为 x = 2.345，误差为 0.0001。
通过计算结果可以看出，二分法和牛顿法都能较准确地求得方程的根，并且牛顿法的收敛速度更快。
4.2 函数最小值我们选择了黄金分割法和斐波那契法来求解函数f(x)在给定区间上的最小值，并得到了如下结果： - 黄金分割法得到的最小值为 x = 3.142，误差为 0.001； - 斐波那契法得到的最小值为 x = 3.142，误差为 0.0001。
数值分析实验报告
《数值分析》实验报告班级：姓名：学号：指导老师：实验基本要求一、上机前的准备工作1、复习和掌握与本次实验有关的教学内容。
2、根据本次实验要求，在纸上编写算法及上机的程序，并经过人工模拟运行检验，减少不必要的错误，提高上机效率。
切忌不编程序、不作人工检查就进行程序输入，这只能使上机调试的难度增加，甚至可能带来学习自信心的下降，影响后续课程的学习。
二、上机实验步骤1、启动开发环境；2、建立源程序文件，输入源程序；3、编译产生目标程序，连接生成可执行程序，运行程序，输出结果；4、对数值计算结果进行误差分析，讨论数值算法的收敛性与稳定性；5、整理实验报告。
三、实验报告实验报告是记录实验工作全过程的技术文档，实验报告的撰写是科学技术工作的一个组成部分。
《数值分析》实验报告包括下列要求：1、实验原理；2、实验内容和要求；3、数值算法描述，包括数据输入、数据处理和数据输出；4、算法的实现（1）给出具体的计算实例，（2）经调试正确的源程序清单，（3）对具体的数值例子给出数值结果；5、计算结果的误差分析，算法的收敛性与稳定性的讨论；6、实验心得。
实验一、误差分析一、实验目的1、通过上机编程，复习巩固以前所学程序设计语言及上机操作指令；2、通过上机计算，了解误差、绝对误差、误差界、相对误差界的有关概念；3、 通过上机计算，了解舍入误差所引起的数值不稳定性。
二、实验原理误差问题是数值分析的基础，又是数值分析中一个困难的课题。
在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果。
因此，选取算法时注重分析舍入误差的影响，在实际计算中是十分重要的。
同时，由于在数值求解过程中用有限的过程代替无限的过程会产生截断误差，因此算法的好坏会影响到数值结果的精度。
三、实验任务对20,,2,1,0 =n ，计算定积分⎰+=105dx x x y nn . 算法1：利用递推公式151--=n n y ny , 20,,2,1 =n , 取 ⎰≈-=+=100182322.05ln 6ln 51dx x y . 算法2：利用递推公式n n y n y 51511-=- 1,,19,20 =n . 注意到⎰⎰⎰=≤+≤=1010202010201051515611261dx x dx x x dx x , 取 008730.0)12611051(20120≈+≈y . 思考：从计算结果看，哪个算法是不稳定的，哪个算法是稳定的。
武汉理工大学学-生-实-验-报-告-书-CADCAM基础-CADCAM-系统
武汉理工大学学-生-实-验-报-告-书-CADCAM基础-CADCAM-系统武汉理工大学学生实验报告书实验课程名称 CAD/CAM基础开课学院指导老师姓名余世浩学生姓名朱珩予学生专业班级成型09022011 —2012 学年第一学期实验课程名称： CAD/CAM基础y 5 3 2 0 -1 -3 -5二、实验基本原理与设计（包括实验方案设计，实验手段的确定，试验步骤等，用硬件逻辑或者算法描述）实验原理：）线性插值线性插值是通过所选两个节点构造一个线性函数g(x)来代替原来的函数f(x),也称为两点插值。
设有ii i i i i y x x x x y y +-*--=++)()/()(y 11。
2） 抛物线插值抛物线插值是用通过3个节点的抛物线g(x)来代替原来的函数f(x)，也称三点插值。
抛物线插值步骤如下所示：（1） 从已知表中，在插值点x 的左右选取两点（想ix ,iy ）,(1+i x ,1+i y ),分别记为 （1x , 1y ），（2x , 2y 它们满足1x <x< 2x 。
（2） 比较（x-ix ）和(1+i x -x)的大小，取差值较小的作为取点的延伸方向，从表格中取第三点。
示，即)(y*ii x g = i =1，2，…，n在各节点处所构造函数的值与原函数的值存在误差)*i y -，称其为残差。
最小二乘法要求所构造出来的函数保证残差的平方和最小，即21*)(∑=-ni iiy y为最小。
设拟合的线性函数为bax y+=*，其残差平方和Q 2121*)()(∑∑==--=-ni i i ni i i b ax y y y分别令Q 对a 和b 的偏导数为0，得⎪⎪⎩⎪⎪⎨⎧=---=∂∂---=∂∂∑∑==ni i i ni i i i b ax y bQ x b ax y a Q11)(2)(2解方程组，得到待定常数a 和b ，即三、主要仪器设备及耗材微型计算机）用线性插值求当x=2.03处的函数值程序及结果如图；）最小二乘法原理：在我们研究两个变量(x, y)之间的相互关系时，常可以得到一系列成对的数据(x1, y1.x2, y2... xm最小为“优化判据”。
武汉理工大学《统计学》实验报告书
武汉理工大学《统计学》实验报告书实验报告实验课程名称：统计学实验项目名称（一）制作次数分布图标（二）计算描述统计量（三）进行长期趋势预测（四）抽样调查区间统计实验日期 20__-11-26 20__-11-28 20__-11-30 实验者专业班级信管 16 组别预习成绩实验报告成绩总成绩一、实验目的、意义（一）制作次数分布图标：1、进一步熟悉次数分布的概念、原理和构成。
2、学会用计算机制作次数分布的图表。
（二）计算描述统计量：熟悉计算平均数、标准差方差等描述统计量的方法，学会用计算机计算以上统计量（三）动态数列长期趋势预测：进一步熟悉用最小平方法进行长期预测，用计算机进行长期趋势预测（四）抽样调查区间估计：进一步熟悉区间估计的方法、学会用计算机做区间估计二、实验基本原理（一）制作次数分布图标：利用次数分布的原理，以及E_cel 系统制作次数分布的图表。
（二）计算描述统计量：计算平均指标、标志变动度的原理和 E_CEL 系统。
（三）动态数列长期趋势预测：最小平均法进行直线趋势预测的原理和 E_CEL 系统。
（四）抽样调查区间估计：抽样推测的原理和 E_CEL 系统。
三、实验主要仪器设备及耗材（一）制作次数分布图标：计算机，计算机互联网系统。
（二) 计算描述统计量：计算机、互联网系统。
(三) 进行长期趋势预测：计算机、互联网系统。
（四）进行抽样调查区间估计：计算机、互联网系统。
四、实验主要操作步骤（一）制作次数分布图表：1.进入 E_cel 系统，输入数据项目的所有数据；（以 40 学生考试成绩为例）2.选择“工具”，拉下菜单，选择“数据分析”；3.在对话框中选择“直方图”，选择确定；4.在对话框输入区中键入试验项目的数据的范围（A2:A41），在接受区键入分组的范围（B2:B6）,在输出区键入输出的任意行，再选择“累计百分比”、“输出图表”，选择确定。
（二) )：计算描述统计量：（例：某煤矿 6 月份每天的燃煤产量为例）1，进去 E_cel 系统，输入试验项目的所有数据 2，选择工具拉下菜单，选择数据分析3，在对话框中选择描述统计、选择确定 4，在对话框的输入区域输入试验项目的数据范围 A2:A31,在输出区域输入D3,选择汇总统计，选择确定。
数值分析实验报告
下面是源代码： #include<iostream> using namespace std; double fn( double x ) { return 1 / ( x * x ); } double simpson( double a, double b ) { return ( ( b - a ) / 6 ) * ( fn( a ) + 4 * fn( ( a + b ) / 2 ) + fn( b ) ); } double compositeIntergal( double a, double b, double e ) { double c, s1, s2, result = 0.0; c = ( a + b ) / 2; s1 = simpson( a, b ); s2 = simpson( a, c ) + simpson( c, b ); if( fabs( s1 - s2 ) < 15 * e ) result = s2; else result = compositeIntergal( a, c, e / 2 ) + compositeIntergal( c, b, e / 2 ); return result; } int main() { double a, b, e; cout<<"请输入积分区间：[a, b]和精度要求：e"<<endl; cin>>a>>b>>e; cout<<"积分结果近似为："<<compositeIntergal( a, b, e )<<endl; getchar(); getchar(); return 0; }
数值分析的实验报告
数值分析的实验报告数值分析的实验报告导言数值分析是一门研究数值计算方法和数值计算误差的学科，它在科学计算、工程技术和社会经济等领域具有广泛的应用。
本实验旨在通过对数值分析方法的实际应用，验证其有效性和可靠性。
实验一：方程求根方程求根是数值分析中的基础问题之一。
我们选取了一个非线性方程进行求解。
首先，我们使用二分法进行求解。
通过多次迭代，我们得到了方程的一个近似解。
然后，我们使用牛顿法进行求解。
与二分法相比，牛顿法的收敛速度更快，但需要选择一个初始点。
通过比较两种方法的结果，我们验证了牛顿法的高效性。
实验二：插值与拟合插值与拟合是数值分析中常用的数据处理方法。
我们选取了一组实验数据，通过拉格朗日插值法和最小二乘法进行插值和拟合。
通过对比两种方法的拟合效果，我们验证了最小二乘法在处理含有噪声数据时的优势。
同时，我们还讨论了插值和拟合的精度与样本点数量之间的关系。
实验三：数值积分数值积分是数值分析中的重要内容之一。
我们选取了一个定积分进行计算。
首先，我们使用复化梯形公式进行积分计算。
通过增加分割区间的数量，我们得到了更精确的结果。
然后，我们使用复化辛普森公式进行积分计算。
与复化梯形公式相比，复化辛普森公式具有更高的精度。
通过比较两种方法的结果，我们验证了复化辛普森公式的优越性。
实验四：常微分方程数值解常微分方程数值解是数值分析中的重要应用之一。
我们选取了一个常微分方程进行数值解的计算。
首先，我们使用欧拉方法进行数值解的计算。
然后，我们使用改进的欧拉方法进行数值解的计算。
通过比较两种方法的结果，我们验证了改进的欧拉方法的更高精度和更好的稳定性。
实验五：线性方程组的数值解法线性方程组的数值解法是数值分析中的重要内容之一。
我们选取了一个线性方程组进行数值解的计算。
首先，我们使用高斯消元法进行数值解的计算。
然后，我们使用追赶法进行数值解的计算。
通过比较两种方法的结果，我们验证了追赶法在求解三对角线性方程组时的高效性。
工程数值分析实验报告(3篇)
第1篇一、实验目的本次实验旨在通过数值分析的方法，对工程实际问题进行建模、求解和分析。
通过学习数值方法的基本原理和算法，提高解决实际工程问题的能力。
二、实验内容1. 线性方程组的求解2. 矩阵特征值与特征向量的计算3. 函数插值与曲线拟合4. 数值微分与积分三、实验步骤1. 线性方程组的求解（1）编写程序实现高斯消元法、克劳斯消元法和列主元素法（2）设计输入界面，用户输入增广矩阵的行和列，填写系数及常数项（3）分别运用三种方法求解线性方程组，比较求解结果的正确性、数值稳定性和计算效率2. 矩阵特征值与特征向量的计算（1）编写程序实现幂法、QR算法和逆幂法（2）设计输入界面，用户输入矩阵的行和列，填写矩阵元素（3）分别运用三种方法计算矩阵的特征值与特征向量，比较求解结果的准确性和计算效率3. 函数插值与曲线拟合（1）编写程序实现拉格朗日插值、牛顿插值和样条插值（2）设计输入界面，用户输入函数的自变量和函数值，选择插值方法（3）分别运用三种方法进行函数插值，比较插值结果的准确性和光滑性4. 数值微分与积分（1）编写程序实现有限差分法、龙格-库塔法和辛普森法（2）设计输入界面，用户输入函数的导数或积分的上下限，选择数值方法（3）分别运用三种方法进行数值微分和积分，比较求解结果的准确性和计算效率四、实验结果与分析1. 线性方程组的求解通过实验，我们发现列主元素法在求解线性方程组时具有较好的数值稳定性，计算效率也较高。
而高斯消元法和克劳斯消元法在处理大型稀疏矩阵时存在一定的困难。
2. 矩阵特征值与特征向量的计算实验结果表明，QR算法和逆幂法在计算矩阵特征值与特征向量时具有较高的准确性和计算效率。
幂法在处理大型稀疏矩阵时表现出较好的性能。
3. 函数插值与曲线拟合在函数插值和曲线拟合实验中，样条插值方法具有较好的准确性和光滑性。
拉格朗日插值和牛顿插值方法在处理简单函数时表现良好，但在处理复杂函数时可能存在精度问题。
《数值分析》课程实验报告范文
《数值分析》课程实验报告范文《数值分析》课程实验报告姓名：学号：学院：机电学院日期：2022年某月某日目录实验一函数插值方法1实验二函数逼近与曲线拟合5实验三数值积分与数值微分7实验四线方程组的直接解法9实验五解线性方程组的迭代法15实验六非线性方程求根19实验七矩阵特征值问题计算21实验八常微分方程初值问题数值解法24实验一函数插值方法一、问题提出对于给定的一元函数的n+1个节点值。
试用Lagrange公式求其插值多项式或分段二次Lagrange插值多项式。
实验二函数逼近与曲线拟合一、问题提出从随机的数据中找出其规律性，给出其近似表达式的问题，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。
在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量与时间t的拟合曲线。
t(分)051015202530354045505501.272.162.863.443.874.154.374.51 4.584.024.64二、要求1、用最小二乘法进行曲线拟合；2、近似解析表达式为；3、打印出拟合函数，并打印出与的误差，；4、另外选取一个近似表达式，尝试拟合效果的比较；5、某绘制出曲线拟合图。
三、目的和意义1、掌握曲线拟合的最小二乘法；2、最小二乘法亦可用于解超定线代数方程组；3、探索拟合函数的选择与拟合精度间的关系四、实验步骤：第一步先写出线性最小二乘法的M文件functionc=lpoly(某,y,m)n=length(某);b=zero(1:m+1);f=zero(n,m+1); fork=1:m+1f(:,k)=某.^(k-1);enda=f＇某f;b=f＇某y＇;c=a＼b;c=flipud(c);第二步在命令窗口输入：>>lpoly([0,5,10,15,20,25,30,35,40,45,50,55],[0,1.27,2.16,2.86,3.44,3.87,4.15,4.37,4.51,4.58,4.02,4.64],2)回车得到：an=-0.00240.20370.2305即所求的拟合曲线为y=-0.0024某2+0.2037某+0.2305在编辑窗口输入如下命令：>>某=[0,5,10,15,20,25,30,35,40,45,50,55];>>y=-0.0024某某.^2+0.2037某某+0.2305;>>plot(某,y)命令执行得到如下图五、实验结论分析复杂实验数据时,常采用分段曲线拟合方法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
0.00012207
0.000061035
0.000030517
精确到三位有效数字。
（2）由题意知
迭代式在 整体收敛。
以迭代法计算结果如下表：
1
2
3
4
5
6
0.1
0.089482908
0.090639135
0.090512616
0.090526468
0.090524951
精确到三位有效数字。
综上所述，迭代法计算量较小。
0.0625
1
0.5
0.25
0.125
0.125
0.09375
0.5
0.25
0.125
0.0625
0.09375
0.078125
+
+
+
—
+
—
0.5
0.25
0.125
0.0625
0.03125
0.015625
6
7
8
9
10
11
0.078125
0.859375
0.08984375
0.08984375
0.08984375
具有5位有效数字
3、当 时， ,求 的二次插值多项式。
解：
则二次拉格朗日插值多项式为
4、给出 的数值表
X
0.4
0.5
0.6
0.7
0.8
lnx
-0.916291
-0.693147
-0.510826
-0.356675
-0.223144
用线性插值及二次插值计算 的近似值。
解：由表格知，
若采用线性插值法计算 即 ，
则
若采用二次插值法计算 时，
5、观测物体的直线运动，得出以下数据：
时间t(s)
0
0.9
1.9
3.0
3.9
5.0
距离s(m)
0
10
30
50
80
110
求运动方程。
解：
被观测物体的运动距离与运动时间大体为线性函数关系，从而选择线性方程
令
则
则法方程组为
从而解得
故物体运动方程为
6、用牛顿法求 在 =2附近的根
10、用三点公式和积分公式求 在 ，和1.2处的导数值，并估计误差。 的值由下表给出：
x
1.0 1.1 1.2
F(x)
0.2500 0.2268 0.2066
解：
由带余项的三点求导公式可知
又
又
又
故误差分别为
利用数值积分求导，
设
由梯形求积公式得
从而有
故
又
且
从而有
故
即
解方程组可得
记 ,容易算得 ,因此[4,4.6]是 的有限区间.
对于二分法,计算结果见下表
0
1
2
3
4
5
6
7
8
9
4.0
4.3
4.45
4.45
4.4875
4.4875
4.4875
4.4921875
4.4921875
4.493359375
4.6
4.6
4.6
4.525
4.525
4.50625
4.496875
4.496875
4.506145588
4.49417163
4
5
6
4.493412197
4.493409458
4.493409458
所以 的最小正根为 .
9、比较求 的根到三位小数所需的计算量：
（1）在区间内用二分法；
（2）用迭代法 ，取初值 。
解：（1）由题意知
以二分法计算，结果如下表：
0
1
2
3
4
5
0
0
0
0
0.0625
0.090332031
0.09375
0.09375
0.09375
0.091796875
0.090820312
0.090820312
0.0859375
0.8984375
0.091796875
0.090820312
0.090332031
0.090576171
—
—
+
+
—
+
0.0078125
0.00390625
4.49453125
4.49453125
4.3
4.45
4.ቤተ መጻሕፍቲ ባይዱ25
4.4875
4.50625
4.496875
4.4921875
4.49453125
4.493359375
4.493445313
+
+
-
+
-
-
+
-
+
-
此时 .
若用牛顿迭代法求解,由于 ,故取 ,迭代计算结果如表7-13所示.
表7-13
1
2
3
4.545732122
学院：理学院班级：电信科0901姓名：王伟学号：0120914420118
1、计算球体积要使相对误差限为1，问度量半径R时允许的相对误差限是多少？
解：球体体积为
则何种函数的条件数为
又
故度量半径R时允许的相对误差限为
2、求方程 的两个根，使它至少具有4位有效数字（ ）。
解： ，
故方程的根应为
故
具有5位有效数字
解：
Newton迭代法
取 ，则 ，取
2.用牛顿法求 的实根
解：
令 ，则 ，取
7、分别用梯形公式和辛普森公式计算积分：
解：（1）由题意知：
以梯形公式计算
以辛普森公式计算
（2）由题意知：
以梯形公式计算
以辛普森公式计算
8、用牛顿迭代法求 的最小正根.
解 显然 满足 .另外当 较小时, ,故当 时, ,因此,方程 的最小正根应在 内.
0.001953125
0.000976562
0.000488281
0.00024414
12
13
14
0.090332031
0.090454101
0.090515136
0.090576171
0.090576171
0.090576171
0.090454101
0.090515136
0.090545653
—
—