This is just an extension of ordinary addition and multiplication with natural numbers. Starting with “2 x 3 = 6” it is inevitable that you’d come up with “(-2) x (-3) = 6”. It’s not some kind of imaginary or weird rule. If we extend our numbers to include negative numbers, and we want to preserve as much behavior as we can from what we observed multiplying positive numbers, then this is the only sensible way of doing things.
When you start with something simple (like positive integers) and extend it, you keep some properties, gain some new properties, and lose some properties. For example, in the transition from rational to real numbers we gain the property that all Cauchy sequences converge. In the transition from real to complex we gain the property that the number of solutions to a nonzero polynomial is equal to the polynomial’s degree, but we lose the property that numbers are ordered.
There are some very deep reasons why complex numbers are a natural choice for doing things in functional analysis. It’s definitely a sweet spot… surprisingly, there’s a concept called “holomorphic functions” which is a very tight constraint on functions, yet simultaneously a right field of study, and it’s the foundation of QM. If you move down the ladder to real numbers, the concept of holomorphic functions does not exist. If you move up the ladder to quaternions or octonions, you lose some critical properties like commutativity.
> …we definitely didn't need them when we adopted them in the 16th century…
They were necessary for solving polynomial equations… even finding real solutions to polynomials with real coefficients.