实用 MATLAB教学资料-7-8.pptVIP

例4　求x=1时的积分 练习被积函数的四种定义方式 q = quad(myfun, a, b) 例5　 quad or quadl 例6　 原理：自适应Gauss-Kronrod积分 四、quadgk 目前MATLAB精度最高的积分函数 匿名函数 myfun=@(x)f(x) 可以计算无穷积分 两种方式定义被积函数 四、quadgk q = quadgk(myfun, a, b) M文件 function y=myfun(x) y=f(x) 四、quadgk q = quadgk(@myfun, a, b) 原理： 自适应积分 五、integral 可以计算反常积分 匿名函数定义 q = integral(myfun, a, b) myfun=@(x)f(x) 例6　 q = quadgk(myfun, a, b) 例7　 quadgk or integral 六、dblquad 求函数f(x,y)在矩形区域 [a,b]×[c,d]上的二重定积分 double quadrature 六、dblquad 格式 q=dblquad(fun, xmin,xmax,ymin,ymax, tol) fun: 被积函数，四种定义方式 xmin,xmax: 被积函数fun(x,y)的x积分区间 ymin,ymax: 被积函数fun(x,y)的y积分区间 tol: 精度要求，缺省值1.e-6 例8　计算二重积分 q=dblquad(fun, xmin,xmax,ymin,ymax) 例9　计算二重积分 q=dblquad(fun, xmin,xmax,ymin,ymax) 例10　计算重积分 七、integral2 double quadrature 求解一般区域（不限于矩形）二重积分 其中 ymin,ymax： y的积分上下限 表达为x的函数 匿名函数定义 q = integral2(fun, xmin, xmax, ymin,ymax) 八、quad2d double quadrature q = quad2d(fun, xmin, xmax, ymin,ymax) 求解一般区域（不限于矩形）二重积分 其中 ymin,ymax： y的积分上下限 表达为x的函数 匿名函数定义 例10　计算重积分 九、triplequad 求函数f(x,y,z)在长方体 [xmin,xmax] ×[ymin,ymax] ×[zmin,zmax]上三重定积分 triple quadrature 九、triplequad triple quadrature 格式 q=triplequad(fun,xmin,xmax,ymin, ymax, zmin,zmax) fun: 被积函数 四种方式 十、integral3 求函数f(x,y,z)在一般区域上的三重定积分 十、integral3 格式 Q=integral3(fun,xmin,xmax,ymin, ymax, zmin,zmax) ymin,ymax：表达为x的函数 zmin,zmax：表达为x,y的函数 匿名函数定义 例11 计算重积分 数值积分功能函数 函数名 功能特点 integral 可计算反常积分 quadgk Gauss求积公式，精度高，可计算反常积分 quad 变步长复化Simpson公式，可计算瑕积分 quadl Gauss求积公式，可计算瑕积分 trapz 复化梯形公式，列表函数求积分 数值积分功能函数 函数名 功能特点 integral2 二重积分，一般区域 quad2d 二重积分，一般区域 dblquad 二重积分，矩形域 integral3 三重积分，一般区域 triplequad 三重积分，积分域为长方体 * * * 实用 MATLAB 已知Lagrange插值算法的命令式m文件，求f(4.5)。改写为函数式m文件，再求f(4.5)。已知节点数据如下表： x -5 -4 -3 -2 -1 0 f(x) 0.038462 0.058824 0.1 0.2 0.5 1 x 1 2 3 4 5 ? f(x) 0.5 0.2 0.1 0.058824 0.038462 ? 作 业 x=[-5 -4 -3 -2 -1 0 1 2 3 4 5];%已知节点 y=[0.038461538