求matlab或mathematica符号运算高手帮忙完成符号运算.
求matlab或mathematica符号运算高手帮忙完成符号运算.
如图,想通过符号运算求两个矩阵的乘积,即jieguo=Tpark*Vg.
![](https://img.yulucn.com/upload/3/70/3706f9dc6fe4dac435f3c069862dda66_thumb.jpg)
得到的结果是除了一次的余弦函数外不含其他三角函数,类似于下图的形式,为简单起见,先只求Tpark的第二行与Vg相乘.矩阵已经写好了.我很需要计算这个,
![](https://img.yulucn.com/upload/d/b4/db4deb178e4951a2c421ddfdd1d44731_thumb.jpg)
syms q b r v1 v5 v7 v11 v13 v17 v19;
Vg = [v1*cos(q)-v5*cos(5*q)+v7*cos(7*q)-v11*cos(11*q)+v13*cos(13*q)-v17*cos(17*q)+v19*cos(19*q);
(1+b)*(v1*cos(q-2*pi/3)-v5*cos(5*(q-2*pi/3))+v7*cos(7*(q-2*pi/3))-v11*cos(11*(q-2*pi/3))+v13*cos(13*(q-2*pi/3))-v17*cos(17*(q-2*pi/3))+v19*cos(19*(q-2*pi/3)));
(1+r)*(v1*cos(q+2*pi/3)-v5*cos(5*(q+2*pi/3))+v7*cos(7*(q+2*pi/3))-v11*cos(11*(q+2*pi/3))+v13*cos(13*(q+2*pi/3))-v17*cos(17*(q+2*pi/3))+v19*cos(19*(q+2*pi/3)));]
T = sqrt(1.5)*[-sin(q) -sin(q-2*pi/3) -sin(q+2*pi/3)];
dq=T*Vg