求matlab或mathematica符号运算高手帮忙完成符号运算.

1个回答

  • 下面的代码,在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')