西南交通大学2016-2017 数值分析报告
序言
数值分析是计算数学的范畴，有时也称它为计算数学、计算方法、数值方法等，其研究对象是各种数学问题的数值方法的设计、分析及其有关的数学理论和 具体实现的一门学科，它是一个数学分支。是科学与工程计算（科学计算）的理论支持。许多科学与工程实际问题（核武器的研制、导弹的发射、气象预报）的解决都离不开科学计算。目前，试验、理论、计算已成为人类进行科学活动的三大方法。
数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现。现在面向数值分析问题的计算机软件有：C,C++,MATLAB,Python,Fortran等。
MATLAB是matrix laboratory的英文缩写，它是由美国Mathwork公司于1967年推出的适合用于不同规格计算机和各种操纵系统的数学软件包，现已发展成为一种功能强大的计算机语言，特别适合用于科学和工程计算。目前，MATLAB应用非常广泛，主要用于算法开发、数据可视化、数值计算和数据分析等，除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。
本实验报告使用了MATLAB软件。对不动点迭代，函数逼近（lagrange插值，三次样条插值，最小二乘拟合），追赶法求解矩阵的解，4RungeKutta方法求解，欧拉法及改进欧拉法等算法做了简单的计算模拟实践。并比较了各种算法的优劣性，得到了对数值分析这们学科良好的理解，对以后的科研数值分析能力有了极大的提高。
西南交通大学2016-2017 数值分析报告
目录
序言 ........................................................................................................................................................ 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 附录
西南交通大学2016-2017 数值分析报告
问题一 非线性方程数值解法
1.1 计算题目
编写不动点迭代法求根程序：把方程x3?4x2?10?0写成至少四种x=g（x）的形式，取初值x0?1.5，进行不动点迭代求根，并比较收敛性及收敛速度。
1.2 迭代法分析
将方程f（x）=0改写成其等价形式
x??(x)
取方程根的某一近似值x0作为初始点，由函数x1??(x0)可计算出x1，如此下去，设当前点为xk，有?（x）计算出xk+1，即
xk?1??(xk) k=0,1，.....
称为迭代公式。
（收敛条件）设?(x)在[a，b]上有连续的一阶导数，并满足：
(1)?x?[a,b],有a??(x)?b;(2)?'(x)?L?1,?x?[a,b].
则有函数φ(x)在区间[a,b]上存在唯一的不动点（方程的根）x*；对任何x0属于[a,b]，可由迭代公式得到序列{xk}均收敛到方程的根x*。设上述条件成立时，算法的中止条件为：
数值分析2016上机实验报告.wps
