源代码见文末
部分源代码：
1 % 符号变量 两种表达方式
2 a=sym('a');
3 class(a);
4
5 syms b;
6 b;
7
8 % 符号常量
9 c=sym('3');
10 c;
11 % 符号表达式 三种表达方式
12 f1='3*x+4'; % char类型，但是可以参与运算。所以这种形式是可以直接用的，比较方便
13
14 syms x
15 f2=3*x+4; % sym类型
16
17 f3=sym('3*y+4'); % sym类型
18
19 % 符号四则运算
20 fadd1 = f1 + f2;
21 fmu1 = f1*f2;
22
23 % 符号表达式的化解
24 syms x y
25 s = (x^2+y^2)^2 + (x^2-y^2)^2;
26 simplify(s);
27
28 % 符号表达式和数值的转换
29 eval(c);
30 syms; % 上面已经涉及到
31
32 % 因式分解、展开与合并同类项
33
34 syms a b x y
35 f1=a^3-b^3;
36 factor(f1);
37
38 f2=(3*x^2+8*y^2)*(-x^2+3*y);
39 expand(f2)
40
41 f3 = 2*x+3*x+4*x+y;
42 collect(f3)
43
44 % 符号矩阵
45 a1=[x x+y;y x+y];
46 a1
47 transpose(a1) % 转置
48
49 % 符号函数值的求解
50 syms x
51 f1 = x^3 -9;
52 subs(f1,3)
53
54 % 符号极限、符号微分、符号积分
55 syms x
56 y=sin(x);
57 limit(y,0)
58
59 y1 = x^2;
60 diff(y1) % 求导、微分
61 diff(y1,2) % 二阶导数
62
63 y3=(3-x^2)^3;
64 int(y3) % 不定积分
65
66 y4=abs(1-x);
67 int(y4,1,2) % 定积分1-2范围内
68
69 %%符号方程求解
70 %代数方程
71 clear
72 syms x
73 solve(x+x*exp(x)-10);
74
75 solve('x+x*exp(x)-10') % 和前面的符号表达式方法联系起来了
76 eval(solve('x+x*exp(x)-10')) % 这是数值的表达结果
77
78 %方程组
79 [x,y] = solve('x+y=10','x-y=2','x,y') % 输出x和y的值
80
81 % 符号的微分方程
82 dsolve('Dy-(x^2+y^2)/x^2/2','x') % 微分方程组合上面的方程组道理一样