求matlab或mathematica符号运算高手帮忙完成符号运算.如图,想通过符号运算求两个矩阵的乘积,即jieguo=Tpark*Vg. 得到的结果是除了一次的余弦函数外不含其他三角函数,类似于下图的形式,为简单起
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/27 17:57:08
![求matlab或mathematica符号运算高手帮忙完成符号运算.如图,想通过符号运算求两个矩阵的乘积,即jieguo=Tpark*Vg. 得到的结果是除了一次的余弦函数外不含其他三角函数,类似于下图的形式,为简单起](/uploads/image/z/6836871-39-1.jpg?t=%E6%B1%82matlab%E6%88%96mathematica%E7%AC%A6%E5%8F%B7%E8%BF%90%E7%AE%97%E9%AB%98%E6%89%8B%E5%B8%AE%E5%BF%99%E5%AE%8C%E6%88%90%E7%AC%A6%E5%8F%B7%E8%BF%90%E7%AE%97.%E5%A6%82%E5%9B%BE%2C%E6%83%B3%E9%80%9A%E8%BF%87%E7%AC%A6%E5%8F%B7%E8%BF%90%E7%AE%97%E6%B1%82%E4%B8%A4%E4%B8%AA%E7%9F%A9%E9%98%B5%E7%9A%84%E4%B9%98%E7%A7%AF%2C%E5%8D%B3jieguo%3DTpark%2AVg.+%E5%BE%97%E5%88%B0%E7%9A%84%E7%BB%93%E6%9E%9C%E6%98%AF%E9%99%A4%E4%BA%86%E4%B8%80%E6%AC%A1%E7%9A%84%E4%BD%99%E5%BC%A6%E5%87%BD%E6%95%B0%E5%A4%96%E4%B8%8D%E5%90%AB%E5%85%B6%E4%BB%96%E4%B8%89%E8%A7%92%E5%87%BD%E6%95%B0%2C%E7%B1%BB%E4%BC%BC%E4%BA%8E%E4%B8%8B%E5%9B%BE%E7%9A%84%E5%BD%A2%E5%BC%8F%2C%E4%B8%BA%E7%AE%80%E5%8D%95%E8%B5%B7)
求matlab或mathematica符号运算高手帮忙完成符号运算.如图,想通过符号运算求两个矩阵的乘积,即jieguo=Tpark*Vg. 得到的结果是除了一次的余弦函数外不含其他三角函数,类似于下图的形式,为简单起
求matlab或mathematica符号运算高手帮忙完成符号运算.
如图,想通过符号运算求两个矩阵的乘积,即jieguo=Tpark*Vg.
得到的结果是除了一次的余弦函数外不含其他三角函数,类似于下图的形式,为简单起见,先只求Tpark的第二行与Vg相乘.矩阵已经写好了.我很需要计算这个,
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
求matlab或mathematica符号运算高手帮忙完成符号运算.如图,想通过符号运算求两个矩阵的乘积,即jieguo=Tpark*Vg. 得到的结果是除了一次的余弦函数外不含其他三角函数,类似于下图的形式,为简单起
下面的代码,在matlab2014a可运行.
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)*[cos(q+pi/2) cos(q-pi/6) cos(q+7*pi/6)]; % 三角转换 -sin(q) = cos(q+pi/2)
dq = T*Vg
dq1 = expand(dq,'ArithmeticOnly',true)
dq2 = combine(dq1,'sincos')