7MATLAB符号计算全解.ppt

第七章 MATLAB符号计算 符号计算基础 符号导数及其应用 符号积分 级数 符号方程求解 所谓符号计算是指在运算时,无须事先对变量赋值,而将所得到结果以标准的符号形式来表示。 MathWorks公司以Maple的内核作为符号计算引擎（Engine），依赖Maple已有的函数库，开发了实现符号计算的两个工具箱：基本符号工具箱和扩展符号工具箱。 7.1 符号计算基础 7.1.1 符号对象 1. 建立符号变量和符号常数 (1)sym函数 sym函数用来建立单个符号量，例如，a=sym(‘a’)建立符号变量a，此后，用户可以在表达式中使用变量a进行各种运算。 符号变量a和在其他过程中建立的非符号变量a是不同的.一个非符号变量在参与运算前必须赋值,变量的运算实际上是该变量所对应值的计算,其运算结果是一个和变量类型对应在的值,而符号变量参与运算前无须赋值,其结果是一个由参与运算的变量名组成的表达式. 另外,与数值运算一样,也可以用+、-、*、/、^运算符实现符号运算。 syms x y z; f=2*x+x^2*x-5*x+x^3+exp(2) f=2*x/(5*x) f=(x*x-y*y)/(x-y) 有时候，MATLAB并未将结果化为最简形式。例如，上面第三个符号表达式的结果不是x+y numden函数在提取各部分之前，将符号表达式有理化后，返回所得的分子和分母。 g=sym(‘(x^2+3)/(2*x-1)+3*x/(x-1)’) [n,d]=numden(g) 如果符号表达式是一个符号数组，numden返回两个新数组n和d，其中n是分子数组，d是分母数组 h=sym(‘[3/2,(2*x+1)/3;a/x+a/y,3*x+4]’) [n,d]=numden(h) 7.2 符号函数及其应用 7.2.1 函数的极限 limit函数的调用格式为： limit(f,x,a)：计算当变量x趋近于常数a时，f(x)函数的极限值。 limit(f,a)：默认变量趋于a时，f(x)函数的极限值。 limit(f)：系统默认变量趋近于0时，f(x)函数的极限值。 limit函数的另一种功能是求单边极限，其调用格式为： limit(f,x,a,right) 或 limit(f,x,a,left) 7.2.2 符号函数求导及其应用 MATLAB中的求导的函数为： diff(f,x,n) diff函数求函数f对变量x的n阶导数。参数x的用法同求极限函数limit，可以缺省，缺省值与limit相同，n的缺省值是1。 例7.3 求函数的导数。 命令如下： syms a b t x y z; f=sqrt(1+exp(x)); diff(f) %求(1)。未指定求导变量和阶数，按缺省规则处理 7.3 符号积分 7.3.1 不定积分 在MATLAB中，求不定积分的函数是int，其调用格式为： int(f,x) int函数求函数f对变量x的不定积分。参数x可以缺省，缺省原则与diff函数相同。 例7.5 求不定积分。 命令如下： x=sym(x); f=(3-x^2)^3; int(f) %求不定积分(1) f=sqrt(x^3+x^4); int(f) %求不定积分(2) g=simple(ans) %调用simple函数对结果化简 f=1/(1+x^2); int(f,-inf,inf) %求定积分(2) int(4*t*x,x,2,sin(t)) %求定积分(3) f=x^3/(x-1)^100; I=int(f,2,3) %用符号积分的方法求定积分(4) double(I) %将上述符号结果转换为数值 7.3.3 积分变换 1. 傅立叶(Fourier)变换 在MATLAB中，进行傅立叶变换的函数是： fourier(fx,x,t) 求函数f(x)的傅立叶像函数F(t)。 ifourier(Fw,t,x) 求傅立叶像函数F(t)的原函数f(x)。 2. 拉普拉斯(Laplace)变换 在MATLAB中，进行拉普拉斯变换的函数是： laplace(fx,x,t) 求函数f(x)的拉普拉斯像函数F(t)。 ilaplace(Fw,t,x) 求拉普拉斯像函数F(t)的原函数f(x)。 例7.10 计算y=x2的拉普拉斯变换及其逆变换. 命令如下： x=sym(x);y=x^2; Ft=laplace(y,x,t)