MATLAB基础(四)

% 符号计算基础
%     符号常量和符号变量
%     符号四则运算
%     符号表达式的化简
%     符号矩阵的运算
% 符号函数
%     极限
%     微分
%     积分
% 符号级数
% 符号方程求解
%     代数方程
%     微分方程
%     
% 符号变量

a = sym('a');
sym b;
% 符号常量
c = sym('3');
% 符号表达式
f = '3*x+6';
syms x;
f2 = 3*x+6;

f3 = sym('3*y+4');
% 符号运算直接可以使用加减乘除
% 符号 四则运算 + - / *
fadd1 = f2+f3;
fmul = f1*f2;

% 符号表达式化简
syms x y
s = (x^2+y^2)^2+(x^2-y^2)^2;
simplify(s)

% 符号表达式和数值的计算
eval(c) %将符号变量转为数值
numeric(c) % 不存在了
% 符号表达式的因式分解和展开
% factor% 因式分解
% expand% 合并同类项
% collect% 展开

% 因式分解,展开和合并同类项
syms a b x y;
f1 = a^3-b^3;
factor(f1)

f2 = (3*x^2+ 8*y^2)*(-x^2+3*y);
expand(f2)

f3 = 3*x^2 + 4*x^2 + 5*x^2*y;
collect(f3)
% 
% 符号矩阵
%     构建方式和数值矩阵一致
%     转置和行列式
%     transpose
%     determ
%     其他使用在数值矩阵中的函数也可以直接用于符号矩阵
%     
a1 = [x x+y; y y^2];a1'
% ans =
% [           conj(x),   conj(y)]
% [ conj(x) + conj(y), conj(y)^2]
% 符号函数的求解
syms x
f1 = x^3-9;
subs(f1,3)
% 符号极限,符号微分,符号积分
syms x a;
y = (sin(x+a)-sin(x-a))/x;
limit(y,0)

y2 = sqrt(1+exp(x));
diff(y2)
diff(y2,2)
diff(y2,3)

% 积分
y3 = (3-x^2)^3;
int(y3)

y4 = abs(1-x);
int(y4, 1,2)

% 符号级数求和
syms n
f = 1/n^2;
s1 = symsum(f,n ,1, inf)

% 泰勒级数
syms x
y = (1 + x + x^2)/(1 - x+ x^2);
taylor(y,6,1)

% 代数方程
clear
syms x
solve(x +x*exp(x) -10)

eval(solve('x+x*exp(x)-10'))

% 方程组
[x y] = solve('x+y -98','x^(1/3)+y^(1/3)','x,y')
[x y] = solve('x^3+y^3-28','1/x+1/y-4','x,y')

% 微分方程
dsolve('Dy - (x^2+y^2)/x^2/2','x')