用Matlab解MTPA的dq轴电流的方程

最大转矩电流比的方程非常难解，于是尝试用matlab做一个解方程的代码，具体的过程就不在各位大神面前班门弄斧了，主要说一下函数。如果使用solve函数是解不出来的，可以尝试用vpasolve函数（不过旧版本可能用不了）

MTPA的公式

Matlab的代码

%计算isq
y= vpasolve( 2 * flux * y *(1 + sqrt( 1 + 4 *( y^2 *( lq - ld ) ^2)/ flux^2)) == te * 2/3,y);
%计算isd
x= vpasolve( 4 * sqrt( x * ( x * ( lq - ld ) - flux)^3 / ( lq - ld ))== te * 2/3 ,x);

本人是一菜鸟，希望各位大神不吝赐教。如有错误，欢迎指正，谢谢！