1楼: Originally posted by _aqing at 2011-11-10 21:09:57:
这个微分方程组，实在解不出。求大家给个意见或参考。能给出解的定重谢。

M=200;c=2;b=2;
odefun=@(t,x)[M*(-3/(8*pi*x(2)))*(17*c*x(1).^3/x(2).^2+3*b);M*(3/(16*pi*x(1)))*(11*c*x(1).^3/x(2).^2-b)];
[t,x]=ode45(odefun,[0 10],[2 3]);
plot(t,x)

5楼: Originally posted by zhaozhh02 at 2011-11-14 16:41:31:
忘了说了，只能数值求解。这个方程组是非线性的，极难有解析解。

9楼: Originally posted by _aqing at 2011-11-17 17:14:38:
给出参数 大家帮忙算算
具体参数
b=34.5086
M=0.0048
c=1.6*10^8*(ln(8.16*L/r)-1)
初值
l=1.35*10^(-4)
r=110*10^(-7)