第 5 章 微积问题的计算机求解.ppt

* MATLAB语言与应用 * 【例5-15】 * MATLAB语言与应用 * 5.2 数值积分问题 5.2.1 由给定数据进行梯形求积 5.2.2 单变量数值积分问题求解 5.2.3 双重积分问题的数值解 5.2.4 三重定积分的数值求解 * MATLAB语言与应用 * 5.2.1 由给定数据进行梯形求积 * MATLAB语言与应用 * * MATLAB语言与应用 * 【例5-16】 x1=[0:pi/30:pi]’ y=[sin(x1) cos(x1) sin(x1/2)]; * MATLAB语言与应用 * 【例5-17】 画图： * MATLAB语言与应用 * 求理论值： 不同步距 * MATLAB语言与应用 * 5.2.2 单变量数值积分问题求解 * MATLAB语言与应用 * 【例5-18】 第三种：匿名函数(MATLAB 7.0) 第二种：inline 函数 第一种: 一般函数方法 * MATLAB语言与应用 * 用inline函数定义: * MATLAB语言与应用 * 【例5-19】 提高求解精度。 y1=quadl(f, 0, 1.5, 1e-16) abs([y1-y0, y-y0]) * MATLAB语言与应用 * 【例5-20】求解 绘制函数: * MATLAB语言与应用 * 调用quadl（）： f =inline(‘exp(x.^2).*(x=2) +’. . . ‘80*(x2) ./(4-sin(16*pi*x) )’ , ‘ x’ ); I1=quadl(f, 0, 4) * MATLAB语言与应用 * 【例5-21】 求两条曲线 y1=x.^2 和 y2=x+2 所围成的二维图形的面积和重心 设两条曲线的交点为（x1, y1), (x2, y2)，则 * MATLAB语言与应用 * 绘制函数: x=linspace(-1,2,200); y1=x.^2; y2=x+2; fill(x,y1,r) area=quadl(inline(‘x+2’,‘x’),-1,2)-,… quadl(inline(x.^2,x),-1,2) xc=quadl(inline(x.*((x+2)-x.^2),x),-1,2)/area yc=quadl(inline(((x+2).^2-x.^4)/2,x),-1,2)/area 求解： * MATLAB语言与应用 * 【例5-22】 求螺旋线的长度。已知螺旋线参数 方程： * MATLAB语言与应用 * 绘制函数: t=linspace(0,6*pi); plot3(sin(t),cos(t),t),grid on length=quadl(inline(sqrt(cos(t).^2+4*,... sin(2*t).^2+1)),0,6*pi) 求解： * MATLAB语言与应用 * 5.2.3 双重积分问题的数值解 * MATLAB语言与应用 * 【例5-23】求解 * MATLAB语言与应用 * Nit工具箱 * MATLAB语言与应用 * 【例5-24】 * MATLAB语言与应用 * 解析解方法： * MATLAB语言与应用 * 5.2.4 三重定积分的数值求解 * MATLAB语言与应用 * 【例5-25】 triplequad(inline(. . . ‘4*x.*z.*exp(-x.*x.*y-z.*z)’, . . . ‘x’,’y’,’z’), 0, 2, 0, pi, 0, pi, 1e-7) * MATLAB语言与应用 * 本章要点简介 本章涉及的函数小结 limit() 极限问题求解 diff（） 导数问题求解 int（） 积分问题求解 trapz（） 梯形法求数值积分 quadl（） 数值积分函数 dblquad（） 矩形区域的二重数值积分 triplequad（）长方形区域的三重数值积分 * MATLAB语言与应用 MATLAB语言与应用 * MATLAB语言与应用 * 第 5 章微积分问题的计算机求解 * MATLAB语言与应用 * 主要内容 微积分问题的解析解 数值积分 * MATLAB语言与应用 * 5.1 微积分问题的解析解 5.1.1 极限问题的解析解 5.1.2 函数导数的解析解 5.1.3 积分问题的解析解 * MATLAB语言与应用 * 5.1.1 极限问题的解析解5.1.1.1单变量函数的极限 * MATLAB语言与应用 * 【例5-1】试求解极限问题 * MATLAB语言与应用 * 【例5-2】求解单边极限问题 * MATLAB语言与应用 * 5.1.1.2 多变量函数的极限 * MATLAB语言