对于下列具有一般形式的时间分数阶非线性偏微分方程：

其中D_t^α为卡普托型微分算子或者Riemann-Liouville型微分算子，$u=u(x, t), n \in \mathbb{N}_{+}, t>0, x \in \mathbb{R}$.方程(1)包含了许许多多非线性时间分数阶数学物理模型.求解形如方程(1)的分数阶非线性偏微分方程比求解整数阶非线性偏微分方程更加困难，因为大多数整数阶偏微分领域的求解方法不适用于分数阶偏微分方程领域，为此需要重新设计一个新的算法.值得一提的是，在文献[10-11]中，芮伟国教授率先提出了分离变量法与齐次平衡原理相结合的新方法.下面介绍这种方法的具体操作步骤：

第1步：对形如类型(1)的时间分数阶非线性偏微分方程，将它们的解假设成下列两种形式：

或者

其中：v=v(x)为空间变量x的待定函数，a₀，a₁，λ，γ₀，γ₁为待定参数且均为常数，这些待定的函数和参数将在下一个步骤中加以确定；E_α，1(λt^α)称为单参数Mittag-leffier函数，其分数阶导数定义为D_t^αE_α，1(λt^α)=λE_α，1(λt^α).由于(2)式带有一个常数a₀，因此(2)式非常适合Caputo型微分算子下的运算，因为只有在Caputo型分数阶导数下，常数的导数才为零.显然，(3)式非常适合Riemann-Liouville型微分算子下的运算.下面仅以(2)式为例来说明第二个步骤的操作要领.

第2步：将(2)式代入到(1)式中必定可以得到下列形式的方程

像文献[10-11]中那样，利用齐次平衡原理，令Mittag-leffier函数E_α(λt^α)的同次项的系数以及常数项分别等于零，会得到下列形式的一些常微分方程组

对常微分方程组进行求解，可以获得这些方程的解v=v(x)的具体表达式以及相关参数a₀，a₁，λ的取值.当这些函数和参数的值被确定之后，通过(2)式，不难获得形如方程(1)的时间分数阶非线性偏微分方程的各种精确解，具体操作见第3步.

第3步：将第2步中所获得的各种解v=v(x)的具体表达式以及相应的参数a₀，a₁，λ的值代入(2)式，就可获得方程(1)的各种精确解.使用类似的操作步骤和方法，利用(3)式，可以获得方程(1)的其他类型的精确解.

求解下列时间分数阶Klein-Gordon型方程^[14]：

的精确解，其中u=u(x，t)，t＞0，x∈R，将方程(5)改写成

如果方程(5)，(6)是Riemann-Liouville型分数阶导数定义下的分数阶偏微分方程，那么我们就假设方程(6)的解具有如下的形式

其中：c₀，c₁，γ₀，γ₁为待定参数，均为常数且c₁≠0，γ₀＞-1，γ₁＞-1；函数v=v(x)为空间变量的待定函数.把(7)式代入方程(6)中可得

在(8)式中，令时间变量t的幂指数全部相等，可得

求解(9)式即得：

为简便起见，记

将(10)式代入(8)式中消去时间项t^-4α后即得：

假设方程(11)具有下列m次多项式形式的解

根据齐次平衡原理可以确定最高次m的值.如果平衡方程(11)中线性项v与非线性项$v \frac{\mathrm{d}^{4} v}{\mathrm{d} x^{4}}$中x的次数，即得m=4.类似地，如果我们平衡高阶线性项$\frac{\mathrm{d}^{2} v}{\mathrm{d} x^{2}}$和最高阶非线性项$v \frac{\mathrm{d}^{4} v}{\mathrm{d} x^{4}}$中x的次数，可得到m=2.

将m=4代入(12)式可得

其中r₀，r₁，r₂，r₃，r₄为待定系数.将(13)式代入(11)式即得

很容易发现方程(14)有两个孤立项-56ac₁²r₄²x⁶和-84ac₁²r₃r₄x⁵.显然，当且仅当r₄=0时，两个项的系数都等于零.因此，方程(11)没有形如(13)式的解.

但如果将m=2代入(12)式中，却可得到

又将(15)式代入(11)式中，不难得到

在方程(16)中，让x的同次幂的所有系数都等于零，不难得到下列方程组：

解方程组(17)，有

将(18)式代入(15)式可以得到待定函数v(x)的具体表达式

将(19)式和(10)式代入到(7)式中，就获得了方程(6)的一个精确解：

其中:c₀,m 是任意常数, $\mathit{\Omega } {\rm{ = }}\frac{{\mathit{\Gamma } \left( {1 - 2a} \right)}}{{\mathit{\Gamma } \left( {1 - 4a} \right)}},\;0 < a < \frac{1}{2}$且$\alpha \neq \frac{1}{4}$.

当方程(6)是卡普托型分数阶导数定义下的分数阶偏微分方程时，不妨假设方程(6)存在其它形式的解，即

其中：b₀，b₁，λ是待定参数，均为常数且b₀≠0；v(x)是关于空间变量x的待定函数.将(21)式代入(6)式即得

在方程(22)中，令E_α(λt^α)和E_α²(λt^α)的系数都等于零可得下列方程组

首先求解(23)式中第一个线性常微分方程，只讨论Δ=a²b₀²-2bb₀λ²＞0的情况，可以得到2种形式的通解：

其中：${\mathit{\Omega } _0} = \sqrt {\frac{{a{b_0} + \sqrt {{a^2}b_0^2 - 2b{b_0}{\lambda ^2}} }}{{2b{b_0}}}} $, ${\mathit{\Omega } _1} = \sqrt {\frac{{a{b_0} - \sqrt {{a^2}b_0^2 - 2b{b_0}{\lambda ^2}} }}{{2b{b_0}}}} $, ${\mathit{\Omega } _2} = \sqrt { - \frac{{a{b_0} + \sqrt {{a^2}b_0^2 - 2b{b_0}{\lambda ^2}} }}{{2b{b_0}}}} $；B₁，B₂，B₃，B₄，C₁，C₂，C₃，C₄为任意常数.

为了计算简便，令B₃=B₄=0，将(24)式代入(23)式的第二个非线性常微分方程中得

令(26)式中$\mathrm{e}^{2 \mathit{\Omega }_{0} x}$和$e^{-2 \mathit{\Omega }_{0} x}$的系数全部等于零，只需要$a^{2} b_{0}+a \sqrt{\Delta}-\frac{4}{3} b \lambda^{2}=0$，解得

将(27)式及B₃=B₄=0代入到(24)式中得

将(28)式代入(21)式中可以得到方程(6)的2个精确解为

其中：b₀，k₁，k₂为任意常数且b₀＞0.同样令B₁=B₂=0，类似可得到和(29)，(30)相同形式的精确解，只需要将(29)和(30)式中的k₁，k₂改为k₃，k₄即可.

运用同样的方法，取C₁=C₂=0，在此条件下就可以得到方程(6)的另外2个精确解

其中b₀，k₅，k₆为任意常数.

解(29)和(31)的3维坐标图形见图 1.