f[lbk]\[lbk]Omega[rbk]_, n_[rbk] :=
Module[lbk]{m = 1, m1 = 1, m2 = 1, k1 = 1, k2 = 0.105,
k3 = 1, \[lbk]Alpha[rbk]1 = 0.5, \[lbk]Alpha[rbk]2 = 0.02,
\[lbk]Alpha[rbk]3 = 0.4, p = 1, \[lbk]Zeta[rbk] = 0.005,
\[lbk]Omega[rbk]n1 = 1.`,
\[lbk]Omega[rbk]n2 = 1.1`,
\[lbk]Phi[rbk]11 = 0.7071067811865476`,
\[lbk]Phi[rbk]12 = 0.7071067811865468`,
\[lbk]Phi[rbk]21 = 0.7071067811865476`,
\[lbk]Phi[rbk]22 = -0.7071067811865482`,
A1 = 0.2250000000000001`,
A2 = 0.07499999999999998`,
A3 = 0.6749999999999999`,
A4 = 0.02499999999999962`,
A5 = 0.02499999999999998`,
A6 = 0.675`,
A7 = 0.07499999999999862`,
A8 = 0.30499999999999994`, sol},
sol = {U1, U2, U3} /.
NSolve[lbk]{(\[lbk]Omega[rbk]n1^2 -\[lbk]Omega[rbk]^2)*U1 +1/(
4 m (\[lbk]Phi[rbk]12 \[lbk]Phi[rbk]21 - \[lbk]Phi[rbk]11 \[lbk]Phi[rbk]22))*(3 A1 U1^3 +A3 \
U1 U2^2 (1 +2 p^2) +(3 A2 U1^2 U2 + 3 A4 U2^3) p) ==
0, (\[lbk]Omega[rbk]n2^2 -\[lbk]Omega[rbk]^2)*U2 +
1/(4 m (\[lbk]Phi[rbk]12 \[lbk]Phi[rbk]21 - \[lbk]Phi[rbk]11 \[lbk]Phi[rbk]22))*(3 A8 U2^3 \
+A6 U1^2 U2 (1 +2 p^2) +(3 A5 U1^3 + 3 A7 U1 U2^2) p ) == 0,
4*m*\[lbk]Zeta[rbk]*\[lbk]Omega[rbk]*((\[lbk]Omega[rbk]n1*U1^2 + \[lbk]Omega[rbk]n2*U2^2)/(
U1 - U2)) - U3 == 0, U1 > 0, U2 > 0, U3 > 0}, {U1, U2, U3},
Reals[rbk]; {\[lbk]Omega[rbk], #[lbk][lbk]n[rbk][rbk]} & /@ sol[rbk]
shuju1 = (f[lbk]#, 1[rbk] &)~ParallelMap~Range[lbk]1.1, 1.5, 0.1[rbk]
shuju2 = (f[lbk]#, 2[rbk] &)~ParallelMap~Range[lbk]1.1, 1.5, 0.1[rbk]
shuju3 = (f[lbk]#, 3[rbk] &)~ParallelMap~Range[lbk]1.1, 1.5, 0.1[rbk]
tu1 = ListPlot[lbk]shuju1, PlotRange -> All[rbk]
tu2 = ListPlot[lbk]shuju2, PlotRange -> All[rbk]
tu3 = ListPlot[lbk]shuju3, PlotRange -> All[rbk]
----
其中U1,U2，U3和omega是变量，我想分别画出以U1,U2，U3为纵坐标，omega为横坐标的图，但是只有第1，2个方程时，它能出结果，如果把关于U3的方程加入，它就会报错没有结果，请大佬帮我看看是什么原因，能有啥解决办法嘛感谢

……跟U3没关系，问题出在你加的那个大于零的限制条件上，以f[1.1, 1]为例，把那几个大于零的限制去掉后会发现方程的解为
{{U1 -> -1.12018, U2 -> 0.041634, U3 -> -0.0237969}, {U1 -> 1.12018, U2 -> -0.041634,
U3 -> 0.0237969}}
明显不符合你给的限制。