20. И не знаю, что сказать (((
Скорее всего, спектр падает. (умножь FFT на комплексно сопряженную, возьми модуль и посмотри).
Если он падает, то ты можешь низкие гармоники перекинуть на высокие, которые имели малые амплитуды. Тогда после ifft высокочастотные колебания в сигнале могут подскочить на порядок и больше.
Это к примеру

21+. Сомневаюсь, что ты работаешь с белым шумом. Тока там спектр ровный

кажись я нашёл ошибку у себя, завтра проверю и отпишусь )

malin, спасибо за наводку!

ещё вопрос, что бы до завтра не оттягивать, после fft нужно же массив в 2 раза обрезать, т.е:

sig_fur=fft(s_out);

s_fur=s_fur(1:n/2);

или не?

23. На водку я люблю ))
Очень подозреваю, что ты мог курочить весь спектр ))) Там две "половинки" и работать с ними нужно "симметрично" (не забывая, что "0"-частота не в счет ))

Удач

*
sig_fur=fft(s_out);

s_fur=sig_fur(1:n/2);

25 а разве симметрия обязательна?

27. Безусловно. Например, если спектр 1000 точек, а ты хочешь убрать 2 самые низкие частоты, то обнулять нужно 2-3 и 999-1000. Иначе чушь получишь.
Для высоких частот- аналогично, только отсчитывай от центральной частоты.
Рекомендую (сам так делаю, чтобы не запутаться) взять сумму разных синусоид и попробать просто отфильтровать, допустим, высокие. После ifft (не забудь делать real(ifft(...)) ) посмотри картинку

Может и мне поможете? 4тоб тему не перекрывать!.....4ё зна4ит это сообшение,постоянно выскакивает? A script on this page may be busy, or it may have stopped responding. You can stop the script now, or you can continue to see if the script will complete.

28 ок, завтра обязательно всё посмотрю!

ещё раз спасибо!!!

29 это откуда? при открытии страницы? браузер какой?

24. !!!
Кстати, если твоя задача не спектральный анализ ряда, то никаких обрезаний - ни себе, ни спектру.

32 задача сделать частотный скремблер )

31 при открытой странице такое вылазит! Mozila

32+ То, что ты написал,- половина симметричного (относительно центральной частоты) результата FFT. В этом случае, когда считают спектры, учитывают, что энергию уменьшили в 2 раза

34 чесно не знаю что это такое (((
35 это понятно, видимо проблема в том что я крутил только половину спектра, а вторую оставлял какой она и была.

если меня не подводит мой склероз - при ФФТ выход в два раза больше,... или это при ДПФ... ыть, склеротик

чета я аж резко задумался, 4 кГц при 8к отсчетов. какая временная длина сигнала? Или шаг дискретизациии? Ты в найквиста вообще укладываешься?

16 бит на отсчёт, 8000 отсчётов в секунду )))
конечно укладываюсь, телефон-то работает )

я просто привык к 2мс и частатам в 250 Гц предел =)