My problem is easy to understand, just how to calculate the matrix of phase operator(or phase gate) acts on multi-qubits so that i can perfrom it in quantum circuit on IBM Quantum Experience
just like the controlled-U gate in phase estimation:
for example:
$$U(y)|φ\rangle=\exp(2πiy)|φ\rangle$$
where $|φ\rangle$ is a two- or three-qubit register.