我用Maple试了一下,这个问题它给不出符号解.我看这个方程应该是是非线性系统的受迫振动,我就试了一下初始条件为x(0)=0,x'(0)=0的数值解.命令是:sol1 := dsolve([diff(diff(x(t), t), t)+x(t)+0.01*x(t)^3 = cos(t), x(0) = 0, (D(x))(0) = 0], numeric); plots[odeplot](sol1, 0 .. 12)
你可以根据你的需要调整初始条件(x(0) = 0, (D(x))(0) = 0就是初始条件),也可以改变作图时x的范围(我给的是0..12).