第3章 符号运算基础.pptVIP

第3章 符号运算基础 数值计算与符号计算是计算机计算科学中各有特点的两类计算方法。 由于计算机所有运算都基于二进制数据，所以数值计算以数值数组作为运算对象给出数值解，具有计算速度快和占用资源少的特点，但受到计算机的存储空间和运算时间的限制，使得任何数值计算都有可能产生累积误差，结果是一个被截断的近似值。 符号计算以符号对象和符号表达式作为运算对象，得到的解析表达式永远是精确的，但符号计算需占用大量的系统资源和消耗较多时间。 数值计算的常用运算单元是数值矩阵，必须事先赋值才能进行计算；参与符号运算的符号表达式，也必须事先对其中的基本对象进行符号定义，才能进行相应的符号计算。 3.1 建立符号对象 3.1.1建立符号量 1、函数sym用来建立单个符号量 一般调用格式是 符号量名=sym(符号字符串) 其中，符号字符串可以是常量、变量、矩阵、函数或表达式。 例如： f1=sym(a*x^2+b*x+c) % 二次三项式 f2=sym(Dy+y^2=1) % 微分方程 A=sym(‘[ ]’) % 符号矩阵 其中，符号矩阵的每一行的两端都有方括号，这是与数值矩阵的一个重要区别。 2、函数syms用来建立多个符号量。一般调用格式是： syms 变量1 变量2 ... 变量n 其中，各个变量之间只能使用空格分隔。 例如： syms a b c x % 定义4个符号变量 f=a*x^2+b*x+c % 组成符号表达式 运算结果： f = a*x^2+b*x+c 3.1.2 建立符号表达式 建立一个符号表达式有3种方法。例如，对于二次多项式，有： 1、利用单引号建立：f=a*x^2+b*x+c； 2、利用函数sym建立：f=sym(a*x^2+b*x+c)； 3、利用函数syms已经定义的符号变量建立： syms a b c x; f=a*x^2+b*x+c 以上3种方式建立的符号表达式结果是相同的。 3.1.3 符号表达式中变量的确定 1、在上面的符号表达式中只有1个变量时，系统会自动将x作为自变量，将a、b、c等作为常量处理； 2、如果符号表达式中有多个变量时，系统会自动将除了i和j之外，在字母位置上最接近x的小写字母作为自变量，例如：…u、v、w、x、y、z。 3、利用函数findsym可以查找符号表达式中变量。例如 syms a b c d x y % 定义6个符号变量 fun=a*x^2+b*x*y+c*y^2+d; % 定义符号表达式 findsym(fun) % 查找符号表达式中变量 运算结果： ans = a, b, c, d, x, y 3.2 符号表达式的基本运算 3.2.1 符号表达式的四则运算 利用运算符+、-、*、/、^实现符号表达式的四则运算，其结果是一个符号表达式。 例3-1 比较符号常数与数值在代数运算时差别。 k1=sym(8);k2=sym(2); % 定义符号变量并赋值 r1=8;r2=2; % 数值变量赋值 fun_k=exp(k1)+sqrt(k2) % 计算符号表达式值 fun_r=exp(r1)+sqrt(r2) % 计算数值表达式值 whos 运算结果： fun_k = exp(8)+2^(1/2) fun_r = 2.9824e+003 Name Size Bytes Class fun_k 1x1 152 sym object fun_r 1x1 8 double array k1 1x1 126 sym object k2 1x1 126 sym object r1 1x1 8 double array r2 1x1 8 double array Grand total is 22 elements using 428 bytes 可见，k1、k2和fun_k是符号量，r1、r2和fun_r是数值量。符号运算是进行数学演算，得到精确的解析表达式；数值计算得到近似的有限数值。 3.2.2 符号表达式分子和分母的提取 如果符号表达式是一个有理分式或是可以展开为有理分式，利用函数numden可以提取其中的分子和分母。其一般调用格式是： [n,d]=numden(s) 其中，s是符号表达有理分式；n是其分子；d是其分母。 例3-2 提取符号量和