07 Matlab符号计算.pptVIP

MATLAB数学实验 第七章 Matlab符号计算 陆利正 * 主要内容 7.1 符号对象 7.2 符号矩阵和符号函数 7.3 符号微积分 7.4 符号方程和符号微分方程 7.5 便捷函数作图 7.6 符号计算局限性和Maple调用 7.1 符号对象 s = sym(str) 将数值或字符串str转化为符号对象, 数值为有理分式表示 s = sym(num, d) 将数值表达式转化为符号表达式, 数值用十进制表示 syms var1 var2 … 定义var1, var2, … 为符号变量? subs(s, old, new) 将符号表达式 s 中的符号变量old用new代替 7.1 符号对象 2. 计算精度和数据类型转换 digits(n) 将数值计算精度设为n位 x = vpa(s) 求 s 的数值结果 x = vpa(s, n) 采用 n 位计算精度求 s 的数值结果 double(s) 将符号对象转化为双精度数值 char(s) 将符号对象转化为字符串 7.2 符号矩阵和符号函数 1. 矩阵 Matlab大部分矩阵和数组运算符及指令都可以用于符号, 如: inv, ./, .\, eig 等 7.2 符号矩阵和符号函数 2. 符号函数计算 factor(expr) 对expr作因式分解 expand(expr 将expr展开 collect(expr, v) 将expr按变量 v 合并同类项 simple(expr) 将expr化简 g = finverse(f, v) 求函数f(v)的反函数g(v) fg = compose(f, g) 求函数f(v)和g(v)的复合f(g(v)) [n, d]=numden(expr) 分式通分, 分子 n, 分母 d funtool 函数分析图形界面 例7.1 已知 f (x, y) = (x－y)3, g(x, y) = (x+y)3, 求 h = f * g 并展开. 解 syms x y; f = (x - y)^3; g = (x + y)^3; h = f * g, s = expand(h), s2 = subs(s, y, x^2+x+1), 7.3 符号微积分 1. 极限和级数 limit(s, x, a) 返回符号表达式 s 当 x ? a 时的极限 limit(s, x, a, right) 返回 s 当 x ? a 时的右极限 limit(s, x, a, left) 返回 s 当 x ? a 时的左极限 symsum(s, n, a, b) 求和 例7.2 syms x n, limit((1+x/n)^n, n, inf) % exp(x) symsum((-1)^n*x^n/n, n, 1, +inf) % -log(1+x) 7.3 符号微积分 2. 微分 diff(s, x) 返回符号表达式 s 对 x 的导函数 diff(s, x, n) 返回 s 对 x 的 n 阶导函数 taylor(s, n, a, x) 返回 s 在 a 点Taylor展开到 n－1 次式, 自变量为 x taylortool Taylor分析图形界面 jacobian(f, x) f 的Jacobian矩阵 例7.3 (1) (2) f (x1, x2) = 的Jacobian矩阵 (3) 函数 g(x) = ln x sin x 在 x = 1 的5次Taylor展开 例7.3 (1) (2) f (x1, x2) = 的Jacobian矩阵 (3) 函数 g(x) = ln x sin x 在 x = 1 的5次Taylor展开 7.3 符号微积分 3. 积分 int(s) 符号表达式 s 的不定积分 int(s, v) s 关于变量 v 的不定积分 int(s, a, b) s 在区间 [a, b] 上的定积分 int(s, v, a, b) s 关于变量 v 从 a 到 b 的定积分 注: 当系统求不出解析解, 会自动求原点附近的一个近似解 例7.4 7.4 符号方程和符号微分方程 1. 代数方程的解 s = solve(表达式1, 表