Landau-Ginzburg thermodynamic formalism is used for the description of the anomalous ferroelectric, ferromagnetic, and magnetoelectric properties of Pb(Fe1/2Ta1/2)x(Zr0.53Ti0.47)1−xO3 and Pb(Fe1/2Nb1/2)x(Zr0.53Ti0.47)1−xO3 micro-ceramics. We calculated temperature, composition, and external field dependences of ferroelectric, ferromagnetic, and antiferromagnetic phases transition temperatures, remanent polarization, magnetization, hysteresis loops, dielectric permittivity, and magnetoelectric coupling. Special attention was paid to the comparison of developed theory with experiments. It appeared possible to describe adequately main experimental results including a reasonable agreement between the shape of calculated and measured hysteresis loops and remnant polarization. Since Landau-Ginzburg thermodynamic formalism appertains to single domain properties of a ferroic, we did not aim to describe quantitatively the coercive field under the presence of realistic poly-domain switching. Information about linear and nonlinear magnetoelectric coupling coefficients was extracted from the experimental data. From the fitting of experimental data with theoretical formula, we obtained the composition dependence of Curie-Weiss constant that is known to be inversely proportional to harmonic (linear) dielectric stiffness, as well as the strong nonlinear dependence of anharmonic parameters in free energy. Keeping in mind the essential influence of these parameters on multiferroic properties, the obtained results open the way to govern practically all the material properties with the help of suitable composition choice. A forecast of the strong enough influence of antiferrodistortive order parameter on the transition temperatures and so on the phase diagrams and properties of multiferroics are made on the basis of the developed theory.