符号对象的生成

最新推荐文章于 2022-07-16 23:14:35 发布

seamanj

最新推荐文章于 2022-07-16 23:14:35 发布

阅读量1.1k

点赞数

分类专栏： MATLAB

本文链接：https://blog.csdn.net/seamanj/article/details/32732685

版权

MATLAB 专栏收录该内容

108 篇文章 19 订阅

订阅专栏

在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为符号矩阵