高等数学实验第章.pptVIP

1、计算机的局限性 2、广义积分? 广义积分数值求解是一个较困难的问题， 通常数值积分方法都不适用。 例5 (无界域积分) 计算积分 I= 例6、(奇点积分) 计算积分 I= 5.5 扩展实验：计算机可信吗 重积分的数值计算可通过单积分组合计算 我们利用梯形法，先将[a,b]区间m等分，hx=(b-a)/m, xi=a+ihx, i=0,1,…,m 3、重积分 再将[c(xi),d(xi)]区间n等分, hy(i)=(d(xi)-c(xi))/n, yij=c(xi)+jhy(i), 则 G(xi) A a b d(x) c(x) 其中 M文件dblquad2.m给出二重积分计算法。 I=dblquad2(fun,a,b,clo,dhi,m,n) fun为被积函数f(x,y)句柄或嵌入表示， x为标量,y为向量,clo和dhi是y的下限和 上限函数;c(x),d(x), 都是x的标量函数； a,b分别为x的下限和上限; m,n 为x和y方向的等分数(缺省值100) 例1 5.6 习题 * 第五章 应用微积分 5.1 预备知识：微积分的基本概念 5.2 数值微积分 MATLAB命令 5.3 计算实验：数值微积分 5.4 建模实验：奶油蛋糕 5.5 计算机可信吗 5.6 习题 5.1 预备知识：微积分的基本概念 1、极限和连续 2、微分和导数 3、多元函数微分学 4、积分 1、极限和连续 数列极限：??0, ? N0 ，使当nN时，有?xn -a??，则 函数极限：如果当x?x0时有f(x) ? A，则 连续：如果当x?x0时，有f(x)? f(x0)，则称 f(x)在x0连续。 闭区间上的连续函数必有最大值和最小值。 函数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点达到局部极大（或局部极小） 函数在x0处切线的斜率 2、 微分与导数 当n=0 得，微分中值定理 f(x) - f(x0) = f’(?) (x- x0) 其中?是x0与x之间某个值 Taylor公式：当f(x)在含有x0某个开区间内具有直到n+1阶的导数， 设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)的偏导数分别定义为： 3、多元函数微分学 函数f(x)在区间[a,b]上的积分定义为 其中a=x0x1…xn=b,?xi=xi-xi-1,?i?(xi-1,xi), i=1,2,…,n 若在[a,b]上, F’(x)=f(x), 则 二重积分定义为 4、积分 5.2 数值微积分MATLAB命令 1、数值差分 2、求导数 3、梯形积分法 4、函数数值积分法 5、矩形区域重积分 1、数值差分 n维向量x=(x1, x2, ?, xn)的差分定义为n-1维向量 diff(x) 如果x是向量，返回向量x的差分 如果x是矩阵，则按各列作差分。 diff(x, k) k阶差分，即差分k次 ? A=[1 3;5 2;6 5;7 7],B=diff(A),C=diff(A,2) A = 1 3 5 2 6 5 7 7 B = 4 -1 1 3 1 2 C = -3 4 0 -1 q=polyder(p) 求得由向量p表示的多项式的导函数的向量表示q； Fx=gradient(F,x) 返回向量F表示的一元函数沿x方向的导函数F’(x)，其中x是与F同维数的向量； [Fx,Fy]=gradient(F,x,y) 返回矩阵F表示的二元函数的数值梯度(F‘x,F’y),当F为m×n矩阵时, x,y 分别为n维和m维的向量； quiver(X, Y, U, V)在(X,Y)平面点上，画(U,V)表示的方向箭头； 2、求导数 例：画函数 的方向导数图 xa=-2:0.2:2;ya=-2:0.5:2;[x,y]=meshgrid(xa,ya);