西南交大数值分析上机实习报告
西南交通大学数值分析上机实验报告
数值分析上机实习报告学号：姓名：专业：联系电话：任课教师：序 (3)一、必做题 (4)1、问题一 (4)1.1 问题重述 (4)1.2 实验方法介绍 (4)1.3 实验结果 (5)2、问题二 (7)2.1 问题重述 (7)2.2 实验原理 (7)雅各比算法：将系数矩阵A分解为：A=L+U+D，则推到的最后迭代公式为： (8)2.3 实验结果 (8)二、选做题 (10)3、问题三 (10)3.1 问题重述 (10)3.2 实验原理 (10)3.3 实验结果 (11)总结 (11)序伴随着计算机技术的飞速发展，所有的学科都走向定量化和准确化，从而产生了一系列的计算性的学科分支，而数值计算方法就是解决计算问题的桥梁和工具。
数值计算方法，是一种研究并解决数学问题的数值近似解方法，是在计算机上使用的解数学问题的方法。
为了提高计算能力，需要结合计算能力与计算效率，因此，用来解决数值计算的软件因为高效率的计算凸显的十分重要。
数值方法是用来解决数值问题的计算公式，而数值方法的有效性需要根据其方法本身的好坏以及数值本身的好坏来综合判断。
数值计算方法计算的结果大多数都是近似值，但是理论的严密性又要求我们不仅要掌握将基本的算法，还要了解必要的误差分析，以验证计算结果的可靠性。
数值计算一般涉及的计算对象是微积分，线性代数，常微分方程中的数学问题，从而对应解决实际中的工程技术问题。
在借助MA TLAB、JA V A、C++ 和VB软件解决数学模型求解过程中，可以极大的提高计算效率。
本实验采用的是MATLAB软件来解决数值计算问题。
MA TLAB是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，其对解决矩阵运算、绘制函数/数据图像等有非常高的效率。
本文采用MATLAB对多项式拟合、雅雅格比法与高斯－赛德尔迭代法求解方程组迭代求解，对Runge-Kutta 4阶算法进行编程，并通过实例求解验证了其可行性，使用不同方法对计算进行比较，得出不同方法的收敛性与迭代次数的多少，比较各种方法的精确度和解的收敛速度。
数值分析2024上机实验报告
数值分析2024上机实验报告数值分析是计算数学的一个重要分支，它研究如何用数值方法来解决数学问题。
在数值分析的学习过程中，学生需要通过上机实验来巩固理论知识，并学会使用相应的数值方法来解决实际问题。
本篇报告将详细介绍2024年度数值分析上机实验的内容和结果。
一、实验内容2024年度数值分析上机实验分为四个部分，分别是：方程求根、插值与拟合、数值积分和常微分方程的数值解。
1.方程求根这部分实验要求使用数值方法求解给定的非线性方程的根。
常见的数值方法有二分法、牛顿法、割线法等。
在实验过程中，我们需要熟悉这些数值方法的原理和实现步骤，并对不同方法的收敛性进行分析和比较。
2.插值与拟合这部分实验要求使用插值和拟合方法对给定的一组数据进行拟合。
插值方法包括拉格朗日插值、牛顿插值等；拟合方法包括最小二乘拟合、多项式拟合等。
在实验中，我们需要熟悉插值和拟合方法的原理和实现步骤，并对不同方法的精度和稳定性进行比较。
3.数值积分这部分实验要求使用数值方法计算给定函数的积分。
常见的数值积分方法有梯形法则、辛普森法则、龙贝格积分等。
在实验过程中，我们需要熟悉这些数值积分方法的原理和实现步骤，并对不同方法的精度和效率进行比较。
4.常微分方程的数值解这部分实验要求使用数值方法求解给定的常微分方程初值问题。
常见的数值方法有欧拉法、改进的欧拉法、四阶龙格-库塔法等。
在实验中，我们需要熟悉这些数值解方法的原理和实现步骤，并对不同方法的精度和稳定性进行比较。
二、实验结果在完成2024年度数值分析上机实验后，我们得到了以下实验结果：1.方程求根我们实现了二分法、牛顿法和割线法，并对比了它们的收敛速度和稳定性。
结果表明，割线法的收敛速度最快，但在一些情况下可能会出现振荡；二分法和牛顿法的收敛速度相对较慢，但稳定性较好。
2.插值与拟合我们实现了拉格朗日插值和最小二乘拟合，并对比了它们的拟合效果和精度。
结果表明，拉格朗日插值在小区间上拟合效果较好，但在大区间上可能出现振荡；最小二乘拟合在整体上拟合效果较好，但可能出现过拟合。
数值分析上机实习报告
指导教师：姓名：学号：专业：联系电话：上海交通大学目录序言 (3)实验课题(一) 雅可比迭代法和高斯-塞得尔迭代法的收敛性和收敛速度 (4)数值分析 (6)实验课题(二) 松弛因子对SOR法收敛速度的影响 (6)数值分析 (12)总结 (13)附录（程序清单） (14)1.雅可比迭代法和高斯-塞得尔迭代法的收敛性和收敛速度 (14)雅可比迭代法： (14)高斯-塞得尔迭代法： (16)2.松弛因子对SOR法收敛速度的影响 (18)松弛法（SOR） (18)序言随着科学技术的发展，提出了大量复杂的数值计算问题，在实际解决这些计算问题的长期过程中，形成了计算方法这门学科，专门研究各种数学问题的数值解法（近似解法），包括方法的构造和求解过程的误差分析，是一门内容丰富，有自身理论体系的实用性很强的学科。
解决工程问题，往往需要处理很多数学模型，这就要花费大量的人力和时间，但是还有不少数学模型无法用解析法得到解。
使用数值方法并利用计算机，就可以克服这些困难。
事实上，科学计算已经与理论分析、科学实验成为平行的研究和解决科技问题的科学手段，经常被科技工作者所采用。
作为科学计算的核心内容——数值分析（数值计算方法），已逐渐成为广大科技工作者必备的基本知识并越来越被人重视。
由于数值方法是解数值问题的系列计算公式，所以数值方法是否有效，不但与方法本身的好坏有关，而且与数值问题本身的好坏也有关，因此，研究数值方法时，不但需要研究数值方法的好坏，即数值稳定性问题，而且还需要研究数值问题本身的好坏，即数值问题的性态，以及它们的判别问题。
数值计算的绝大部分方法都具有近似性，而其理论又具有严密的科学性，方法的近似值正是建立在理论的严密性基础上，根据计算方法的这一特点。
因此不仅要求掌握和使用算法，还要重视必要的误差分析，以保证计算结果的可靠性。
数值计算还具有应用性强的特点，计算方法的绝大部分方法如求微分方程近似解，求积分近似值，求解超越方程，解线性方程组等都具有较强的实用性，而插值法，最小二乘法，样条函数等也都是工程技术领域中常用的，有实际应用价值的方法。
数值分析第一次上机练习实验报告
数值分析第一次上机练习实验报告一、实验目的本次实验旨在通过上机练习，加深对数值分析方法的理解，并掌握实际应用中的数值计算方法。
二、实验内容1. 数值计算的基本概念和方法在本次实验中，我们首先回顾了数值计算的基本概念和方法。
数值计算是一种通过计算机进行数值近似的方法，其包括近似解的计算、误差分析和稳定性分析等内容。
2. 方程求解的数值方法接下来，我们学习了方程求解的数值方法。
方程求解是数值分析中非常重要的一部分，其目的是找到方程的实数或复数解。
我们学习了二分法、牛顿法和割线法等常用的数值求解方法，并对它们的原理和步骤进行了理论学习。
3. 插值和拟合插值和拟合是数值分析中常用的数值逼近方法。
在本次实验中，我们学习了插值和拟合的基本原理，并介绍了常见的插值方法，如拉格朗日插值和牛顿插值。
我们还学习了最小二乘拟合方法，如线性拟合和多项式拟合方法。
4. 数值积分和数值微分数值积分和数值微分是数值分析中的两个重要内容。
在本次实验中，我们学习了数值积分和数值微分的基本原理，并介绍了常用的数值积分方法，如梯形法和辛卜生公式。
我们还学习了数值微分的数值方法，如差商法和牛顿插值法。
5. 常微分方程的数值解法常微分方程是物理和工程问题中常见的数学模型，在本次实验中，我们学习了常微分方程的数值解法，包括欧拉法和四阶龙格-库塔法。
我们学习了这些方法的步骤和原理，并通过具体的实例进行了演示。
三、实验结果及分析通过本次实验，我们深入理解了数值分析的基本原理和方法。
我们通过实际操作，掌握了方程求解、插值和拟合、数值积分和数值微分以及常微分方程的数值解法等数值计算方法。
实验结果表明，在使用数值计算方法时，我们要注意误差的控制和结果的稳定性。
根据实验结果，我们可以对计算结果进行误差分析，并选择适当的数值方法和参数来提高计算的精度和稳定性。
此外，在实际应用中，我们还需要根据具体问题的特点和条件选择合适的数值方法和算法。
四、实验总结通过本次实验，我们对数值分析的基本原理和方法有了更加深入的了解。
西南交通大学-数值分析--报告
数值分析课程上机作业计算报告班级：学号：姓名：专业：大地测量学及测绘工程指导老师：联系电话：西南交通大学-数值分析--报告序言通过数值分析的理论知识的学习，此次实验将我们学过的理论知识运用于实践之中。
本次实验，我选用的计算机语言为MATLAB，其主要有一下几个特点。
1．编程效率高MATLAB是一种面向科学与工程计算的高级语言，允许使用数学形式的语言编写程序，且比BASIC、FORTRAN和C等语言更加接近我们书写计算公式的思维方式，用MATLAB编写程序犹如在演算纸上排列出公式与求解问题。
因此，MATLAB语言也可通俗地称为演算纸式科学算法语言。
由于它编写简单，所以编程效率高，易学易懂2. 用户使用方便MATLAB语言与其他语言相比，较好的解决了上述问题，把编辑、编译、链接和执行融为一体。
它能在同一画面上进行灵活操作，快速排除输入程序中的书写错误、语法错误以至语义错误，从而加快了用户编写、修改和调试程序的速度，可以说在编程和调试过程中它是一种比VB还要简单的语言。
3. 方便的绘图功能MATLAB的绘图是十分方便的，它有一系列绘图函数（命令），例如线性坐标、对数坐标、半对数坐标及极坐标，均只需调用不同的绘图函数（命令），在图上标出图题、XY轴标注，格（栅）绘制也只需调用相应的命令，简单易行。
另外，在调用绘图函数时调整自变量可绘出不变颜色的点、线、复线或多重线。
这种为科学研究着想的设计是通用的编程语言所不能及的。
目录1.实验一 (1)1.1题目 (1)1.2计算思路 (1)1.3计算结果 (1)1.4总结 (6)2.第二题 (7)2.1题目 (7)2.2 松弛思想分析 (7)2.3问题的求解 (7)2.4总结 (10)3.第三题 (11)3.1题目 (11)3.2 Runge-Kutta法的基本思想 (11)3.3 问题的求解 (12)3.4问题的总结 (14)总结 (15)附件 (16)实验一程序设计 (16)实验二程序设计 (16)实验三程序设计 (17)实验一：插值问题1.1题目已知：a=-5,b=5, 以下是某函数f(x)的一些点（x k,y k）, 其中x k=a+0.1(k-1) ,k=1,..,101；（数据略）。
数值分析2016上机实验报告
序言数值分析是计算数学的范畴，有时也称它为计算数学、计算方法、数值方法等，其研究对象是各种数学问题的数值方法的设计、分析及其有关的数学理论和具体实现的一门学科，它是一个数学分支。
是科学与工程计算（科学计算）的理论支持。
许多科学与工程实际问题（核武器的研制、导弹的发射、气象预报）的解决都离不开科学计算。
目前，试验、理论、计算已成为人类进行科学活动的三大方法。
数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现。
现在面向数值分析问题的计算机软件有：C,C++,MATLAB,Python,Fortran等。
MATLAB是matrix laboratory的英文缩写，它是由美国Mathwork公司于1967年推出的适合用于不同规格计算机和各种操纵系统的数学软件包，现已发展成为一种功能强大的计算机语言，特别适合用于科学和工程计算。
目前，MATLAB应用非常广泛，主要用于算法开发、数据可视化、数值计算和数据分析等，除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。
本实验报告使用了MATLAB软件。
对不动点迭代，函数逼近（lagrange插值，三次样条插值，最小二乘拟合），追赶法求解矩阵的解，4RungeKutta方法求解，欧拉法及改进欧拉法等算法做了简单的计算模拟实践。
并比较了各种算法的优劣性，得到了对数值分析这们学科良好的理解，对以后的科研数值分析能力有了极大的提高。
目录序言 (1)问题一非线性方程数值解法 (3)1.1 计算题目 (3)1.2 迭代法分析 (3)1.3计算结果分析及结论 (4)问题二追赶法解三对角矩阵 (5)2.1 问题 (5)2.2 问题分析（追赶法） (6)2.3 计算结果 (7)问题三函数拟合 (7)3.1 计算题目 (7)3.2 题目分析 (7)3.3 结果比较 (12)问题四欧拉法解微分方程 (14)4.1 计算题目 (14)4.2.1 方程的准确解 (14)4.2.2 Euler方法求解 (14)4.2.3改进欧拉方法 (16)问题五四阶龙格-库塔计算常微分方程初值问题 (17)5.1 计算题目 (17)5.2 四阶龙格-库塔方法分析 (18)5.3 程序流程图 (18)5.4 标准四阶Runge-Kutta法Matlab实现 (19)5.5 计算结果及比较 (20)问题六舍入误差观察 (22)6.1 计算题目 (22)6.2 计算结果 (22)6.3 结论 (23)7 总结 (24)附录问题一 非线性方程数值解法1.1 计算题目编写不动点迭代法求根程序：把方程010423=-+x x 写成至少四种x=g （x ）的形式，取初值5.1x 0=，进行不动点迭代求根，并比较收敛性及收敛速度。
数值分析上机实践报告
数值分析上机实践报告一、实验目的本次实验主要目的是通过上机操作，加深对数值分析算法的理解，并熟悉使用Matlab进行数值计算的基本方法。
在具体实验中，我们将实现三种常见的数值分析算法：二分法、牛顿法和追赶法，分别应用于解决非线性方程、方程组和线性方程组的求解问题。
二、实验原理与方法1.二分法二分法是一种常见的求解非线性方程的数值方法。
根据函数在给定区间端点处的函数值的符号，不断缩小区间的长度，直到满足精度要求。
2.牛顿法牛顿法是求解方程的一种迭代方法，通过构造方程的泰勒展开式进行近似求解。
根据泰勒展式可以得到迭代公式，利用迭代公式不断逼近方程的解。
3.追赶法追赶法是用于求解三对角线性方程组的一种直接求解方法。
通过构造追赶矩阵，采用较为简便的向前追赶和向后追赶的方法进行计算。
本次实验中，我们选择了一组非线性方程、方程组和线性方程组进行求解。
具体的实验步骤如下：1.调用二分法函数，通过输入给定区间的上下界、截止误差和最大迭代次数，得到非线性方程的数值解。
2.调用牛顿法函数，通过输入初始迭代点、截止误差和最大迭代次数，得到方程组的数值解。
3.调用追赶法函数，通过输入追赶矩阵的三个向量与结果向量，得到线性方程组的数值解。
三、实验结果与分析在进行实验过程中，我们分别给定了不同的参数，通过调用相应的函数得到了实验结果。
下面是实验结果的汇总及分析。
1.非线性方程的数值解我们通过使用二分法对非线性方程进行求解，给定了区间的上下界、截止误差和最大迭代次数。
实验结果显示，根据给定的输入，我们得到了方程的数值解。
通过与解析解进行比较，可以发现二分法得到的数值解与解析解的误差在可接受范围内，说明二分法是有效的。
2.方程组的数值解我们通过使用牛顿法对方程组进行求解，给定了初始迭代点、截止误差和最大迭代次数。
实验结果显示，根据给定的输入，我们得到了方程组的数值解。
与解析解进行比较，同样可以发现牛顿法得到的数值解与解析解的误差在可接受范围内，说明牛顿法是有效的。
数值分析上机实践报告
数值分析上机实践报告一、实验目的本实验的目的是通过编写数值分析程序，掌握解决数学问题的数值计算方法，并通过实际应用来检验其有效性和准确性。
具体包括以下几个方面的内容：1.掌握二分法和牛顿迭代法的基本原理和实现方法；2.熟悉利用矩阵的LU分解和追赶法解线性方程组的过程；3.通过具体的实例应用，比较不同方法的计算效果和精度。
二、实验内容本实验分为三个部分，每个部分包括一个具体的数学问题和相应的数值计算方法。
1.问题一：求方程f(x)=x^3-5x^2+10x-80=0的近似解。
在问题一中，我们通过二分法和牛顿迭代法来求解方程的近似解，并比较两种方法的精度和收敛速度。
2.问题二：用LU分解解线性方程组。
问题二中，我们通过矩阵的LU分解方法解线性方程组Ax=b，然后和直接用追赶法解线性方程组进行对比，验证LU分解的有效性和准确性。
三、实验结果及分析1.问题一的结果分析：通过二分法和牛顿迭代法求解方程f(x)=x^3-5x^2+10x-80=0的近似解，得到的结果如下：从结果来看，两种方法得到的近似解均与真实解x≈5非常接近。
但是，通过比较可以发现，牛顿迭代法的计算速度比二分法更快，迭代的次数更少。
因此，在需要高精度近似解的情况下，牛顿迭代法是一个更好的选择。
2.问题二的结果分析：通过LU分解和追赶法解线性方程组Ax=b，得到的结果如下：-用LU分解解线性方程组得到的结果为x1≈1.0，x2≈2.0，x3≈3.0；-用追赶法解线性方程组得到的结果为x1≈1.0，x2≈2.0，x3≈3.0。
从结果来看，两种方法得到的结果完全一致，而且与真实解非常接近。
这表明LU分解方法和追赶法均可以有效地解决线性方程组问题。
但是，在实际应用中，当方程组规模较大时，LU分解方法的计算复杂度较高，因此追赶法更加适用。
四、实验总结通过本实验，我掌握了二分法和牛顿迭代法以及LU分解和追赶法的基本原理和实现方法。
通过具体的数学问题实例应用，我比较了不同方法的计算效果和精度，得出以下结论：1.在求解函数的近似解时，牛顿迭代法相对于二分法具有更快的收敛速度和更高的计算精度；2.在解决线性方程组问题时，LU分解方法在计算准确性方面与追赶法相当，但在处理较大规模的问题时，计算复杂度较高，追赶法更适合。
西南交通大学数值分析上机实习
目录解题： (1)题目一： (1)1.1计算结果 (1)1.2结果分析 (1)题目二： (2)2.1计算结果 (2)2.2结果分析 (3)题目三： (4)3.1计算结果 (4)3.2结果分析 (5)总结 (5)附录 (6)Matlab程序： (6)题目一： (6)第一问Newton法： (6)第二问Newton法： (6)第一问Steffensen加速法： (7)第二问Steffensen加速法： (7)题目二 (8)1、Jacobi迭代法 (8)2、Causs-Seidel迭代法 (8)题目三： (9)题目一：分别用牛顿法，及基于牛顿算法下的Steffensen 加速法(1)求ln(x +sin x )=0的根。
初值x0分别取0.1, 1,1.5, 2, 4进行计算。
(2)求sin x =0的根。
初值x0分别取1，1.4，1.6, 1.8，3进行计算。
分析其中遇到的现象与问题。
1.1计算结果求ln(x +sin x )=0的根，可变行为求解x-sinx-1=0的根。
1.2结果分析从结果对比我们可发现牛顿—Steffensen 加速法比牛顿法要收敛的快，牛顿法对于初值的选取特别重要，比如第（1）问中的初值为4的情况，100次内没有迭代出来收敛解，而用Steffensen 加速法，7次迭代可得；在第（2）问中的初值为1.6的情况，收敛解得31.4159，分析其原因应该是x x f cos )('=，x0=1.62π≈,0)('≈x f ；迭代式在迭代过程中会出现分母趋近于0，程序自动停止迭代的情况，此时得到的x 往往非常大，而在第一问中我们如果转化为用x+sinx=1,则可以收敛到结果。
用雅格比法与高斯－赛德尔迭代法解下列方程组Ax=b，研究其收敛性，上机验证理论分析是否正确，比较它们的收敛速度，观察右端项对迭代收敛有无影响。
(1)A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4]；b1=[-3,2,4]T,b2=[100,-200,345]T,(2) A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1]；b1=[3,2,1]T,b2=[5,0,-10]T,(3)A行分别为A1=[1,3],A2=[-7,1]；b=[4,6]T2.1计算结果初值均为0矩阵带入(1)A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4]；b1=[-3,2,4]T,b2=[100,-200,345]T2) A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1]；b1=[3,2,1]T,b2=[5,0,-10]TT2.2结果分析ρ小于1,故方程组雅可比迭代收第一小题的经计算谱半径为5427B(=).0敛。
西南交通大学《数值分析报告》上机报告
目录目录 (1)序言（1） (2)1.1 C语言简介及结构 (2)1.2 C语言特点及优点 (2)1.3 选用原因 (3)第一题 (3)用雅格比法与高斯－赛德尔迭代法解下列方程组 (3)1.1题目 (4)1.2原理和思路 (4)1.3计算结果与分析 (7)第二题 (10)松弛因子对SOR法收敛速度的影响 (10)2.1题目 (10)2.2原理和思路 (10)2.3计算结果与分析 (12)序言（2） (16)第三题 (17)利用四阶Runge-Kutta算法求解微分方程的初值问题 (17)3.1题目 (17)3.2原理和思路 (17)3.3计算结果与分析 (18)附录1 Jacobi迭代法C语言源程序 (21)附录2 Gauss-Seidel迭代法程序代码 (23)附录3 SOR迭代法C语言源程序 (25)附录4 四阶Runge-Kutta算法程序代码 (27)总结与体会 (29)本报告系西南交通大学2011级硕士研究生《数值分析》课程的上机实习报告，由本人严格按照实习要求独立完成。
序言（1）在第一次给定的四道上机题中，我选择的是第三题（雅格比迭代法、高斯—赛德尔迭代法求解方程组的问题）和第四题（松弛因子对SOR法收敛速度的影响），本次上机实习基于Microsoft Visual 平台进行程序建立，采用C语言面向对象语言，从界面设计到结果输出，完成一个具有针对性的可视化Windows 应用程序的编制。
在此序言中，我将阐述C语言的基本结构、优特点以及选用这种语言进行编程的主要原因。
1.1C语言简介及结构C语言是一种计算机程序设计语言，由美国贝尔研究所的D.M.Ritchie于1972年推出。
1978后，C语言已先后被移植到大、中、小及微型机上，是目前世界上流行、使用最广泛的高级程序设计语言之一。
C语言既有高级语言的特点，又具有汇编语言的特点。
它不仅可以作为工作系统设计语言，编写系统应用程序；也可以作为应用程序设计语言，编写不依赖计算机硬件的应用程序。
数值分析实习报告总结
一、实习背景数值分析是数学的一个重要分支，它研究如何用数值方法求解数学问题。
随着计算机技术的飞速发展，数值分析在各个领域得到了广泛的应用。
为了提高自己的实践能力，我选择了数值分析作为实习课题，希望通过这次实习，能够掌握数值分析的基本方法，并将其应用于实际问题中。
二、实习过程1. 实习初期在实习初期，我首先了解了数值分析的基本概念、理论和方法。
通过阅读相关教材和文献，我对数值分析有了初步的认识。
接着，我学习了数值分析的基本方法，如泰勒展开、牛顿法、高斯消元法等。
2. 实习中期在实习中期，我选择了几个实际问题进行数值计算。
首先，我使用泰勒展开法求解一个简单的微分方程。
通过编写程序，我得到了微分方程的近似解。
然后，我运用牛顿法求解一个非线性方程组。
在实际计算过程中，我遇到了一些问题，如收敛性、迭代次数过多等。
通过查阅资料和请教导师，我找到了解决方法，成功求解了方程组。
3. 实习后期在实习后期，我进一步学习了数值分析的高级方法，如复化梯形公式、复化Simpson公式、自适应梯形法等。
这些方法在解决实际问题中具有更高的精度和效率。
我选择了一个具体的工程问题，运用复化梯形公式求解定积分。
在计算过程中，我遇到了区间细分、精度控制等问题。
通过不断尝试和调整，我得到了较为精确的积分值。
三、实习收获与体会1. 理论与实践相结合通过这次实习，我深刻体会到理论与实践相结合的重要性。
在实习过程中，我不仅学习了数值分析的理论知识，还将其应用于实际问题中。
这使我更加深刻地理解了数值分析的基本方法，提高了自己的实践能力。
2. 严谨的学术态度在实习过程中，我养成了严谨的学术态度。
在编写程序、进行数值计算时，我注重细节，力求精确。
这使我更加注重学术规范，提高了自己的学术素养。
3. 团队合作精神实习过程中，我与其他同学进行了交流与合作。
在解决实际问题时，我们互相学习、互相帮助，共同完成了实习任务。
这使我更加懂得团队合作的重要性，提高了自己的团队协作能力。
数值分析上机实习报告
数值分析上机实习报告随着现代科学技术的迅猛发展，计算机科学的应用日益广泛，数值分析作为计算机科学中重要的分支之一，其在工程、物理、生物学等领域的应用也越来越受到重视。
本学期，我们在数值分析课程的学习中，进行了多次上机实习，通过实习，我们对数值分析的基本方法和算法有了更深入的理解和掌握。
在实习过程中，我们使用了MATLAB软件作为主要的工具，MATLAB是一种功能强大的数学软件，它提供了丰富的数值计算函数和图形显示功能，使我们能够更加方便地进行数值计算和分析。
第一次实习是线性插值和函数逼近。
我们学习了利用已知数据点构造插值函数的方法，并通过MATLAB软件实现了线性插值和拉格朗日插值。
通过实习，我们了解了插值的基本原理，掌握了插值的计算方法，并能够利用MATLAB软件进行插值计算。
第二次实习是解线性方程组。
我们学习了高斯消元法、列主元高斯消元法和克莱姆法则等解线性方程组的方法，并通过MATLAB软件实现了这些算法。
在实习过程中，我们通过实际例子了解了这些算法的应用，掌握了它们的计算步骤，并能够利用MATLAB软件准确地求解线性方程组。
第三次实习是求解非线性方程和方程组。
我们学习了二分法、牛顿法、弦截法和迭代法等求解非线性方程的方法，以及雅可比法和高斯-赛德尔法等求解非线性方程组的方法。
通过实习，我们了解了非线性方程和方程组的求解方法，掌握了它们的计算步骤，并能够利用MATLAB软件求解实际问题。
通过这次上机实习，我们不仅深入学习了数值分析的基本方法和算法，而且锻炼了利用MATLAB软件进行数值计算和分析的能力。
同时，我们也认识到了数值分析在实际问题中的应用价值，增强了解决实际问题的能力。
总之，这次上机实习使我们受益匪浅，对我们学习数值分析课程起到了很好的辅助作用。
数值分析上机实习报告(西南交通大学)
数值分析上机实习报告姓名：学号：专业：大地测量学与测量工程电话：序言1.所用程序语言：本次数值分析上机实习采用Visual c#作为程序设计语言，利用Visual c#可视化的编程实现方法，采用对话框形式进行设计计算程序界面，并将结果用表格或文档的格式给出。
2.程序概述：（1）第一题是采用牛顿法和steffensen法分别对两个题进行分析，编好程序后分别带入不同的初值，观察与真实值的差别，分析出初值对结果的影响，分析两种方法的收敛速度。
（2）第二题使用Visual c#程序设计语言完成了“松弛因子对SOR法收敛速度的影响”，通过在可视化界面下输入不同的n和w值，点击按钮直接可看到迭代次数及计算结果，观察了不同的松弛因子w对收敛速度的影响。
目录一．用牛顿法，及牛顿-Steffensen法............ 错误！未定义书签。
1. 计算结果.................................... 错误！未定义书签。
2. 结果分析 (5)3. 程序清单 (5)二．松弛因子对SOR法收敛速度的影响 (8)1. 迭代次数计算结果 (8)2. 计算X()结果 (10)3. 对比分析 (12)4. 程序清单： (12)三．实习总结 (14)实验课题（一）用牛顿法，及牛顿-Steffensen法题目：分别用牛顿法，及牛顿-Steffensen法(1)求ln(x+sin x)=0的根。
初值x0分别取0.1, 1,1.5, 2, 4进行计算。
(2)求sin x=0的根。
初值x0分别取1，1.4，1.6, 1.8，3进行计算。
分析其中遇到的现象与问题。
1、计算结果由于比较多每种方法中只选取了其中两个的图片例在下面：2、结果分析通过对以上的牛顿法和steffensen法的练习，我发现在初值的选取很重要，好的初值选出后可以很快的达到预定的精度，要是选的不好就很慢，而且在有的时候得出的还是非数字，所以初始值的选取很重要。
西南交大实习报告(共8篇)
西南交大实习报告（共8篇）第1篇：工程造价实习报告(西南交大)一、毕业实习是我们毕业前的重要环节之一，是学生对以后社会生活的适应起着决定性的作用。
具有较强的实践性、社会性，其目的是通过实习更好的实现理论与实践相结合，使学生毕业前巩固和校验学过的理论知识、锻炼和提高学生的实践技能、提高分析和解决实际问题的能力，提高学生对专业现状及未来的认识，增强个人在社会上的生存能力，为毕业后参加工作，能更好的完成各项实际工作，实现自我价值，打下一个良好的基础。
二、2021年1月1日至今，我在重庆城建控股（集团）有限责任公司市政部门做预算员。
这是第一次正式接触预算工作，每天在规定的时间上下班，上班期间要认真准时地完成自己的工作任务，不能草率敷衍了事。
我们的肩上开始扛着责任，凡事得谨慎小心，否则随时可能要为一个小小的错误承担严重的后果付出巨大的代价，再也不是一句对不起和一纸道歉书所能解决。
三、实习期间我也下工地进行了实践，看了开挖沟槽土方时当开挖沟槽深度大于5m时施工单位必须向监理单位及建设单位重新提出专项施工方案，同时在施工过程中放坡系数也与沟槽深度小于5m的段落不同。
管道安装、沟槽回填、检查井钢筋绑扎、检查井模板安装、检查井砼的浇筑、道路路基和路面施工。
对道路工程的施工顺序和施工工艺有了一定的了解后，在计算工程量的时候有很大的帮助。
四、绑扎钢筋专门看了一下，以前只是老师说钢筋在一个工程中占据的费用很大。
现在亲眼所见，果不其然，在施工现场，放眼望去整条道路上的检查井基本都在施工钢筋，工人们正在忙碌的绑钢筋，大家分工明确，都很认真。
我也前去绑扎了点点钢筋，加工处理了一些箍筋，但是最后还是有很多不合格，最后还是现场的师傅给我讲解了一些原理后我才自己动手加工出来看几个合格箍筋。
其实我觉得我们能动手操作就是我们最大的进步，在学校我根本就不能把书本上的东西用于实践中去，这次我就好好的实践了一番，可以说是受益匪浅。
五、我们在学校里学的不是知识，而是一种叫做自学的能力。
数值分析上机实践报告
数值分析上机实践报告班级：计算机1002姓名：陈斯琪学号：20102686课题三A ． 实验题目：线性方程组的迭代法B ． 实验要求（1） 应用迭代法求解线性方程组，并与直接法作比较；（2） 分别对不同精度要求，如5-4-3-10，10，10=ε，利用所需迭代次数体会该迭代法的收敛快慢；（3） 对方程组（2），（3）使用SOR 方法时，选取松弛因子=0.8,0.9,1,1.1,1.2等，试观察对算法收敛性的影响，并找出你所选用松弛因子的最佳值； （4） 编制出各种迭代法的程序并给出计算结果。
C ． 目的和意义（1） 通过上机了解迭代法求解线性方程组的特点；掌握求解线性方程组的各类迭代法；（2） 体会上机计算时，终止准则‖X^(k+1)-X^k ‖∞<ε，对控制迭代精度的有效性； （3） 体会初始值和松弛因子的选择，对迭代收敛速度的影响 D ． 实验方程组（1）线性方程组⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡1-421534100368-24-3-81-012029137-2621-234179-11-1003524-31-23-6217758-6233-761-62911-31-512-301-231-2-2010563-5-6000121-3-20416084-0484⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡10987654321x x x x x x x x x x =⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-2119381346323125精确解Tx )2,1,1,3,0,2,1,0,1,1(*--=.(2) 对称正定线性方程组⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡---=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡45152211236601924-3-360024-3-36014110-3-5211144-3-310-4221-8-13-4-1-612-53-8-1141-2312-1-204204-2004204-2487654321x x x x x x x x精确解T*)2,0,1,1,2,0,1,1(--=x .（3）三对角线性方程组⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡554141262135741-000000001-0041-0000001-41-0000001-41-0000001-41-0000001-41-0000001-41-0000001-41-0000001-400000001-000000001-410987654321x x xx x x x x x x精确解Tx )1,1,0,3,2,1,0,3,1,2(*---=.E ． 实验程序代码及截图（1） 应用Jacobi 迭代法求解方程组代码如下： #include<iostream.h> #include<math.h>#define N 10 //十阶矩阵 staticdoubleA[N][N]={4,2,-3,-1,2,1,0,0,0,0,8,6,-5,-3,6,5,0,1,0,0,4,2,-2,-1,3,2,-1,0,3,1,0,-2,1,5,-1,3,-1,1,9,4,-4,2,6,-1,6,7,-3,3,2,3,8,6,-8,5,7,17,2,6,-3,5,0,2,-1,3,-4,2,5,3,0,1,16,10,-11,-9,17,34,2,-1,2,2,4,6,2,-7,13,9,2,0,12,4,0,0,-1,8,-3,-24,-8,6,3,-1};//方程组左侧系数矩阵 static double B[N]={5,12,3,2,3,46,13,38,19,-21}; //右侧值static double Y[N]; //输出比较项static double Y[N];static double X[N]; //输出项static double G[N]; //X = BX' + G的G矩阵int i,j,k; //计数器double eps;int M=100;bool distance(){ //求两输出项的差的范数是否满足精度要求double temp=0;for (i=0;i<N;i++){temp=temp+fabs(X[i]-Y[i]);}if (temp>eps)return false;elsereturn true; //满足精度要求则结束程序}void main(){cout<<"最大迭代次数为100次"<<endl;cout<<"你希望的精度是多少？"<<endl;cout<<"eps=";cin>>eps;//形成迭代矩阵B，存放到A中for (i=0;i<N;i++){if (fabs(A[i][i])<eps){cout <<"打印失败"<<endl;return;}double T=A[i][i];for (j=0;j<N;j++){A[i][j]=-A[i][j]/T;}A[i][i] = 0;G[i]=B[i]/T;}int counter=0;while (counter<M){//迭代for (i=0;i<N;i++){double temp=0;for (j=0;j<N;j++){temp=temp+A[i][j]*Y[j];}X[i]=G[i]+temp;}if (distance()==true)break;else{//交换X，Y向量；for(i=0;i<N;i++){Y[i]=X[i];}}counter++;}//打印Xcout << "迭代次数为："<<counter<<"次。
西南交通大学研究 数值分析上机实习报告2012
数值分析上机实习报告要求1．应提交一份完整的实习报告。
具体要求如下：（1）要有封面，封面上要标明姓名、学号、专业和联系电话；（2）要有序言，说明所用语言及简要优、特点，说明选用的考量；（3）要有目录，指明题目、程序、计算结果，图标和分析等内容所在位置，作到信息简明而完全；（4）要有总结，全方位总结机编程计算的心得体会；（5）尽量使报告清晰明了，一般可将计算结果、图表及对比分析放在前面，程序清单作为附录放在后面，程序中关键部分要有中文说明或标注，指明该部分的功能和作用。
2．程序需完好保存到期末考试后的一个星期，以便老师索取用于验证、询问或质疑部分内容。
3．认真完成实验内容，可以达到既学习计算方法又提高计算能力的目的，还可以切身体会书本内容之精妙所在，期间可以得到很多乐趣。
4．拷贝或抄袭他人结果是不良行为，将视为不合格。
5．报告打印后按要求的时间提交给任课老师。
上机实习必须在规定的时间内完成，可要求在考前或考后一个星期内提交。
不合格者和不交者不通过数值分析2012年上机试题1. 已知：a=-5,b=5, 以下是某函数f(x)的一些点（x k,y k）, 其中x k=a+0.1(k-1) ,k=1,..,101x k=a+0.1k,请用插值类方法给出函数f(x)的一个解决方案和具体结果。
并通过实验考虑下列问题实验前分析：所给的节点一共有101个，用Lagrange插值发最高可以做次数100的插值多项式做低次插值是需要进行分组，间插值区域分成互不重叠的区间。
分组时须按区间位置高低的顺序依次划分，有可能最后余下的区间节点数不足够做n次的Lagrange插值多项式，需要特殊处理。
分完组后在分得的区间套用插值公式得到：(1)(2)(),,,m n n n L L L ,那么最后的插值多项式是分段的多项式：(1)1(2)2()(),(),()(),n nn m nmL x x I L x x I L x L x x I ⎧∈⎪∈⎪=⎨⎪⎪∈⎩实验结果：下面列出其中某些次的插值多项式5()L x 的系数(1) Ln(x)的次数n 越高，逼近f(x)的程度越好？ 答：不是，比如2()1/(1)f x x =+时1201()1()1()()nn j j j n jx Ln x x x x x ωω+=+='+-∑当n →∞时只在 3.63x ≤内收敛，在这区间外是发散的。
西南交通大学计算机上机实习报告
计算机实习报告一、实习目的计算机办公软件的使用是交通运输、交通工程等专业学生的基本功，本实习的目的是使学生能够通过短时间的集训式训练，熟练掌握办公软件的使用方法、CAD的基本使用方法以及SPSS的使用。
本实习是交通运输、交通工程等专业学生非常重要的一个学习环节。
本次实习的目的是：将在第一学年中课堂所学习到的计算机知识在实践中综合运用，以期能够熟练掌握办公室自动化软件、CAD、SPSS的使用方法，为以后的进一步学习和工作打下文档编辑和制作的基础。
二、实习内容1、办公室自动化1、能运用办公室自动化软件的文字编辑软件实现基本的图文编辑、制表、格式修饰以及排版打印操作。
2、利用电子表格软件实现工作表的建立、编辑与格式化，利用电子表格软件中的函数与公式对工作表的数据进行处理，实现数据管理和分析。
3、利用演示文稿软件制作演示文稿，掌握不同视图方式在演示文稿的制作和显示中的不同作用和优势，掌握演示文稿编排的基本操作、动画效果制作和母版的基本操作。
2、制图软件使用1、熟悉CAD的操作界面：标题栏、菜单栏与快捷菜单、工具栏、绘图窗口、命令行与文本窗口、状态行。
2、掌握CAD的基本功能：绘制与编辑图形、标注图形尺寸、输出与打印图形。
三、实习环境与方式1、实习环境（1）软件winXP操作系统环境，以及办公自动化软件。
（2）硬件能够安装需要的办公自动化软件的计算机。
2、实习方式由带队老师统一组织，在校内机房进行。
四、实习心得和体会在时间不长的几天实习中，确确实实感觉受益良多，这次实习的内容实际上趋于简单，而且我一开始认为这几个软件已经被用得熟透了，没有什么好学习的，因此一开始抱着很随便的态度，感觉随便糊弄一下就好了。
但一开始接触实习指导书上的内容，我就知道，井底之蛙四字何谓，看上去很简单的办公软件却还有许多我未曾掌握的功能和用途，word的使用还算顺利，因为大部分的功能和我经常使用的WPS基本相同，只是有些功能的位置不太一样，还有操作的界面也有不同，但上手还不算太难，我只是有些困惑，如果在学校里练习的是微软公司的办公软件，但是以后工作又经常用WPS系列软件岂不是会很难使用，因为位置不同，做一个操作都需要去找，很麻烦。
数值分析实习报告总结
数值分析实习报告总结首先，我想对我所参加的数值分析实习课程表示由衷的感谢。
这次实习让我对数值分析这门学科有了更深入的理解，并且让我在实际操作中掌握了许多有用的技能和知识。
在这篇实习报告总结中，我将回顾我在实习过程中的学习经历，总结我在实习中学到的主要内容，并分享我的一些感悟。
实习的第一周，我主要学习了数值分析的基本概念和方法。
通过阅读教材和参加课堂讨论，我了解了数值分析的重要性以及在工程、科学和商业领域中的应用。
我学习了插值、线性代数、微分方程等数值方法的原理和实现方式。
此外，我还通过实际编程练习，掌握了使用数值分析方法解决实际问题的基本技能。
在实习的第二周，我深入学习了Lagrange插值和数值线性代数。
我了解到Lagrange插值是一种构造多项式以通过一组给定的点的方法，它在插值和逼近方面有广泛的应用。
通过编写代码实现Lagrange插值算法，我学会了如何利用已知的数据点来预测未知的点。
此外，我还学习了数值线性代数中的矩阵运算、特征值问题和线性方程组的求解方法，这些方法对于解决实际问题非常重要。
在实习的第三周，我学习了数值微积分和数值求解微分方程的方法。
我了解到数值微积分是利用数值方法近似计算积分和导数的过程，它在信号处理和物理模拟等领域有广泛应用。
通过编写代码实现数值积分和数值导数算法，我学会了如何近似计算函数的积分和导数。
此外，我还学习了如何使用数值方法求解常微分方程和偏微分方程，这些方法对于解决工程和科学领域中的问题非常重要。
在实习的过程中，我也遇到了一些困难和挑战。
例如，在实现数值算法时，我常常会遇到编程错误和数值误差的问题。
通过与同学和老师的讨论和交流，我学会了如何调试代码和减小数值误差的方法。
这些经验让我更加熟悉编程和数值分析的方法，并且提高了我的问题解决能力。
通过这次数值分析实习，我不仅学到了许多关于数值分析的知识和技能，还提高了自己的编程能力和问题解决能力。
我相信这些知识和技能将在我未来的学习和工作中发挥重要作用。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数值分析上机实习报告要求
1．应提交一份完整的实习报告。
具体要求如下：
（1）要有封面，封面上要标明姓名、学号、专业和联系电话；
（2）要有序言，说明所用语言及简要优、特点，说明选用的考量；
（3）要有目录，指明题目、程序、计算结果，图标和分析等内容所在位置，作到信息简明而完全；
（4）要有总结，全方位总结机编程计算的心得体会；
（5）尽量使报告清晰明了，一般可将计算结果、图表及对比分析放在前面，程序清单作为附录放在后面，程序中关键部分要有中文说明或标注，指明该部分
的功能和作用。
2．程序需完好保存到期末考试后的一个星期，以便老师索取用于验证、询问或质疑部分内容。
3．认真完成实验内容，可以达到既学习计算方法又提高计算能力的目的，还可以切身体会书本内容之精妙所在，期间可以得到很多乐趣。
4．拷贝或抄袭他人结果是不良行为，将视为不合格。
5．报告打印后按要求的时间提交给任课老师。
数值分析上机试题10 （选择其中两个题目）
1．给定三个n阶线性方程组Ax=b,其中A的元素a ij (i,j,=1,...,n)与阶数分别为（1）a ij=(i+j-1)2, n=3,4,5, (9)
（2）a ij=(i+j)2, n=3,4,5
（3）a ij=1/(i+j-1), n=3,4,5,6;
b的元素b i=a i1+a i2+…,a in, (i=1,2,…,n)
已知其准确解为x=(1,1,...,1)T, 用列主元素法分别求解上列方程组。
输出各步主元，解释所遇到的现象。
2．用雅格比法与高斯－赛德尔迭代法解下列方程组Ax=b，研究其收敛性，上机验证理论分析是否正确，比较它们的收敛速度，观察右端项对迭代收敛有无影响。
(1)A行分别为A1=[6,2,-1],A2=[1,4,-2],A3=[-3,1,4]；b1=[-3,2,4]T, b2=[100,-200,345]T,
(2) A行分别为A1=[1,0,8,0.8],A2=[0.8,1,0.8],A3=[0.8,0.8,1]；b1=[3,2,1]T, b2=[5,0,-10]T,
(3)A 行分别为A 1=[1,3],A 2=[-7,1]；b =[4,6]T ,
3. 松弛因子对SOR 法收敛速度的影响。
用SOR 法求解方程组Ax =b ,其中
⎪⎪⎪⎪⎪⎪⎪⎪⎭
⎫ ⎝⎛=⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----3-2-.2-2-3-B 4114114114,...... 要求程序中不存系数矩阵A ,分别对不同的阶数取w=1.1, 1.2, ...,1.9进行迭代，记录近似解x (k)达到||x (k)-x (k-1)||<10-6时所用的迭代次数k ，观察松弛因子对收敛速度的影响，并观察当w ≤0或w ≥2会有什么影响？
4. 某实际问题中两次过程记录测得(x ,y )数据点如下：其中y i =f (x i ),i =1,…,21
第一次的数据为：
x =[ -5.0000 -4.5000 -4.0000 -3.5000 -3.0000 -2.5000 -2.0000
-1.5000 -1.0000 -0.5000 0 0.5000 1.0000 1.5000
2.0000 2.5000
3.0000 3.5000
4.0000 4.5000
5.0000];
y =[-0.0000 -0.0001 -0.0002 -0.0003 -0.0004 -0.0048 -0.0366
-0.1581 -0.3679 -0.3894 0 0.3894 0.3679 0.1581
0.0366 0.0048 0.0004 0.0003 0.0002 0.0001 0.0000]
第二次的数据为：
x =[-5.0000 -4.7000 -4.4000 -4.1000 -3.8000 -3.5000 -3.2000
-2.9000 -2.6000 -2.3000 -2.0000 -1.7000 -1.4000 -1.1000
-0.8000 -0.5000 -0.2000 0.1000 0.4000 0.7000 1.0000
1.3000 1.6000 1.9000
2.2000 2.5000 2.8000
3.1000
3.4000 3.7000
4.0000 4.3000 4.6000 4.9000 ]
y =[ -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0001
-0.0006 -0.0030 -0.0116 -0.0366 -0.0945 -0.1972 -0.3280
-0.4218 -0.3894 -0.1922 0.0990 0.3409 0.4288 0.3679
0.2399 0.1237 0.0514 0.0174 0.0048 0.0011 0.0002
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000]
现考虑分段低次或样条插值，或直接插值，总之用不超过7次的分段或不分段插值多项式近似该过程的对应规律y =f (x )，请给出一种解决方案使得误差尽量小表示尽量简单。
5. 用Ru n ge-Kutt a 4阶算法对初值问题y /=-20*y ，y (0)=1按不同步长求解，用于观察稳定区间的
作用，推荐两种步长h=0.1,0.2。
注：此方程的精确解为：y =e -20x。