matlab中符号运算的学习

1.符号变量定义

%符号变量
syms a b c 
% 符号方程
y=str2sym('a*x+x^2')
%符号矩阵
syms alpha
M=[cos(alpha)  -sin(alpha)]
% 运算 加减乘除

%符号表达式的化简
simplify(y)

%因式分解,可以对常数，也可以对表达式
factor(12)
syms m n 
y=m^ 3  - n^3
factor(y)

%多项式展开
expend(y)
%多项式合并
collect(z,x)  %合并的式子，和变量指定

%对常数计算分子和分母
[z1,z2]=numden(sym(2,5))
%还可以得到符号变量的分子和分母

同时为了使得结果显示更加自然，可以采用mupad工具箱和实时脚本功能，对于版本比较高的matlab。

2.求导和差分运算

%一元函数求导
syms x 
y=x^4-3*x^2
dy=diff(y)
dy=diff(y,2)
simplify(y)

%多元函数的导数
syms x1 x2 x3
y1=x1^5*x2+x2*x3
py1=diff(y1,x1,1)
py2=diff(y1,x1,2)
py3=diff(y1,x1,x2)
py4=diff(y1,x2,x1)


%如果diff函数作用的对象不是函数是矩阵就是做差分。
%对于行向量是后面减前面做差
A=[ 1 2 6]
diff(A)  %[ 1 4]
diff(A,2) % [3]

%对于矩阵是下一行减去上一行求二阶差分
A=[ 1 2 ; 3  4]
diff(A)  %[ 2 2]
diff(A,2) %同上
diff(A,2,1)  %最后的1表示在上面的基础上再进行差分
diff(A,1,2)   %下一列减去上一列求一阶差分，2表示在列上进行差分

3.不定积分和定积分

syms x
y=x^2
int(y,x) 

y=1/x
int(y,x)  %区分，这里不会有绝对值符号

%定积分的求解，加入上下界范围
syms x a b
y = exp(x)
int(y,x,a,b)

%注意不是所有函数都可以求解出一个精确值
%这里就可以用计算数值积分,计算近似解
syms x
y=@(x) 1./exp(x).*log(x+2.*x+sin(x))
integral(y,0.4)

4.求解方程和方程组

注意不同版本可能会有不兼容的问题，所以如果发现问题，可以取查找对应帮助文档的内容

%求解单变量
syms x
eqn =(sin(x)==1);
answ=solve(eqn,x)
[answ,params,condions]=solove(eqn,x,'ReturnConditions',true)
因为三角函数是周期函数，如果要得到所有解，需要加上条件

%多变量方程求解
clear;clc
syms a b c x
eqn =(a*x^2+b*x+c==0);  %注意这里在写等式的时候一定要是两个等于号，一个是赋值是不可以的
answ1=solve(eqn,x)
answ2=solve(eqn,a) %分别把x和a当成自变量

%方程组求解,用向量表示
eqn=[2*u+v==a,u-v==1];
answ=solve(eqn,[u,v])
answ.u
answ.v
[answ_u,answ_v] = solve(eqn,[u,v]);

%产生警告的原因
%求解符号表达式是比较困难的，对于很多方程，只能近似求解。
%更加强大的功能  vpasolve 函数求解
输入变量，方程，变量，开始点，搜索区间
syms x 
eqn = sin(x)==x^2-1;
vpasolve(eqn,x,[0,2])
vpasolve(eqn,x,[-1 0])
vpasolve(eqn,x,'random',true)
vpasolve(eqn,x,-5)

%复杂函数用fsolve求解
x0=[0,0]
result_x=fsolve(@my_fun,x0)