La compréhension de la supraconductivité à haute température critique exhibée par certains oxydes de métaux de transition demeure une question centrale en physique théorique de la matière condensée. Dans ce contexte et depuis la proposition historique de P. W. Anderson, le modèle de Hubbard répulsif en géométrie bidimensionnelle est devenu un paradigme pour tenter de capturer les propriétés essentielles des matériaux supraconducteurs non conventionnels. Cependant, la détermination de son état fondamental exact se heurte à la complexité exponentielle du problème quantique à N-corps. Dans cette thèse, nous poursuivons l’objectif de développer un schéma variationnel d’approximation s’affranchissant de toute hypothèse concernant les ordres magnétiques, de charge ou supraconducteurs susceptibles d’émerger de l’Hamiltonien à basse énergie. L’approche trouve son originalité dans l’introduction de corrélations par la restauration, avant variation, de symétries délibérément brisées dans un état d’essai en forme d’une superposition de fonctions d’onde versatiles de type Hartree-Fock et Bogoliubov-de Gennes. Pour des petits amas à deux et quatre sites, nous montrons analytiquement que cette méthode de champ moyen enchevêtré par les symétries de l’Hamiltonien permet de retrouver l’état fondamental exact quelle que soit l’intensité des interactions. Pour de plus grands réseaux dopés en trous et dans les régimes fortement corrélés, nous mettons en exergue un arrangement des moments magnétiques en spirale ou suivant une onde de densité de spin qui est alors accompagnée d’inhomogénéités sous forme de rayures périodiquement réparties. De tels ordres sont de plus entrelacés avec des corrélations d’appariement en onde d à longue distance qui, à la limite thermodynamique, signent la supraconductivité. Ces résultats sont obtenus grâce à des simulations systématiques dans une géométrie de tubes de plaquettes pouvant être réalisée expérimentalement à l’aide d’atomes froids piégés dans des réseaux optiques.