Matlab的学习资料.ppt

MATLAB 7.X程序设计 第5章 符号计算 (1) 主要内容： ①创建符号常量和变量； ②符号多项式的运算； ③符号微积分的计算； ④符号积分变换的计算； ⑤符号微分方程和代数方程的求解。 5.1 符号计算基础 5.1.1 概述 符号计算的特点： 1）符号计算定义在符号变量的基础上，符号表达式计算前必须定义符号变量。 2）符号计算是精确计算。 3）符号计算的计算速度较慢。 4）符号计算的运算符和基本数学函数与数值计算中的运算符和基本数学函数几乎完全相同。 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.2 符号表达式的基本操作 5.2.1 符号表达式的基本运算 1 算术运算、关系运算 2 常用的函数 1) 三角函数和双曲函数 2) 指数和对数函数 3) 矩阵分析函数 例5-3 创建符号矩阵，并求逆、对角阵、行 列式和秩。 4) 取整函数 5）复数函数 5.2.2 符号数值精度控制 三种不同的算术运算： 1) 数值型 MATLAB的浮点运算。 2) 有理数型 Maple的精确符号运算。 3) VPA型 Maple的任意精度运算。 第3种运算的精度可以任意设置,由digits和vpa函数实现。 digits(n) 功能：设定计算精度和改变默认的有效位数函数。 其中n为所期望的有效位数，默认值为32位。 Vpa 的格式为： S=vpa(s,n) 功能：将s表示为n位有效位数的符号对象。 例5-6 应用digits和vpa函数设置运算精度。 5.2.3 符号对象转换为数值对象 double(S) 功能：把符号矩阵S转换为双精度浮点数矩阵。 single(S) 功能：把符号矩阵S转换为单精度浮点数矩阵。 例5-7 建立符号常数矩阵，并转换 为数值矩阵。 s=sym([1/3 5/7;sqrt(3) 8/9]) %建立符号常数矩阵s single(s) %把s转换为数值矩阵 5.2.4 自由变量的确定 符号表达式允许使用自由变量。确定自由变量的原则： 1）小写字母i和j不能作为自由变量。 2）符号表达式中如果有多个字符变量，则按照以下顺序选择自由变量： 首先选择x作为自由变量；如果没有x，则选择在字母顺序中最接近x的字符变量；如果与x相同距离，则在x后面的优先。 3）大写字母比所有的小写字母都靠后。 findsym的格式为： findsym(EXPR,n) 功能：确定EXPR中的自由变量。 其中EXPR可以是符号表达式或符号矩阵；n为按顺序得出符号变量的个数，当n省略时，则不按顺序给出EXPR中所有的符号变量。 例5-8 创建符号表达式，然后确定符 号变量。 f1=sym(a*x^2+B*x+c) %创建符号表达式a*x^2+ B*x+c findsym(f1) %确定B, a, c, x为符号变量 findsym(f1,1) %确定x为符号变量 findsym(f1,2) %确定x,c为符号变量 findsym(f1,3) %确定x,c,a为符号变量 findsym(f1,4) %确定x,c,a,B为符号变量 5.2.5 变量置换 subs (S) 功能：用MATLAB工作空间中的变量替换S符号表达式中的所有变量。 例5-9 根据subs的置换规则，分析下 列语句的功能。 syms a b x; %创建符号变量a， b，和 x f=a*sin(x)+b; %创建符号表达式 f1=sub