In a method to generate magnetic resonance (MR) images of an examination subject, MR signals are detected simultaneously with multiple coils, each coil having its own coil characteristic. In the detection of the MR signals, raw data space is incompletely filled with MR signals; with raw data space being undersampled in a central raw data region with a coherent acquisition pattern that is composed of a spatially repeating set of raw data points; and raw data space outside of the central raw data region is sampled with an incoherent acquisition pattern. The MR image is reconstructed from the detected MR signals, step-by-step in an iterative reconstruction procedure using a reconstruction matrix A, starting from an initial estimate; wherein the reconstruction matrix has continuing information about the coil characteristics with which the MR signals were detected.