Projet matlab IRM

Madame, Monsieur, je suis actuellement en 1ere année d'ingénieur spécialité à l'ESIR de Rennes. Je suis confronté à un projet Matlab(logiciel de programmation/simulation) concernant le contraste en IRM. Sachant que vous avez des connaissances de ce domaine j'aimerais a voir des informations sur le fonctionnement de ce procédé. Aussi, je voudrais savoir si vous pouviez m'aider à solutionner le projet décrit ce dessous : " Avec un certain type de séquence, on peut montrer que le signal délivré par une IRM est proportionnel à p[-exp(-TR/T1)]exp(-TE/T2), où p est la densité de proton, T1 le temps de relaxation spin-réseau et T2 le temps de relaxation spin-spin. Ces trois quantités sont des fonctions spatiales dépendant des tissus dans lesquels se trouvent l'eau. TE et TR sont deux paramètres fixés par l'opérateur. Ecrire un programme permettant de simuler la dépendance du contraste en fonction des différents paramètres (images pondérées en densité, en T1 ou en T2)". Mes connaissances en IRM et Matlab sont aux points mort et ce sont des domaines assez complexes et difficiles à comprendre. En espérant une réponse de votre part, Cordialement, NOYER Germain