[2018年最新整理]2符号计算与数值计算.docVIP

第二章 符号计算 数学计算有数值计算与符号计算之分。这两者的根本区别是：数值计算的表达式、矩阵变量中不允许有未定义的自由变量，而符号计算可以含有未定义的符号变量。对于一般的程序设计软件如C,C++等语言实现数值计算还可以，但是实现符号计算并不是一件容易的事。而Matlab自带有符号工具箱Symbolic Math Tooibox,再借助于软件Maple，就具有了强大的符号运算功能。 在数值计算（包括输入、输出及中间变量在内的）过程中，所运作的变量都是被赋了值的数值变量。而在符号计算的整个过程中，所运作的是符号变量。注意：在符号计算中所出现的数字也都是当作符号处理的。 符号变量的创建 在Matlab的数据类型中，字符型与符号型是两种重要而又容易混淆的数据类型。一定要注意这两种类型的不同之处。在Matlab中利用sym命令来创建符号变量，sym的使用格式为： S=sym(‘A’) 定义单个符号变量S syms a b c …. 定义多个符号变量a,b,c… syms命令的使用要比sym简便，它一次可以定义多个符号变量，而且格式简练。因此一般用syms 来创建符号变量。注意各符号变量之间是空格。 例1：比较字符型与符号型 s=123456↙ s = 123456 S1=sym(123456)↙ S1 = 123456 whos↙ Name Size Bytes Class S1 1x1 136 sym object s 1x6 12 char array 说明：从例1可以看到符号变量的保存是不同于矩阵形式的单独保存方式。 符号表达式的创建 创建符号表达式的目的就是把表达式赋值给一符号变量，方便表达式的使用。创建符号表达式由两种方法。 直接创建 例2：fn=sym(2*x^2+b*x+6)↙ fn = 2*x^2+b*x+6 b↙ ?? Undefined function or variable b 说明：符号表达式2*x^2+b*x+6赋给了符号变量fn，但对于符号表达式中的变量x,b并未创建，因此系统并不认识符号b. 间接创建 在创建符号表达式之前，先把符号表达式中的所有变量定义为符号变量，然后直接创建表达式。 例3：syms a b c x↙ f=a*x^2+b*x-c↙ f = a*x^2+b*x-c f+c↙ a*x^2+b*x 符号方程的创建 符号方程与符号表达式不同，表达式只是一个由数字和变量组成的代数式，而方程则是由表达式和等号组成的等式。符号方程的创建方法只有一种。 equ=sym(‘eqution’) 例4：eq1=sym(‘5*x=6+a’)↙ eq1 = 5*x=6+a eq2=sym(‘x^2+y^2+z^2=1’)↙ eq2 = x^2+y^2+z^2=1 符号矩阵的创建 直接创建 利用sym 命令，矩阵元素可以是任何的符号变量、符号表达式及方程，且元素的长度可以不同。 例5：A=sym([sin(x),x+y=2,a;2,6+cos(tan(x)),1/x;abc,y^2=x,x])↙ A = [ sin(x), x+y=2, a] [ 2, 6+cos(tan(x)), 1/x] [ abc, y^2=x, x] 间接创建 在创建符号矩阵之前，先把符号矩阵的所有变量定义为符号变量，而在创建符号矩阵时，只需按创建普通矩阵的格式输入即可完成。 例6：创建例5 的矩阵A. 解:Matlab命令: syms x y a abc↙ e1=sym(x+y=2);↙ e2=sym(y^2=x);↙ A=[sin(x),e1,a;2,6+cos(tan(x)),1/x;abc,y^2,x]↙ A = [ sin(x), x+y=2, a] [ 2, 6+cos(tan(x)), 1/x] [ abc, y^2, x] 由数值矩阵转化为符号矩阵 将一个数值矩阵M转化为符号矩阵S的命令为：S=sym(M) 例7：M=[1,2;3,4];↙ S=sym(M)↙ S = [ 1, 2] [ 3, 4] whos↙ Name Size Bytes Class M 2x2