本试验取材于中南大学《科学计算与MATLAB语言》
- 符号函数的极限
- 符号函数的导数
- 符号函数的积分
符号函数的极限
- 求符号函数极限的命令为limit,其调用格式为:
limit(f,x,a)
即求函数f关于变量x在a点的极限,若x省略,则采用系统默认的自变量,a的默认值为0.
- limit函数的另一种功能是单边极限,其调用格式为:
limit(f,x,a'right')
limit(f,x,a,'left')
>> syms a m x n;
>> f=(x^(1/m)-a^(1/m))/(x-a);
>> limit(f,x,a)
ans =
a^(1/m)/a/m
>> g=(1+1/n)^n;
>> limit(g,n,inf)
ans =
exp(1)
>>
符号函数的导数
MATLAB中的求导函数为:
diff(f,x,n)
即求函数f关于变量x的n阶导数。参数x的用法同求极限函数limit,可以缺省,默认值与limit相同,n的默认值是1.
求下列函数的导数。
第一题:
>> syms x y;
>> f=sqrt(1+exp(x));
>> diff(f)
ans =
1/2/(1+exp(x))^(1/2)*exp(x)
第二题
>> g=x*exp(y)/y^2;
>> diff(g,x)
ans =
exp(y)/y^2
>> diff(g,y)
ans =
x*exp(y)/y^2-2*x*exp(y)/y^3
>>
极限、导数、微分的概念是紧密关联的。有极限是可导的前提,而导数是微分之商,因此导数也称为微商。
符号函数的积分
不定积分
在MATLAB中,求不定积分的函数是int,其常用的调用格式为:
int(f,x)
求函数f对变量x的不定积分。
求下列不定积分
第一题
>> syms x t;
>> f=(3-x^2)^3;
>> int(f)
ans =
27*x-1/7*x^7+9/5*x^5-9*x^3
第二题
>> g=5*x*t/(1+x^2);
>> int(g,t)
ans =
5/2*x*t^2/(1+x^2)
>>
定积分
在MATLAB中,定积分的计算也使用int命令,但调用格式有区别:
int(f,x,a,b)
其中,a,b分别表示定积分的下限和上限.
- 当函数f关于变量x在闭区间[a,b]可积时,函数返回一个定积分结果。
- 当a、b中有一个是inf时,函数返回一个广义积分.
- 当a、b中有一个符号表达式时,函数返回一个符号函数.
求下列定积分
>> syms x t;
>> int(abs(1-x),1,2)
ans =
1/2
>> int(1/(1+x^2),-inf,inf)
ans =
pi
>> int(4*x/t,t,2,sin(x))
ans =
4*log(sin(x))*x-4*log(2)*x
>>