求matlab或mathematica符号运算高手帮忙完成符号运算.如图,想通过符号运算求两个矩阵的乘积,即jieguo=Tpark*Vg. 得到的结果是除了一次的余弦函数外不含其他三角函数,类似于下图的形式,为简单起

来源:学生作业帮助网 编辑:作业帮 时间:2024/07/19 20:27:21
求matlab或mathematica符号运算高手帮忙完成符号运算.如图,想通过符号运算求两个矩阵的乘积,即jieguo=Tpark*Vg. 得到的结果是除了一次的余弦函数外不含其他三角函数,类似于下图的形式,为简单起
x[oGǿ*R\vgw;Rվ% ]{!>ZK4i!J-{m W` !P=g993YZ^/#1꿬W[k[͠Ph=YISu?-'?o7w~ gGqpx7ݑhr8v[K/wҼ|c~aRTԩN⤗M+|~jĻRd_IDK<5) j: |9'go3MיQ`}n^:Ǯ l K)k.tULrt\e:נ0榥@À8q QX4Bz,e: & rd4·dye$#!h[aC]hk8cBl F4cN UAM 0Y^ @tiC3#Tm:$KgWz͝kqq*Kqt.npX jwgҍG5)DqwhY\HD8ך{Rdͳ5[=;s,"rMH51T Ð1CbyIiqdFI#k珑qst ȦXCP@ʐ,f}0TRѕ! iA,?mۚΤr<#cb$V1:m0D>3@hsVY A1jNhSq_>~ T_]~7_ĥ> ````9lyZ_,aBD/M^N>{뺽ųnbjZ\%TB~l>|vQ/!_jPf>l) TR*#PD +mҵcFN*Jx'M'җ+79\A:hJ^9u]KNo!2/XEdHRe_

求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')