数学建模 符号计算.ppt

通知明天（4.15）上午的培训课因有部分同学要参加考试，所以暂停一次。下周六下午正常。望相互转告！;符号计算;本节目标;本节重点;本节内容;Matlab的符号计算功能;符号运算与数值运算的区别：;2.1 符号对象和符号表达式;2.1.1 符号对象的创建和衍生;2 符号数字的定义;2 符号数字的定义;3. 基本符号变量：;4. 自由符号变量;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.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);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 r=a*theta; x=r*cos(theta); y=r*sin(theta); dLdth=sqrt(diff(x,theta)^2+diff(y,theta)^2); L=simple(int(dLdth,theta,0,phi));例2.3-12 画阿基米德螺线r = a*θ及其曲线长度图形。;2.4 微分方程的符号解法;说明： （1）输入量包括三部分：微分方程、初始条件、指定独立变量（不指定时，默认为t）