This is one example using Matlab.
http://www.mathworks.com/access/help...G=Search&meta=
You can also do this in IRIS.
http://www.astrosurf.com/buil/iris/new533/new533_us.htm
Ideally you could somehow capture the Impulse response of the imaging system and use that as the PSF. The image of a star perhaps.
Regards
Fahim