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

MATLAB语言与应用 第 5 章微积分问题的计算机求解 现代设计与分析研究所 王 雷 主要内容 微积分问题的解析解 数值积分 5.1 微积分问题的解析解 5.1.1 极限问题的解析解 5.1.2 函数导数的解析解 5.1.3 积分问题的解析解 5.1.1 极限问题的解析解5.1.1.1单变量函数的极限 【例5-1】试求解极限问题 【例5-2】求解单边极限问题 5.1.1.2 多变量函数的极限 【例5-3】求出二元函数极限值 5.1.2 函数导数的解析解5.1.2.1 函数的导数和高阶导数 【例5-4】  5.1.2.2 多元函数的偏导 【例5-5】 【例5-6】 5.1.2.3 隐函数的偏导数 【例5-7】 5.1.2.4 参数方程的导数 已知参数方程 ,求 5.1.3 积分问题的解析解5.1.3.1 不定积分的推导 【例5-9】  用diff() 函数求其一阶导数，再积分，检验是否可以得出一致的结果。  对原函数求4 阶导数，再对结果进行4 次积分 【例5-10】证明 【例5-11】不可积问题 的积分问题求解。  5.1.3.2 定积分与无穷积分计算 【例5-13】 【例5-14】  【例5-15】 5.2 数值积分问题 5.2.1 由给定数据进行梯形求积 5.2.2 单变量数值积分问题求解 5.2.3 双重积分问题的数值解 5.2.4 三重定积分的数值求解 5.2.1 由给定数据进行梯形求积 【例5-16】 【例5-17】 5.2.2 单变量数值积分问题求解 【例5-18】 用inline函数定义: 【例5-19】高精度数值解（Lobbato算法） 【例5-20】求解 5.2.3 双重积分问题的数值解 【例5-23】求解 【例5-24】 解析解方法： 5.2.4 三重定积分的数值求解 【例5-25】 本章要点简介 本章涉及的函数小结 x1=[0:pi/30:pi]’ y=[sin(x1) cos(x1) sin(x1/2)]; 画图： 求理论值： 不同步距 第三种：匿名函数(MATLAB 7.0) 第二种：inline 函数 第一种: 一般函数方法 提高求解精度。 y1=quadl(f, 0, 1.5, 1e-16) abs([y1-y0, y-y0]) 绘制函数: 调用quadl（）： f =inline(‘exp(x.^2).*(x=2) +’. . . ‘80*(x2) ./(4-sin(16*pi*x) )’ , ‘ x’ ); I1=quadl(f, 0, 4) 【例5-21】 求两条曲线 y1=x.^2 和 y2=x+2 所围成的二维图形的面积和重心 设两条曲线的交点为（x1, y1), (x2, y2)，则 绘制函数: 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 求解： 【例5-22】 求螺旋线的长度。已知螺旋线参数 方程： 绘制函数: 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) 求解： Nit工具箱 * *  函数及其一阶导数图形： f4=diff(f, x, 4); pretty(f4) 求 和 三维曲面： 引力线： syms x y z; f=sin(x^2*y)*exp(-x^2*y-z^2); df=diff (diff (diff (f, x, 2), y), z); df=simple(df) 【例5-8】 syms t; y=sin(t)/(t+1)^3; x=cos(t)/(t+1)^3; d4=diff(y,t,4)/diff(x,t,4) syms x; y=sin(x)/(x^2+4*x+3); y1=diff(y); y0=int(y1); pretty(y0) y4=diff(y,4); y0=int( int( int( int( y4) ) ) ); simplify