工程中的数值分析Word下载.docx

1、（3）迭代输入初值x，输入迭代格式，并往下复制下去（4）在输入f的计算公式，往下复制下去，通过观察数值是否收敛，若收敛，则取收敛到后面的数值；若发散，则更改定义迭代格式，再重新重复以上步骤进行计算。x3-x+1区间端点a=-1b=xf（x）-0.9-0.629-0.8-0.312-0.7-0.043-0.60.184-0.50.375-0.40.536-0.30.673-0.20.792-0.10.899迭代式:xk+1=（xk-1）1/311-0.4999938 1.812-0.4999979 1.313-0.4999993 14-0.4999998 15-0.4999999 1.116-0

2、.5000000 1.417181.9191.3752021f（x19）=1.375不同迭代格式的收敛性：假定迭代函数（1）对任意（2）存在正数L1，使对任意则迭代过程对于任意初值（3）若方程有根，。1.3 Newton迭代法的原理和算法及Excel实现。原理：Newton迭代法的基本思想是“以直代曲”，将f（x）=0在每一步近似为线性方程来求解，具体方法如下：将f（x）在xk作Taylor一阶展开f（x）=f（xk）+f（xk）（x-xk）+1/2!f（）（x-xk）2,介于x和xk之间.略去上式中的二次项，得到线性方程，解出x，作为新的近似根xk+1：xk+1=xk-f（xk）/f（xk）

3、,k=0,1,2,3称为Newton迭代法 先假定方程的有根区间为a,b，计算a,b区间内各个点（整数点）的函数值，当函数值出现f（a0）0时，a0,b0即为方程的有根区间。将有根区间的长度若干等分，求出对应的点的函数值。将此数据绘图，并根据所绘的图求得初始值。求得方程f（x）的一次求导公式f（x），得到迭代公式xk+1=xk-f（xk）/f（xk），将初始值代入迭代公式中计算出下一项的x值，并计算对应的函数值，新的x值代入迭代公式中继续计算出下一项的x值，重复步骤，直到x的值相同不再变化，此x值即为方程的近似解。迭代法求方程x3-x-1确定初值方程化为等价方程，并定义迭代公式为x-（x3-x

4、-1）/3x2-1上图知迭代初值1.412作图数据区-0.7691.2-0.472-0.1030.3441.50.8751.61.4961.72.2133.0323.9595迭代公式为x-（x3-x-1）/3x2-1不动点迭代kxkf（xk）0.09.06038E-0531.79368E-094F（x4）=0,方程解为1.72.1 线性方程组的数值求解的原理和算法及Excel实现。Gauss消去法原理: 设有线性方程组，将其增广矩阵（A丨b）通过初等行变化为（A（n）丨b（n），A（n）为上三角阵，在经过回代解除与原方程组同解的三角形方程组A（n）x=b（n）的解，得到方程组的解。把方程组化为

5、上三角形方程组，做消元的步骤，再做回带的步骤，解上三角形方程组A（n）x=b（n）。Excel实现：x1+x2-4x4=1-x1+4x2+x3+3x4=-2x1+3x2+5x3-4x4=-42x2+2x3-3x4=-2Ab-4-2-36-50.74.3-4.30.3-3.0三角分解法原理：将系数矩阵A分解为两个三角形矩阵的乘积A=LU，进而将原方程组的求解转化为两个三角形方程组的求解。若有三角阵LU,使A=LU,则方程组Ax=b与方程组LUx=b等价，而后者等价于两个三角形线性方程组：Ly=b，Ux=y。将线性方程组的系数矩阵A分解为三角形方程组的乘积LU，称为矩阵A的LU分解；再将线性方程组

6、的求解转换为三角形方程组的求解。A稠密-LU分解法A对称-LDL分解法A正定-LL分解法A三对角线-追赶法新建Excel表格,依次按顺序输入矩阵数据一句矩阵与逆矩阵相乘为单位矩阵原理,依次从A-D列数据从下至上依照公式计算逆矩阵数据上三角形矩阵求逆UU-10.25-0.750.4375-0.253.1 Lagrange插值的原理和算法及Excel实现；将待求的n次多项式插值函数pn（x）改写成另一种表示方式，再利用插值条件确定其中的待定函数，从而求出插值多项式。n=1时，设.作直线方程：令，称为两点式插值或线性插值.时，设令：称为三点式插值或抛物插值.先建立一个Excle数据表:插值节点xiB

7、CDyiEFGH插值点与函数计算值L0L1L2L3L3（x）a在单元格中输入插值点a求基函数L0=（a-B）*（a-C）*（a-E）/（E-F）/（E-G）/（E-H） L1=（a-A）*（a-C）*（a-D）/（F-E）/（F-G）/（F-H）以此类推求至L3,再求出L3（x）.再输入最后一个基函数L3（x）的计算公式：=SUMPRODUCT公式得到f（x）的近似值2.5-0.06250.562517.5点数:1001.030.94589550.0877635-0.04321350.009554518.2956131.060.8935640.171108-0.0829080.01823618

8、.5727041.090.84297850.2501145-0.11916450.026071518.8316511.120.7941120.324864-0.1520640.03308819.0728321.150.74693750.3954375-0.18168750.039312519.2966251.180.7014280.461916-0.2081160.04477219.5034083.2 Newton插值的原理和算法及Excel实现。牛顿插值通过求各阶差商，递推得到的一个公式：f（x）=fx0+fx0,x1（x-x0）+fx0,x1,x2（x-x0）（x-x1）+.fx0,.xn

