matlab定积分的近似计算.pptVIP

* * * * Matlab演示 * * * * * Matlab演示 * * 用Matlab演示 * 用Matlab演示 * 用Matlab演示 * 例：计算二重积分 dblquad(@(x,y)4*x*y+3*x.^2 , -1, 1, 0, 2) 指定 x、y 分别是第一和第二积分变量 dblquad(inline(4*x*y+3*x.^2) ,-1, 1, 0, 2) 被积函数 f (x,y) 的另一种定义方法：匿名函数 dblquad(@(y,x)4*x*y+3*x.^2 , -1, 1 , 0, 2 ) 下面的命令运行结果和上面的一样吗？ dblquad 举例 * int(f,a,b) 计算 f 关于默认自变量 的定积分，积分区间为[a,b]。 int(f) 计算 f 关于默认自变量 的不定积分。 int(f,v,a,b) 计算函数 f 关于自变量 v 的定积分，积分区间为 [a, b] int(f,v) 计算函数 f 关于自变量 v 的不定积分 findsym(f,1) int 符号积分： int * 例：用 int 函数计算定积分： 解： syms x; f=1/(1+x^2); int(f,x,0,1) f=sym(1/(1+x^2)); int(f,x,0,1) int(1/(1+x^2),x,0,1) 或 int(1/(1+x^2),0,1) 或 或 int 举例 * double(a) 将 a 转化为双精度型，若 a 是字符，则取对应的 ASCII 码 a=3; double(a) double(a) 例： ans = 3 ans = 97 相关函数 * x=1:0.001:2; y=exp(x.^(-2)); trapz(x,y) 梯形法： 抛物线法： quad(exp(x.^(-2)),1,2,10e-10) 符号积分法： syms x int(exp(x^(-2)),x,1,2) 例：用 Matlab 函数近似计算定积分 数值实验 * 抛物线法： dblquad(inline(x+y^2),0,2,-1,1) 符号积分法： f=int(x+y^2,y,-1,1); int(f,x,0,2) 数值实验 例：用 Matlab 函数近似计算二重积分 1. 导数、单调性与极值 当f ’(x0)0, 函数在x0点附近是上升的， f ’(x0)0,函数在x0点附近是下降的； 函数在x0点达到局部极大(或局部极小)的 充分条件是f ’(x0)=0 且f ’’(x0)0(或f ’’(x0)0) 考虑函数 f(x)=x2cos(x2+3x-4) 在 [-2, 2]内的图象特征。 建模实验：奶油蛋糕 ? 2 奶油蛋糕 某数学家的学生要送一个特大的蛋糕来庆贺他90岁生日。为了纪念他提出的口腔医学的悬链线模型，学生们要求蛋糕店老板将蛋糕边缘半径作成下列悬链线函数 r = 2-(exp(2h)+exp(-2h))/5, 0h1 (单位：米) 。问如何计算重量？ 解 设高为H，半径 r, 比重为k 若蛋糕是单层圆盘的，则蛋糕的重量为： W=k?Hr2 r H r1 r2 若蛋糕是双层的，每层高H/2 ，下层半径r1, 上层半径r2，则 W=k?H(r12+ r22)/2 如果蛋糕是n层的，每层高H/n，半径分别r1, …, rn, 则 若蛋糕边缘是曲线r = r(h), 0hH，各层半径近似为ri = r((i-1/2)H/n), i=1, …, n, 那么 当n??， 3 一半径为5m的球形水罐充满了水，底部有一半径为b=0.1m的小孔漏水， 问多少时间以后，水面下降至离底部0.5m? 解 水从孔漏出的速度由下列能量方程决定 g(z+R)=u2/2 ， u是速度, z表示从球心 测量的水面高度, g为重力加速度。 考虑在时间dt内水面变化dz，漏水的体积为 uAdt= - ?x2dz 其中x为高度z水面的半径, A=?b2 由于R2=z2+x2 得 dt = 在顶部z=R水降到0.5m时， z=0.5-R，从而 t = 5m 0.5m 0 * 上机作业 1.分别用梯形法与抛物线法，计算 取n = 120，并常识直接使用函数trapz()、quad()进行计算求解，比较结果的差异。 2.试计算定积分 注意：可以运用trapz()、quad() * * * * * * * * * * * * * * * * * * * MATLAB数学建模与仿真 定积分的近似计算 * 定积分计算的基本公式是牛顿－莱