符號计算.pptVIP

第2章 符号计算;教学目标;教学重点;教学内容;Matlab的符号计算功能;符号运算与数值运算的区别：;2.1 符号对象和符号表达式;2.1.1 符号对象的创建和衍生;2 符号数字的定义;2 符号数字的定义;3. 基本符号变量：;4. 自由符号变量;syms u v w z a5 f=sym(3); Eq=sin(f)*u*z^2+v*z+f*w-a5;;syms a b x X Y k=sym(3); z=sym(c*sqrt(d)+y*sin(t)); EXPR=a*z*X+(b*x^2+k)*Y;;4. 自由符号变量;2.1.2 符号计算中的算符;2.1.2 符号计算中的算符;2.1.3 符号计算中的函数指令;(2) 指数,对数函数 sqrt, exp, expm在两者中用法相同．符号计算中只有自然对数log，而没有数值计算中的log2, log10;(4) 矩阵代数指令 在符号计算中, matlab提供的常用矩阵代数指令有：diag, tril, inv, det, rank, eig, svd（ Singular value decomposition奇异值分解）等;数值计算对象,符号计算对象,字符串是MATALB中最常用的数据对象.他们遵循各自不同的运算法则,但有时在外形上却十分相似.MATLAB提供了一些识别不同数据对象的指令,常用的有class, isa, whos ;SizeMn=size(Mn), SizeMc=size(Mc), SizeMs=size(Ms);(5) 利用whos观察内存变量的类别和其他属性;2.1.5 符号运算机理和变量假设;2.1.5 符号运算机理和变量假设;2.1.5 符号运算机理和变量假设;【例2.1-6】syms 对符号变量限定性假设的影响;2.1.6 符号帮助体系 ;2.1.6 符号帮助体系 ;2.1.6 符号帮助体系 ;2.2 符号数字及表达式的操作;2.2.1双精度数字与符号数字之间的转换;2.2.1双精度数字与符号数字之间的转换;digits 显示当前环境下符号数字“十进制浮点” 表示的有效数字位数； digits(n) 设定 “十进制浮点”表示的有效数字位数； xs=vpa(x) 据表达式x得到digits指定精度下的符号数字xs xs=vpa(x,n) 据表达式x得到n位有效数字的符号数字xs;reset(symengine) sa=sym(1/3+sqrt(2)) sa =2^(1/2) + 1/3 digits Digits = 32 format long a=1/3+sqrt(2) sa_Plus_a=vpa(sa+a,20) sa_Minus_a=vpa(sa-a,20) a = 1.747546895706428 sa_Plus_a = 3.4950937914128567869 sa_Minus_a = -0.000000000000000022658064826339973669 ;2.2.3 符号表达式的基本操作;2.2.4 表达式中的置换操作;2.2.4 表达式中的置换操作;(5); 符号限定假设对方程根不起作用。;通用置换指令 RES＝subs(ES,old,new) RES＝subs(ES,new);dfdvn=diff(f,v,n) 求 fjac=jacobian(f,v) 求多元向量函数f(v)的jacobian矩阵 r=taylor(f,n,v,a) 把f(v)在v=a处进行泰勒展开;【例2.3-1】试求; syms a t x f=[a, t^3; t*cos(x), log(x)]; df=diff(f) dfdt2=diff(f,t,2) dfdxdt=diff( diff(f,x) ,t) ;例2.3-5:设cos(x+siny)=siny, 求dy/dx(隐函数求导).;sym x r=taylor(x*exp(x),9,x,0) pretty(r) ;MATLAB求解通式求和 问题的指令为： s=symsum(f,v,a,b) 求通式f在指定变量v取遍 [a,b]中所有整数时的和。 ;syms n k f1=1/(k*(k+1)); s1=symsum(f1,k,1,n) s1 =1 - 1/(n + 1) ;2.3.3 符号积分;2.3.3 符号积分;例2.3-11求积分;syms a r theta phi