第2章 matlab的符号运算.pptVIP

Matlab程序设计 Matlab程序设计 第2章 MATLAB的符号运算 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 介绍教材第2章内容 Matlab程序设计 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号矩阵分析和代数方程解 符号微积分 微分方程符号解法 Matlab程序设计 §1 符号运算的基本操作 什么是符号运算 与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算。 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 Matlab程序设计 特点： ? 运算对象可以是没赋值的符号变量 ? 可以获得任意精度的解 Symbolic Math Toolbox——符号运算工具包通过调用Maple软件实现符号计算。 Maple软件——主要功能是符号运算， 它占据符号软件的主导地位。 Matlab程序设计 §2 符号对象与符号表达式 例: f1 sym ‘[1/3,pi/7;sqrt 5 ,pi+sqrt 5 ]’ %符号矩阵 2 f2 sym Dy+y^2 1 %符号微分方程 3 syms a b c x; f3 a?x^2+b?x+c %二次三项式 2.1 生成符号对象的基本规则 1 用sym函数或syms函数定义基本符号对象 2 含符号对象的表达式或方程一定是符号对象 Matlab程序设计 例2.1-5: 区分数值矩阵、字符矩阵、符号矩阵 a 1;b 2;c 3;d 4; Mn [a,b;c,d] %2*2double型矩阵 Mc ‘[a,b;c,d]’ %1*9char矩阵 Ms sym Mc %2*2sym矩阵 Matlab程序设计 2.2 符号数字 sc sym ‘Num’ %符号常数sc的值精确等于Num 例：a pi + sqrt 5 %a为数值类常量 sa sym ‘pi + sqrt 5 ’ %sa为符号数字常量 % sa pi + sqrt 5 , sym型; eval sa 为5.3777, double型 2.3 符号参数 如表达式sin bx 中的b称为参数. syms Para1 Para2 ParaN %定义符号参数 syms Para1 Para2 ParaN Flag %定义Flag指定属性的符号参数 Flag: positive 正实数 , real 实数 , unreal 复数 例： syms k positive %等价于下式 k sym k,positive ; Matlab程序设计 2.4 符号变量 符号变量与符号参数的创建方法相同，但表达式或方程中作用不同. 确定自由符号变量： findsym EXPR , N %确认EXPR中距离x最近的N个自由符号变量, 略去N表示全部 例2.1-1 用符号计算研究方程uz2+vz+w 0的解 syms u v w z Eq u*z^2+v*z+w; %符号方程 r_1 solve Eq %一个方程只能解一个未知数w 离x最近 findsym Eq,1 %只找一个自由符号变量,则找到w r_2 solve Eq,z Matlab程序设计 2.5 符号运算基本运算符 矩阵运算： + , - , * , / , \ , ^ , 数组运算： + , - , .* , ./ , .\ , .^, . 关系运算: , ~ 2.6 符号运算中函数指令 表2.1-1 三角、双曲函数：sin、cosh等 指数、对数函数：exp、expm、log 即ln 复数函数：conj 共轭 、real、abs 模 矩阵代数：diag、inv、det、eig 方程求解：solve 微积分函数：diff、int 绘图函数：ezplot Matlab程序设计 3.1 数值数字与符号数字间转换 sym Num,’a’ —— 将数值类数字转换为a指定的符号数字输出， a取值: r: 广义有理表示 默认项 , d: 十进制浮点近似 e: 带eps误差的理性近似, f: 十六进制浮点近似 double x —— 将符号常量转为16位双精度数值量 3.2 符号数字的任意精度计算 digits n —— 设置今后数值计算精度，缺省32位 vpa x,n —— n位有效数字的数值型符号结果,略去n由digits决定运算精度 §3 符号数字及表达式操作 Matlab程序设计 例2.2-1： digits sym