另一方面它们在实际问题中也有着许多直接的应用函数的积.docVIP

积分是非常重要的数学工具（微分方程、概率论等的基础），另一方面它们在实际问题中也有着许多直接的应用.函数的积分计算可分为数值积分和符号积分两类方法.数值积分是求积分的近似值的一类近似计算的方法.对于定积分d，如果对应的不定积分d不易求出，或者根本不能表示为初等函数时，那么就只能用数值积分的方法求其近似值了.例如，因为不定积分ed，d，ed，d无法由初等函数表示，所以计算这种类型的定积分只能用数值方法.至于由离散数据或图形表示的函数的定积分，理所当然地属于数值积分的范畴.符号积分是指求积分的精确值.本章重点介绍数值积分及其用MATLAB计算的方法.为了将数值积分的值与精确值比较，在第一节中简单介绍用MATLAB 作符号积分的方法. 9.1 积分的MATLAB符号计算 计算机软件MATLAB的系统提供了符号计算定积分的函数int.下面主要介绍用MATLAB软件进行符号计算定积分、变上（下）限积分及其应用. 9.1.1 定积分的MATLAB符号计算 如果函数在区间上连续，且是的一个原函数，则函数在区间上的定积分为d.我们可以用MATLAB软件中的函数int对定积分进行符号计算，它的调用格式和功能如表9–1所示. 表9–1 函数int求定积分的调用格式和功能 函数int求定积分的调用格式 功 能 F=int(S,a,b) 输入量：S可以是被积函数，也可以是由几个被积函数，构成的矩阵，积分上、下限a 和 b可以是数，也可以是符号标量，积分变量是. 输出量：F是S关于积分变量的从a 到 b的定积分. F=int(S,t,a,b) 输入量：S、a 和 b同上，是积分变量. 输出量：F是S关于积分变量从a 到 b的定积分. 例9.1.5 由，所围成的平面区域D.求平面区域D的面积S. 解 （1） 画出由，所围成的封闭平面区域的图形，输入作函数图形的程序 x=-1:0.001:2; F1= sin(x); F2=cos(x); plot(x ,F1,b-,x ,F2,g-), axis([-1,pi/4+1,-1.3,1.3]), xlabel(x), ylabel(y), title(y=sinx , y=cosx 和x=-0.5及x=1.5所围成的平面区域的图形) 运行后屏幕显示图9–1. 图9–1 由，所围成的平面区域图形 （2）求平面区域D的面积S. 解方程组得在内的两条曲线交点的横坐标，所求面积为 dd. 输入计算面积S的程序 syms x f1= cos(x)-sin(x); f2=-f1; S1 =int(f1,x,-0.5,pi/4); S2=int(f2,x, pi/4,1.5); S=S1+S2,Sj= double (S) 运行后屏幕显示计算面积的值 S 及其近似值 Sj 如下 S = 2*2^(1/2)+sin(1/2)-cos(1/2)-sin(3/2)-cos(3/2) Sj = 1.36203791318826 9.1.2 变限积分的MATLAB符号计算 例9.1.7 已知ed，求. 解 因为 edteded 所以输入程序 syms x t F1=int(exp(t)*sin(2+sqrt(t^3)),x,0); F2=int(exp(t)*sin(2+sqrt(t^3)),0,x^2); Fi= F1+ F2; dF=diff(Fi) 运行后屏幕显示计算变限积分的导数如下 dF = -exp(x)*sin(2+(x^3)^(1/2))+2*x*exp(x^2)*sin(2+(x^6)^(1/2)) 9.2 数值积分的思想及其MATLAB程序 9.2.3 矩形公式的MATLAB程序 用矩形公式计算数值积分可用下面函数sum和cumsum实现，调用格式如下： （一） 函数sum的调用格式 函数sum的两种调用格式如下： 调用格式一：sum(X) 如果输入向量X，则输出向量X的所有元素的和，可用于按矩形公式（9.3）和（9.4）计算积分；如果输入矩阵X，则输出矩阵X的每列向量的和. 调用格式二：sum (X,DIM) 输入N-D 数组X，输出为X的每个DIM向量的和. 例9.2.2 用MATLAB和矩形公式（9.3）、（9.4）计算ed，并与精确值比较. 解 将分成20等分，步长为，输入程序 h=pi/40; x=0:h:pi/2; y=exp(sin(x)); z1=sum(y(1:20))*h, z2=sum(y(2:21))*h, 运行后屏幕显示矩形公式计算结果分别如下 z1 = z2 = 3.0364 3.1713 即矩形公式（9.3）计算结果3.036 4,