第4章 符号运算概要1.ppt

* 也与通常的算术表达式一样， * 以找到对应的最简形式。 1. 用sym命令直接创建符号矩阵 这时sym命令的使用方法与前面创建符号表达式及方程的用法类似。所创建的符号矩阵的元素可以是任何符号对象，且元素的长度允许不同。在输入格式上，矩阵行之间以“；”分割，各矩阵元素之间用“，”或空格分隔。 【例4-15】 用sym函数创建符号矩阵。 A=sym(‘[a,b;c,d]) A = [ a, b] [ c, d] B=sym(‘[x+3*y,5*z+6*x;y-x,z/y]’) B = [ x+3*y, 5*z+6*x] [ y-x, z/y] 2. 以数值矩阵生成方法创建符号矩阵（即用syms函数创建） 用这种方法创建符号矩阵之前，需要预先定义所有需要的符号变量。 【例4-16】用生成数值矩阵的方法创建符号矩阵 syms x y z B=[x+3*x,5*z+6*z;y-y,z/z] B = [ 4*x, 11*z] [ 0, 1] 3. 由数值矩阵转换为符号矩阵 由于数值型对象和符号型对象分属于两个不同的数据类型，它们之间不能直接运算，但却可以相互转换。将数值对象M转化为符号对象S时，可以应用sym函数，格式为： S=sym(M)。 【例4-17】使用sym函数将3阶Hilbert矩阵（对称正定矩阵）转换为符号矩阵。 h=hilb(3) h = 1.0000 0.5000 0.3333 0.5000 0.3333 0.2500 0.3333 0.2500 0.2000 h1=sym(h) h1 = [ 1, 1/2, 1/3] [ 1/2, 1/3, 1/4] [ 1/3, 1/4, 1/5] 从本例可以看出，不管原来数值矩阵M是以分数还是浮点数形式赋值的，但当它被转化为符号矩阵后，都将以最接近原数的精确有理式给出。 4.4.2 符号矩阵的运算 符号矩阵的基本代数运算包括矩阵的四则运算、乘方、转置等，这些运算与数值矩阵的运算相同，这里仅举一例介绍其用法。 【例4-18】符号矩阵的基本运算。 syms t R=[cos(t),sin(t);-sin(t),cos(t)] R = [ cos(t), sin(t)] [ -sin(t), cos(t)] R.‘ %符号矩阵R的转置 ans = [ cos(t), -sin(t)] [ sin(t), cos(t)] D=det(R) %求矩阵R的行列式值 D = cos(t)^2+sin(t)^2 simplify(D) ans = 1 A=R.*R A = [ cos(t)^2+sin(t)^2, 0] [ 0, cos(t)^2+sin(t)^2] simplify(A) ans = [ 1, 0] [ 0, 1] 4.5 符号微积分 微积分是高等数学的基础，MATLAB的符号数学工具箱提供了许多关于微积分计算的功能。 4.5.1符号极限 limit函数用来求符号函数的极限。其格式如下： limit(F,x,a) 计算符号表达式F在x→a条件下的极限； limit(F,a) 计算符号表达式F中由默认自变量趋向于a条件下的极限； limit(F) 计算符号表达式F在默认自变量趋向于0条件下的极限； limit(F,x,a,‘right’) 和limit(F,x,a,’left’) 计算符号表达式F在x→a条件下的右极限和左极限。 【例4-19】计算如下表达式 、 、 、 分别输入下列语句： syms x a limit(sin(x)/x) limit(1/x,x,0,right) limit(1/x,x,0,left) v=[(1+a/x)^x, exp(-x)]; %v是由两个符号表达式为元素的符号矩阵 limit(v,x,inf,‘left’) %inf为无穷大 返回结果依次为： ans = 1 ans = inf ans = -inf ans = [ exp(a), 0] 4.5.2符号微分 diff函数用来求符号微分，其格式如下： diff(S)，求符号表达式S对于默认自变量的微分； diff(S,‘v’)，求符号表达式S对于自变量v的微分； diff(S,n)，求符号表达式S对于默认自变量的n次微分； diff(S, ‘v’, n)，求符号表达式S对自变量v的n次微分。 【例4-20】 符号表达式的微分运算。 S1 = sym(6*x^3-4*x^2+b