You are welcome- I decided to expand on the analysis to handle the case of variable energy-at-infinity, which I call E
The case of E=1 corresponds to "at rest at infinity".
We can write the radial ingoing geodesic using the usual formula:
http://www.fourmilab.ch/gravitation/orbits/
\frac{dt}{d\tau} = \frac{E}{1-2m/r}
\frac{dr}{d\tau} = -\sqrt{E^2 - 1 + 2m/r}
We choose the negative sign on the sqrt because for infalling objects dr/d\tau must be negative.
To convert to EF coordinates we define the tortise coordinate r*, and use v=t + r*.
http://en.wikipedia.org/w/index.php?title=Eddington–Finkelstein_coordinates&oldid=538982166
r* = r + 2m\;ln \left| r/2m -1 \right|
for both r>2m and r<2m we can write
\frac{dr\!*}{dr} = \frac{1}{1-2m/r}
As we agrued previously, v is constant for any infalling light-beam, and can be thought of as the "time of emission at infinity".
Then the doppler shift is the rate of change of v, the time of emission at infinity, with \tau, the proper time.
Thus, knowing that v = t+r*, we can write
\frac{dv}{d\tau} = \frac{dt}{d\tau} + \left( \frac{dr*}{dr} \right) \left( \frac{dr}{d\tau}\right) = \frac{E}{1-2m/r} - \frac{\sqrt{E^2-1+2m/r}}{1-2m/r}
We can see that we have the difference of two infinite quantites at r=2m :-(.
To get the doppler shift at r=2m, it's convenient to do a series expansion of the above expression around r=2m
This gives:
\frac{dv}{d\tau} \approx \frac{1}{2E} + \frac{1}{16\,m\,E^3} \left(r -2m \right) + \frac{2E^2-1}{64\,m^2\,E^5} \left(r -2m \right)^2 + ...
For values at or near the horizon we can use the series expansion, for values far away from the horizon we should use the original expression (the one which becomes singular at the horizon).
For E<1, the object can never reach infinity. For the limit of low E, one has a blueshift at the horizon, which makes sense (the object didn't have enough energy to reach infinity, and can be considered to be dropped in from some low value of r). It started out experiencing blueshift, and continues to experience blue shift as it falls in
For E=1, we see the doppler shift I mentioned of 1/2 at the event horizon. E>1 implies starting out with some velocity towards the black hole at infinity which gives more redshift.
I believe the computed value for redshift also makes sense if we take the limit as r->infinity, but I haven't looked at that as closely as I might. Clearly the doppler shift is 1 (no doppler shift) for E=1 at r=infinity, which is correct.