代码如下：
eql = {\[EmptySet]1'[t] ==
w1 + L*I1[t] +
k*(2*Sqrt[I1[t]*I2[t]]*Sin[\[EmptySet]2[t] - \[EmptySet]1[t]] -
Sqrt[I2[t]/I1[t]] (I2[t] -
I1[t]) Sin[\[EmptySet]2[t] - \[EmptySet]1[t]]) +
k*(2*Sqrt[I1[t]*I3[t]]*Sin[\[EmptySet]3[t] - \[EmptySet]1[t]] -
Sqrt[I3[t]/I1[t]] (I3[t] -
I1[t]) Sin[\[EmptySet]3[t] - \[EmptySet]1[
t]]), \[EmptySet]2'[t] ==
w2 + L*I2[t] +
k*(2*Sqrt[I1[t]*I2[t]]*Sin[\[EmptySet]1[t] - \[EmptySet]2[t]] -
Sqrt[I1[t]/I2[t]] (I1[t] -
I2[t]) Sin[\[EmptySet]1[t] - \[EmptySet]2[t]]) +
k*(2*Sqrt[I2[t]*I3[t]]*Sin[\[EmptySet]3[t] - \[EmptySet]2[t]] -
Sqrt[I3[t]/I2[t]] (I3[t] -
I2[t]) Sin[\[EmptySet]3[t] - \[EmptySet]2[
t]]), \[EmptySet]3'[t] ==
w3 + L*I3[t] +
k*(2*Sqrt[I1[t]*I3[t]]*Sin[\[EmptySet]1[t] - \[EmptySet]3[t]] -
Sqrt[I1[t]/I3[t]] (I1[t] -
I3[t]) Sin[\[EmptySet]1[t] - \[EmptySet]3[t]]) +
k*(2*Sqrt[I2[t]*I3[t]]*Sin[\[EmptySet]2[t] - \[EmptySet]3[t]] -
Sqrt[I2[t]/I3[t]] (I2[t] -
I3[t]) Sin[\[EmptySet]2[t] - \[EmptySet]3[t]]),
I1'[t] == -2 (k*
Sqrt[I2[t]*I1[t]] (I2[t] -
I1[t]) Cos[\[EmptySet]2[t] - \[EmptySet]1[t]] +
k*Sqrt[I3[t]*I1[t]] (I3[t] -
I1[t]) Cos[\[EmptySet]3[t] - \[EmptySet]1[t]]),
I2'[t] == -2 (k*
Sqrt[I2[t]*I1[t]] (I1[t] -
I2[t]) Cos[\[EmptySet]1[t] - \[EmptySet]2[t]] +
k*Sqrt[I3[t]*I2[t]] (I3[t] -
I2[t]) Cos[\[EmptySet]3[t] - \[EmptySet]2[t]]),
I3'[t] == -2 (k*
Sqrt[I3[t]*I1[t]] (I1[t] -
I3[t]) Cos[\[EmptySet]1[t] - \[EmptySet]3[t]] +
k*Sqrt[I3[t]*I2[t]] (I2[t] -
I3[t]) Cos[\[EmptySet]2[t] - \[EmptySet]3[t]])};
ics = {I1[0] == 0.5, I2[0] == 0.5,
I3[0] == 0.5, \[EmptySet]1[0] == 0.1, \[EmptySet]2[0] ==
0.22, \[EmptySet]3[0] == 0.04};
params = {w1 -> -2, w2 -> -1, w3 -> 3, k -> 2.25, L -> 0.5};
soll = NDSolve[{eql, ics} /. params, {I1, I2,
I3, \[EmptySet]1, \[EmptySet]2, \[EmptySet]3}, {t, 0, 1000},
MaxSteps -> \[Infinity]];
m[t] = Mod[
Evaluate[{\[EmptySet]2[t] - \[EmptySet]1[t], \[EmptySet]3[
t] - \[EmptySet]1[t]}] /. soll, 2 Pi]
ParametricPlot3D[Evaluate[{m[t], I2[t]} /. soll], {t, 0, 1000}]
这是运行结果：
求大神指出哪里有问题，万分感激

你执行下
Evaluate[{m[t], I2[t]} /. soll] /. t -> 1
就能看出哪不对了。然后，这方程明显还有别的问题：
Plot[Evaluate[{m[t], I2[t]} /. soll // Flatten], {t, 0, 1000}]