第三章符号运算-201410202215298讲述.pptVIP

一、字符串与字符串函数 关于字符串的约定 （1）在MATLAB中所有的字符串都用单引号括起来。 （2）字符串的每个字符（包括空格）都是矩阵相应的一个元素。 （3）在MATLAB中，字符串的字符数组（或矩阵）基本上是等价的。 字符串和数组之间的转换： ? 将字符串转换为数值：double() 把符号矩阵或任意精度表示的矩阵A转换成为双精度矩阵。 符号表达式包括符号函数和符号方程，两者的区别是：前者不包括等号，而后者必须带有等号。但两者的创建方式是相同的。 符号运算的最终目的是为了求得数值解，这就要对解析解进行数值转换。 digits 函数设置返回有效数字个数 vpa 函数返回符号表达式在digits函数设置下的精度的数值解。 问： 请给出sin(pi/8)保留6位有效数字的结果。 3、subs(S,OLD,NEW) 将符号表达式中的OLD变量替换为NEW变量。 4、numeric(S) 将不含自由变量的符号表达式转换为数值形式，其效果同double(sym(S)) 四、符号函数的运算 1. 复合函数的运算 compose(f,g) 返回当f=f(x)和g=g(y)时的复合函数f(g(y))。 compose(f,g,z) 返回自变量为z 的复合函数 compose(f,g,x,z) 返回复合函数f(g(z))且使得x为f的独立变量，y为g的独立变量。 compose(f,g,x,y,z) 返回复合函数f(g(z))并使得x 为f的独立变量，y为g的独立变量。 例：syms x y z t u; f=1/(1+x^2); g=sin(y); h=x^t; p=exp(-y/u); compose(f,g) 五、 符号矩阵的创立 1．使用sym函数直接创建符号矩阵 矩阵元素可以是任何不带等号的符号表达式，各符号表达式的长度可以不同，矩阵元素之间 可用空格或逗号分隔。 例：a=sym([1/s+x,sin(x),cos(x)^2/(b+x);9,exp(x^2+y^2),log(tanh(y))]) 2．用创建子阵的方法创建符号矩阵 这种方法不需要调用sym命令，但要保证同一列的各字符串具有相同的长度。 3．将数值矩阵转化为符号矩阵 不论数值矩阵的元素原先是用分数还是用浮点数表示，转化后的符号矩阵将都以最接近的精确有理数形式给出。 a = 0.6667 1.4142 0.2220 1.4000 4.3478 1.0986 b=sym(a) b = [ 2/3, sqrt(2), 111/500] [ 7/5, 100/23, 709893870346*2^(-52)] 六、 符号矩阵的运算 1．基本运算 （1） 符号矩阵的四则运算：和数值矩阵的四则运算完全相同。+，-，*，\，/ （2） 符号矩阵的其他一些基本运算：和数值矩阵的运算格式相同，包括转置（’）、行列式（det）、 逆(inv)、秩(rank)、幂（^）和指数(exp 和expm)等运算。 3、符号矩阵的简化 （1）factor() 因式分解符号矩阵S的各个元素，如果S包含的所有元素为整数，则最佳因式分解式将被计算。 例：syms x factor(x^9-1) ans = (x-1)*(x^2+x+1)*(x^6+x^3+1) expand(sin(x+y)) ans = sin(x)*cos(y)+cos(x)*sin(y) 2. expand() 对符号矩阵的各个元素的符号表达式进行展开。 例：syms x y expand((x+1)^3) ans = x^3+3*x^2+3*x+1 符号表达式的化简 （3） collect() 函数实现功能为将符号表达式中同类项合并 R=collect(s ,v) 将表达式S 中v 的相同次幂的项进行合并。如果 v 没有指定，则缺省地将含有x 的相同次幂的项进行合并。 例子 collect(x^2*y + y*x - x^2 - 2*x) ans = (y-1)*x^2+(y-2)*x f = -1/4*x*exp(-2*x)+3/16*exp(-2*x); collect(f,exp(-2*x)) ans = (-1/4*x+3/16)*exp(-2*x) (4).simple() 对表达式S尝试多种不同的算法简化，以显示S表达式的长度最短