第六讲MATLAB符号计算.doc

第六讲 MATLAB符号运算 MATLAB 不仅具有数值运算功能，还开发了在MATLAB环境下实现符号计算的工具包Symbolic Math Toolbox和大量的符号运算函数，符号运算的主要功能有 (1)符号表达式、符号矩阵的创建； (2)符号线性代数； (3)因式分解、展开和简化； (4)符号代数方程求解； (5)符号微积分； (6)符号微分方程。 第一节 符号运算的基本操作 一、概述 1、什么是符号运算 符号运算与数值运算的区别： (1)数值运算中必须先对变量赋值，然后才能参与运算； (2)符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。?? 2、特点 (1)运算对象可以是没赋值的符号变量； (2)可以获得任意精度的解析解。 二、符号矩阵的创建 数值矩阵A=[1,2;3,4]可以，A=[a,b;c,d]不可以。 MATLAB提供了两个建立符号对象的函数：sym和syms，两个函数的用法不同。 1、sym函数 sym函数用来建立单个符号量，一般调用格式为： 符号量名=sym(符号字符串) 该函数可以建立一个符号量，符号字符串可以是常量、变量、函数或表达式。 应用sym函数还可以定义符号常量，使用符号常量进行代数运算时和数值常量进行的运算不同。 例： A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [ 3*a, 0] 这就完成了一个符号矩阵的创建。 注意：符号矩阵的每一行的两端都有方括号，这是与 MATLAB数值矩阵的一个重要区别。 把字符表达式转换为符号变量。 例： y=sym(2*sin(x)*cos(x)) y = 2*cos(x)*sin(x) y=simple(y) y = sin(2*x) 2、syms函数 函数sym一次只能定义一个符号变量，使用不方便。MATLAB提供了另一个函数syms，一次可以定义多个符号变量。syms函数的一般调用格式为： syms 符号变量名1 符号变量名2 … 符号变量名n 用这种格式定义符号变量时不要在变量名上加字符串分界符(‘)，变量间用空格而不要用逗号分隔。 例： ①用符号计算验证三角等式。 syms fai1 fai2 ; y=simple(sin(fai1)*cos(fai2)-cos(fai1)*sin(fai2)) y = sin(fai1 - fai2) ②求矩阵的行列式值、逆和特征根。 syms a11 a12 a21 a22;A=[a11,a12;a21,a22] A = [ a11, a12] [ a21, a22] DA=det(A),IA=inv(A),EA=eig(A) DA = a11*a22 - a12*a21 IA = [ a22/(a11*a22 - a12*a21), -a12/(a11*a22 - a12*a21)] [ -a21/(a11*a22 - a12*a21), a11/(a11*a22 - a12*a21)] EA = a11/2 + a22/2 - (a11^2 - 2*a11*a22 + a22^2 + 4*a12*a21)^(1/2)/2 a11/2 + a22/2 + (a11^2 - 2*a11*a22 + a22^2 + 4*a12*a21)^(1/2)/2 符号常数形成中的差异 a1=[1/3,pi/7,sqrt(5),pi+sqrt(5)] %a1是数值常数 1 a2=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) %最接近的有理表示 2 a3=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)],e) %带估计误差的有理表示 3 a4=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) %绝对准确的符号数值表示 4 a24=a2-a4 a1 = 0.3333 0.4488 2.2361 5.3777 a2 = [ 1/3, pi/7, 5^(1/2), 189209612611719/35184372088832] a3 = [ 1/3 - eps/12, pi/7 - (13*eps)/165, (137*eps)/280 + 5^(1/2), 189209612611719/35184372088832] a4 = [ 1/3, pi/7, 5^(1/2), pi + 5^(1/2)] a24 = [ 0, 0, 0, 189209612611719/35184372088832 - 5^(1/2) - pi] 第二节 因式分解、展开和简化 一、因式分解 1、factor