|
|
Вы так и не сказали толком, что же Вам нужно.
Дело то ведь будет определяться и окном и типом "колбасни".
А колбасню Вы скромно опускаете, оставляете одно окно.
И какие проблемы в этой ситуации с объединением?
Вам даже обратное Фурье не надо делать, запихивайте в OUT входной массив, да и дело с концом.
Если уж очень хочется сделать обратное БПФ и затем сшить куски, не колбася спектр, как во всех Ваших текстах.
Y = fft( W .* X );
Y = ifft( Y );
После IFFT Y=W*X с точностью до ошибок округления.
Ну и делите Y на W, получите X.
А в инетовском примере делается следующее, на примере 8 точечных Фурье, 8-ми точечного окна и 4-х точечных кусков сигнала.
1. X={0,0,0,0,IN0,IN1,IN2,IN3}
1a. Y (после IFFT)={0,0,0,0,W4*IN0+W5*IN1+W6*IN2+W7*IN3}
1b. OUT={0,0,0,0,0,0,0,0,}
1c. Y0=Y={0,0,0,0,W4*IN0+W5*IN1+W6*IN2+W7*IN3}
2. X={IN0,IN1,IN2,IN3,IN4,IN5,IN6,IN7}
2a. Y={W0*IN0,W1*IN1,W2*IN2,W3*IN3,W4*IN4,W5*IN5,W6*IN6,W7*IN7}
2b. OUT={0,0,0,0,W4*IN0+W0*IN0,W5*IN1+W1*IN1,W6*IN2+W2*IN2,W7*IN3+W3*IN3}
2c. Y0=Y={W0*IN0,W1*IN1,W2*IN2,W3*IN3,W4*IN4,W5*IN5,W6*IN6,W7*IN7}
В OUT: W4+W0=W7+W3~=1, W5+W1=W6+W2~=1. Но W4+W0!=W5+W1 (для окна Хэннинга, для треугольного равны и все равенства точные). Так что OUT после второй ступени:
OUT~={0,0,0,0,IN0,IN1,IN2,IN3}
Равенство, кстати, можно сделать точным, если точный OUT поэлементно поделить на следующий массив:
{1,1,1,1,W4+W0,W5+W1,W6+W2,W7+W3}.
Ну и для следующих ступеней аналогично.
Какие проблемы, если со спектром ничего не делается, хоть в Ваших вариантах, хоть в найденном?
Но Вы же колбасите как-то спектр. Вот колбаса вместе с окном и будет результат определять. Что Вы хотите получить, Вы так и не сказали.
E-mail: 
info@telesys.ru
