Matlab符号运算PPT课件

《Matlab符号运算PPT课件》由会员分享，可在线阅读，更多相关《Matlab符号运算PPT课件（26页珍藏版）》请在装配图网上搜索。

1、西南交通大学数学建模西南交通大学数学建模Matlab 符号运算Matlab 符号运算介绍符号运算介绍q Matlab 符号运算是通过符号运算是通过符号数学工具箱符号数学工具箱（Symbolic Math Toolbox）来实现的。来实现的。Matlab 符号数学工具箱是建立在功能符号数学工具箱是建立在功能强大的强大的 Maple 软件的基础上的，当软件的基础上的，当 Matlab 进行符号运算时，进行符号运算时，它就请求它就请求 Maple 软件去计算并将结果返回给软件去计算并将结果返回给 Matlab。q Matlab 的符号数学工具箱可以完成几乎所有得符号运的符号数学工具箱可以完成几乎所有

2、得符号运算功能。主要包括：算功能。主要包括：符号表达式的运算符号表达式的运算，符号表达式的复符号表达式的复合合、化简化简，符号矩阵的运算符号矩阵的运算，符号微积分符号微积分、符号作图符号作图，符符号代数方程求解号代数方程求解，符号微分方程求解符号微分方程求解等。此外，该工具箱等。此外，该工具箱还支持还支持可变精度运算可变精度运算，即支持以指定的精度返回结果。，即支持以指定的精度返回结果。Matlab 符号运算特点符号运算特点u 计算以计算以推理方式推理方式进行，因此不受计算误差累积所带来的进行，因此不受计算误差累积所带来的困扰。困扰。u 符号计算指令的调用比较简单，与数学教科书上的公式符号计算

3、指令的调用比较简单，与数学教科书上的公式相近。相近。u 符号计算可以给出完全正确的封闭解，或任意精度的数符号计算可以给出完全正确的封闭解，或任意精度的数值解（封闭解不存在时）。值解（封闭解不存在时）。u 符号计算所需的运行时间相对较长。符号计算所需的运行时间相对较长。Matlab 符号运算举例符号运算举例u 求一元二次方程 ax2+bx+c=0 的根 solve(a*x2+b*x+c)u 求的根 f(x)=(cos x)2 的一次导数 x=sym(x);diff(cos(x)2)u 计算 f(x)=x2 在区间 a,b 上的定积分 syms a b x;int(x2,a,b)q 在进行符号运算

4、时，必须先定义基本的在进行符号运算时，必须先定义基本的符号对象符号对象，可以是，可以是符号常量、符号变量、符号表达式等。符号对象是一种数据符号常量、符号变量、符号表达式等。符号对象是一种数据结构。结构。符号对象与符号表达式符号对象与符号表达式q 含有符号对象的表达式称为含有符号对象的表达式称为符号表达式符号表达式，Matlab 在内在内部把符号表达式表示成字符串，以与数字变量或运算相区部把符号表达式表示成字符串，以与数字变量或运算相区别。别。q 符号矩阵符号矩阵/数组：元素为符号表达式的矩阵数组：元素为符号表达式的矩阵/数组。数组。u sym 函数用来建立单个符号变量，一般调用格式为：q 符号

5、对象的建立：sym 和 syms符号对象符号对象的建立的建立例：a=sym(a)符号变量=sym(A)参数 A 可以是一个数或数值矩阵，也可以是字符串a 是符号变量b 是符号常量 b=sym(1/3)c 是符号矩阵 c=sym(1 ab;c d)q 符号对象的建立：sym 和 syms符号对象符号对象的建立的建立u syms 命令命令用来建立多个符号变量，一般调用格式为：syms 符号变量符号变量1 符号变量符号变量2.符号变量符号变量n 例：syms a b c a=sym(a);b=sym(b);c=sym(c);q 符号表达式的建立：例：建立符号表达式通常有以下2种方法：(1)用 sym

6、 函数直接建立符号表达式。(2)使用已经定义的符号变量组成符号表达式。y=sym(sin(x)+cos(x)x=sym(x);y=sin(x)+cos(x)符号符号表达式表达式的建立的建立Matlab 符号运算采用的运算符和基本函数，在形状、名称和使用上，都与数值计算中的运算符和基本函数完全相同符号对象符号对象的的基本基本运算运算q 基本运算符u 普通运算：+、-、*、/、u 数组运算：.*、.、./、.u 矩阵转置：、.例：X=sym(x11,x12;x21,x22;x31,x32);Y=sym(y11,y12,y13;y21,y22,y23);Z1=X*Y;Z2=X.*Y;符号对象符号对象

7、的的基本基本运算运算sin、cos、tan、cot、sec、csc、asin、acos、atan、acot、asec、acsc、exp、log、log2、log10、sqrtabs、conj、real、imagrank、det、inv、eigdiag、triu、tril三角函数与反三角函数、指数函数、对数函数等q 基本函数符号表达式的替换符号表达式的替换subs(f,x,a)用用 a 替换字符函数替换字符函数 f 中的字符变量中的字符变量 x a 是可以是是可以是 数数/数值变量数值变量/表达式表达式 或或 字符变量字符变量/表达式表达式若 x 是一个由多个字符变量组成的数组或矩阵，则 a 应

8、该具有与 x 相同的形状的数组或矩阵。q 用给定的数据替换符号表达式中的指定的符号变量subs 举例举例 f=sym(2*u);subs(f,u,2)f2=subs(f,u,u+2)a=3;subs(f2,u,a+2)subs(f2,u,a+2)syms x y f3=subs(f,u,x+y)subs(f3,x,y,1,2)ans=4f2=2*(u+2)ans=14ans=2*(a+2)+2)f3=2*x+2*yans=6u 例：指出下面各条语句的输出结果f=2*u六六类常见符号运算类常见符号运算q 因式分解、展开、合并、简化及通分等q 计算极限q 计算导数q 计算积分q 符号求和q 代数方

