matlab第三讲符号计算.ppt

符号矩阵的基本运算与数值矩阵的基本运算相类似。 1) 基本运算符：+、-、*、\、/、、.*、.\、./、.^、’ 、.’ 2) 三角函数与反三角函数：sin、cos、tan、… … 3) 指数、对数函数：sqrt、exp、log、… … 4) 复数函数：real、imag、conj、 abs 5) 矩阵函数：det、inv、rank、 … … （没有norm） 6) 矩阵元素的抽取：diag、tril、triu 因式分解、展开、合并、简化及通分等 因式分解：factor syms x f=x^6 +1 s=factor(f) s =(1+x^2)*(x^4-x^2+1) factor 也可用于正整数的分解 问题： factor(sym(300000000000000000000000000000000000000)) ans = (2)^76*(23)*(172628997402371) factor(sym(300000000000000000000000000000000000000)) ans = (2)^38*(3)*(5)^38 多项式展开 三角函数展开 该函数经常用于多项式展开，也常用于三角函数、指数函数和对数函数的展开中。 collect(f,v): 按指定变量 v 的次数合并系数； collect(f): 合并 f 中的默认自变量的各项系数。 findsym(f,1) simple(f): 对 f 尝试多种不同的算法简化， 返回其中最短的简化形式； [R,HOW]=simple(f): R为f的最短简化形式， HOW中记录的为简化过程中使用的主要方法。 f R HOW 2*cos(x)^2-sin(x)^2 3*cos(x)^2-1 simplify (x+1)*x*(x-1) x^3-x combine(trig) x^3+3*x^2+3*x+1 (x+1)^3 factor cos(3*acos(x)) 4*x^3-3*x expand simple函数示例 simplify(f): 简化函数 注：多次 使用 simple 可以达到最简表达。 例：化简 [N,D]=numden(f): N为通分后的分子，D为通分后的分母 例： 因式分解、展开、合并、简化及通分等 计算极限 limit(f,x,a): 计算 limit(f,a): 计算默认自变量趋向于a时f的极限 limit(f): 计算 a=0 时的极限 limit(f,x,a,’right’)：右极限 limit(f,x,a,’left’)：左极限 syms h n x L=limit((log(x+h)-log(x))/h,h,0) M=limit((1-x/n)^n,n,inf) L=1/x M=exp(-x) syms x L=limit(abs(x)/x,x,0,’left) R=limit(abs(x)/x,x,0,’right) L=-1 M=1 因式分解、展开、合并、简化及通分等 计算极限 计算导数 diff(f): 计算 f 关于默认自变量的导数 diff(f,’v’): 计算 f 关于变量 v 的导数 diff(f,n),diff(f,’v’,n),diff(f,n,’v’): n次求导 syms a x y=sin(a*x) A=diff(y,x) B=diff(y,a) C=diff(y,x,2) D=diff(y,a,2) A=cos(a*x)*a B=cos(a*x)*x C=-sin(a*x)*a^2 D=-sin(a*x)*x^2 因式分解、展开、合并、简化及通分等 计算极限 计算导数 计算积分 int(f,v,a,b):计算定积分 int(f,a,b): 计算 f 关于默认自变量 的定积分 int(f,v)：计算不定积分 int(f)：计算 f 关于默认自变量 的不定积分 syms x f=(x^2+1)/(x^2-2*x+2)^2; I=int(f) g=cos(x)/(sin(x)+cos(x)); J=int(g,x,0,pi/2) h=exp(-x^2); K=int(h,x,0,inf) I=3/2*atan(x-1)+1/4*(2*x-6)/(x^2-2*x+2) J=1/4*pi K=1/2*pi^(1/2) 因式分解、展开、合并、简化及通分等 计算极限 计算导数 计算积分 符号求和 symsum(f,v,a,b): 求和 symsum(f,a,b): 关于默认自变量 求和。 syms n S=symsum(1/n^2,n,1,inf)