第4章MATLAB符号计算论述.ppt

第四章 MATLAB符号计算 4.1符号计算基础 4.2微分运算 4.3积分运算 4.4级数 4.5求解方程 4.6积分变换 4.1符号计算基础 符号对象的创建 sym函数：用来创建单个符号变量 syms函数：用来创建多个符号变量 格式: a=sym(‘a’); syms a b c ：符号变量a和非符号变量a是不同的,非符号变量运算前必须赋值,变量的运算实际上是对应值的运算.符号变量运算前不需要赋值,结果 为一个由参与运算的变量名组成的表达式. 例: a=sym(a);b=sym(b);c=sym(c);d=sym(d); A=[a,b;c,d]; A = [ a, b] [ c, d] det(A) ans = a*d-b*c 注意: sym还可定义符号常数,但与数值常数有些差别 a=sym(‘8’); b=8; sqrt(a); ans = 2*2^(1/2) sqrt(b); ans = 2.8284 创建符号表达式 方式： 1、直接使用sym函数创建。 2 、使用已经定义的符号变量来创建。 例: f=sym(‘a*x^2+b*x+c’); f = a*x^2+b*x+c syms a b c x f= a*x^2+b*x+c 创建符号矩阵 与创建符号表达式第二种方法相近 syms a11 a12 a13 a21 a22 a23 a31 a32 a33 A=[a11 a12 a13;a21 a22 a23;a31 a32 a33] A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] 基本符号运算 factor(f) 对符号表达式或符号矩阵进行因式分解 expand(f) 对符号表达式或符号矩阵进行展开 collect(f) 对符号表达式或符号矩阵合并同类项 collect(f,v) 对符号表达式或符号矩阵按变量v进行合并同类项 例 syms x,y; f=(2*x+3*y)*(x+1); a=expand(f) a = 2*x^2+2*x+3*y*x+3*y a=collect(a,x) a = 2*x^2+(2+3*y)*x+3*y factor(a) ans = (2*x+3*y)*(x+1) 4.2微分运算 符号运算的极限 格式：limit(f,x,a) 说明：参数x,a可以省略,当x省略时按默认原则确定极限变量;当a省略时,matlab默认变量趋近于0. 例： 求极限 syms x; limit(sin(x)/x,x,0); ans = 1 limit函数求单边极限 格式： limit(f,x,a,’right’); limit(f,x,a,’lift’); 例：求 syms x;limit(1/x,x,0) ans= NaN limit(1/x,x,0,’left’) ans= - inf limit(1/x,x,0,’right’) ans= inf 符号函数的微分 格式: diff(f,x,n) syms x; diff(x*cos(x),x) ans= cos(x)-x*sin(x) diff(x*cos(x),x,2) ans= -2*sin(x)-x*cos(x) 4.3积分运算 符号函数的不定积分 格式: int(f,x) 例如:求不定积分 syms x int(1/(x^2+1),x) ans = atan(x) syms x y f=int(x*y,x); f = 1/2*x^2*y int(f,y) ans = 1/4*x^2*y^2 有时积分过于复杂，其积分未必存在，或存在但计算机无法找到他，MATLAB会简单的返回命令窗口，包含积分变量但未进行任何计算。当积分含有其他参变量时，当常数处理。 符号函数的定积分 int(f,x,a,b) 例: syms x t f1=x^2/(x+1) ; a1=int(f1,x,1,2) a1 = log(3)+1/2-log(2) double(a1) ans = 0.9055 f2=1/(x^2+1) ; a2=int(f2,x,-inf,inf) a2= pi f3=2*t*x; a3=int(f3,x,1,sin(t)); a3 = t*(sin(t)^2-1) 4.4级