第三讲MATLAB的符号运算.ppt

第三讲 MATLAB的符号运算 —— matlab 不仅具有数值运算功能，还开发了在matlab环境下实现符号计算的工具包Symbolic Math Toolbox 符号运算的功能 符号表达式、符号矩阵的创建 符号线性代数 因式分解、展开和简化 符号代数方程求解 符号微积分 符号微分方程 一、符号运算的基本操作 1. 什么是符号运算 与数值运算的区别 ※ 数值运算中必须先对变量赋值， 然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 2. 符号变量与符号表达式 MATLAB提供了两个建立符号对象的函数：sym和 syms，两个函数的用法不同。 sym函数 sym函数用来建立单个符号量，一般调用格式为： 符号量名=sym(符号字符串) 该函数可以建立一个符号量，符号字符串可以是常 量、变量、函数或表达式。 数值矩阵A=[1,2;3,4] A=[a,b;c,d] —— 不识别 1） 用matlab函数sym创建矩阵命令格式： A=sym([ ]) ※ 符号矩阵内容同数值矩阵 ※ 需用sym指令定义 ※ 需用‘ ’ 标识 例：A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] 这就完成了一个符号矩阵的创建。 2）用字符串直接创建矩阵 二、符号运算 数值运算中，所有矩阵运算操作指令都比 较直观、简单。例如：a=b+c; a=a*b ； A=2*a^2+3*a-5等。 而符号运算就不同了，所有涉及符号运算 的操作都可以用专用函数来进行。这些函数作用于 符号矩阵时，是分别作用于矩阵的每一个元素。 符号矩阵运算的函数： symadd(a,b) —— 符号矩阵的加 symsub(a,b) —— 符号矩阵的减 symmul(a,b) —— 符号矩阵的乘 symdiv(a,b) —— 符号矩阵的除 sympow(a,b) —— 符号矩阵的幂运算 symop(a,b) —— 符号矩阵的综合运算 例1： f= 2*x^2+3*x-5; g= x^2+x-7; 法①： h= symadd(f,g) % 符号矩阵的加法运算 h= 3*x^2+4*x-12 符号一些特殊运算函数： symsize —— 求符号矩阵维数 charploy —— 特征多项式 determ —— 符号矩阵行列式的值 eigensys —— 特征值和特征向量 inverse —— 逆矩阵 transpose —— 矩阵的转置 jordan —— 约当标准型 simple —— 符号矩阵简化 2. 符号微积分与积分变换 diff(f) —对f表达式的缺省变量求微分 diff(f,v) —对f表达式的指定变量v求微分 diff(f,v,n) —对f表达式的指定变量v求n阶微分 int(f) — 对f表达式的缺省变量求积分 int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在（a,b）区间求定积分 int（被积表达式，积分变量，积分上限， 积分下限）—— 定积分 例2、计算二重不定积分 例3、 例4、 3. 符号代数方程求解 f=a*x^2+b*x+c; ％利用单引号来生成符号表达式。 solve(f) ％对缺省变量x求解 ans = [ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))] solve(f , b ) %对指定变量b求解 ans = -(a*x^2+c)/x f1=solve(cos(x)=sin(x)), f1 = 1/4*pi f2=solve(tan(2*x)=sin(x)) f2 = [ 0] [