Matlab符号计算(1)

一、建立符号对象

1. 建立符号变量和符号常量

sym()函数用来建立单个符号量，调用格式为：
符号量名 = sym('符号字符串')；
该函数可以建立一个符号量，符号字符串可以是常量、变量、函数或表达式，注意：该方法只能建立一个符号变量。

一次性创建多个符号变量的方法：
syms 符号变量明1 符号变量名2 ··· 符号变量明n
注意：变量名之间用空格分割，不要用逗号等其他字符。正确使用方法如下：
syms a b c d
该方法可以一次性创建多个变量名。
PS:符号常量的计算结果是精确的数学表达式，而数值常量计算结果是约分后的常数。

2.建立符号表达式

1. 单引号生成符号表达式

y = '1/sqrt(2*x)'; 
g = 'cos(x^2)-sin(x)=0';%符号方程

2. sym()函数生成符号表达式

Y = sym('3*x');
G = sym('[a,b;c,d]');%矩阵表达式

3. 使用已经定义的符号变量组成符号表达式

syms x y
V = 3*x^2-5*y+2*x*y+6;

二、符号表达式的计算

1. 四则运算

• 使用函数实现

matlab
%其中f,g为符号表达式
symadd(f,g);    %加法
symsub(f,g);    %减法
symmul(f,g);    %乘法
symdiv(f,g);    %除法
sympow(f,g);    %幂运算

• 直接使用"+、-、*、/、^"运算符实现。
  PS：Matlab对运算之后的结果不一定会化到最简形式，需要额外增添一步化简操作。

三、因式分解与展开

factor(s);   				%对符号表示式分解分解因式
expands(s);  				%对s进行展开
collect(s);  				%对s合并同类项
collect(s,v);				%对s按变量v合并同类项。

3. 符号表达式的化简

simplify(s); 				%应用函数规则对s进行化简
simple(s);   				%调用MATLAB的其他函数对s进行综合化简，并显示化简过程。

4. 符号表达式与数值之间的转换

sym(1.5);						%数值表达式转换为符号表达式
numeric('sqrt(5)');	%符号表达式转换为数值表达式
eval('sqrt(5)');		%符号表达式转换为数值表达式

四、符号表达式变量的确定

findsym(s,n); %返回符号表达式s中的n个符号量，若没有指定n，则返回s中全部的符号变量。

##四、符号矩阵
使用sym()函数可以建立符号矩阵并简化
m = sym('[1/(a+x),1;2,1/(b+x)]');
对一般矩阵操作的函数同样适用于符号矩阵

transpose(s);       %返回s矩阵的转置矩阵
determ(s);          %返回s矩阵的行列式值
diag(s);            %以矩阵s的元素作为矩阵X的主对角线元素
triu(s);            %返回矩阵s上三角矩阵
tril(s);            %返回矩阵s下三角矩阵
inv(s);             %返回矩阵s的逆矩阵
det(s);             %返回矩阵s的行列式的值
rank(s);            %返回矩阵的秩
eig(s);             %返回矩阵的特征值和特征向量

五、符号表达式转化为数值

通常采用subs()函数实现符号表达式中变量的替换，用法如下：

subs（f,var,value）;

将符号表达式f中的变量组var替换成数值value，但需要注意，替换之后变量类型仍为sym。通常需要将subs函数和vpa,eval,double函数结合使用，如果单独使用subs()函数，只是将变量替换成了符号常量，但类型仍旧是sym，因此需要外加vpa、eval或double函数，进行类型之间的转换。

syms x y
z = x^2+y^2-2*x*y;
subs(z,[x y],[1,2]);
type1 = class(subs(z,[x y],[1,2]));
double(subs(z,[x y],[1,2]));
type2 = class(double(subs(z,[x y],[1,2])));

运行后如下：