We present 3D models of the Galactic magnetic field including regular and turbulent components, and of the distribution of matter in the Galaxy including relativistic electrons and dust grains. By integrating along the line of sight, we construct maps of the polarized Galactic synchrotron and thermal dust emissions for each of these models. We perform a likelihood analysis to compare the maps of the Ka, Q, V and W bands of the Wilkinson Microwave Anisotropy Probe (Wmap) and the 353 GHz Archeops data to the models obtained by varying the pitch angle of the regular magnetic field, the relative amplitude of the turbulent magnetic field and the extrapolation spectral indices of the synchrotron and thermal dust emissions. The best-fit parameters obtained for the different frequency bands are very similar and globally the data seem to favor a negligible isotropic turbulent magnetic field component at large angular scales (an anisotropic line-of-sight ordered component can not be studied using these data). From this study, we conclude that we are able to propose a consistent model of the polarized diffuse Galac- tic synchrotron and thermal dust emissions in the frequency range from 33 to 353 GHz, where most of the CMB studies are performed and where we expect a mixture of these two main foreground emissions. This model can be very helpful to estimate the contamination by foregrounds of the polarized CMB anisotropies, for experiments like the Planck satellite.