[第3章 matlab符号计算.pptxVIP

第三章 符号计算 符号计算是指：解算数学表达式、方程时，不是在离散化的数值点上进行，而是凭借一系列恒等式和数学定理，通过推理和演绎，获得解析结果。符号运算和数值运算的区别符号运算数值运算运算方式推理、演绎数字化运算运算前变量是否需要赋值可以不赋值必须赋值运算结果结果是完全准确地解析解具有一个误差的数值解3.1 符号计算基础3.2 符号导数及应用3.3 符号积分3.4 级数3.5 代数方程符号求解3.6 常微分方程符号求解3.1 符号计算基础3.1.1 符号对象Ⅰ、建立符号常数语法：sc=sym(Num)作用：创建符号常量sc,其值为Num说明：1.Num代表一个具体的数字 2.Num必须用‘’括起来例3.1 创建符号常量，查看其与数值常量的区别。a=sym(pi+sqrt(5))%定义符号常量b=pi+sqrt(5)%定义数值常量Ca=class(a)%通过class查看数据类别Cb=class(b)%在32位精度意义上，计算两类数值的差值vpa(a-b)Ⅱ、建立符号参量与变量语法：1.sc=sym(para, flag) 2.syms para1 para2 …paraN flag说明：1.flag表示变量的数域，其取值可以为： 2.语法1中para与flag必须用括起来 3.语法2中多个para间仅能以空格隔开positive 正实数real实数Unreal复数缺省复数Ⅲ 、运算符1.基本运算符加+减-乘*左除\右除/幂^点乘.*点左除.\点右除./点幂.^共轭转置非共轭转置. 2.关系运算符 在符号对象的逻辑运算中，仅有等于（==），不等于（~=）Ⅳ 、建立符号表达式1.定义：由常量，参数，变量以及运算符构成的式子。2.创建方法：U=sym(‘表达式’)首先申明表达式中的符号常量、符号变量等，然后再定义表达式 例3.2 ：建立表达式3x2+5y+2xy+3方法一：U=sym(3*x^2+5*y+2*x*y+3)方法二：syms x y V=3*x^2+5*y+2*x*y+3Ⅵ、建立符号矩阵?用函数sym创建矩阵 命令格式：A=sym([ ]) 例如：A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0]注意：符号矩阵的每一行的两端都有方括号，这是与 matlab数值矩阵的一个重要区别。?用字符串直接创建矩阵 模仿matlab数值矩阵的创建方法,需保证同一列中各元素字符串有相同的长度。例如：A =[[ a,2*b]; [3*a, 0]] A = [ a, 2*b] [3*a, 0]? 符号矩阵的修改a.直接修改 可用?、 ?键找到所要修改的矩阵，直接修改b.指令修改 用A(索引)=‘new’来修改； 用A1=subs(A, ‘new’, ‘old’)来修改例如：A =[ a, 2*b] [3*a, 0]A(2,2)=4*bA = [ a, 2*b] [3*a, 4*b]A2=subs(A1, c, b) A2 =[ a, 2*c] [3*a, 4*c] ? 符号矩阵与数值矩阵的转换符号矩阵转化为数值矩阵调用格式： numeric(A)A =[ 1/3, 5/2][10/7, 2/5]numeric(A)ans = 0.3333 2.5000 1.4286 0.4000数值矩阵转化为符号矩阵A=[1/3,2.5;1/0.7,2/5]A = 0.3333 2.5000 1.4286 0.4000sym(A)ans =[ 1/3, 5/2][10/7, 2/5]Ⅶ 、表达式中的变量1.默认变量原则：1）与x的ASCII码值差的绝对值小的字母优先;2）差绝对值相同时，ASCII码值大的优先。例3.3：默认自变量情况a*x^2+b*x+c %默认自变量为x1/(4+cos(t))%默认自变量为t4*z/y%默认自变量为y2*z+w %默认自变量为w2.符号表达式中变量的认定指令格式：findsym(EXPR,n)说明：1.n缺省时，返回EXPR中所有自由变量； 2.n取值N时，返回EXPR中离x最近的N 个自由符号变量； 3.当EXPR为符号矩阵时，该指令的确认是针对整个矩阵的。3.1.2 基本的符号运算Ⅰ、符号表达式运算(1)符号表达式的四则运算例3.4符号表达式的四则运算示例。syms x y z;f=2*x+x^2*x-5*x+x^3 %符号表达式不是最简形式f=2*x/(5*x) %符号表达式为最简形式f=(x+y)*(x-y) %符号表达式的结果不是x^2-y^2 而是(x+y)*(x-y)(2)因式分解与展开1)factor(S) ：%