Case 1: Invertible Matrices
For each $A \in M_{n\times n}(\mathbb{C})$, let us consider the characteristic polynomial
\begin{align}
p_A(z) = \det(zI-A).
\end{align}
Then it is clear that $A$ is not invertible if and only if zero is a root of of $p_A(z)$. However, if we perturb $A$ by a diagonal matrix, i.e.
\begin{align}
B= A+\epsilon I
\end{align}
then we see that $B$ is invertible since the characteristic polynomial of $B$ is $p_B(z) = p_A(z-\epsilon)$ no longer has zero as its root for any $\epsilon \neq 0$. Hence the class of invertible matrices are dense in $M_n(\mathbb{C})$. (In fact, the set of invertible matrices is actually open in $M_n(\mathbb{C})$. )
Case 2: Unitary Matrices
A quick way to see that this class is not dense is to consider the following matrix
\begin{align}
A =
\begin{pmatrix}
2 & 0\\
0 & 2
\end{pmatrix}
\end{align}
which is not unitary. In particular, for all unitary matrix $U$ we have
\begin{align}
\operatorname{Tr}|A-U|\geq \big|\operatorname{Tr}|A|-\operatorname{Tr}|U| \big|= 1.
\end{align}
Since all norms are equivalent on finite dimensional vector spaces, we have that the class of unitary matrices can't be dense in $M_n(\mathbb{C})$.
Case 3: Symmetric Matrices
A quick way to see that this class is not dense is to consider the matrix
\begin{align}
C =
\begin{pmatrix}
1 & 100\\
0 & 1
\end{pmatrix}
\end{align}
which is clearly not symmetric. Now, for any symmetric matrix $S$, we have
\begin{align}
C-S =
\begin{pmatrix}
1-a & 100-b\\
-b & 1-c
\end{pmatrix}.
\end{align}
If we take the Frobenius norm of $C-S$, we get
\begin{align}
\|C-S\|_F = \sqrt{|1-a|^2+|1-c|^2+b^2+|100-b|^2}\geq \sqrt{b^2+|100-b|^2} \geq \frac{100}{\sqrt{2}}.
\end{align}
Hence the class of symmetric matrices is not dense in $M_n(\mathbb{C})$.
Case 4: Diagonalizable Matrices
Let $A$ be an arbitrary matrix in $M_n(\mathbb{C})$. Let us show that $A$ can be approximated by diagonal matrix.
Using the fact that a matrix is diagonalizable if it has $n$ distinct eigenvalues, we shall construct a $B$ matrix arbitrarily close to $A$ in norm and $B$ has distinct $n$-eigenvalues.
Suppose $p_A(z) =\det(zI-A) =(z-\lambda_1)(z-\lambda_2)\cdots (z-\lambda_n)$ is the characteristic polynomial of $A$ where $|\lambda_1\leq |\lambda_2|\leq \ldots \leq |\lambda_n|$, then consider a diagonal matrix
\begin{align}
D=
\begin{pmatrix}
d_1 & 0 & \ldots & 0\\
0 & d_2 & \ldots & \vdots\\
\vdots & \ldots &\ddots & 0\\
0 & \ldots & 0 & d_n
\end{pmatrix}
\end{align}
where the $|d_1|<|d_2|<\ldots <|d_n|$. Then by considering the Jordan canocial form of $A$, we have that
\begin{align}
A+\epsilon D
\end{align}
has distinct eigenvalues, i.e. diagonaliable for any $\epsilon\neq 0$. Hence the set of diagonalizable matrices are dense in $M_n(\mathbb{C})$. (My explanation of this part is flimsy. Essentially, I'm choosing $d_1, \ldots, d_n$ so that $|\lambda_1+d_1|<|\lambda_2+d_2|<\ldots<|\lambda_n+d_n|$.)