Usted podría pensar de los números complejos como especiales, transformaciones lineales, de la forma $\begin{bmatrix} a & -b \\ b & a\end{bmatrix}$ donde esta matriz representa el número complejo a $z = a + ib$, una transformación entre los puntos en $\mathbb{R}^2$. Esta matriz representa una rotación y escala uniforme de transformación, donde primero girar su punto y, a continuación, modificar la escala. Esta tal vez puede ser mejor observados a partir de la exponencial de la representación, si $z = re^\theta$, entonces la corepsonding matriz de transformación sería
$$
\begin{bmatrix} a & -b \\ b & a\end{bmatrix} = \begin{bmatrix} r & 0 \\ 0 & r\end{bmatrix}\begin{bmatrix} \cos(\theta) & -\sin(\theta) \\ \sin(\theta) & \cos(\theta)\end{bmatrix} = \begin{bmatrix} r\cos(\theta) & -r\sin(\theta) \\ r\sin(\theta) & r\cos(\theta)\end{bmatrix}
$$
Multiplicar dos números complejos mediante la multiplicación de sus coresponding matrices y multiplicación de matrices de esta forma es comutative. A continuación, $0,1$ $i$ se definen
$$
0 = \begin{bmatrix} 0 & 0 \\ 0 & 0\end{bmatrix} \\
1 = \begin{bmatrix} 1 & 0 \\ 0 & 1\end{bmatrix} \\
i = \begin{bmatrix} 0 & -1 \\1 & 0\end{bmatrix}
$$
Como se puede ver cuando los números complejos se definen de esta manera no hay nada místico en el número de $i$. Es sólo otra transformación, $90^{\circ}$ rotación en sentido antihorario.
EDITAR
También esta forma de pensar acerca de los números complejos intuitivamente explica por qué el Cauchy-Riemann ecuaciones deben poseer para funciones diferenciables.
Una función compleja $f(z) = v(x,y) + iu(x,y)$ puede considerarse como un simple $\mathbb{R}^2 \to \mathbb{R}^2$ transformación (de "xy" coordinación de sys. "uv" ). Así que la mejor aproximación lineal de la función de $f$ cerca del punto de $z = \begin{bmatrix} x \\ y \end{bmatrix}$ (si la función es derivable en z) es el Jacobiano de la matriz $\begin{bmatrix} u_x & v_y \\ u_y & v_y \end{bmatrix}$ que es una generalización de la derivada en más de una dimensión.
$$f(z + \begin{bmatrix} \Delta x \\ \Delta y \end{bmatrix}) \approx f(z) + \begin{bmatrix} u_x & v_y \\ u_y & v_y \end{bmatrix} \begin{bmatrix} \Delta x \\ \Delta y \end{bmatrix}$$
Pero debido a que $f$ es una función compleja, la diferenciabilidad en $z$ significa
$$f( z + \Delta z ) = f(z) + f'(z)\Delta z +o(\Delta z)$$
Donde $\Delta z = \Delta x + i\Delta y$. Esto significa que la función es lineal cerca de $z$ ($\Delta z \to 0$), y por lo tanto
$$f( z + \Delta z ) \approx f(z) + f'(z)\Delta z$$
donde $f'(z)$ es un número complejo que cuando representados por una transformación lineal es de hecho el Jacobiano de esta función. Pero debido a que todos los números complejos son especiales transformaciones lineales de la forma $\begin{bmatrix} a & -b \\ b & a\end{bmatrix}$ la de Cauchy-Riemann ecuaciones seguir, $u_x = v_y, v_x = -u_y$.