数学实验第六次(数值积分与微分)重点.pptVIP

数学实验 之六：数值积分与微分 数值积分与微分 积分和微分是我们在高等数学中学过的两个最基本的运算。它们在实际问题中有着许多直接的应用。 微分运算比较简单，但是遇到由离散数据或图形表示的函数就只能求助于数值方法。 积分运算存在“积不出来”，计算这种类型的定积分只能用数值方法，至于由离散数据或图形表示的函数的定积分，理所当然地属于数值积分的范畴。 一、提出几个用数值积分和微分求解的实际问题 二、介绍数值积分的梯形公式、辛普森公式 三、介绍数值微分的基本方法 四、提出的几个问题的求解过程 一、卫星轨道长度、人口增长率 一、卫星轨道长度、人口增长率 二、数值积分 二、数值积分 二、数值积分 三、数值微分 四、问题求解 * 2008-11 桂林电子科技大学信息科技学院 桂林电子科技大学信息科技学院 问题1：人造地球卫星轨道可视为平面上的椭圆。我国第一颗人造地球卫星近地点距地球表面439km，远地点距地球表面2384km，地球半径为6317km，求该卫星的轨道长度。 问题2：已知20世纪美国人口统计数据如表1，试计算表中这些年份的人口增长率。 又已知某地区20世纪70年代的人口增长率如表2，且1970年人口为210（百万），试估计1980年的人口。 数值积分可用下面几种命令实现： sum(x) 输入数组x，输出为x的和，可用于按矩形公式(1)、(2)计算积分。 cumsum(x)输入数组x，输出数组为x的依次累加和，也可用于按矩形公式计算积分。 trapz(x)输入数组x，输出为按梯形公式(3)计算的x的积分(h=1). trapz(x,y)输入x，y为同长度的数组，输出y对x的积分(按梯形公式计算，但步长不一定相等)。 quad(‘fun’,a,b,tol)用辛普森公式(4)计算以fun.m文件命名的函数（或库函数如’sin’）在区间(a,b)上的积分，相对误差tol（缺省时为10^-6），输出积分值。 quadl(‘fun’,a,b,tol)用自适应Gauss-Lobatto公式计算。 对于表1中给出的数据，若记时刻t的人口为x(t)，则人口（相对）增长率为r(t)=(dx/dt)/x(t),表示每年人口增长的比例。记1900年为k=0，1910至2000依次为k=1…10,k年人口记为xk，年增长率为rk，用数值微分的三点公式计算（将每10年的增长率变为每年的增长率） 人口增长率（续）