第九讲应用微积分.ppt

第九讲 应用微积分 第五章 应用微积分 5.1 预备知识：微积分的基本概念 5.2 数值微积分MATLAB指令 5.3 作业 5.1 预备知识：微积分 多元函数Taylor公式 多元函数极值 梯度：记为?f 。 f x,y 在 x0,y0 取得局部极大或极小的必要条件是 ?f x0,y0 0，充分条件是 ?f x0,y0 0且下列Hesse矩阵负定 局部极大 或正定 局部极小 曲线 曲面 平面曲线 x t , y t , a t b的长度为 空间曲线 x t , y t , z t , a t b的长度为 曲面 z z x,y , x,y in G的面积为 5.2 数值微积分MATLAB指令 例： x [1 2 4 7]; diff x ans 1 2 3 diff x,2 ans 1 1 例1 1、y quad 'exp -x.^2 ',-1,1 y 1.4936 2、y quadl 'exp -x.^2 ',-1,1 y 1.4936 例1 程序： y dblquad 'x.*exp x.^2+y.^2 ',0,2,-2,2 y 881.8304 例2 程序： y triplequad 'y.*sin t +z.*cos t ',0,pi,0,1,-1,1 y 2.0000 作业 1、计算下列多项式的根 2、计算函数 在[0,1]区间上的积分。 3、计算 在 上的二重积分。 * * 1.极限和连续 数列极限: ?? 0, ? N 0 ，使当n N时 有?xn -a? ?，则 函数极限: 如果当x?x0时有f x ? A， 则 连续: 如果当x?x0时，有f x ? f x0 则称 f x 在x0连续。 闭区间上连续函数必有最大值和最小值。 2. 微分与导数 函数f x 在点x x0的导数为 若f x 在x0可导则在x0可微，dy Adx 当f’ x0 0,函数在x0点附近是上升的； 当f’ x0 0,函数在x0点附近是下降的； 当f’ x0 0, x0为驻点, 若x0为驻点且f” x0 0 或f” x0 0 ，则 f x 在x0点达到局部极大（或局部极小） 当n 0 得，微分中值定理 f x - f x0 f’ ? x- x0 其中?是x0 与 x 之间某个值 Taylor公式:当f x 在含有x0某个开区间内具有直到n+1阶的导数， 3.多元函数微分学 设f x,y 在点 x0,y0 附近有定义，当 x,y 以任何方式趋向于 x0,y0 时，f x,y 趋向于一个确定的常数A，则 若 A f x0,y0 , 称f x,y 在 x0,y0 点连续 f x,y 在点 x0,y0 的偏导数分别定义为 二元函数在 x0,y0 附近局部线性化: f x,y ? f x0,y0 + f x’ x0,y0 x-x0 + f y’ x0,y0 y-y0 4 . 积分? 函数f x 在区间[a,b]上的积分定义为 其中 a x0 x1 … xn b, ?xi xi-xi-1, ?i? xi-1, xi , i 1,2,…,n 若在[a,b]上, F’ x f x , 则 二重积分定义为 1.数值微分 若f x 在x a可导, 设h 0且足够小 称为向前差商 向后差商 中心差商 数值差分?: n维向量x x1, x2, ?, xn 的差分定义为 n-1维向量?x x2-x1, x3-x2, ?, xn- xn-1 。 diff x 如果x是向量，返回向量x的差分 如果x是矩阵，则按各列作差分。 可得到列差分矩阵：[x 2:m,: -x 1:m-1,: ] diff x, k k阶差分，即差分k次。 例：diff x,2 就等同于diff diff x 。 2.多项式的导数 q polyder p 求得由向量p表示的多项式的 导函数，并用q向量表示. 例： p [1 0 -2 -5]; q polyder p 3.梯形积分法 z trapz x,y 返回积分的近似值，其中x 表 示积分区间的离散化向量; y是与x同维 数的向量,表示被积函数 。 例1 解 ? x -1:0.1:1; y exp -x.^2 ; ? trapz x,y 4.高精度数值积分 quad Fun,a, b 自适应步长Simpson积分法 求得Fun在区间[a, b]上的定积分z quadl Fun,a,b 高精度Lobatto积分法. 格式同quad. 注：trapz, quad, quad1都不能用于求广义积分，对于一些假奇异积分也不能直接求解 。