MR signals are acquired with a method for phase contrast magnetic resonance (MR) imaging with speed encoding, in order to acquire raw data for multiple MR images. The multiple MR images are reconstructed. For this purpose, matrix elements are determined for numerous matrices, wherein the sum of the numerous matrices results in a pixel matrix. The pixel matrix has matrix elements that represent the pixel values for a reference MR image with flow compensation. The pixel matrix has further matrix elements that represent the pixel values for the at least one MR image with speed encoding. The matrix elements of the numerous matrices are determined such that a first matrix of the numerous matrices fulfills a first condition.