matlab求解复杂隐函数,MATLAB复杂隐函数求解!求大神指导!

问题描述：

MATLAB复杂隐函数求解!求大神指导!

syms D1 D2 G1 G2 al E1 E2 Km w6

G1=10^(-5);G2=10^(-5); E1=2*(10^11);E2=2*(10^11);Km=1.2;w6=2*(10^9);for D1=1.3:0.1:1.8;

for D2=1.3:0.1:1.8;

ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1)));

ac2=(G2^2)/((Km*0.01/2)^(2/(D2-1)));

w=2*(pi^(1/2))*E1*(G1^(D1-1))*D1*(al^(D1/2))/3;

r=(2/(3-2D1))*(al^((3-2D1)/2)-ac1^((3-2D1)/2))-((al^((2-D2)/2))*2/(5-2D1-D2)*(al^((5-2D1-D2)/2)-ac1^((5-2D1-D2)/2)));

t=Km*w6*D1*(al^(D1/2))/2;

u=(2/(2-D1)*(ac1^((2-D1)/2))-((al^((2-D2)/2)))*2/(4-D2-D1)*(ac1^((4-D2-D1)/2)));

P2=w*r+t*u;

w1=2*(pi^(1/2))*E2*(G2^(D2-1))*D1*(al^((D1+D2-2)/2))/3;

r1=2/(5-2D1-D2)*(al^((5-2D1-D2)/2)-ac2^((5-2D1-D2)/2));

t1=Km*w6*D1*(al^((D1+D2-2)/2))/2;

u1=2/(4-D2-D1)*(ac2^((4-D2-D1)/2));

P1=w1*r1+t