DSP engineer here:- the typical way to do this for better performance would be to do the FM modulation of each channel at a lower rate, say 100kHz, then use an FFT-based channeliser to mix each sub-band to the right carrier frequency and upsample in one step. That helps to break the n^2 problem mentioned in the article.