MATLAB符号运算功能.ppt

《MATLAB符号运算功能.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第三章符号运算功能——matlab不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包SymbolicMathToolbox符号运算的功能符号表达式、符号矩阵的创建符号线性代数因式分解、展开和简化符号代数方程求解符号微积分符号微分方程一、符号运算的基本操作什么是符号运算与数值运算的区别※数值运算中必须先对变量赋值，然后才能参与运算。※符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。2.符号表达式与方程（1）单引号法f1='ax^2+bx+c'——二次三项式f2='ax^2+bx+c=0'——方程f3='Dy+y^2=1'——微分方程符号表达式一定要用''

2、单引号括起来matlab才能识别。（2）sym或syms法格式:sym(‘变量’,参数)syms(‘变量1’,‘变量2’,···,参数)或syms变量1变量2···参数(其中参数为positive,real,unreal)y1=sym(‘a*exp(x)+2*cos(2*x)’)或：symsaxy1=a*exp(x)+2*cos(2*x)3.符号矩阵的创建函数sym创建矩阵（symbolic的缩写)命令格式：A=sym('[]')※符号矩阵内容同数值矩阵例如：A=sym('[a,2*b;3*a,0]')A=[a,2*b][3*a,0]这就完成了一个符号矩阵的创建。4.自由变量的确定除了或指定符

3、号变量以外，还允许使用自由变量，其确定原则：（a）小写字母不能作为自由变量；（b）有多个变量时，优先作为自由变量，若没有，则在词典序中最靠近的则优先作为自由变量；（c）小写字母的顺序优先于大写字母；（d）findsym函数可以确定自由变量的个数。格式：findsym(S,n)5.符号变量的修改a.直接修改可用鼠标修改矩阵，直接修改b.指令修改用A=subs(S,‘old’,‘new’)来修改或A=subs(S,‘new’)——new替代S中的自由变量注：多个参量需替代时，应放入{}中例如：A=sym('[a,2*b;3*a,0]')A=[a,2*b][3*a,0]A1=subs(A,0,’4*

4、b’)A(2,2)='4*b'A=[a,2*b][3*a,4*b]A2=subs(A1,‘b',‘c')A2=[a,2*c][3*a,4*c]A3=subs(A,{‘b',0}，{‘a’,‘c’})A3=[a,2*a][3*a,c]将数值矩阵转化为符号矩阵函数调用格式：sym(A)A=[1/3,2.5;1/0.7,2/5]A=0.33332.50001.42860.4000sym(A)ans=[1/3,5/2][10/7,2/5]6.符号对象转换为数值对象将符号矩阵转化为数值矩阵函数调用格式：double(A)、single(A)A=[1/3,5/2][10/7,2/5]single(A)an

5、s=0.33332.50001.42860.4000符号矩阵运算符号运算与数值运算中的运算符和基本函数几乎完全相同。二、符号运算例1：f=2*x^2+3*x-5;g=x^2+x-7;symsxf=2*x^2+3*x-5;g=x^2+x-7;h=f+gh=3*x^2+4*x-12例2：f=cos(x);g=sin(2*x);symsxf=cos(x);g=sin(2*x);f/g+f*gans=cos(x)/sin(x)+cos(x)*sin(x)2.任意精度的数学运算在symbolic中有三种不同的算术运算：数值类型matlab的浮点算术运算有理数类型maple的精确符号运算vpa类型mapl

6、e的任意精度算术运算任意精度算术运算digits(n)——设置可变精度，n位有效数字,缺省32位vpa(x,n)——将x表示为n位有效数字digits(25)vpa(1/2+1/3)ans=.8333333333333333333333333diff(f)—对缺省变量求微分diff(f,v)—对指定变量v求微分diff(f,v,n)—对指定变量v求n阶微分int(f)—对f表达式的缺省变量求积分int(f,v)—对f表达式的v变量求积分int(f,v,a,b)—对f表达式的v变量在（a,b)区间求定积分3.符号微积分与积分变换limit(f)—对自变量趋于0时的极限limit(f,x,a)—对

7、自变量x趋于a时的极限limit(f,x,a,‘right’)—自变量x从右边趋于a时的极限limit(f,x,a,‘left’)—自变量x从左边趋于a时的极限symsum(s,x,a,b)—对s表达式的级数求和taylor(f,n)——泰勒级数展开n-1次幂多项式taylor(f,a)——在x=a点泰勒级数展开例1.计算二重不定积分F=int(int('x*exp(-x*y)','x'),'y'