What is the general unitary matrix for two- and three-qubit states?

As pointed out in the QisKit tutorial here, for one qubit there exists a general unitary (see the expression for it in the previous link). I wonder if there exists equally unambiguous expressions for two and three qubits that allow transformation from separable states to all the interesting states, like entangled states in the case of two qubits or GHZ states in the case of three qubits. I was looking at how plot_bloch_multivector works as shown here. Is this in the right direction? An educated guess could be using the SU(N) matrices? It is a 4x4 matrix for two qubits and a 8x8 matrix for three qubits; and then decompose such SU(N) matrices in terms of Pauli and Clifford gates? I don't know if this is possible because the SU(N) matrices are supposed to be given in a parametric form. Another alternative that seems compelling is using Givens rotations, which are a general technique used to parametrize unitary matrices. However, I have no idea on how to implement this.