(no subject)
May. 18th, 2014 01:02 amМне подбросили задачку, я от нее в некой растерянности.
Есть цель сделать красивую рентгенограмму из плохой, негодной, шумной рентгенограммы. В идеале стоило бы усреднить несколько измерений, но такой возможности по какой-то причине нет (например, измеряемые процессы настолько быстры, что не получается накопить измерений на хороший сигнал).
Человек предлагает сделать следующую имитацию нескольких измерений:
Считается, что интенсивность Y в точке X - это истинная величина плюс-минус некая ошибка, зависящая от значения Y. Ошибка принимает значение (-y0,+y0). Из этого интервала выбирается некое случайное значение (распределение ошибки считается гауссовым), делается поправка для интенсивности (получается "как бы" еще одно измерение), для новой интенсивности процедура повторяется, и так несколько раз. Random walk, короче.
Затем для полученных значений интенсивности считается среднее арифметическое. С ним дальше и работаем.
Мне кажется, что так нельзя. Может появиться ошибка, которая будет накапливаться.
Но насколько нельзя? Как будет выглядеть эта ошибка?
Если чуть-чуть нельзя, то, может быть, стоит пожертвовать ошибкой в пользу такого метода фильтрации шума?
Или в силу случайности процесса ошибка накапливаться не будет? Все же это не последовательные приближения.
Update: при достаточном количестве повторений среднее арифметическое, кажется, просто сведется к исходной интенсивности (ну ок, я не считала, я не люблю нестационарные случайные процессы), поэтому повторений делается немного. не породит ли это еще больше шума?
Есть цель сделать красивую рентгенограмму из плохой, негодной, шумной рентгенограммы. В идеале стоило бы усреднить несколько измерений, но такой возможности по какой-то причине нет (например, измеряемые процессы настолько быстры, что не получается накопить измерений на хороший сигнал).
Человек предлагает сделать следующую имитацию нескольких измерений:
Считается, что интенсивность Y в точке X - это истинная величина плюс-минус некая ошибка, зависящая от значения Y. Ошибка принимает значение (-y0,+y0). Из этого интервала выбирается некое случайное значение (распределение ошибки считается гауссовым), делается поправка для интенсивности (получается "как бы" еще одно измерение), для новой интенсивности процедура повторяется, и так несколько раз. Random walk, короче.
Затем для полученных значений интенсивности считается среднее арифметическое. С ним дальше и работаем.
Мне кажется, что так нельзя. Может появиться ошибка, которая будет накапливаться.
Но насколько нельзя? Как будет выглядеть эта ошибка?
Если чуть-чуть нельзя, то, может быть, стоит пожертвовать ошибкой в пользу такого метода фильтрации шума?
Или в силу случайности процесса ошибка накапливаться не будет? Все же это не последовательные приближения.
Update: при достаточном количестве повторений среднее арифметическое, кажется, просто сведется к исходной интенсивности (ну ок, я не считала, я не люблю нестационарные случайные процессы), поэтому повторений делается немного. не породит ли это еще больше шума?