QM Basis Transformation Through Unitary Operator

I’ve a short question about basis transformations in QM. Suppose I have two bases ${|{phi_n}rangle}$ and ${|{phi_n’}rangle}$. For brevity, we can make them orthonormal. I know that any state vector can be expanded in terms of both bases:

$$|psirangle = sum_n |phi_nrangle langle phi_n | psi rangle = sum_n |phi_n’rangle langle phi_n’ | psi rangle$$

I also understand that the operator mapping $|phi_nrangle$ to $|phi_n’rangle$, $hat{U}$, is a unitary operator.

Now I have always been told that a change of basis does not change the state vector. If I have expressed my state vector in the basis ${|{phi_n}rangle}$ and want to express it in terms of ${|{phi_n’}rangle}$ I should just apply the identity transformation $hat{1} = sum_n |phi_n’ranglelanglephi_n’|$ and I will have changed the basis.

In this process, the corresponding column matrix will change as if a unitary matrix was applied. From this observation, Zettili’s book page 115-116 and the lecture notes of a course I’m following conclude that in order to express a state vector in a new basis, we should actually apply a unitary operator like this:

$$|psi_text{new}rangle = hat{U}|psi_text{old}rangle$$

I’ve seen this conclusion in multiple places now, but cannot quite follow. Why do we change the state $|psi rangle$? Isn’t the state invariant of the basis chosen? That is the approach of for instance Sakurai in section 1.5.

Or are we keeping the basis fixed while applying the unitary operator? An analogy I see often is a 2D-rotation. Rotating the $(x,y)$ coordinate system through an angle $theta$ to become $(x’, y’)$ while keeping the vector fixed is a transformation of basis (from the $(x,y)$ to the $(x’,y’)$-system). But mathematically, I can also rotate the state vector through an angle $-theta$ and treat the old $(x,y)$-axes as the new $(x’, y’)$-axes. Is this what the sources that confuse me are doing (in a QM context)?