9、程和微分方程求解因式分解因式分解u 因式分解factor(f)syms x;f=x6+1;factor(f)l factor 也可用于正整数的分解 s=factor(100)factor(syml 大整数的分解要转化成符号常量大整数的分解要转化成符号常量函数展开函数展开u 函数展开expand(f)syms x;f=(x+1)6;expand(f)l 多项式展开l 三角函数展开 syms x y;f=sin(x+y);expand(f)合并同类项合并同类项u 合并同类项collect(f,v):按指定变量按指定变量 v 进行进行合并合并collect(f):按按默认默认变量变量进行进行合并合并

10、 syms x y;f=x2*y+y*x-x2+2*x;collect(f)collect(f,y)函数简化函数简化u 函数简化y=simple(f):对对 f 尝试多种不同的算法尝试多种不同的算法进行进行简化简化，返回其中最简短的形式，返回其中最简短的形式How,y=simple(f):y 为为 f 的的最简短形式最简短形式，How 中记录的为简化过程中使用的方法。中记录的为简化过程中使用的方法。fRHOW2*cos(x)2-sin(x)23*cos(x)2-1simplify(x+1)*x*(x-1)x3-xcombine(trig)x3+3*x2+3*x+1(x+1)3factorcos

11、(3*acos(x)4*x3-3*xexpand函数简化函数简化u 函数简化y=simplify(f):对对 f 进行简化进行简化 syms x;f=sin(x)2+cos(x)2;simplify(f)syms c alpha beta;f=exp(c*log(sqrt(alpha+beta);simplify(f)函数简化函数简化举例举例 syms c alpha beta;f=(1/x3+6/x2+12/x+8)(1/3);y1=simplify(f)g1=simple(f)g2=simple(g1)l 多次使用多次使用 simple 可以达到最简表达。可以达到最简表达。例：简化例：简化

12、分式通分分式通分u 函数简化N,D=numden(f):N 为通分后的分子，为通分后的分子，D 为通分后的分母为通分后的分母 syms x y;f=x/y+y/x;N,D=numden(f)计算极限计算极限limit(f,x,a):计算计算limit(f,a):当当默认变量默认变量趋向于趋向于 a 时的极限时的极限limit(f):计算计算 a=0 时的极限时的极限limit(f,x,a,right):计算右极限计算右极限limit(f,x,a,left):计算左极限计算左极限例：计算例：计算 ，syms x h n;L=limit(log(x+h)-log(x)/h,h,0)M=limit(

13、1-x/n)n,n,inf)计算计算导数导数g=diff(f,v)：求符号表达式求符号表达式 f 关于关于 v 的导数的导数g=diff(f)：求符号表达式求符号表达式 f 关于关于默认变量默认变量的导数的导数g=diff(f,v,n)：求求 f 关于关于 v 的的 n 阶导数阶导数q diff syms x;f=sin(x)+3*x2;g=diff(f,x)计算计算积分积分int(f,v,a,b):计算定积分计算定积分int(f,a,b):计算关于计算关于默认变量默认变量的定积分的定积分int(f,v):计算不定积分计算不定积分int(f):计算关于计算关于默认变量默认变量的不定积分的不定积

14、分 syms x;f=(x2+1)/(x2-2*x+2)2;I=int(f,x)K=int(exp(-x2),x,0,inf)例：计算例：计算 和和符号求和符号求和 syms n;f=1/n2;S=symsum(f,n,1,inf)S100=symsum(f,n,1,100)symsum(f,v,a,b):求和求和symsum(f,a,b):关于关于默认变量默认变量求和求和例：计算级数例：计算级数 及其前及其前100项的部分和项的部分和例：计算函数级数例：计算函数级数 syms n x;f=x/n2;S=symsum(f,n,1,inf)代数方程和微分方程代数方程和微分方程求解求解u 代数方程

15、求解solve(f,v)：求方程关于指定自变量的解，求方程关于指定自变量的解，f 可以是可以是用字符串表示的方程用字符串表示的方程、符号表达式符号表达式或或符号方程符号方程；l solve 也可解方程组也可解方程组(包含非线性包含非线性)；l 得不到解析解时，给出数值解。得不到解析解时，给出数值解。u 微分方程求解其它运算其它运算u 反函数finverse(f,v)：求求 f 关于指定变量关于指定变量 v 的反函数的反函数finverse(f)：求求 f 关于默认变量的反函数关于默认变量的反函数 syms x t;f=x2+2*t;g1=finverse(f,x)g2=finverse(f,t)例：计算函数例：计算函数 的反函数的反函数