matlab符号函数

一、符号函数简介

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

创建符号函数的方法:

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

1.1 syms方法

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

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

1.2 利用sym和@方法

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

matlab 复制代码
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的方法实现)

matlab 复制代码
function output = xpower(n)
    syms x
    output = x^n;
end

二、符号函数操作

2.1 符号函数的数值运算

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

matlab 复制代码
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)函数来合并多项式中相同的项。

    matlab 复制代码
    syms x t;
    f=(1+x)*t+x*t;
    collect(f)
  • 因式分解
    使用factor(f)对多项式进行因式分解,如果无法在有理数的范围内作分解, 那么返回的结果还是输入值。

    matlab 复制代码
    syms f;
    f=x^3-6*x^2+11*x-6;
    factor(f);
  • 多项式展开
    使用expand(f)实现多项式展开。

    matlab 复制代码
    syms x
    f=x*(x*(x-1)+3)+2;
    expand(f);
  • 符号函数化简
    使用simplify(f)函数和simple(f)函数对符号函数进行化简

    matlab 复制代码
    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()函数可以求解代数方程(多项式方程)的符号解析解。

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

2.4 符号函数求导

2.4.1 一元符号函数求导

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

matlab 复制代码
syms x;
y = x^4 - 5*x^2 +6;

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

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

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

2.4.2 多元符号函数求导

  • 第一个参数为待求导符号函数
  • 第二个参数为要对哪个自变量求偏导
  • 第三个参数为要求几阶导数
    当连续对自变量求偏导时,放在前面的自变量先求偏导,后面的自变量后求偏导。
matlab 复制代码
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()函数,可以对符号函数进行定积分计算。

matlab 复制代码
% 使用形式: 
int(被积函数,积分变量,积分下限,积分上限); 
% 如果只有一个参量,可以省略第2项
matlab 复制代码
syms x b;
f = b + 2*x;
int_f = int(f,x,-6,6);
相关推荐
缓风浪起38 分钟前
【力扣】2011. 执行操作后的变量值
算法·leetcode·职场和发展
gsfl42 分钟前
双指针算法
算法·双指针
郝学胜-神的一滴42 分钟前
矩阵的奇异值分解(SVD)及其在计算机图形学中的应用
程序人生·线性代数·算法·矩阵·图形渲染
万粉变现经纪人3 小时前
如何解决 pip install -r requirements.txt 私有索引未设为 trusted-host 导致拒绝 问题
开发语言·python·scrapy·flask·beautifulsoup·pandas·pip
qq_479875433 小时前
C++ std::Set<std::pair>
开发语言·c++
云知谷5 小时前
【C++基本功】C++适合做什么,哪些领域适合哪些领域不适合?
c语言·开发语言·c++·人工智能·团队开发
电子_咸鱼6 小时前
LeetCode——Hot 100【电话号码的字母组合】
数据结构·算法·leetcode·链表·职场和发展·贪心算法·深度优先
仰泳的熊猫6 小时前
LeetCode:785. 判断二分图
数据结构·c++·算法·leetcode
rit84324996 小时前
基于MATLAB实现基于距离的离群点检测算法
人工智能·算法·matlab
l1t6 小时前
DeepSeek辅助利用搬移底层xml实现快速编辑xlsx文件的python程序
xml·开发语言·python·xlsx