工程数学(例题).ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 解 做差商表,由于是等距离节点, 由第二类边界条件得 解方程得 将Mi代入式4.4.14)得 由于 故 例题 例5.1 设 求f(x)在区间[0,1]上的一次最佳平方逼近多项式. 解 设 由于 故法方程为 解得 例题 例5.2 设函数y=f(x)的离散数据如下表所示 试用二次多项式拟和上述数据,并求平方误差. 0 1 2 3 4 5 0 0.2 0.4 0.6 0.8 1 1.000 1.221 1.492 1.822 2.226 2.718 解 可得 解方程组得 所以拟合二次函数为 例5.3 地球温室效应问题 下表统计了近100年内地球大气气温上升的数据.试根据表中数据建立一数学模型即拟和曲线,并根据这一模型,预报地球气温何年会比1860年的平均温度高 年份N 1860年后地球气温增加值 年份N 1860年后地球气温增加值 1880 0.01 1940 0.10 1890 0.02 1950 0.13 1900 0.03 1960 0.18 1910 0.04 1970 0.24 1920 0.06 1980 0.32 1930 0.08 解 为简化数据,从1880年起年份记N,其变换n=(N-1870)/10.将地球气温增加值改记为t=1,2,3,4,6,8,10,13,18,24,32,也就是将原气温增加值扩大100倍,根据新数据绘制图 可以看出,气温t与变换n大致服从指数函数增长过程,因此,可以假设t与n满足指数函数关系 为决定参数α,β将上式改写成 记 则有 这是已知数据相应地变为如下表所示 n 1 2 3 4 5 6 7 8 9 10 11 ln1 ln2 ln3 ln4 ln6 ln8 ln10 ln13 ln19 ln24 ln32 取n=1,m=10,并将上表已知数据带入得 解方程组得: 相应的t 与 n 的指数型拟合曲线关系为 就是所求地球温室效应的指数函数的数学模型,以此进行预报,即已知t值求 以地球气温比1860年上升 为例,即以t=700代入上式可得: N(7)=2078(年) 常用的几个积分公式 梯形公式(n=1) Simpson公式(n=2) Newton公式(n=3) Cotes公式(n=4) 例题 例题 例题 * * * * * 例题 例题 Gauss-Seidel迭代法的算法 例题 例题 例题 例题 例题 抛物线插值的精度与正弦函数表完全一样。 （3）相应的误差估计： Newton插值计算 例题 例4.2.1 试用Newton插值公式计算sinx在x=π/12处的近似值。 解 先列差商表如表4.2.2所示，所以得 例题 例题 例题 例题 例题 例4.4.1 已知函数y=f(x)的数表如下表所示。 求满足边界条件 x 0 0.15 0.30 0.45 0.60 f(x) 1 0.97800 0.91743 0.83160 0.73529 Gauss全主元消元算法 选列主元的Gauss-Jordan消元法 改进的cholesky分解算法 改进的cholesky分解算法 例题 例题 例题 例题 例题 Doolittle分解 K=1时 n-1次除法 K=2,…, n 时 (k-1)(n-k+1)次乘法,加减法 k(n-k)次乘除法, (k-1)(n-k)加减法 第二步 次乘除法 例题 例题 例题 例题 算法（原地求逆法） 例题 例题 例题 例题 例题 例题 例题 对称矩阵范数 例题 3.5 病态方程组与矩阵的条件数 3.5.1 病态方程组与扰动方程组的误差分析 病态方程组与扰动方程组的误差分析 病态方程组与扰动方程组的误差分析 病态方程组与扰动方程组的误差分析 病态方程组与扰动方程组的误差分析 例题 例题 例 1.1 求 在区间[1, 2]之间的零点。 解 in 误差不超过10