matlab符号函数

一、符号函数简介

符号函数是matlab的一个重要功能，可以用来表示数学上的函数，同时也可以用来进行数值计算。
创建符号函数的方法：

• 利用syms定义符号变量生成符号函数
• 利用sym和@生成符号函数
• 利用函数文件生成符号函数

1.1 syms方法

首先利用syms命令定义自变量，然后再生成符号函数

syms sigamn m ft sigam3;
sigamn = (12*m*ft + (12*m+5)*sigam3) / (24*m -1); % sigamn 就是一个符号函数

1.2 利用sym和@方法

首先使用@声明变量，在其后写符号函数表达式，之后使用sym转换成符号函数，就可以进行其他符号运算。

y = @(t) sin(t) - cos(t) + exp(t);
ys = sym(y);

y1 = @(a,b) cos(a) + cos(b);
ys1 = sym(y1)；

1.3 函数文件法

利用函数文件来生成符号函数。(实际上，也是用1.1或者1.2的方法实现)

function output = xpower(n)
    syms x
    output = x^n;
end

二、符号函数操作

2.1 符号函数的数值运算

对于符号函数，可以计算符号函数在某些点处的函数值，此时需要使用matlabFunction命令做一次转换。

syms x;
y = sin(x) - cos(x) + exp(x);
y_f = matlabFunction(y1);
x0 = [0, pi];
y_v = y_f(x0); % 计算符号函数y_f在x0处的函数值，x0可以是单个值，也可以是数组 

2.2 符号函数化简

• 合并同类项
  使用collect(f)函数来合并多项式中相同的项。

  syms x t;
  f=(1+x)*t+x*t;
  collect(f)
  

• 因式分解
  使用factor(f)对多项式进行因式分解，如果无法在有理数的范围内作分解, 那么返回的结果还是输入值。

  syms f;
  f=x^3-6*x^2+11*x-6;
  factor(f);
  

• 多项式展开
  使用expand(f)实现多项式展开。

  syms x
  f=x*(x*(x-1)+3)+2;
  expand(f);
  

• 符号函数化简
  使用simplify(f)函数和simple(f)函数对符号函数进行化简

  syms a b c;   %定义为符号变量
  r = 2;
  ss = a^r/((a-b)*(a-c))+b^r/((b-c)*(b-a))+c^r/((c-a)*(c-b)); %多项式公式
  simplify(ss); %%简化多项式
  simple(ss);
  % simple是通过对表达式尝试多种不同的方法（包括simplify）进行化简
  % 以寻求符号表达式的最简形式。
  

2.3 方程求解

solve()函数可以求解代数方程(多项式方程)的符号解析解。

syms x;
y = x*sin(x)-x; 
solve(y, x);  % 只有一个自变量，可以省略自变量，同solve(y)
solve(y==0);

2.4 符号函数求导

2.4.1 一元符号函数求导

可以使用命令diff对符号函数进行求导，并且当diff作用到向量或者矩阵使可以求差分。

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

diff(y,1); % 求1阶导，同diff(y)
diff(y,2); % 求2阶导

当diff作用于符号函数上时，可以对符号函数进行求导。

• 第一个参数传入符号函数
• 第二个参数传入要求几阶导，默认求1阶导。

2.4.2 多元符号函数求导

• 第一个参数为待求导符号函数
• 第二个参数为要对哪个自变量求偏导
• 第三个参数为要求几阶导数
  当连续对自变量求偏导时，放在前面的自变量先求偏导，后面的自变量后求偏导。

syms x1 x2 x3;
y = x1^5*x2 + x2*x3 - x1^2*x3;
py1 = diff(y,x1,1); % 对x1求一阶偏导
py2 = diff(y,x1,2); % 对x1求二阶偏导
py3 = diff(y,x1,x2); % 先对x1求偏导，再对x2求偏导
py4 = diff(y,x2,x1); % 先对x2求偏导，再对x1求偏导

如果diff函数作用的对象不是符号函数，而是矩阵，那么对应的功能就是求差分。

2.5 符号函数定积分

使用int()函数，可以对符号函数进行定积分计算。

% 使用形式： 
int(被积函数,积分变量，积分下限，积分上限); 
% 如果只有一个参量，可以省略第2项

syms x b;
f = b + 2*x;
int_f = int(f,x,-6,6);