Well - ok! I do admit that I was simply not willig to even consider a bug in a function which belongs to the prominent Fourier-family ...
At least PeriodogramArray[]
seems to give a correct result:
n = 100;
data = Cos[2 \[Pi] Range[0, n - 1]*0.25] // N;
pdg = Periodogram[data, PlotRange -> All, ScalingFunctions -> {"Linear", "Linear"}, PlotStyle -> Dotted];
pdata = MapIndexed[{(First[#2] - 1)/n, #1} &, PeriodogramArray[data]];
Show[ListLinePlot[pdata[[;; 50]], PlotRange -> All, GridLines -> {{.25}, None}], pdg]
