Nous utilisons une fonction de minimisation par moindres carrés fournie par Matlab, appelée « lsqlin » dans laquelle on résout le même système matriciel expliqué précédemment mais où on impose deux contraintes : que les nombres de la matrice solution soient compris entre 0 et 255.

 

Voici les résultats par simulation dans le cas de fortes diffusions :

​

​

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

On voit ici que la diffusion est tellement forte que le texte disparaît (image en haut à droite).

 

On résout le système matriciel à doubles contraintes et on obtient une image X. Lorsqu’on multiplie cette image X avec la matrice de transfert, on peut observer que l’on retrouve le texte. On ne retrouve évidemment pas la totalité de l’information (granulosité) car X est une solution partielle.