第二节符号运算.pptVIP

第二章 符号计算 2.1 符号对象和符号表达式 2.1.2 符号计算中的算符和基本函数 2.1.3 对象类别的识别 2.1.4 符号表达式中自由变量的确定 说明： ● 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:deltx:2*pi; y=subs(df_p, x, xx); plot(xx, y); 2) 符号序列求和 对于数学上的求和 可用matlab的求和指令解决.其指令为: s=symsum(f,v,a,b)---求通式f字指定变量v取遍[a,b]中所有整数时的和. 说明: ● f是矩阵时,对矩阵的元素逐个求和,但自变量定义在整个矩阵上 ● v缺省时 f中的自变量由findsym自动辨认;b可以取有限整数, 也可以取无穷大inf ● a,b可同时缺省, 此时默认求和的自变量区间为[0,v-1] 例2.3-2 求 syms n; f=[n^2, 1/n^2]; s1=symsu