The internal bit growth should not be an issue for the field engineer.
Your target is to get unity gain.
Your input sinusoid is at 2^15-1 amplitude(orthogonal)
Your fft output says it is sqrt((16068^2 + (-3196)^2)) = 16383
Hence all you need a scale factor of 2.
I expect your exponent to be -1, however you got -9. This issue was raised before in another recent thread. I think it is still not clear if altera got the exponent issue right or the designers got it wrong.
Your post and another post indicate independantly that the altera exponent must be made as exponent + 8 i.e. -9+8 = -1 in order to achieve unity gain.
Edit: Matlab scaling Vs altera scaling
In the case of Matlab scaling, if you want to get unity gain from its fft then you should divide by the fft resolution(by 2048). This leads to unity gain of amplitude across positive and negative frequencies.
Thus in your case:
y = 1/2048 * fft(complex(sink_real,sink_imag));
should give you 32767 as amplitude (one sided).
However if the frequency is two sided then Matlab scaling implies that the sum of two sides = input amplitude. In this case the comparison with hardware may not be useful if you want to use full resolution of both the real and imag bitwidth. It follows that we must be clear about our target scaling. We may target unity amplitude gain or target full use of bit resolution of each of real and imaginary widths