MATLAB数学实验.ppt

《MATLAB数学实验.ppt》由会员分享，可在线阅读，更多相关《MATLAB数学实验.ppt（28页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、MATLAB数学实验 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望第五章第五章 应用微积分应用微积分5.1 预备知识：微积分的基本概念预备知识：微积分的基本概念5.2 数值微积分数值微积分MATLAB指令指令5.3 计算实验：数值微积分计算实验：数值微积分5.4 建模实验：奶油蛋糕建模实验：奶油蛋糕1.极限和连续数列极限:0,N0，使当nN时有xn-a0,函数在x0点附近是上升的；当f(x0)0,函数在x0点附近是下降的；当f(x0)=0,x0为驻点,当n=0得

2、微分中值定理f(x)-f(x0)=f()(x-x0)其中是x0与x之间某个值Taylor公式:当f(x)在含有x0某个开区间内具有直到n+1阶的导数，Taylor展开的几何解释展开的几何解释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)的偏导数多元函数Taylor公式n n二元函数在二元函数在二元函数在二元函数在(x(x0 0,y,y0 0)附近局部线性化附近局部线性化附近局部线性化附近局部线性化(

3、平面近似曲面平面近似曲面平面近似曲面平面近似曲面):):f f(x,yx,y)f f(x x0 0,y y0 0)+)+f f x x(x x0 0,y y0 0)()(x-xx-x0 0)+)+f f y y(x x0 0,y y0 0)(y-yy-y0 0)函数极值n n一元可微一元可微f(x)在内点在内点x0取得局部极大或极小的必取得局部极大或极小的必要条件是要条件是f(x0)=0.局部极大局部极大(或极小或极小)充分条充分条件是件是f(x0)=0,f(x0)0).n n二元可微二元可微f(x,y)在内点在内点(x0,y0)取得局部极大或极取得局部极大或极小的必要条件是梯度小的必要条件是

4、梯度 f(x0,y0)=0,充分条件是充分条件是 f(x0,y0)=0且下列且下列Hesse矩阵负定矩阵负定(局部极大局部极大)或正定或正定(局部极小局部极小)n n梯度梯度 f=(),Hesse矩阵矩阵4.积分函数f(x)在区间a,b上的积分定义为其中a=x0 x1xn=b,xi=xi-xi-1,i(xi-1,xi),i=1,2,n若在a,b上,F(x)=f(x),则牛顿莱布尼兹公式二重积分定义为曲线曲线&曲面曲面n n平面曲线平面曲线(x(t),y(t),atb的长度为的长度为n n空间曲线空间曲线(x(t),y(t),z(t),atz=quadl(x)exp(-x.2),-1,1)z=i

5、ntegral(x)exp(-x.2),-1,1)注意：积分中注意：积分中fun里用数组运算里用数组运算!fun=(x)exp(-x.2).*log(x).2;q=integral(fun,0,Inf)%integral能求反常积分5.重积分z=dblquad(Fun,a,b,c,d)求得二元函数Fun(x,y)在矩形区域的重积分.z=triplequad(Fun,a,b,c,d,e,f)求得三元函数Fun(x,y,z)在长方体区域上的三重积分。z=quad2d(Fun,a,b,cx,dx)求得二元函数Fun(x,y)的重积分。a,b为变量x的下、上限；cx,dx为变量y的下、上限函数(自变量

6、为x)z=integral2(Fun,a,b,cx,dx)类似quad2dz=integral3(Fun,a,b,cx,dx,exy,fxy)求得三元函数Fun(x,y,z)的三重积分例5.2计算重积分fun=(x,y)1./(sqrt(x+y).*(1+x+y).2);q=integral2(fun,0,1,0,(x)1-x)fun=(z,x,y)log(3+x+y+z)+z.2.*cos(x);xmin=(z)-sqrt(1-z.2);xmax=(z)sqrt(1-z.2);ymin=(z,x)-sqrt(1-x.2-z.2);ymax=(z,x)sqrt(1-x.2-z.2);q=int

7、egral3(fun,0,1,xmin,xmax,ymin,ymax)1.1.数值微分数值微分若f(x)在x=a可导,设h0且足够小称为向前差商向后差商中心差商5.3 计算实验：数值微积分计算实验：数值微积分自动步长的中心差商程序：自动步长的中心差商程序：deriv.mfunction d=deriv(fname,a,h0,e)h=h0;d=(fname(a+h)-fname(a-h)/2/h;d0=d+2*e;%为了保证为了保证abs(d-d0)e成立成立while abs(d-d0)e d0=d;h0=h;h=h0/2 d=(fname(a+h)-fname(a-h)/2/h;endder

8、iv(x)sin(x),pi/6,0.1,1e-4)2.奶油蛋糕例5某数学家的学生要送一个特大的蛋糕来庆贺他90岁生日。为了纪念他提出的口腔医学的悬链线模型，学生们要求蛋糕店老板将蛋糕边缘半径作成下列悬链线函数r=2-(exp(2h)+exp(-2h)/5,0hfun=(h)(2-(exp(2*h)+exp(-2*h)/5).2;pi*integral(fun,0,1)即旋转体体积公式解：t=0时z=R.考虑在时间dt内水面变化dz，漏水的体积为uAdt=-x2dz其中u为出水速度，A为小孔面积，x为高度z水面的半径5m0.5m0例6一半径为R=5m的球形水罐充满了水，底部有一半径为b=0.1

9、m的小孔漏水，问多少时间以后，水面下降至离底部0.5m?分析：求时间分析：求时间t与水面高度与水面高度z的对应关的对应关系。系。z表示从球心表示从球心 测量的水面高度测量的水面高度A=b2,x2=R2-z2注意：水从孔漏出的速度注意：水从孔漏出的速度u不是常数！不是常数！zxt=0时z=R.在顶部水降到0.5m时，z=0.5-R，从而t=0+5m0.5m0u由下列能量方程决定mg(z+R)=mu2/2m是质量，g为重力加速度。%M脚本脚本eg5_6.mclear;R=5;b=0.1;g=9.81;z1=0.5-R;z2=R;n=100;h=(z2-z1)/n;z=z1:h:z2;f=(R2-z.2)./(b2*sqrt(2*g*(z+R);I=trapz(z,f)/60/60n n结果为结果为0.5144小时。小时。Page94 习题习题n nex4n nex5(2)(4)(6)(7)(8)n nex6n nex10n nex12 n nex13