MATLAB应用基础第五章.docVIP

第5章 符号运算 MATLAB的工具包Symbolic Math Toolbox以MAPLE为内核，实现了符号计算，如符号表达式与矩阵的操作、微积分、线性代数、方程求解、化简与展开、特殊数学函数等。 5.1符号表达式和符号矩阵的创建 1、用字符串创建符号表达式 如： f=sin(x)^2； ％创建函数sin2(x)赋给f q1=a*x^2+b*x+c = 0； ％创建方程赋给q1 de=Dy+y^2=1； ％创建微分方程y+y2 =1并赋给de 2、创建符号矩阵的方法 1）直接使用Sym命令矩阵 如： k=sym([a,b;c,d]) 运行结果如下： k = [ a, b] [ c, d] msy=sym([1/(a+x),sin(x),(b-x)/(a+x);2*x,1/(8+x^2),exp(x)]) 运行结果如下： msy = [ 1/(a+x), sin(x), (b-x)/(a+x)] [ 2*x, 1/(8+x^2), exp(x)] 使用Sym命令创建符号变量，再利用数组运算生成符号矩阵 如，下面的命令生成5阶符号矩阵： t=sym(t); ％创建符号变量t [I,J]=meshgrid(1:5) H=1./(I+J-t) 运行结果如下： I = 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 J = 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 5 5 5 5 5 H = [1/(2-t),1/(3-t),1/(4-t),1/(5-t),1/(6-t)] [1/(3-t),1/(4-t),1/(5-t),1/(6-t),1/(7-t)] [1/(4-t),1/(5-t),1/(6-t),1/(7-t),1/(8-t)] [1/(5-t),1/(6-t),1/(7-t),1/(8-t),1/(9-t)] [1/(6-t),1/(7-t),1/(8-t),1/(9-t),1/(10-t)] 3）使用Sym命令把数值矩阵化成符号矩阵 设A为数值矩阵，则: msy = Sym(A) 将A转化为符号矩阵msy。 如: A=[1 2 3;7 8 9;4 5 6] M=sym(A) 运行结果如下： A = 1 2 3 7 8 9 4 5 6 M = [ 1, 2, 3] [ 7, 8, 9] [ 4, 5, 6] 3、符号矩阵元素的引用和修改 Sym(S,i,j) ％引用符号矩阵元素Sij A=[1 2 3;7 8 9;4 5 6] M=sym(A) disp(第3行2列的元素为:) sym(M,3,2) 运行结果如下： A = 1 2 3 7 8 9 4 5 6 M = [ 1, 2, 3] [ 7, 8, 9] [ 4, 5, 6] 第3行2列的元素为: ans = 5 5.2 符号表达式与矩阵的基本运算 1、加、减、乘、除 Symadd(A,B) %A+B 与sym(A)+sym(B)等价 SymSub(A,B) %A–B 与sym(A)-sym(B)等价 Symmul(A,B) %A×B 与sym(A)*sym(B)等价 Symdiv(A,B) %A/B 与sym(A)/sym(B)等价 2、幂 sympow(S,p) %Sp 与sym(S)^sym(p)等价 若S为标量符号表达式，p可为标量符号或数值表达式； 若S为符号矩阵，p必须为整数 5.3因式分解、展开和简化及微积分 设s是一个符号表达式。 1、因式分解 factor(s) 如： f=sym(x^4-5*x^3+5*x^2+5*x-6) factor(f) 运行结果如下： f = x^4-5*x^3+5*x^2+5*x-6 ans = (-1+x)*(x-2)*(x-3)*(x+1) 2、展开 expand(s) 如： f=sym(x^4-5*x^3+5*x^2+5*x-6) g=fact