matlab积分的计算及其简单应用论文.docVIP

积分的计算及其简单应用 摘要：本文简要的概述了MATLAB在高等数学中积分的计算及应用；利用MATLAB中符号积分和数值积分的命令，计算定积分和不定积分。同时，也可以通过这些命令来解决一些实际问题，如求解面积，弧长，体积等等一些问题。 关键词：MATLAB，定积分，不定积分，面积，弧长。 一、概述 计算机的广泛应用和迅速普及，促成了它在数学方面的发展，例如，MATLAB强大的计算功能为以实验的方式学习和研究数学理论创造了良好的条件。在高等数学的教学中,严格的推理证明,容易使学生陷入僵化,没有感性认识。利用计算机进行仿真计算是一种行之有效的教学辅助方法。在此,应用MATLAB数学软件,以定积分，不定积分为例进行数学实验设计,通过课堂教学实践证明,对帮助学生理解定积分，不定积分概念和几何意义,具有良好的效果。 定积分，不定积分是高等数学中最基本的运算之一，它的运算思路很简单，就是求一条曲线，一个空间曲面及空间曲体在一定坐标系下对应的面积或体积。在MATLAB中，它的功能十分强大。以下是用MATLAB来解决一些定积分与不定积分的问题。 二、应用实例 1、用MATLAB计算不定积分 在MATLAB中，int函数用于求符号函数的不定积分，其调用格式为：int（f,v）,而对于定积分，只要加上其上限和下限就可以了，如：int(f,v,a,b).其中，f为被积函数，v为自变量，a和b分别为定积分的下限和上限。例如以下几个例子。 例1：用MATLAB求f= 的不定积分。 程序运行如下： syms x; %定义符号变量x f=(x+3)/(x^2-5*x+6) ;%给出被积函数f int(f) %求定积分 运行结果： ans = -5*log(x-2)+6*log(x-3) 2、符号积分和数值积分的计算 例2：计算f=在[0,1]上的定积分ans = 1/2*erf(1)*pi^(1/2) 当积分限由某一具体数值变为正负无穷时，定积分就转变为广义积分，也只需将积分限变为无穷，就可得到相应函数的广义积分值。如： 例3：求函数f=,g=在负无穷到正无穷的积分。 程序运行如下： syms x;%定义符号变量x f=1/(x^2+3*x+4);%被积函数f g=1/(x^2+3*x-4);%被积函数g subplot(2,1,1);%选择1*2个区中的1号区 ezplot(f,-10,10);%画出函数f的图形 subplot(2,1,2);%选择1*2个区中的2号区 ezplot(g,-10,10);%画出函数g的图形 f1=int(f,-inf,inf)%求函数f的积分 g1=int(g,-inf,inf)%求函数g的积分 结果如下： f1 = 2/7*pi*7^(1/2) ??? Error using == sym/maple Error, (in int/definite) wrong number (or type) of parameters in function has Error in == D:\matlab\MATLAB\toolbox\symbolic\@sym\int.m On line 51 == r = maple(map,int,f,[x.s = a.s .. b.s]); 图1 函数f和g的图像 有图1可以得出，函数f的积分结果为曲线f和x轴的面积，函数g的积分为无穷大，由图1的曲线也可以看出。 在求解定积分，不定积分时，还可以用matlab的数值积分命令：quad,quadl, trapz, cumtrapz, 等。其中，trapz(x,y)为梯形法沿列方向求函数Y关于自变量X的积分cumtrapz(x,y)为梯形法沿列方向求函数Y关于自变量X的累计积分f=cos(x)*cos(2*x)在[, pi/2]上的定积分f1 = 2/3 ans = 0.6667 解法二： x=-pi/2:pi/1000:pi/2;%数值计算适合于有限区间上,取有限个采样点 f=cos(x).*cos(2*x);%给出被积函数f trapz(x,f) %求定积分 结果如下： ans = 0.6667 解法三： syms x%定义符号变量x f=inline(cos(x).*cos(2*x));%定义一个语句函数f f1=quad(f,-pi/2, pi/2) %