matlab符号运算结果的分数转化为小数vpa/digits

1、我们以下面的代码为例：

clc

syms a1 a2 a3 real

Rx=@(x)[1 0 0;

    0 cos(x) sin(x);

    0 -sin(x) cos(x)];

Ry=@(x)[cos(x) 0 -sin(x);

    0 1 0;

    sin(x) 0 cos(x)];

Rz=@(x)[cos(x) sin(x) 0;

    -sin(x) cos(x) 0;

    0 0 1];

R121=Rx(a3)*Ry(a2)*Rx(a1);

R121_0=subs(R121,a2,pi)

其执行结果如图。可以看出结果中的一部分运算结果以分数的形式显示了。

2、下面用vpa函数对R121_0进行处理，在命令窗口中输入：

vpa(R121_0)

运行结果如图。可以看出已经之前的分数被转换成了小数，但是小数位数过于冗长，那么怎么控制显示位数呢？

3、下面我们在执行代码：

vpa(R121_0,4) %4就是表示要保留的位数

可以看到结果的小数保留了4位，看起来明显好多了。

4、除了使用vpa(R121_0,4)来控制小数位数外，还可以使用digits函数来修改所有符号运算结果的小数保留位数，执行下面的代码：

d_pre=digits(4) %设置显示位数，返回设置之前的位数

d_cur=digits %查询当前的显示位数

vpa(R121_0)

运算结果如图。可知，第二步中之所以显示的那么长的小数，是因为符号运算默认的显示小数位数为32位，通过digits(4)后把符号运算的显示小数位数改为了4位，再执行vpa(R121_0)就与第三步的结果一样了。