I've never done it in Perl, but did do it in C a few million years ago using algorithms from Numerical Recipes. A search on "Levenberg Marquardt perl" turns up this:

It's not a totally canned tool for what you want, but you should be able to use Levmar to do what you want, at least fit-wise