声明:本文章中数据来自清风老师数学建模课程
文章目录
MATLAB符号变量的创建和简单运算1、符号变量1. 1 符号变量的创建1.2 符号方程的创建3 符号矩阵的创建 2、符号运算2.1 简单运算2.2 表达式的整理2.3 因式分解2.4 多项式展开2.5 合并2.6 计算分子与分母2.7 让结果显示的更加自然MATLAB符号变量的创建和简单运算
1、符号变量
1. 1 符号变量的创建
使用syms创建符号变量
syms x % 创建一个符号变量,sym是symbolic的缩写syms a b c% 创建多个符号变量
1.2 符号方程的创建
方法一:先创建符号变量后创建符号方程
syms a x% 创建符号变量y = a*x+x^2% 创建符号方程
方法二:str2sym命令创建,在MATLAB b后可以使用
y = str2sym('a*x+x^2')
创建符号方程后,y的类型在工作区中观察到为1x1sym
3 符号矩阵的创建
%% 符号矩阵的创建syms alphaM = [cos(alpha) -sin(alpha);sin(alpha) cos(alpha)]
2、符号运算
了解如何符号变量后,我们可以使用符号变量进行一些简单的运算、多项式运算、分式运算等等
2.1 简单运算
%% 简单运算syms a b c d ey = a + b% a + bx = c - d% c - dy1 = x*y% (a + b)*(c - d)y2 = y1/y% c - dy3 = y1^3% (a + b)^3*(c - d)^3y4 = sqrt(y3) % ((a + b)^3*(c - d)^3)^(1/2)y5 = exp(y4)% exp(((a + b)^3*(c - d)^3)^(1/2))
2.2 表达式的整理
有时候我们的表达式是可以继续通分化简的,这时我们可以使用simplify命令来进行通分化简。如:
%% 符号表达式的整理clear;clc% 化简syms ay=(cot(a/2)-tan(a/2))*(1+tan(a)*tan(a/2))simplify(y)% 2/sin(a)
2.3 因式分解
使用MATLAB中的factor命令可以对整数或者多项式进行因式分解
%% 因式分解factor(12) % 对常数进行因式分解 %223syms m n xy = -24*m^2*x-16*n^2*xfactor(y)% [ -8, x, 3*m^2 + 2*n^2]y1=m^3-n^3factor(y1)% [ m - n, m^2 + m*n + n^2]% 立方差公式的因式分解 y2 = m^2-n^2factor(y2)% [ m - n, m + n]% 平方差公式的因式分解
2.4 多项式展开
使用expand将因式展开为多项式
%% 多项式展开syms a xy = a*(x^2-a)^2+(x-2)expand(y)% a^3 - 2*a^2*x^2 + a*x^4 + x - 2
2.5 合并
使用collect命令可以进行类似合并同类项的操作
%% 合并syms x yz = (x+y)^2*y+5*y*x-2*x^3% expand(z) - 2*x^3 + x^2*y + 2*x*y^2 + 5*x*y + y^3collect(z,x)% y*x^2 - 2*x^3 + (2*y^2 + 5*y)*x + y^3collect(z,y)% y^3 + 2*x*y^2 + (x^2 + 5*x)*y - 2*x^3
2.6 计算分子与分母
使用numden来计算分式的分子和分母
返回值为两个[z1,z2],其中z1为分子,z2为分母
%% 计算分子与分母% [z1,z2] = numden(2.5) % 会报错,因为numden的输入变量不能是数值,只能是符号变量% ans = sym(2.5); % sym函数可以将数值2.5转换为符号[z1,z2] = numden(sym(2.5)) % 对常数计算分子与分母% z1 = 5% z2 = 2syms x yz = 1/x*y+x/(x^2-2*y)[z1,z2] = numden(z) %z1分子,z2分母% z1 = - x^2*y - x^2 + 2*y^2% z2 = x*(- x^2 + 2*y)
2.7 让结果显示的更加自然
看下面命令
syms x yM = (1/x*y+x/(x^2-2*y)-x^2/(3+y)^2)^2;expand(M)
通过expand命令对M进行展开,结果为:
y^2/x^2 + x^4/(y^4 + 12*y^3 + 54*y^2 + 108*y + 81) + (2*x^3)/(- x^2*y^2 - 6*x^2*y - 9*x^2 + 2*y^3 + 12*y^2 + 18*y) - (2*y)/(- x^2 + 2*y) + x^2/(x^4 - 4*x^2*y + 4*y^2) - (2*x*y)/(y^2 + 6*y + 9)
看到这样的结果你是不是非常绝望呢??????????????????
如果我们要把这样的结果转化为人们习惯的书写方式放到论文中,是不是会崩溃呢??????
其实我们可以利用MATLAB提供的mupad工具箱(一个商用计算机代数系统,“Multi Processing Algebra Data Tool”的缩写)来帮助我们
方法:在MATLAB中的命令行中输入mupad,回车打开工具箱,然后将上面M的展开结果复制到里面,回车就OK了
当然你可能会遇到这样的情况:
这是说明在后序的发行版本中可能会将该工具箱移除,建议我们使用通过新建实时脚本来解决。
创建实时脚本
将M的展开结果粘贴过去,运行
如果觉得《MATLAB符号变量的创建和简单运算》对你有帮助,请点赞、收藏,并留下你的观点哦!