在Matlab中,符号计算可以用推理解析的方式进行,避免数值计算带来的截断误差,同时符号计算可以得到正确的封闭解或者正确的数值解。
首先搞清楚符号对象的相关概念。
符号对象包括符号变量、符号常量、符号表达式、符号向量、符号矩阵。
符号变量:符号变量就是用于符号运算的变量,例如分解一个关于x的多项式,那么x就是一个符号变量(因为他没有具体的数值),或者一个类似于[x,1,8,7]这样的向量也是符号变量的一种,符号变量一般都用于符号运算,如分解因式,微积分一类的运算
符号常量:符号常量就是可以确定值的符号变量
符号表达式:含有符号常量或者符号变量的表达式就称符号表达式。符号表达式包括符号函数和符号方程。其中符号函数没有等号,而符号方程必须带有符号。可见符号表达式包含了符号常量和符号变量的范围。
符号向量、符号矩阵:向量(矩阵)中的元素都是符号表达式即为符号向量(符号矩阵)。
符号对象的生成
1.符号变量的生成
可以通过sym和syms函数来生成符号变量
通过sym来生成符号变量的格式为:
符号变量名=sym('变量名')
syms函数的一般调用格式为:
syms var_1 var_2 ... var_n
用这种格式定义符号变量时不要在变量名上加字符分界符('),变量间用空格而不要用逗号分隔。
>> x = sym('x')
x =
x
>> syms y z
>> y
y =
y
>> z
z =
z
>> help syms
SYMS Short-cut for constructing symbolic objects.
SYMS arg1 arg2 ...
is short-hand notation for
arg1 = sym('arg1');
arg2 = sym('arg2'); ...
SYMS arg1 arg2 ... real
is short-hand notation for
arg1 = sym('arg1','real');
arg2 = sym('arg2','real'); ...
SYMS arg1 arg2 ... positive
is short-hand notation for
arg1 = sym('arg1','positive');
arg2 = sym('arg2','positive'); ...
SYMS arg1 arg2 ... unreal
is short-hand notation for
arg1 = sym('arg1','unreal');
arg2 = sym('arg2','unreal'); ...
Each input argument must begin with a letter and must contain only
alphanumeric characters.
By itself, SYMS lists the symbolic objects in the workspace.
Examples:
syms x beta real
is equivalent to:
x = sym('x','real');
beta = sym('beta','real');
syms k positive
is equivalent to:
k = sym('k','positive');
To clear the symbolic objects x and beta of 'real' or 'positive' status, type
syms x beta unreal
See also sym.
Reference page in Help browser
doc syms
>>
考察一下符号变量与数值变量的区别
>> a = sym('a');b = sym('b');c = sym('c');d = sym('d'); %define 4 symbolic variables,can also write it as 'syms a b c d' .
>> w = 10; x = 5; y = -8; z = 11;%define 4 numerical variables
>> A = [a,b;c,d];%construct symbolic matrix A
>> B = [w,x;y,z];%construct numerical matrix B
>> det(A)
ans =
a*d-b*c
>> det(B)
ans =
150
>>
2.符号常量的生成
符号常量的生成通过sym函数生成,调用格式如下:
符号常量名=sym('常量') 或者 符号常量名=sym('常量表达式')
>> a = sym('pi')
a =
pi
>> f = sym('2-3*6')
f =
2-3*6
>>
>> pi1 = sym('pi');k1 = sym('8');k2 = sym('2');k3 = sym('3'); %define symbolic constants
>> pi2 = pi; r1 = 8; r2 = 2; r3 = 3;%define numerical constants
>> sin(pi1/3)
ans =
1/2*3^(1/2)
>> sin(pi2/3)
ans =
0.8660
>> sqrt(k1)
ans =
2*2^(1/2)
>> sqrt(r1)
ans =
2.8284
>> sqrt(k3+sqrt(k2))
ans =
(3+2^(1/2))^(1/2)
>> sqrt(r3+sqrt(r2))
ans =
2.1010
>>
3.符号表达式的生成
符号表达式的生成有三种方法:
(1)用sym函数建立符号表达式
(2)使用已经定义的符号量(变量、常量)组成符号表达式。
>> syms x
>> y1 = '1/sqrt(2*x)'
y1 =
1/sqrt(2*x)
>> y2 = sym('1/sqrt(2*t)')
y2 =
1/sqrt(2*t)
>> y3 = 1/sqrt(2*x)
y3 =
1/2*2^(1/2)/x^(1/2)
注意y1并不是符号表达式,而是字符串
4.符号矩阵的生成
>> syms x
>> A = [x 1;1 1]
A =
[ x, 1]
[ 1, 1]
>> B = sym('[y 1;1 1]')
B =
[ y, 1]
[ 1, 1]
>> y
??? Undefined function or variable 'y'.
>> C = [y 1;1 1]
??? Undefined function or variable 'y'.
>>
这里A,B为符号矩阵