第10章工程数值计算方法实验指导.doc

《第10章工程数值计算方法实验指导.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第10章工程数值计算方法实验指导10.1引言数值计算主要用于解决工程中常见的基本数学问题，研究其相关的数值解法，包含了线性方程组、非线性方程组的解法，矩阵求逆、数值逼近、数值微积分、常微分方程及偏微分方程的数值解法等。数值计算的基本理论和研究方法建立在数学建模的基础上，与计算机科学密切相关。采用计算机解决这些问题的主要步骤是：分析实际问题、建立数学模型、设计算法、编写程序代码并上机实现。本章精选前面章节中重点内容组成了7个实验项目，要求学生上机实现各个实验项目并编写实验报告。实验一线性方程组的直接解——列主元消去法解线性方程组；实验二线性方程组的迭代解——雅可比法、高斯—赛德尔迭代法解线

2、性方程组；实验三非线性方程的近似解——二分法、牛顿法求非线性方程的根；实验四插值问题——拉格朗日插值、牛顿插值；实验五曲线拟合问题——最小二乘法；实验六数值积分——复化辛甫生公式；实验七求解常微分方程的初值问题——改进欧拉方法、四阶龙格库塔方法。10.2实验一线性方程组的直接解——列主元消去法解线性方程组一、实验目的掌握列主元消去法解线性方程组的理论。二、实验环境PC机一台，C语言、MATLAB任选。三、实验内容用Gauss列主元消去法求解方程组。四、实验原理Gauss列主元消去法第一步：消元过程对于,将进行变换为，其中是上三角矩阵。即(1)选列主元选取第列中绝对值最大元素作为主元。(2

3、)换行(3)归一化(4)消元第二步：回代过程由解出。五、实验步骤1．要求上机实验前先编写出程序代码；2．编辑录入程序；3．调试程序并记录调试过程中出现的问题并修改程序；4．记录运行时输入数据和输出结果；5．撰写实验报告。六、思考题进一步思考全主元消去法的实现步骤。七、参考代码#include"stdio.h"#include"math.h"#definen3main(){inti,j,k;intmi;floatmv,tmp;floata[n][n]={{0.01,2,-0.5},{-1,-0.5,2},{5,-4,0.5}};floatb[n]={-5,5,9},x[n];for(k=0

4、;kmv){mi=i;mv=fabs(a[i][k]);}if(mi>k){tmp=b[k];b[k]=b[mi];b[mi]=tmp;for(j=k;j

5、mp;}}x[n-1]=b[n-1]/a[n-1][n-1];for(i=n-2;i>=0;i--){x[i]=b[i];for(j=i+1;j

6、TLAB任选。三、实验内容分别采用雅可比迭代法、高斯—赛德尔迭代法求解线性方程组。四、实验原理1．雅可比迭代法设方程组的系数矩阵的对角线元素，为迭代次数容许的最大值，为容许误差。(1)取初始向量，令。(2)对，计算(3)如果并且，则输出方程的解，结束；如果，则说明方程组不收敛，输出方程组无解，终止程序；否则，转(2)。2．高斯-赛德尔迭代法(1)判断线性方程组是否主对角占优，即判断是否满足。(2)取初始向量，令。(3)直接分离，，并建立高斯-赛德尔迭代格式五、实验步骤1．要求上机实验前先编写出程序代码；2．编辑录入程序；3．调试程序并记录调试过程中出现的问题并修改程序；4．记录运行时输入

7、数据和输出结果；5．撰写实验报告。六、实验注意事项雅可比迭代法和高斯—赛德尔迭代法解线性方程组应注意迭代初值的选取和迭代终止的条件。七、思考题为节省计算时间，提高精度，如何选择合适的数据结构使得迭代收敛的速度更快？八、参考代码1．雅可比迭代法#include#include#include#include#defineEPS1e-6#defineMAX100fl