工程优化设计--MALAB的符号运算.ppt

2、MATLAB的符号运算 ;符号运算的功能;? 运算对象可以是没赋值的符号变量 ? 可以获得任意精度的解 Symbolic Math Toolbox——符号运算工具包通过调用Maple软件实现符号计算的。 maple软件——主要功能是符号运算，它占据符号软件的主导地位。　　 ;1.符号运算的基本操作;2. 符号变量与符号表达式; ;用字符串直接创建矩阵 模仿matlab数值矩阵的创建方法 需保证同一列中各元素字符串有相同的长度。 例：A =[[ a,2*b]; [3*a, 0]] A = [ a, 2*b] [3*a, 0] ;用函数sym创建矩阵（symbolic) 命令格式：A=sym([ ]) ※ 符号矩阵内容同数值矩阵 ※ 需用sym指令定义 ※ 需用 标识 ;例如：A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] 这就完成了一个符号矩阵的创建。 注意：符号矩阵的每一行的两端都有方括号，这是与 matlab数值矩阵的一个重要区别。;符号矩阵的修改 a.直接修改 可用?、 ?键找到所要修改的矩阵，直接修改 b.指令修改 用A1=subs(A, new, old)来修改 ;例如： A =[ a, 2*b] [3*a, 0];符号矩阵与数值矩阵的转换 将数值矩阵转化为符号矩阵 函数调用格式：sym(A) A=[1/3,2.5;1/0.7,2/5] A = 0.3333 2.5000 1.4286 0.4000 sym(A) ans = [ 1/3, 5/2] [10/7, 2/5];将符号矩阵转化为数值矩阵 函数调用格式： numeric(A) A = [ 1/3, 5/2] [10/7, 2/5] numeric(A) ans = 0.3333 2.5000 1.4286 0.4000;符号矩阵运算 数值运算中，所有矩阵运算操作指 令都比较直观、简单。例如：a=b+c; a=a*b ；A=2*a^2+3*a-5等。 而符号运算就不同了，所有涉及符 号运算的操作都有专用函数来进行。;例1：f= 2*x^2+3*x-5; g= x^2+x-7; syms x f=2*x^2+3*x-5; g= x^2+x-7; h=f+g h = 3*x^2+4*x-12 例2：f=cos(x);g= sin(2*x); syms x f=cos(x);g=sin(2*x); f/g+f*g ans = cos(x)/sin(2*x)+cos(x)*sin(2*x) ;2.因式的分解、展开、化简;例3：syms a x;f1=x^5*x^3+5*x^2+5*x-6;factor(f1) ans = (x-1)*(x-2)*(x-3)*(x+1) 例4： syms x factor(x^9-1) ans = (x-1)*(x^2+x+1)*(x^6+x^3+1) ;例5：syms x y; expand(cos(x+y)) ans = cos(x)*cos(y) - sin(x)*sin(y) 例6： f=sym(cos(x)^2+sqrt(x^2+2*x+1)+sin(x)^2); F=simple(f); F=simple(F);2.符号极限;(3) limit(f)：求符号函数f(x)的极限值。符号函数f(x)的变量为函数findsym(f)确定的默认变量；没有