Использую сигналы после преобразования Гильберта для построения диаграмм мгновенной фазы, частоты и амплитуды. На различных модуляциях (простые амплитудная, частотная, фазовая) получается примерно одинаковая картина. Действую по формулам:
Phase_out[i]=atan(_in[i].imag()/_in[i].real());
Amp_out[i]=sqrt(pow(_in[i].real(),2)+pow(_in[i].imag(),2));
Freq_out[i]=((_in[i].imag()*_in[i-1].real()) - (_in[i-1].imag()*_in[i].real()))/ (2*pow(IPP_PI,2)*(pow(_in[i].real(),2)+pow(_in[i].imag(),2)));
Что получается - можно увидеть в прикрепленной картинке.
Что я делаю не так?