考虑定义在$\mathbb{R}^d$有界区域上的含记忆热方程，该方程由$m$个作用于区域内部的控制输入驱动。我们的目标是数值构造一个状态反馈控制器，使得对于每个初始状态，闭环系统的解以大于给定速率$\omega>0$的指数衰减率趋于零，即解决含记忆热方程的$\omega$-镇定问题。我们首先证明了与此方程相关的状态算子$A$的谱在$-\omega_0<0$处具有聚点。给定$\omega\in(0,\omega_0)$，我们证明了若满足某些可验证的关于控制算子$B$的条件，则含记忆热方程的$\omega$-镇定问题是可解的。随后，我们研究了含记忆热方程的一个适当LQR问题。对于每个$n\in\mathbb{N}$，我们构造了$A$和$B$的有限维近似$A_n$和$B_n$，并证明通过求解相应的LQR问题的近似问题，可以计算出反馈算子$K_n$，使得$A_n + B_n K_n$的所有特征值实部均小于$-\omega$。我们进一步证明当$n$充分大时，$K_n$解决了含记忆热方程的$\omega$-镇定问题。我们的证明中一个关键且非平凡的步骤是建立$(A_n+\omega I, B_n)$的一致（关于$n$）可镇定性。我们通过两个数值例子验证了理论结果：一维单位区间上的例子以及二维方形区域上的例子。

We consider a heat equation with memory which is defined on a bounded domain in $\mathbb{R}^d$ and is driven by $m$ control inputs acting on the interior of the domain. Our objective is to numerically construct a state feedback controller for this equation such that, for each initial state, the solution of the closed-loop system decays exponentially to zero with a decay rate larger than a given rate $\omega>0$, i.e. we want to solve the $\omega$-stabilization problem for the heat equation with memory. We first show that the spectrum of the state operator $A$ associated with this equation has an accumulation point at $-\omega_0<0$. Given a $\omega\in(0,\omega_0)$, we show that the $\omega$-stabilization problem for the heat equation with memory is solvable provided certain verifiable conditions on the control operator $B$ associated with this equation hold. We then consider an appropriate LQR problem for the heat equation with memory. For each $n\in\mathbb{N}$, we construct finite-dimensional approximations $A_n$ and $B_n$ of $A$ and $B$, respectively, and then show that by solving a corresponding approximation of the LQR problem a feedback operator $K_n$ can be computed such that all the eigenvalues of $A_n + B_n K_n$ have real part less than $-\omega$. We prove that $K_n$ for $n$ sufficiently large solves the $\omega$-stabilization problem for the heat equation with memory. A crucial and nontrivial step in our proof is establishing the uniform (in $n$) stabilizability of the pair $(A_n+\omega I, B_n)$. We have validated our theoretical results numerically using two examples: an 1D example on a unit interval and a 2D example on a square domain.