Question

为什么运行程序时MATLAB一直处于busy状态？？？求高手指点

syms G1 E1 G2 E2 a b Wa1 Wa2 Wb1 Wb2 Wc1 Wc2 t1 t2 u1 u2 Ec Ac X1 X2 q G Geb1 Geb2 Ec Ac k l p;
f1=((8/15-k+2/3*k^3-k^5/5)+1.7625*(G2*(1-k^3)+Geb2*k^3)/(E2*X1^2*a^2))*Wa1-0.8059*(G2*(1-k^3)+Geb2*k^3)/(E2*X1^2*a^2)*Wa2-(2/3-k+k^3/3)*Wc1;
f2=((8/15-k+2/3*k^3-k^5/5)*0.51953+1.019337*(G2*(1-k^3)+Geb2*k^3)/(E2*X2^2*a^2))*Wa2-0.8059*(G2*(1-k^3)+Geb2*k^3)/(E2*X2^2*a^2)*Wa1-0.51953*(2/3-k+k^3/3)*Wc2;
f3=((8/105-l^3/3+2/5*l^5-l^7/7)+1.32187*G1*(4/5-l+2*l^3-9*l^5/5)/(E1*X1^2*b^2))*Wb1-0.60435*(G1*(4/5-l+2*l^3-9*l^5/5)+Geb1*(l-2*l^3+9*l^5/5))/(E1*X1^2*b^2)*Wb2-0.60435*(G1-Geb1)*(l-l^3)/(E1*X1^2*b^2)*Wc2+1.32187*(l-l^3)*(G1-Geb1)/(E1*X1*b)*u1-0.60435*(G1-Geb1)*(l-l^3)/(E1*X1*b)*u2-(2/15+l^3/3-l^5/5)*Wc1;
f4=((8/105-l^3/3+2/5*l^5-l^7/7)*0.51953+0.764498*G1*(4/5-l+2*l^3-9*l^5/5)/(E1*X2^2*b^2))*Wb2-0.60435*(G1*(4/5-l+2*l^3-9*l^5/5)+Geb1*(l-2*l^3+9*l^5/5))/(E1*X2^2*b^2)*Wb1+1.32187*(G1-Geb1)*(l-l^3)/(E1*X2^2*b^2)*Wc1-0.60435*(l-l^3)*(G1-Geb1)/(E1*X1*b)*u1+0.60435*(G1-Geb1)*(l-l^3)/(E1*X1*b)*u2-0.51953*(2/15+l^3/3-l^5/5)*Wc2;
f5=-1.8557*(2/3-k+k^3/3)*Wa1-(1.8557*(2/15-l^3/3+l^5/5)*E1*t1*b/(E2*t2*a)+2.4530*t1*(G1-Geb1)*(l^3-l)/(E2*t2*a*b*X1^2))*Wb1+1.1215*t1*(G1-Geb1)*(l^3-l)/(E2*t2*a*b*X1^2)*Wb2+(1.8557*(1-k)+1.8557*(1/3-l^3/3)*E1*t1*b/(E2*t2*a)+1.8557*Ec*Ac/(E2*t2*a)+2.4530*(G1*t1*(1-l)/(E2*t2*a*b*X1^2)+Geb1*t1*l/(E2*t2*a*b*X1^2)))*Wc1-1.1215*t1*(G1*(1-l)+Geb1*l)/(E2*t2*a*b*X1^2)*Wc2+2.4530*t1*(G1*(1-l)+Geb1*l)/(E2*t2*a*X1)*u1-2.80755*t1*(G1*(1-l)+Geb1*l)/(E2*t2*a*X1)*u2;
f6=-0.9641*(2/3-k+k^3/3)*Wa2+1.1215*t1*(G1-Geb1)*(l^3-l)/(E2*t2*a*b*X2^2)*Wb1-(0.9641*(2/15-l^3/3+l^5/5)*E1*t1*b/(E2*t2*a)+1.4185*t1*(G1-Geb1)*(l^3-l)/(E2*t2*a*b*X2^2))*Wb2-1.1215*(t1*G1*t1*(1-l)/(E2*t2*a*b*X2^2)+Geb1*t1*l/(E2*t2*a*b*X2^2))*Wc1+(0.9641*(1-k)+0.9641*(1/3-l^3/3)*E1*t1*b/(E2*t2*a)+0.9641*Ec*Ac/(E2*t2*a)+1.4185*(G1*t1*(1-l)/(E2*t2*a*b*X2^2)+Geb1*t1*l/(E2*t2*a*b*X2^2)))*Wc2-0.448*t1*(G1*(1-l)+Geb1*l)/(E2*t2*a*X2)*u1+1.4185*t1*(G1*(1-l)+Geb1*l)/(E2*t2*a*X2)*u2;
f7=2.4530*(G1*(1-l)+Geb1*l)*u1*X1*b-1.1215*(G1*(1-l)+Geb1*l)*u2*X2*b+2.4530*(G1*(1-l)+Geb1*l)*Wc1-1.1215*(G1*(1-l)+Geb1*l)*Wc2-2.4530*(G1*(l^3-l)+Geb1*(l-l^3))*Wb1+1.1215*(G1*(l^3-l)+Geb1*(l-l^3))*Wb2-1.8556*(0.5478*q+0.4176*p)/(4*t1*X1);
f8=-1.1215*(G1*(1-l)+Geb1*l)*u1*X1*a+1.4185*(G1*(1-l)+Geb1*l)*u2*X2*a-1.1215*(G1*(1-l)+Geb1*l)*Wc1+1.4185*(G1*(1-l)+Geb1*l)*Wc2+1.1215*(G1*(l^3-l)+Geb1*(l-l^3))*Wb1-1.4185*(G1*(l^3-l)+Geb1*(l-l^3))*Wb2-0.9146*(0.4419*q+0.09245*p)/(4*t1*X2);
S=solve(f1,f2,f3,f4,f5,f6,f7,f8,'Wa1,Wa2,Wb1,Wb2,Wc1,Wc2,u1,u2');

Answer 1

这么多指数幂级数运算。。里面如果再有矩阵，那是要跑一会

Answer 2

符号运算是相当慢的