第四讲、微分、积分程序.pptVIP

第四讲、微分、积分程序 求极限、导数和极值 级数求和、Taylor展开 微分、数值微分 若f(x)在x=a可导，设h0且足够小，向前差商、向后差商、中心差商分别为 二阶导数公式 求极限、导数和极值的MATLAB命令 limit 符号求极限； taylor 符号Taylor展开； symsum符号级数求和； polyder 多项式求导； diff 数值微分或符号微分；fmin一元函数极值； gradient 数值梯度； fmins 多元函数极值。 1.数值求导 2. 函数极值 3. 解析运算 limit(s,x,a) 表达式 s当x?a时的极限。 symsum(s,n,a,b) 表达式 s的通项当自变量 n由a到b的和。 diff(s,x,n) 表达式 s对x的n阶导数。 taylor(s,n,a) 表达式 s在a点Taylor展开到n次式。 积分 求积分的MATLAB命令 trapz 梯形法积分； dblquad 矩形区域积分； quad 变步长数值积分； int 符号积分； quad8 高精度数值积分。 梯形积分法 例 1 2. 变步长数值积分 z=quad8(‘Fun’,A,B,Tol) Fun—被积函数的M函数名； A,B—下限，上限； Tol—精度，缺省值为1e-4. 例 3. 重积分 (1) dblquad只能求矩形区域二重积分。 （2） M文件dblquad2.m给出的二重计算方法 用法 I=dblquad2(‘f_name’,a,b,’c_lo’,’d_hi’.m,n) ‘f_name’—为被积函数f(x,y)的字符串，x为标量，y为向量；a,b为x的下上限，‘c_lo’和‘d_hi’是y的下上上限函数c(x),d(x),都是x的标量函数，m,n分别为 x、y方向的等分数（缺省值为100)。 4. 符号积分 int(s) 符号表达式s的不定积分； int(s,v)符号表达式s关于变量v的不定积分; int(s,a,b)符号表达式s的定积分,a,b为下上限； int(s,v,a,b)符号表达式s关于变量v从a到b的定积分。 例3 注意 例1 人口增长率 中国人口统计数据如下(数据来源：国家统计局《2003年中国人口统计数据年鉴》，中国人口与发展研究中心整理)。 1、试计算这些年份的人口增长率； 2、预测未来若干年的人口数据。 表1　中国人口统计数据 例2 教堂顶部曲面面积 某阿拉伯国家有一座著名巨大拱形圆顶的伊斯兰教堂，年久失修，国王下令将教堂顶部重新贴金箔装饰。据记载，大厅的顶部形状为半球面，其半径为30m。考虑到可能的损耗和其他技术因素，实际用量将会比教堂顶部面积多1.55%。据此，国王的财政大臣拨出了可制造5800m2有规定厚度金箔的黄金。建筑商人哈桑略通数学，他计算了一下，觉得黄金会有盈余。于是，他以较低的承包价得到了这项装饰工程。但在施工前的测量中，工程师发现教堂顶部实际上并非是一个精确的半球面而是半椭球面，其半立轴恰是30m，而半长轴和半短轴分别是30.6m和29.6m。这一来哈桑犯了愁，他担心黄金是否还有盈余？甚至可能短缺。最后的结果究竟如何呢？ 例1求解:预测 三种简单的模型： xk=x0(1+r)k （指数） x(t)=x0ert (Malthus) (Logistic) 增长率：可用平均增长率。亦即 0.0148 0.0209 0.0196 0.0159 0.0122 0.0085 的平均。 可把2000年的人口作为x0. 使用Logistic模型要 给定最大的人口容量xm 求解例2，关键在于计算教堂顶部表面积，取椭球面中心为坐标原点建立直角坐标系，则教堂顶部半椭球面的方程可写为 其中R=30，a=30.6，b=29.6，引入广义极坐标变换， , ，不难得到表面积为 这个积分式十分复杂，必须用数值方法 来处理，由于被积表达式在 r=1处的奇 性，不能直接对此式用上述积分方法，再 作一个变换 上式变为： x=0:0.25:1; y=0:pi./2:2*pi; A=[1 2 2 2 1;2 4 4 4 2;2 4 4 4 2;2 4 4 4 2;1 2 2 2 1]; w=0; for j