matlab符号运算(一)

基础篇

f = sym('x^2-1');% 定义符号函数
disp(f);
syms x y;        % 定义符号
g = x*2+1-y;
disp(g);
p = g + f;       % 符号函数的运算
disp(p);
disp(factor(f)); % 因式分解

f = sym('(x + 2)^2');
disp(f);
disp(expand(f)); % 将式子展开

f =sym('x + x + x');
disp(f);         % 默认完成了合并同类项的工作

disp(limit(f,'x',3,'right'));
disp(limit(f,x,3,'right'));  % 两个算出来结果是一样的

disp(diff(f,x,1));  % 符号函数实体关于变量X求1阶导数

例题展示：

syms x a n; % 定义符号变量 a x n;
y =  1/(x * log(x)^2) - 1/(x -1)^2;
f1 = limit(y,x,1,'right');  % 取极限
f = a*n^3+(a-1)*n^2;
f2 = symsum(f,n,0,50); % 级数求和 

clc;
clear all;
syms a x;
f = sin(a*x);
dfx = diff(f,x); % 偏导数
dfa = diff(f,a); % 
f1 = x*log(1+x);
int1 = int(f1,x);   % f1 关于x的不定积分
int2 = int(f1,x,0,1); % f1 关于x积分在[0,1]上
% 以下为输出操作
disp(int1);
disp(int2);

常微方程的求解

% 常微方程的求解
f1 = dsolve('Dy + 3*x*y = x * exp(-x/2)','x');
f2 = dsolve('x*Dy + 2*y-exp(x) =0','y(1) = 2*exp(1)','x');
f3 = dsolve('D2y+2*Dy + exp(x) = 0','x');
disp(f1);
disp(f2);
disp(f3);