第3章节 MATLAB符号计算3.pptVIP

第3部分 MATLAB符号计算 2011-8-31 本部分内容 3.1 符号计算基础 3.2 微分运算 3.3 积分运算 3.4 泰勒级数 3.5 积分变换 3.6 方程求解 3.7 符号函数画图 3.1 符号计算基础 3.1.1 符号对象 符号表达式是代表数字、函数、算子和变量的MATLAB字符串，或者字符串组； 符号运算是使用已知的规则和给定符号恒等式求解符号方程的实践； 符号变量和非符号变量区别：非符号变量参与计算之间必须赋值；而符号变量无须赋值。 1 建立符号变量和符号常数 (1) sym函数 S=sym(A) sym函数用来建立单个符号变量，例如，a=sym(‘a’)建立符号变量a，此后，用户可以在表达式中使用变量a进行各种运算。 如果A是一个数值、数值矩阵或数值表达式，可不加‘’号，其输出是将数值对象转化为符号对象 如 pi=sym(2*pi)；a=sym(3/4) 例3.1考察符号变量和数值变量的差别。 在 MATLAB命令窗口，输入命令： a=sym(a);b=sym(b);c=sym(c);d=sym(d); %定义4个符号变量 w=10;x=5;y=-8;z=11; %定义4个数值变量 A=[a,b;c,d] %建立符号矩阵A B=[w,x;y,z] %建立数值矩阵B det(A) %计算符号矩阵A的行列式 det(B) %计算数值矩阵B的行列式 例3.2 比较符号常数与数值变量在代数运算时的差别。 在 MATLAB命令窗口，输入命令： pi1=sym(‘pi’);k1=sym(‘8’);k2=sym(‘2’);k3=sym(‘3’); % 定义符号常数 pi2=pi;r1=8;r2=2;r3=3; % 定义数值变量 pi1/3 % 计算符号表达式值 pi2/3 % 计算数值表达式值 sqrt(k1) % 计算符号表达式值 sqrt(r1) % 计算数值表达式值 sqrt(k3+sqrt(k2)) % 计算符号表达式值 sqrt(r3+sqrt(r2)) % 计算数值表达式值 (2)syms函数 syms用于定义多个符号变量。 syms函数的一般调用格式为： syms var1 var2 … varn 函数定义符号变量var1,var2,…,varn等。用这种格式定义符号变量时不要在变量名上加字符分界符(‘)，变量间用空格而不要用逗号分隔。 (3)符号常量 数值常量如果作为函数命令sym()的输入参量，即建立了一个符号对象---符号常量 可以用class()函数来检测数据的类似。 2. 建立符号表达式 例3.3 用2种方法建立符号表达式。 在MATLAB窗口，输入命令： U=sym(3*x^2+5*y+2*x*y+6) %定义符号表达式U syms x y; %建立符号变量x、y V=3*x^2+5*y+2*x*y+6 %定义符号表达式V 2*U-V+6 %求符号表达式的值 例3.4 计算3阶范得蒙矩阵行列式的值。设A是一个由符号变量a,b,c确定的范得蒙矩阵。 命令如下： syms a b c; U=[a,b,c]; A=[[1,1,1];U;U.^2] %建立范得蒙符号矩阵 det(A) %计算A的行列式值 (1)符号表达式的代数运算 与矩阵运算基本一致 M=sym（‘a+b’）;N=sym(‘b+c’)；M+N (2)符号矩阵的基本运算命令 transpose(A) 符号矩阵的转置 determ(A) 符号矩阵的行列式计算 inv(A) 符号矩阵的求逆 rank(A) 符号矩阵的求秩 eig 符号矩阵的特征值、特征向量求解 (3)符号与数值的转换（★） Numeric，eval：将符号常数转换为数字或计算表达式的函数； Sym2poly：将符号多项式变换为等价系数向量。 Poly2sym：要求用户指定用于表达式的变量（x，y等） 例3.6 f=sym(‘2*x^2+7*x+9’) n=sym2poly(f) f1=poly2sym([1 2 3],’t’) (4)因式分解与展开 factor(S): 对S分解因式，S是符号表达式