matlab分子分母同时扩大或缩小,MATLAB语言 微积分

解 :format long %将MATLAB 的结果精度扩大

quad (‘1./(1+x.^2)’，0,1) (抛物线法求数值积分) ans=0.78539814924326 例10 判别广义积分?

+∞1

1dx x p

、?+∞∞--πdx e x 22

21与?-2

02)1(1dx x 的敛散性,收敛时计算积分值。

解 对第一个积分输入命令

syms p real;int(1/x^p,x,1,inf)

得结果ans =limit(-1/(p-1)*x^(-p+1)+1/(p-1),x = inf).由结果看出当1

p 时,ans=1/(p-1),这与课本例题是一致的。

对第二个积分输入命令

int(1/(2*pi)^(1/2)*exp(-x^2/2),-inf,inf) 得结果:

ans=7186705221432913/18014398509481984*2^(1/2)*pi^(1/2) 由输出结果看出这两个积分收敛。

对后一个积分输入命令 int(1/(1-x)^2,0,2)

结果得ans=inf .说明这个积分是无穷大不收敛。 4、MATLAB 在微分方程中的应用

在科学技术和现代经济管理中有些实际问题，常常需要通过未知函数及其导数所满足

的关系式去求未知函数，这种关系式就叫做微分方程。其理论产生于17世纪末，并在力学、天文学、生物学以及现代经济管理中已得到较为广泛的应用，并且取得了巨大成就。微分方程是高等数学的重要组成部分。而解微分方程是一个繁杂且困难的工作，MATLAB 提供了非常强大的求解微分方程的功能。求解微分方程的调用格式如下：

函数 dsolve

格式 r = dsolve('eq1,eq2,…','cond1,cond2,…','v')

说明 对给定的常微分方程(组)eq1,eq 2,…中指定的符号自变量v ，与给定的边界条件和初始条件cond1,cond2,….求符号解(即解析解)r ；若没有指定变量v ，则缺省变量为t ；在微分方程(组)的表达式eq 中，大写字母D 表示对自变量(设为x)的微分算子：D=d/dx ，D2=d2/dx2，…。微分算子D 后面的字母则表示为因变量，即待求解的未知函数。初始和边界条件由字符串表示：y(a)=b ，Dy(c)=d ，D2y(e)=f ，等等，分别表示b )x (y a x ==，

d )x (y c x ='=，f )x (y

e x =''=；若边界条件少于方程(组)的阶数，则返回的结果r 中会出现任意常数C1，C2，…；dsolve 命令最多可以接受12个输入参量(包括方程组与定解条件个数，当然我们可以做到输入的方程个数多于12个，只要将多个方程置于一字符串内即可)。若没有给定输出参量，则在命令窗口显示解列表。。但是，我们知道，有大量的常微分方程虽然从理论上讲，其解是存在的，但我们却无法求出其解析解，此时，我们需要寻求方程的数值解，在求常微分方程数值解方面，MATLAB 具有丰富的函数，我们将其统称为solver ，其一般格式为：