Power Spectral Density of Digital Modulations
At first I found this formula promising, because it does not require a bit sequence,


That's the point about PSDs: They are an expectation value of a stochastic process! That is, the randomness (in this case, the data you're transmitting) is being "averaged away".

Deriving the spectral shape of any modulation (including the pulse shaping, don't forget that – you don't want to end up with a result like the formula you cite, as that is hardly band-limited!) is actually not that hard: You find a formula for the autocorrelation of the signal, and fourier transform that, and then you find the expectation. The trick is in finding a suitable form of the autocorrelation formula!


  but I do not understand the last addend.


That's a speciality of ASK: Most constellations are designed so that the mean of the information symbols is 0; that's not the case for ASK. 

To make it short: The PSD $\Phi_{vv}$ of the baseband signal $v$, consisting of a sequence of information symbols $I_n$ with variance $\sigma_i^2$ and mean $\mu_i$, modulated at a symbol period $T$ with a pulse $g(t)$ (the pulse's spectrum being $G(f)$) is (from Proakis "Digital Communications", 3rd ed., eq. 4-4-17, on page 207):

$$\Phi_{ii}(f) = \frac{\sigma_i^2}{T} \left\lvert G(f) \right\rvert^2 + \frac{\mu_i^2}{T^2} \sum\limits_{m=-\infty}^{\infty} \left\lvert G\left(\frac mT\right) \right\rvert^2\delta\left(f-\frac mT\right)\tag{4-4-17}$$

where $\delta$ is the dirac delta distribution. 

This applies to all linear modulations – PSK, QAM, ASK, OOK, but not to non-linear ones (FSK etc).

In the case of the example you cite, a rectangular pulse shaper $g(t)$ was chosen, which is $1$ for the period $-\frac T2, \frac T2$ around a "symbol time", and $0$ else. That one has a $\DeclareMathOperator{\sinc}{sinc}\sinc$ Fourier transform shape. That's where your $\sinc^2$ comes from! 

The second term, the $\frac{A^2}4 \delta$ is what the sum in eq. $(4-4-17)$ collapses to, because the $\sinc$ is zero in every $mT$ but for $m=0$.

As said above, most constellations are designed to have zero mean ($\mu_i = 0$), so that the sum term just disappears, and your spectral shape is only defined by the spectrum of your pulse. (Only ASK and OOK differ, as far as commonly used modulations are concerned, as far as I can tell.)

Your formula uses $(f-f_c)$ instead of $f$, because it's already the bandpass signal's shape. Don't do that to you! Just consider the signal spectra in baseband, because we can mathematically shift those anywhere we like (especially by an $f_c$).

If you don't have access to good literature: Proakis 3rd ed can be had only for less than a fancy coffee at an American coffee house chain at online used book sellers. It's not my favourite edition of the book, as it cuts the math parts where the second edition really excelled, but it's still a very solid book, and as far as I can see, the cheapest edition of the "very established" digital comms books (being Proakis, Sklar, Tse/Viswanath, and, if you're German, Kammeyer). I like Lapidoth's book, too, and it can be had for free on his website. (update: I just found out that you can get Tse online, too, which is awesome! Might be a bit too advanced, though.)