实验二积分初步.ppt

实验二 MATLAB中的积分运算初步 一、实验目的 熟悉MATLAB软件中关于积分运算的基本命令，掌握利用MATLAB软件进行求不定积分、定积分等积分运算的方法。 二、相关知识 在高等数学中，我们曾经学习了求函数的不定积分和定积分的运算，那时我们根据微积分学的原理，学习了一整套各种各样的方法，其中包括了许多技巧，现在我们尝试用软件来解决这样的问题。 MATLAB提供的命令函数int()可以完成积分运算，其调用格式有如下几种： int(fun) 计算函数fun关于默认变量的不定积分 int(fun,x) 计算函数fun关于变量x的不定积分 int(fun,x, a, b) 计算函数fun关于变量x从a到b的定积分 我们通过例子来学习具体的用法： 例1：计算不定积分 。 解：可以用下面的程序完成： clear int(fun) 计算函数fun关于默认变量的不定积分 int(fun,x) 计算函数fun关于变量x的不定积分 int(fun,x, a, b) 计算函数fun关于变量x从a到b的定积分 我们通过例子来学习具体的用法： 例1：计算不定积分 。 解：可以用下面的程序完成： clear y=sym(‘x^5+x^3-sqrt(x)/4’) int(y) pretty(ans) 计算结果为：1/6*x^6+1/4*x^4-1/6*x^(3/2) 例2：计算定积分 。 解：可以用下面的程序实现计算： clear y=sym(‘x^5+x^3-sqrt(x)/4’) int(y) pretty(ans) 例2：计算定积分 。 解：可以用下面的程序实现计算： clear syms x y=(x*exp(x))/(1+x)^2; int(y,0,1) 计算结果为：1/2*exp(1)-1 MATLAB的函数int同样可以用于计算广义积分 例3 计算广义积分 解：可以用下面的程序实现计算： clear syms x y=(x*exp(x))/(1+x)^2; int(y,0,1) MATLAB的函数int同样可以用于计算广义积分 例3 计算广义积分 解：可以用下面的程序实现计算： clear syms x int(1/x/sqrt(x*x-1),x,-2,-1) 计算结果为：-1/3*pi。 例4 计算广义积分 解：可以用下面的程序实现计算： clear syms x syms x int(1/x/sqrt(x*x-1),x,-2,-1) 计算结果为：-1/3*pi。 例4 计算广义积分 解：可以用下面的程序实现计算： clear syms x int(1/(x*x+1)/(x*x+4),x,-inf,inf) 计算结果为：1/6*pi 并非所有的积分都可以用MATLAB软件通过函数int( )来完成运算，函数int( )完成的是符号运算，而不是数值运算，因此，当积分不存在初等形式的原函数时，函数int( )便不能完成积分运算。 例5计算定积分 。 解：先用语句： smys x int(x*sin(x)/(1+cos(x)*cos(x)),x,0,pi) 运行后，系统会返回： 并非所有的积分都可以用MATLAB软件通过函数int( )来完成运算，函数int( )完成的是符号运算，而不是数值运算，因此，当积分不存在初等形式的原函数时，函数int( )便不能完成积分运算。 例5 计算定积分 。 解：先用语句： smys x int(x*sin(x)/(1+cos(x)*cos(x)),x,0,pi) 运行后，系统会返回： ans = int(x*sin(x)/(1+cos(x)^2),x = 0 .. pi)， 这表示系统无法完成运算，为了解决这个问题，我们可以用另一个MATLAB函数quad( )，这个函数用于计算函数的数值积分，我们通过例子来说明quad( )的使用用法和步骤： 1．在磁盘上建立一个函数文件f.m，内容如下： function y=f(x) y=x.*sin(x)./(1+cos(x).^2) 2．在另一个文件或命令窗口中输入： quad(@f,0,pi) ans = int(x*sin(x)/(1+cos(x)^2),x = 0 .. pi)， 这表示系统无法完成运算，为了解决这个问题，我们可以用另一个MATLAB函数quad( )，这个函数用于计算函