第3讲 符号运算功能.pptVIP

* 武汉理工大学数学系 尹强 数学实验三MATLAB的符号运算功能 函数的简化 微分、积分 、Taylor展开等 代数方程求解 微分方程求解 行列式、矩阵求逆、特征值、特征向量 可变精度计算 基本 概念 符 号 变 量：预先不要求有确定值的量 符号表达式：数字、函数、算子和变量的 Matlab字符串 符号方程式：含有等号的符号表达式 符 号 矩 阵：数组（其元素为符号表达式） 符号变量的定义方式： syms x y z c 定义符号表达式的方式： 1、隐含方式，如‘sin(x)’ 2、显式方式，如 sym(‘[a b;c d]’) 通常在Matlab可以确定变量类型的场合下，可以不用显示方式,如： diff sin(x) diff(‘sin(x)’) 但是在大多数场合，显式说明是必需的。 findsym(‘s^2+y^2’) 符号变量的确定： 在符号表达式中，含有多于一个符号变量时，当对表达式进行操作时，只有一个符号变量是独立变量，具体规则如下： 如果有x字母，则x为独立变量 如果没有x字母，则选取在字母表顺序中最接近x的字母（后面优先）不过在很多函数中，可对变量具体指定 简单表达式的产生： syms x,y f=‘a*x+x^2*y+sqrt(x)’ g=‘c*y-x’ A=sym(‘[a,b;c,d]’) add=symadd(f,g) f+g ; sub=symsub(f,g) f-g mul=symmul(f,g) f*g ; div=symdiv(f,g) f/g pow=sympow(f,g) 分子和分母提取 [n,d]=Numden(f) 注： numerator（分子） denominator（分母） 如： [n,d] = numden(x/y + y/x) n = x^2+y^2 , d = y*x 函数的复合 compose(f,g) ?f(g(y)) ，其中：f = f(x) and g = g(y); compose(f,g,u,v)…含有不同独立变量函 数的 复合 f = cos(x/t); COMPOSE(f,g,x,z); returns cos(g(z)/t) COMPOSE(f,g,t,z); returns cos(x/g(z)) COMPOSE(f,g,x,y,z) syms x y z t u; f = 1/(1 + x^2); g = sin(y); h = x^t; p = exp(-y/u) compose(f,g) returns 1/(1+sin(y)^2) compose(f,g,t) returns 1/(1+sin(t)^2) compose(h,g,x,z) returns sin(z)^t compose(h,g,t,z) returns x^sin(z) compose(h,p,x,y,z) returns exp(-z/u)^t compose(h,p,t,u,z) returns x^exp(-y/z) 变量替换函数 Subs(f) subs(f,new,old) syms a c1 t a = ‘980’; c1 =‘3’ y = exp(-a*t)*c1 subs(y,{‘a’,’c1’},{a,c1}) ? 3*exp(-980*t) subs(cos(a)+sin(b),{a,b},{sym(alpha),2}) ? cos(alpha)+sin(2) subs(exp(a*t),a,-magic(2)) returns [ exp(-t), exp(-3*t)] [ exp(-4*t), exp(-2*t)] subs(x*y,{x,y},{[0 1;-1 0],[1