[an error occurred while processing this directive]
|
Может не понял чего, но ...
В спектре действительного сигнала (N точек, N четное) комплексно-сопряженной симметрией обладают (N-2)/2 спектральных отсчета. Спектр. отсчеты на 0 и на N/2 сами по себе, как кошки.
В Вашем случае комплексного сигнала отсчет на половине новой ч-ты дискретиз. неизвестен. Т.е. дополнять нужно так:
Re(X[N+1])=Re(X[N-1]), Im(X[N+1]=-Im(X[N-1]);
Re(X[N+2])=Re(X[N-2]), Im(X[N+2]=-Im(X[N-2]);
...
X[N] особая точка, к тому же неизвестная, к 0 ее можно.
N - длина исходного FFT.
У Вас вроде немного по другому, может просто описка, конечно:
Re(N+1)=Re(N),Im(N+1)=-Im(N).
Ну и для всех частот, кроме 0 и N, амплитуду в 2 раза надо бы уменьшить.
E-mail: info@telesys.ru