The article proposes a method for analyzing the Cauchy problem for a wide class of evolutionary linear partial differential equations with variable coefficients. By applying the (inverse) Fourier transform, the original equation is reduced to an integro-differential equation, which can be considered as an ordinary differential equation in the corresponding Banach space. The selection of this space is carried out in such a way that the principle of contraction mappings can be used. To carry out the corresponding estimates for the operators generated by the transformed equation, we impose the conditions of finiteness in the space variable for the inverse Fourier transform of the coefficients, and the spaces of the coefficients of the original equation are determined from the Paley-Wiener Fourier transform theorems. In this case, the apparatus of the theory of the Bochner integral in pseudo-normed spaces, countably-normed spaces and Sobolev spaces is used. Classes of functions are distinguished in which the existence and uniqueness of solutions are proved. For equations with coefficients with separated variables, exact solutions are obtained in the form of a Fourier transform of finite sums for operator exponentials.