第四讲 符号运算与微积分.ppt

数学实验 第5章 符号计算 主要内容： ①创建符号常量和变量； ②符号多项式的运算； ③符号微积分的计算； ④符号积分变换的计算； ⑤符号微分方程和代数方程的求解。 5.1.1 Matlab 符号运算介绍 Matlab 符号运算特点 符号对象与符号表达式 5.1.2 创建符号常量 sym函数创建符号常量，格式1为： sym(常量) sym函数的格式2为： sym(常量,参数) 功能：把常量按参数的格式要求转换为符号常量。其中参数的含义： d 返回最接近的十进制数(默认位数为32位) f 返回该符号值最接近的浮点表示 r 返回该符号值最接近的有理数型(为系统默认方式)，可表示为p/q、p*q、10^q、p/q、2^q和sqrt(p)形式之一 e 返回最接近的带有机器浮点误差的有理值 5.1.3 创建符号变量 sym(‘变量’,参数) 功能：把变量定义为符号对象。其中参数用来设置限定符号变量的数学特性，有三种选择： ’positive’ 表示为“正、实”符号变量， ’real’ 表示为“实”符号变量， ’unreal’ 表示为“非实”符号变量。如果不限定则参数可省略。 syms函数的格式为： syms(‘arg1’, ‘arg2’, …,参数) syms arg1 arg2 …参数 功能：创建多个符号变量。 符号表达式的建立： 5.1.5符号数值精度控制 第3种运算的精度可以任意设置,由digits和vpa函数实现。 digits(n) 功能：设定计算精度和改变默认的有效位数函数。 其中n为所期望的有效位数，默认值为32位。 Vpa 的格式为： S=vpa(s,n) 功能：将s表示为n位有效位数的符号对象。 例5-6 应用digits和vpa函数设置运算精度。 5.1.6 符号对象转换为数值对象 double(S) 功能：把符号矩阵S转换为双精度浮点数矩阵。 single(S) 功能：把符号矩阵S转换为单精度浮点数矩阵。 例5-7 建立符号常数矩阵，并转换 为数值矩阵。 s=sym([1/3 5/7;sqrt(3) 8/9]) %建立符号常数矩阵s single(s) %把s转换为数值矩阵 5.1.7 变量置换 用给定的数据替换符号表达式中的指定的符号变量 例5-9 根据subs的置换规则，分析下 列语句的功能。 syms a b x; %创建符号变量a， b，和 x f=a*sin(x)+b; %创建符号表达式 f1=subs(f,a,sym(2)) %以符号常量2替换符号变量a f2=subs(f,{a,b},{sym(2),sym(3)}) %分别以符号常量2、3替换符号变量a、b f3=subs(f,{a,b},{3,2}) %分别以标量2、3替换符号变量a、b f4=subs(f,sin(x),sym(y)) %以符号变量y替换符号表达式sin(x) f5=subs(f,{a,b,x},{2,2,sym(pi/3)}) %分别以标量2、3替换符号变量a、b，以符号常量pi/3替换符号变量x f6=subs(f,{a,b,x},{2,2,pi/3}) %分别以标量2、3 、pi/3替换符号变量a、b、x 5.3 符号多项式函数 5.3.1 合并同类项 collect (S,v) 功能：合并多项式S的同类项。其中S是符号表达式，v是合并同类项的变量。 例5-14 合并多项式 的同类项。 5.3.2 因式分解 factor(s) 功能：把符号表达式S转换为多个因式相乘的形式，各多项式的系数均为有理数。 例5-16 分解因式 syms x a;f=factor(x^3-a^3) ; pretty(f,x) %显示符号多项式的易读表示形式 执行结果如下： ans = 2 2 -(a - x) (a + a x + x ) 例5-19 将数字1025分解为正整数之积。 factor(1025) 执行结果如下： ans = 5 5 41 分析：由执行结果可以看出1025=5×5×41，恰巧是素数之积。 例5-18 将多项式表示为嵌套形式。 syms x; horner([x^6-x^4-5*x^3+5*x^2+5*x-6]) 例5-20问入取何值时，齐次方程组有非0解？ 5.3.4 化简 simplify（s）