第8章 符运算.pptVIP

第8章 符号运算 以前的Matlab版本中，符号运算的功能比较薄弱，用户在解决比较复杂的符号运算时，单纯使用matlab往往无法求解 5.3开始，Matlab在符号运算方面有了很大提高，使符号计算和数值计算形式上统一。 7.0之后提供了更为强大的符号运算功能，提供了专门工具箱Symbolic Math Toolbox Matlab本身开发了许多专用于符号运算的函数 Matlab根据不同的模块，将这些函数分成如下几类 符号表达式和符号矩阵的操作 符号微积分 符号线性方程 符号微分方程 8.1 符号变量的生成和使用 sym函数和syms函数 8.1.1符号变量、符号表达式和符号方程 sqrt(2) sqrt(sym(2)) 2/5+1/3 sym(2)/sym(5)+sym(1)/sym(3) 使用sym函数定义符号表达式（两种方法） Eg.定义符号表达式 a=sym(‘a’); b=sym(‘b’); c=sym(‘c’); x=sym(‘x’); f=a*x^2+b*x+c 也可以整体定义 f=sym(‘a*x^2+b*x+c’) g=f^2+4*f-2 2使用syms定义符号变量和符号表达式 Syms a b c x f=sym(‘a*x^2+b*x+c’) 3符号方程的生成 Equation1=sym(‘sin(x)+cos(x)=1’) 8.1.2符号变量的基本操作 1 findsym函数:寻找符号变量 syms a alpha b x1 y findsym(alpha+a+b) findsym(cos(alpha)*b*x1+14*y,2) findsym(y*(4+3*j)+6*j) 2任意精度的符号表达式 digits函数 digits查询当前精度， digits(a)讲精度设置为a vpa函数 vpa(s)显示s在当前精度下的值 3 数值型与符号型变量的转换形式 使用sym可以转换成4种形式的符号变量 sym(t)或sym(t,’r’)有理数形式 sym(t,’f’)浮点数形式 sym(t,’e’)指数形式 sym(t,’d’)数值精度形式 A=hilb(4) sym(A) 8.1.3符号表达式的操作 1符号表达式的四则运算 syms x y a b fun1=sin(x)+cos(y) fun2=a+b fun1+fun2 fun1*fun2 2 合并符号表达式的同类项 collect函数 collect(S,v),collect(S) syms x y collect(x^2*y+y*x-x^2-2*x) 3符号多项式的因式分解 horner函数 syms x fun1=2*x^3+2*x^2-32*x+40 horner(fun1) 4 符号表达式的简化 simple函数 s=2*cos(s)^2-sin(x)^2; simple(s) 5 替换求值 subs(S),subs(S,new), subs(S,old,new) syms x y f=x^2*y+5*x*sqrt(y) subs(f,x,3) subs(f,y,3) subs(f,{x,y},{sym(‘alpha’),2}) 6反函数的运算 finverse()函数 syms x y f=x^2+y finverse(f,y) finverse(f) 7符合函数的运算 compose()函数 compose(f,g):f=f(x),g=g(y),f(g(y)) compose(f,g,z):f=f(x),g=g(y),f(g(z)) compose(f,g,x,z):f=cos(x/t),cos(g(z)/t) compose(f,g,x,y,z):f=cos(x/t),g=sin(y/u), cos(sin(z/u)/t) syms x y z t u f=x^t g=exp(-y/u) compose(f,g) compose(f,g,t,u,z) 9.2 符号矩阵的生成和运算 9.2.1符号矩阵的生成 1使用sym函数直接生成符号矩阵 a1=sym(‘[1/3 2/3 5/7;9/11 11/13 13/17;17/19 19/23 23/29]’) 2用生成子矩阵的方法生成符号矩阵 a=[‘[100,cos(x)]’;’[1/s,x ]’] 3由数值矩阵转换为符号矩阵 a=magic(3) s=sym(a) 9.2.2符号矩阵及符号数组的运算 1符号矩阵的四则运算 m=sym(‘[