syms r P o q t v u b a e m
solve('b^2/(a^2-b^2)*(1+a^2/r^2)*q-a^2/(a^2-b^2)*(1+b^2/r^2)*o-u/(1-u)*P*(r/t)^v-2/(a^2-b^2)*(b^2*q-a^2*o)-(1-2u)/(1-u)*P*(r/t)^v-e*(m*(2/(a^2-b^2)*(b^2*q-a^2*o)+(1-2u)/(1-u)*P*(r/t)^v)/e+1)^0.5')

用的是matlab2012b，除了r是自变量以外，其他字母都是已知量

matlab解复杂方程最好自己编程

修正版
syms r P o q t u b a e m
r=solve('b^2/(a^2-b^2)*(1+a^2/r^2)*q-a^2/(a^2-b^2)*(1+b^2/r^2)*o-u/(1-u)*P*(r/t)^（2+u/(1-u))=2/(a^2-b^2)*(b^2*q-a^2*o)+(1-2*u)/(1-u)*P*(r/t)^（2+u/(1-u))+e*(m*(2/(a^2-b^2)*(b^2*q-a^2*o)+(1-2*u)/(1-u)*P*(r/t)^（2+u/(1-u)))/e+1)^0.5','r')
结果还是一样的