matlab符号函数

简介: 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);
相关文章
|
1月前
|
索引
matlab--------矩阵重构,重新排列的相关函数说明
matlab--------矩阵重构,重新排列的相关函数说明
77 0
matlab--------矩阵重构,重新排列的相关函数说明
|
1月前
|
Serverless
基于Logistic函数的负荷需求响应(matlab代码)
基于Logistic函数的负荷需求响应(matlab代码)
|
1月前
|
数据安全/隐私保护
地震波功率谱密度函数、功率谱密度曲线,反应谱转功率谱,matlab代码
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
|
1月前
|
调度
知识分享|分段函数线性化及matlab测试
知识分享|分段函数线性化及matlab测试
|
1月前
雷达模糊函数及MATLAB仿真(三)
雷达模糊函数及MATLAB仿真
64 0
|
1月前
雷达模糊函数及MATLAB仿真(二)
雷达模糊函数及MATLAB仿真
80 0
|
1月前
雷达模糊函数及MATLAB仿真(一)
雷达模糊函数及MATLAB仿真
72 0
|
1月前
|
机器学习/深度学习 算法
基于自适应支持向量机的matlab建模与仿真,不使用matlab的SVM工具箱函数
基于自适应支持向量机的matlab建模与仿真,不使用matlab的SVM工具箱函数
|
1月前
|
算法
【数学建模竞赛】Matlab逻辑规则,结构基础及函数
【数学建模竞赛】Matlab逻辑规则,结构基础及函数
34 0