A systematic and general method for obtaining orthogonal far-field radiation patterns using any electronically steerable parasitic array radiator (ESPAR) embedded in a propagation environment with an arbitrary power angular spectrum (PAS) is described. The method is useful in designing single-RF multiple-input multiple-output (MIMO) systems where the orthogonal patterns can be utilized as a basis set for beamspace MIMO. The method utilizes the correlation matrix of the open-circuit radiation patterns and is termed the pattern correlation decomposition method (PCDM). PCDM is applicable to any antenna type (including realistic antennas with losses), array configuration, PAS and overcomes the limitations of previously proposed techniques. PCDM also provides closed-form solutions for the antenna currents to generate the orthogonal patterns. It also provides an estimate of the effective aerial degrees-of-freedom of the ESPAR. In addition, efficient optimization approaches to finding the optimal reactive loads that excite the required orthogonal radiation patterns are provided. Simulation examples using a 5-element and a 8-element planar inverted-F antenna array including material losses in various PAS are provided. These show that the proposed approach increases system capacity by up to 7 bits/s/Hz or equivalently reduces the required transmit power by up to 73% compared with previous approaches.