第二章符号运算96563课件.pptVIP

第二章 符号计算 2.1 符号对象和符号表达式 2.1.2 符号计算中的算符和基本函数 2.1.3 对象类别的识别 2.1.4 符号表达式中自由变量的确定 2.2.3 符号数值精度控制和任意精度的计算 数值计算受计算机字长的限制,每次数值操作都带有截断误差,因此任何数值计算不管采用什么算法都将产生积累误差.在matlab中,每个算术操作结果的相对精度约为16位数字.但是符号计算的结果是绝对准确的,不包含任何计算误差.符号计算中与数值精度计算有关的指令有： double(x)---把符号常数转化为16位相对精度的浮点数值对象 digits(n)---设置以后的数值计算以n位相对精度进行 xs=vpa(x)---在digits指定精度下，给出x的数值型符号结果xs xs=vpa(x,n)---在n位相对精度下，给出x的数值型符号结果xs 说明： ● double指令运作所得结果一定是双精度数值对象 ● 除了vpa(x,n)对特定符号对象指定具体精度外,所有vpa(x) 的精度都受到其前面的digits指令控制,digits指令缺省精 度为32位 ● x可以是符号对象或数值对象,但运行后所得结果xs一定是符 号对象 2.2.4 符号对象与其他数据对象间的转换 数值,符号,字符是matlab中三种不同的数据类型,matlab为每种数据类型提供了各自特定的生成指令和操作指令.为实现不同数据类型间的交互, matlab 提供了一系列的转换指令． 数值型 符号结果 符号常数 数值 符号量(表达式) 字符串(表达式) ASCII码 vpa double vpa sym double str2num,str2double, sscanf int2str, num2str, mat2str char sym char abs double 数值数字和符号数字之间的相互转换 1) 数值数字?符号数字 sk=sym(num, flag) flag: r----数值类数字的广义有理表达 d---十进制浮点近似表达 e---带误差的理性近似表达 f---十六进制浮点近似表达 PI=sym(pi, r); PI = pi 注意：在符号运算中，数值类数字会自动转换为符号数字 syms t; y=1/3*sin(t)+2/4 y = 1/3*sin(t)+1/2 sk=sym(num) sk=sym(num) PI=sym(pi, d); PI = 3.1415926535897931159979634685442 2) 符号数字?数值数字 double(num_sym) PI=sym(pi); double(PI) ans = 3.1416 2.2.5 符号表达式的复合函数和反函数 1) 复合函数运算 ● compose(f, g)---返回复合函数f(g(y)), 这里f=f(x), g=g(y) ● compose(f,g,z)---返回自变量为z的复合函数f(g(z)), 这里 2) 反函数 ● g=finverse(f)---返回符号函数f的反函数g,满足g(f(x))=x ● g=finverse(f,v)---返回自变量为v的符号函数f的反函数使得 f=f(x), g=g(y) g(f(v))=v,当f包含不止一个符号变量时,使用这种格式. 例2.2.5-1 复合函数和反函数示例 syms x y; f=1/x^3; g=tan(y); fg=compose(f,g) fin=finverse(fg) fg = 1/tan(y)^3 fin = atan(1/y^(1/3)) 2.3 符号微积分 1) 符号极限 limit(F,x,a)---计算F在x?a时的极限 limit(F,x,a, left/right)---计算左/右极限 例5.3-1 极限示例 syms x; f=3/(x^2+6); val1=limit(f,x,inf) val2=limit(f,x,0) val1 = 0 val2 = 1/2 例 已知f(x)=sin(|x|), 求fx(0), fx(x) f(x)=sin(|x|)曲线图 syms x; syms dx positive; x0=sym(0); f_p=sin(x); df_p=limit((subs(f_p, x, x+dx)-f_p)/dx, dx,0) df_p0=limit((sin(x0+dx)-sin(x0))/dx, dx, 0) 数学分析： f(x)=sin(|x|) 绘图 要绘制原函数在x0时的导函数曲线，应先将导函数离散化，采用通用置换指令进行离散化 deltx=pi/200; xx=0:de