The following matrix cropped up in a model I am building of a dynamical system:
$$A= \begin{bmatrix} 1 - \alpha & \alpha/2 & 0 & 0 &\cdots & 0 & 0 & \alpha/2\\ \alpha/2 & 1-\alpha & \alpha/2 & 0 &\cdots & 0 & 0 & 0\\ 0 & \alpha/2 & 1-\alpha & \alpha/2 &\cdots & 0 & 0 & 0\\ \vdots & \vdots & \vdots & \vdots & \ddots & \vdots & \vdots & \vdots\\ 0 & 0 & 0 & 0 &\cdots & \alpha/2 & 1-\alpha & \alpha/2\\ \alpha/2 & 0 & 0 & 0 &\cdots & 0 & \alpha/2 & 1-\alpha\\ \end{bmatrix}$$
It is a stochastic matrix and a circulant matrix, and has equal values in the diagonal.
I am interested in the eigenvalues of this matrix, and it was easy to derive them from the properties listed here. It turns out that for size $n$,
$$ \lambda_k = 1 - \alpha \left(1 - \cos\frac{\pi k (n-2)}{n}\right), \qquad k\in\{0,1,\dots,n-1\}, $$
and in the limiting case,
$$ \lim_{n\rightarrow\infty} \lambda_k = 1 - \alpha(1 + (-1)^k) = \begin{cases}1-2\alpha & k \textrm{ even}\\ 1 & k \textrm{ odd}\end{cases} $$
This is interesting for my study, because an eigenvalue of $1$ that is independent of $\alpha$ implies a marginally stable system that cannot be fully stabilized.
Now, I am interested in a slightly modified system, represented by the matrix below. This matrix is exactly like the one above save for the first and last rows, and is still a stochastic matrix with equal values in the diagonal.
I am wondering whether it is possible to derive the eigenvalues of this matrix, even if only for the limiting case.
$$A^\prime= \begin{bmatrix} 1 - \alpha & \color{red} \alpha & 0 & 0 &\cdots & 0 & 0 & \color{red} 0\\ \alpha/2 & 1-\alpha & \alpha/2 & 0 &\cdots & 0 & 0 & 0\\ 0 & \alpha/2 & 1-\alpha & \alpha/2 &\cdots & 0 & 0 & 0\\ \vdots & \vdots & \vdots & \vdots & \ddots & \vdots & \vdots & \vdots\\ 0 & 0 & 0 & 0 &\cdots & \alpha/2 & 1-\alpha & \alpha/2\\ \color{red} 0 & 0 & 0 & 0 &\cdots & 0 & \color{red} \alpha & 1-\alpha\\ \end{bmatrix}$$
Oddly enough, spectral graph theory has the analytic solution for this problem.
Let $L$ be the graph Laplacian for the path graph. Then $L$ is tri-diagonal with $[1,2,2,...,2,1]$ on the diagonal and $-1$'s on the super and sub diagonals. Then the matrix $A'$ is given by $$A'=I-\frac{\alpha}2 L.$$ Now since every vector in $\mathbb{C}^N$ is an eigenfunction of $I$, every eigenfunction $\chi_j$ of $L$ with eigenvalue $\lambda_j$ will be an eigenfunction of $A'$ with eigenvalue $1-\frac{\alpha}{2}\lambda_j$. Thus, to solve the problem stated hereinbefore, one needs only to solve eigenvalue problem of the path graph. The derivation of which can be found here.
To summarize the results, the eigenvalues of $L$ are $\lambda_j = 2-2\cos(\frac{\pi j}{N})$ and may be obtained by mapping the cycle graph with $2N$ vertices to the path graph with $N$ vertices.
Interestingly, the cycle graph with $N$ vertices also yields a method for finding the eigenvalues and eigenvectors of $A$ as $$A=I-\frac{\alpha}2 L_c,$$ where $L_c$ is the graph Laplacian for the cycle graph with $N$ vertices. If you let $v_j=\frac{1}{\sqrt{N}}[\omega^{0j}, \omega^{1j}, ..., \omega^{(N-1)j}]^T$ for $\omega=e^{2\pi i/N}$ (i.e. $v_j$ is a column of the DFT matrix) then $$L_cv_j = (2\omega^{0}-\omega^{j}-\omega^{-j})v_j = (2-2\cos(2\pi j/N))v_j.$$ Hence, $Av_j = (1-\alpha(1-\cos(2\pi j/N)))v_j$ (so I think the indexing in the question might be a bit off) and the DFT matrix forms an orthonormal eigenbasis that diagonalizes $A$. To get a real orthonormal eigenbasis for $A$ (which is guaranteed to exist because $A$ is real and symmetric), taking combinations $w_j=1/2(v_j+v_{N-j})$ and $w_{N-j} = i/2 (v_j-v_{N-j})$ will do the trick.