第四章 matlab的符号运算.pptVIP

第四章 MATLAB的符号运算 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号代数方程求解 符号微积分 符号微分方程 §1 符号运算的基本操作 什么是符号运算 与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 特点： ? 运算对象可以是没赋值的符号变量 ? 可以获得任意精度的解 Symbolic Math Toolbox——符号运算工具包通过调用Maple软件实现符号计算的。 Maple软件——主要功能是符号运算， 它占据符号软件的主导地位。　　 2. 符号变量与符号表达式 例: (1) f1=sym([1/3,pi/7;sqrt(5),pi+sqrt(5)]) (2) f2=sym(Dy+y^2=1) %微分方程 (3) syms a b c x; f3= a?x^2+b?x+c %二次三项式 3.符号矩阵的创建 数值矩阵A=[1,2;3,4] A=[a,b;c,d] %a,b,c,d未定义,则不识别 （1）用函数sym创建 (symbolic的缩写) 命令格式：A=sym([ ]) A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] 4. 符号对象的修改 a.直接修改 可用?、 ?键找到所要修改的数组，直接修改 b.指令修改(通用转换指令subs) 用A1=subs(A, old, new)来修改 5. 符号矩阵与数值矩阵的转换 将数值矩阵转化为符号矩阵 函数调用格式：sym(A) A=[1/3,2.5;1/0.7,2/5] A = 0.3333 2.5000 1.4286 0.4000 §2 符号运算 符号矩阵运算 数值运算中，所有矩阵运算操作指 令都比较直观、简单。例如：a=b+c; a=a*b ；A=2*a^2+3*a-5等。 涉及符号运算的操作增加了一些专用函数 符号运算函数: (运算规则与数值运算基本相同) 矩阵运算： + , - , * , / , \ , ^ ,diag , tril, triu inv, det, rank, rref, null, poly, expm, eig, svd 数组运算： + , - , .* , ./ , .\ , .^ 关系运算； = = , ~ = 三角、双曲函数：sin、cosh等 指数、对数函数：sqrt、exp、log(即ln) 复数函数：conj(共轭)、real、abs (模) 矩阵代数：diag、inv、det(特征值) 例1：计算：S=(-7x2-8y2)*(-x2+3y2) 解：syms x y; s= (-7*x.^2-8*y.^2)*(-x.*x+3*y.*y) expand(s) %对x展开 ans = 7*x^4-13*x^2*y^2-24*y^4 collect(s,x) %对S按x合并同类项 ans = 7*x^4-13*x^2*y^2-24*y^4 factor(ans) %分解因式 ans = (8*y^2+7*x^2)*(x^2-3*y^2) findsym(ans) %找自由符号变量 ans = x, y 例2：化简 S=(x2+y2)2+(x2-y2)2 syms x y; S=(x^2+y^2)^2+(x^2-y^2)^2 simple(S) %系统自动试探各种函数化简 simple(ans) %使用多次找到最少字母的简化式 2. 符号数值精度控制 *符号运算的结果是绝对准确的，用有理数表示。 例： sym(1/2)+(1/3)