第3章Matlab的符号运算功能.pptVIP

第三章 MATLAB的符号运算 除了数值运算以外，在数学、工程和其他应用科学中还经常用到符号运算。MATLAB和著名的符号运算语言MAPLE相结合，为用户提供了集符号运算与符号可视化为一体的符号运算功能. 符号运算的内容 符号线性代数 微分、积分、泰勒展开 符号代数方程求解 符号微分方程 符号运算与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算. 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达. 特点： ? 运算对象可以是没赋值的符号变量 ? 可以获得任意精度的解 2. 符号变量与符号表达式 在使用符号变量之前，要声明哪些是符号变量，定义方法如下： x=sym(‘x’); % symbolic的缩写 y=sym(‘y’); 或： syms x y %以空格间隔，不能用“，” 符号表达式 f=x*y+2*x 3.符号矩阵的定义及代数运算 符号矩阵定义，与普通矩阵类似。 syms a b c d x A= [a b;c d]; =sym(‘[a b;c d]’) B=[2*a+b,3*b;5*c+d,2*d] =sym(‘[2*a+b,3*b;5*c+d,2*d]’) 符号矩阵的加、减、乘、除运算 A+B; symadd(A,B) sym(A)+sym(B) A-B; symsub(A,B) sym(A)-sym(B) A*B; symmul(A,B) sym(A)*sym(B) A/B; symdiv(A,B) sym(A)/sym(B) 显然： 参与加减运算的符号矩阵必须具有相同的维数. 参与乘除运算的矩阵符合矩阵乘法定义. 4.符号变量替换 当得到一个符号解或者在一个符号表达式中，需要将一些符号变量替换成数字或其他符号. subs(s,new)：用新变量new 替换s中的默认变量。 subs(s,new,old):用新变量new 替换s中的指定变量old。 5. 符号微积分与积分变换 diff(f) — 对缺省变量求微分 diff(f,v) — 对指定变量v求微分 diff(f,v,n) —对指定变量v求n阶微分 int(f) — 对f表达式的缺省变量求积分 int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在（a,b) 区间求定积分 例3-1 syms x t f1=sin(x^2); f2=t^6; r1=diff(f1,x); r2=diff(f2,t,6); 结果 r1=2*x*cos(x^2) r2=720 int（‘被积表达式’，‘积分变量’，‘积分下限’， ‘积分上限）—— 定积分 6.符号代数方程求解 例3-2. f = ax2+bx+c 求解 syms a b c x; f=a*x^2+b*x+c; solve(f) —— 对缺省变量x求解 ans = [1/2/a*(-b+(b^2-4*a*c)^(1/2))] [1/2/a*(-b-(b^2-4*a*c)^(1/2))] 例3-3. 符号方程cos(x)=sin(x) 求解 f1=solve(cos(x)=sin(x)), f1 = 1/4*pi numeric(f1) ans = 0.7854 例3-4. 解方程组 x+y+z=1 x-y+z=2 2x-y-z=1 g1=x+y+z=1,g2=x-y+z=2,g3=2*x-y-z=1 f=solve(g1,g2,g3) f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1) f = x: [1x1 sym] f.x ans =2/3 y: [1x1 sym] f.y ans =-1/2 z: [1x1 sym] f.z ans =5/6 [x,y,z]=solve(x+y+z=1,x-y+z=2,2*x-y-z=1) x = 2/3 y =-1/2 z =5/6 或者 7. 符号微分方程求解 —— 用一个函数可以方便地得到微 分方程的符号解 符号微分方程求解指令：dsolve 命令格式：dsolve(f,g) f —— 微分方程，可多至12个微分方程的求