我有一个pde $$ u_ {tt} + u_t = u_ {xx} - \ sin（u）+ f $$，它是在无界样本（无边界）上定义的。我可以通过使用u（t，x）矩阵形式的解决方案的方法在matlab中解决它。问题是，因为我的解决方案是一个传播的解决方案，我的计算域是有限的，很快，解决方案的感兴趣部分就会离开域，我不再能够看到它。这是我的代码和建议的一部分，我从来没有把它放在代码中。
x = [-10:0.1:10];
u0 = asin(f) + 4*atan(exp(x));
[t,u] = ode113(@pde,[t0 tf],u0, options);
plot(x,u(end,:));
这只是我的代码中的几行，让你知道解决方案是一个在力的作用下移动的扭结＆＃34; f＆＃34;向右或向左。我怎样才能在整合中心处纠结？ 我想出了这个想法，以获得“活动”的帮助。作为ode求解器的选项，当它从中心离开特定距离时检测我的解决方案，但不知道我还应该做什么。