9、（x-x0）.（x-xn-1）+Rn（x）。改写记则两点公式可改为：三点公式可改为：这种插值形式的基函数为，.，系数称为差商（均差）. l0一阶二阶三阶四阶五阶（1）计算差商表假设n次输入一阶差商的计算公式“=（B-A）/（2-1）”以此类推往下拉输入二阶差商的计算公式用一阶的值相隔两数相减除以x对应相隔两数相减的值,以此类推往下拉三阶,四阶,N阶如此算下去（2）计算插值点处的函数值输入插值点；分别输入Newdon插值函数N1，N2N-1的计算公式；分别得到插值点处的1阶至n-1阶插值函数值.24差商表fi8-3.5-0.975-5.55.7-3.7-1011.5-7.5-11-9N1N2N3

10、N43.733.617.53515.39313.8668251.0512.412.5662512.504512.0718640612.813.11513.00112.21582513.213.6462513.48912.13.614.1613.96812.69121.2514.6562514.437513.0087890614.415.13514.89713.3708251.3514.815.5962515.34613.4.1 数据拟合的最小二乘法的原理和算法；当实验提供了大量数据时,由于观测数据往往不准确,因此不能要求y=f（x）通过所有点,只要求i=f（xi）-yi（i=1,2,m）严格为

11、零,使近似曲线尽量反映所给数据点的变化趋势同时偏差平方和最小，常采用欧式范数作为误差度量的标准,此即称为最小二乘法原理。关于最小二乘法的一般提法是：对给定的一组数据（xi,yi）（i=0,1，.，m），要求在函数类=Span0（x）,1（x）,2（x）,n（x）中求函数 （） 1使误差平方和 2为了使问题的提法更有一般性，通常在最小二乘法中考虑加权平方和 （4.3）处的数据比重不同，称为权系数，例如可表示在点处重复观测的次数。 按条件式（4.3）求函数的方法称为数据拟合的最小二乘法，用几何语言，即称为曲线拟合的最小二乘法。为最小二乘解，S（x）为拟合。函数。4.2 直线拟合最小二乘法的Exce

12、l实现建立Excle数据表,输入实验数据输入拟合多项式的次数列出法方程组在B6:F9中并输入计算公式计算出结果.之后分解方程组再回代入方程中,并且计算平方误差,作图X2.85.4Y2.128.141.9W次数法方程组13.783.656.93381.81解法方程组6.8541.8-11.3.230.9.7P（X）3.117.点数1.2421.2840.11.3260.44.3 曲线拟合最小二乘法的Excel实现。建立Excle数据表,输入实验数据,依照数据变化趋势设想y=f（x）的方程,再用线性函数S（u）来拟合数据.将数据取倒数变换到下方,再有法方程组输入公式计算,进行矩阵分解以及回代结果.

13、计算平方误差最后确定初值输出作图数据.实验数据t7910y6.58.018.799.39.59.8610.210.3210.4210.5110.5810.6210.7变化数据0.50.20.1250.0625w0.80.095969293.30.6平方误差0.0003289674.74.64.1.455.61.755.36.02.056.2.26.22.356.95.1 数值积分的原理和算法；将函数图形与x轴形成的图形等分求面积即求其积分.从不同角度出发，通过各种途径来构造数值求积公式，常用的一个方法是，利用插值多项式来构造数值求积公式，具体做法如下：在积分区间a,b上取一组点：a=x0x1x

14、n=b,做f（x）的n次插值多项式： 其中lk（x）（k=0,1，n）为n次Lagrange插值基函数，用Ln（x）近似代替被基函数f（x），则有：若记得数值求积公式：xk称为求积节点，Ak称为求积系数例如把图形分成n份,n=1时用梯形公式,n=2时用Sinmpson公式,n=4时用Cotes公式计算代入将每一小块求和5.2 数值积分的的Excel实现；建立一个Excle数据表,在节点区输入节点值于B列,之后计算积分精确值最后运用梯形公式,Sinmpson公式与Cotes公式计算核对节点-1.5函数值积分值f（-2）f（-1.5）f（-1）f（-0.5）f（0）精确值梯形值Simpson值Co

15、tes值x22.252.7x3-8-3.375-0.125x45.06256.46.7ex0.比较.6.常微分方程的数值解法的原理和算法；采取“进步式”和“离散化”。“进步式”是指求解过程依节点排列的次序一步一步地向前推进。描述这类算法，只需给出用已知信息yn,yn-1,yn-2,计算yn+1的递推公式.“离散化”是指通过一定的方法将连续的问题转化为关于离散变量的相应问题。“离散化”的常见方法有：直接用磋商代替微商发、Taylor级数展开法、数值积分法等。一阶方程的初值问题y=f（x，y），x属于a,b，y（a）=y0只要函数f（x，y）在axb，|y|+上连续，且关于y满足Lipschitz

16、条件：|f（x,y1）-f（x,y2）|L|y1-y2|,则方程存在唯一解y=y（x）。所谓微分方程数值解法，就是需求解函数y（x）在一系列离散节点上的近似值：yi60;y（xi）,ax1x2Xn=b.通常采用等距节点Xi=a+ih，i=0,1,2,n,其中h=（b-a）/n称为步长。常微分方程的数值解法的的Excel实现建立Excel数据表,在基本数据区域输入常微分方程的初步数据和步长值,计算节点A列输入序数值B列求出节点dy/dx=f（x,y）,y（x0）=y0,先计算节点之后用Euler法写出求解公式计算值并用改进Euler求解公式计算值各自复制后面,最后作图基本数据x0y0h数值解Euler法改进Euler法精确解iy（xi）