我的程序如下:有13个变量,13个方程
syms id iq udc igd igq ed eq wk mdc mid miq r a
f1=('((1/0.0054)-udc/19.44)*ed-725*id/27+(wk-udc*pi/36)*iq+(udc/19.44)*(3600-udc+10*mdc-2*id+50*mid)=0');
f2=('((1/0.0054)-udc/19.44)*eq-725*iq/27+(udc*pi/36-wk)*id-(udc/19.44)*(2*iq-50*miq)=0');
f3=('(id/64.8)*(ed+iq*0.54*pi-(3600-udc+10*mdc-2*id+50*mid))+(iq/64.8)*(eq-id*0.54*pi+2*iq-50*miq)-(udc/0.225)=0');
f4=('1376529-(cos(a-r))*ed*500-(sin(a-r))*eq*500-0.5*igd+100*pi*igq=0');
f5=('(-ed*500*sin(r-a))-(cos(a-r))*eq*500-0.5*igq+100*pi*igd=0');
f6=('(igd*5000*cos(r-a))+(sin(r-a))*igq*5000-5000*id+wk*eq=0');
f7=('(sin(a-r))*igd*5000+(cos(r-a))*igq*5000-5000*iq-wk*ed=0');
f8=('-3500*(sin(a-r))*igd-3500*(cos(r-a))*i