Некоторые методы статистического обнаружения и оценивания сигналов в астрономии Многие методы обработки сигналов оказываются неэффективными при наблюдениях. - презентация

1 Некоторые методы статистического обнаружения и оценивания сигналов в астрономии Многие методы обработки сигналов оказываются неэффективными при наблюдениях слабых звезд в условиях дефицита квантов. Излагаются новейшие техники статистического обнаружения и оценивания сигналов в разреженных потоках квантов. ****************** Many techniques fail to work in the condition of light quantum deficiency, in particular, on short time scale or if we observe faint stars. Modern techniques make possible the detection of events comparable to the limit of the counting statistic. Detection and evaluation the signal in sparse quantum fluxes is presented using modern statistical approaches.

2 Содержание Обнаружение сигналов в разреженных пуассоновских потоках квантов с использованием критериев на основе энтропии Фурье и вейвлет спектральный анализ

3 Detection of Signal in Sparse Poissonian Streams of Quanta Using Entropy - Based Criteria Definition Simulation Empirical quantilies for entropy function Deterministic signal Sinusoidal signal with random phase Chaotic signal See UHF&Entropy.doc

4 Энтропия В физике где Г, E – число квантовых состояний подсистемы в интервале энергий E (микроканонический ансамбль), w(E) – распределение вероятностей энергии. Таким образом, энтропию можно определить как среднее значение логарифма функции распределения энергии подсистемы. В теории информации p i 0 – вероятность какого-либо испытания (опыта). Энтропия равна нулю, когда одно из p i равно единице, т.е. исход опыта достоверен. Энтропия достигает максимума, когда все исходы равновероятны.

5 Entropy is a common concept in many fields, mainly in signal processing. Classical entropy-based criteria describe information- related properties for an accurate representation of a given signal. The entropy E must be an additive cost function ( ) such that E(0) = 0. The following example lists different entropy criteria, many others are available and can be easily integrated. In the following expressions, s is the signal and (s i ) i the coefficients of s in an orthonormal basis. ( ) cost function - функция стоимости

6 The (nonnormalized) Shannon entropy where s is the signal and (s i ) i the coefficients of s in an orthonormal basis. The threshold entropy E4(s i ) = 1 if |s i | > p and 0 elsewhere, so E4(s) = #{i such that |s i | > p} is the number of time instants when the signal is greater than a threshold p.

7 Распределение Пуассона Распределение Пуассона - дискретное распределение с одним параметром. Параметр является средним и дисперсией распределения. Гистограмма – результат моделирования Пуассоновского распределения с интенсивностью = 0.1 (1000 испытаний). Красные кружки – теоретическое распределение с той же интенсивностью.

8 Sparse Poissonian stream: ( Poissonian intensity =.1, i.e. 90% of zeros in average ) L=1000; Lambda =.1*(1+cos (2*pi/50*(1:L))); for i =1:L; v(i) = poissrnd(Lambda(i),1,1); end; Sparse signal modulated by harmonic. **************************************************** L=1000; Lambda0 =.1; for i =1:L; v0(i) = poissrnd(Lambda0,1,1);end; Sparse signal of the same mean intensity without modulation (No signal) Codes: MATLAB

9 Однако, хотя мы имеем разреженный модулированный поток квантов гармонический сигнал детектируется методами спектрального анализа Визуально ни кривые блеска, ни гармонический сигнал не распознаются Signal No signal

10 Simulation (quantile estimator - оценка квантилей) No signal: A=0; signal amplitude = 0; 100 trials for j=1:100; L=1000; v = zeros(1,L); Lamb =.1*(1+A*cos(2*pi/50*(1:L))); for i =1:L;v(i)=poissrnd(Lamb(i),1,1);end;V=[V;v];e = entropy(v,'shannon');E=[E;e];end; The same with a signal: A=1 mean(E) = ; No signal (top) mean(E) = ; With signal (bottom) The signal reduces entropy COMMENTS: Есть обнаружение примерно в 25% случаев (при значении энтропии Шеннона < ~ 60) Необходим более строгий критерий

11 Simulation of the Shannon entropy: Poisson intensity = 0.1 per sampling time The number of measurements N = Three trials in hundreds, no signal n = 100; N = ; Lambda =.1; for j=1: n; v = poissrnd(Lambda,1,N); e = wentropy(v,'shannon'); E = [E;e]; end; done triply : E1, E2, E3 – three arrays shown in the above Figure mean([E1 E2 E3]) : – the mean values of arrays std([E1 E2 E3]): – their standard deviations, abs values std([E1 E2 E3])./mean([E1 E2 E3]) : – relative values Когда теоретическая функция распределения (ф.р.) какой-либо случайной величины не известна ее можно заменить эмпирической ф.р., полученной в результате моделирования. Ф.р. позволяет вычислять квантили – такие значения случайной величины, когда отклонения можно считать не случайными.

12 Empirical quantilies for the entropy function quantile: квантиль порядка p, 0 < p = 0.95)) % answer = 3050, 3100, 3150, …3300 K 0.95 = 3050 – the 95% empirical quantile for the entropy function Shown are the individual and averaged integral probability distributions for the Shannon entropy function resulted from three simulation trials. The empirical quantilies: 3050 – 95% 3150 – 99%

13 Now we shall take an harmonic signal A=1; N=100000; vs=zeros(1,N); Lambda =.1*(1+A*cos(2*pi/50*(1:N))); for i =1:N;vs(i)=poissrnd(Lambda(i),1,1);end;e = wentropy(vs,'shannon');% with signal 3 trials: e1 = 4585 ; e2 = 4608 ; e3 = 4158 THUS, e1, e2, e3 > 3150 – DETECTION PRESENTS! The same, A = 0.5: e1= 3179, e2 = 3436, e3 = 3184 > 3150 Thus, detection is present (marginally detected at 99% confidence; A=0.5; N=100000) sum(vs) = (quanta) Вывод: зарегистрировав всего лишь квантов в измерений, можно с вероятностью 99% обнаружить гармонический сигнал с глубиной модуляции, равной ½, используя функцию энтропии Шеннона. Codes: MATLAB

14 An algorithm for the threshold entropy function 1) Let mean intensity per sampling time dt is equal to LAMBDA ( ). 2) Let N is the number of measurements. 3) Then the measured value of threshold entropy function ENTROPY_T may be defined as ENTROPY_T = sum(x > THRESHOLD), where x are the counting in a series of measurements. Let THRESHOLD = 1 for sparse fluxes of quanta. 4) The mean theoretical value ENTROPY_TT = (1-poisscdf(THRESHOLD, LAMBDA))*N, where poisscdf - Poisson cumulative distribution function. floor n наименьшее целое число, не превосходящее n. 5) The two-sigma value corresponding to the 95% confidence level is defined as ENTROPY_TT + 2*sqrt(ENTROPY_TT). Let a signal S = Lambda*(1+A*cos(2*pi/P*(1:N)+rand(1,N)*2*pi)), where A, P are the amplitude and period of a harmonic signal with random phase (uniformly distributed on the interval 0, 2*pi).

15 The empirical quantile for threshold entropy function n=100; N=100000; Lambda =.1 for j = 1: n; v = poissrnd (Lambda,1,N); V = [V; sum(v>1)]; end; min(V) = 421; max(V) = 520 Итак, по результатам моделирования для заданных N и Lambda : K 0.95 = 510 – 95% эмпирический квантиль K 0.99 = 520 – 99% эмпирический квантиль

16 A sinusoidal signal with random phase A = …; N = ; n = 3 (three times repeating) for j =1 : n; L=.1*(1+A*cos(2*pi/50*(1:N)+rand(1,N)*2*pi)); for i = 1 : N; vsp(i) = poissrnd (L(i),1,1); end; e = wentropy(vsp,'threshold',1);E=[E;e]; end; % (MATLAB Codes) A = 1.0; E = ; DETECTS WITH CERTAINTY (99% = 520) A = 0.5; E = ; DETECTS MARGINALLY (99% = 520) A = 0.25; E = ; THERE IS NO DETECTION (95% = 510) For certain detection of a signal with small amplitude A we must increase N – the number of measurements Further amplitude decreases with ½ power of N (Почему?)

17 Now we shall take a chaotic signal A=1; N=100000; Lambda =.1*(1+A*cos(2*pi/50*(1:N))); L = eralash(Lambda); (Destroyed harmonic signal) for i =1: N; vs (i) = poissrnd (L(i),1,1); end; E = wentropy(vs,'shannon'); 3 trials: E = 3998, 4544, 4043 > 3150 where 3150 – the 99% empirical quantile THUS, DETECTION PRESENTS! Lambda (top) and L (bottom) signals COMMENT FFT detects harmonic modulation easily (left Fig) Chaotic signal reveals the Shannon entropy detector

18 1.Естественные источники излучения (пуассоновские и гауссовские потоки квантов) характеризуются минимальным информационным содержанием (только первым и вторым статистическими моментами – средним значением и дисперсией) и, следовательно, максимальной энтропией. Любое информационное сообщение или полезный сигнал уменьшают энтропию потока квантов. Измеряя энтропию, можно обнаружить сигнал. 2.Можно сконструировать различные обнаружители сигналов, например, основанные на энтропии Шеннона, энергетических, пороговых или определенных пользователем функциях энтропии. 3.Классические энтропийные критерии сравнительно просто описывают информационные свойства стационарных случайных и периодических сигналов. 4.Пороговый статистический обнаружитель может измерять энтропию потока данных и выдавать сообщение о наличии полезного сигнала. Таким образом происходит статистическое обнаружение сигнала. 5.Последующий детальный статистический анализ направлен на оценивание амплитудно-частотных характеристик обнаруженного сигнала. Резюме

19 Detection and Evaluation of Signals Фурье и вейвлет спектральный анализ Наличие в любых измерениях ошибок, случайных по своей природе, требует применения статистических методов при обработке данных. Ни один статистический вывод не может быть достоверным. Вероятностная природа статистических выводов приводит к тому, что можно отвергнуть верную гипотезу (ошибка 1-го рода) и принять ложную гипотезу (ошибка 2-го рода). Это плата за ту неопределенность, которую вносят в нашу жизнь шумы измерений. «Искусство предположений» - так вначале называл то, что со временем стало теорией вероятности и математической статистикой, ее основатель Даниил Бернулли.

20 Фурье спектр Спектр Фурье по своему физическому смыслу пропорционален корреляции сигнала с монохромной волной на заданной частоте (в случае комплексного представления волны спектр - комплексная величина). Поскольку монохромная волна не локализована в координатном пространстве, то при интегрировании по пространственной координате остается зависимость спектра только от частоты (переход от координатного представления сигнала к частотному представлению). Сигнал f(t) можно разложить в ряд Фурье:

21 Фурье спектр Прямое и обратное преобразование Фурье (Discrete Fourier Transform, DFT) позволяют оперировать с сигналом х как во временной, так и в частотной области. Оба представления сигнала равноправны и связаны взаимно однозначным соответствием.

22 Рассмотрим гармонический сигнал Его спектр на частотах n = ± l будет равен В случае нормального «белого шума» с дисперсией отношение сигнал / шум в спектре мощности имеет вид

23 Где - распределение хи-квадрат с двумя степенями свободы, - процентная точка распределения хи-квадрат. При = критическое значение отношения сигнал / шум в спектре мощности на уровне доверительной вероятности 99.5% будет равно Q = Естественно возникает вопрос о выборе доверительной вероятности. Можно предложить следующее правило: Если выбрать величину = 1 – 2 / N, то в среднем один шумовой пик будет иметь величину больше, чем. Такой уровень доверительной вероятности можно называть one noise peak level. Можно установить one tenth/hundredth…level. Это будет означать, что только в одном случае из 10/100 пик в спектре может быть шумовым. Итак, вероятность пика может быть оценена.

24 Вейвлет спектр В вейвлет анализе рассматривается свертка сигнала с пространственно ограниченной волной, поэтому вейвлет спектр - двумерная величина, зависящая от двух координат (обычно временных величин): протяженности ограниченной волны и несущей частоты волны (ее периода). Поскольку частота строго определена только для монохромной волна (бесконечной протяженности), то говорят не о периоде несущей частоты, а о шкале, или эквивалентном Фурье-периоде.

25 Базисные вейвлет функции 4 различных базисных вейвлета: Морли, Поля и производные гауссианы 2-го и 6-го порядка во временной области. Эти функции имеют нулевое среднее и локализованы как во временной, так и в частотной области. Функции имеют действительную и мнимую части, или амплитуду и фазу. Они определяют вейвлет спектр мощности. Узкие функции дают хорошее разрешение по времени, широкие – по частоте.

26 Фурье и вейвлет спектры используются далее для анализа быстрой гамма вспышки. Данные получены на борту космического аппарата BATSE под номером

27 Sample plot of time-tagged event data, BATSE trigger number binned to 2.5-millisecond resolution. Four channels from 20 to 300 keV.

28 A flare itself has an apparent duration of about 30 milliseconds. Two lower Figures show the power spectral density estimate of the photon flux in an energy channel keV using a periodogram with a rectangular and a Tukey window, respectively. The power spectra indicate clearly an oscillation feature at 175 Hz during the outburst phase. The upper dashed line is the 99 % confidence spectrum.

29 Sample plot of time-tagged event data BATSE trigger number (left plot). The wavelet power spectra. Fourier period vs. time (right plot). Good algorithm what sees more than eye.

30 Harmonic at ~ 175 Hz in the Fourier power spectrum. Fourier period of ~ 1/175 = sec in the wavelet power spectrum. Good coincidence between the Fourier and wavelet.

31 Как и в случае Фурье спектров можно найти уровень доверительной вероятности и для вейвлет спектра. Для белых шумов локальный вейвлет спектр имеет распределение Отсюда легко установить, например, 99.5% границу для пиков в вейвлет спектрах. Так квантиль для хи- квадрат распределения с двумя степенями свободы будет равен K 1-α = - 2 ln (1- α) Его значение для доверительной вероятности 99.5% равно 10.6.

32 Белые контуры соответствуют 95% доверительной вероятности на вейвлет спектрах. Одновременные вспышки колебаний в трех спектральных каналах доказывают их реальность.

33 Резюме Статистические обнаружители сигналов, основанные на энтропии Шеннона или других функциях энтропии способны выдавать сообщение о наличии полезного сигнала. Энтропийные критерии могут измерять энтропию потока данных в реальном времени и оценивать информационные свойства стационарных случайных и периодических сигналов. Классический анализ Фурье – непревзойденный инструмент для обнаружения стационарных сигналов во временных рядах данных, оценивания их амплитудных и частотных характеристик. Вейвлет анализ – мощный инструмент для анализа временных рядов, содержащих вариации на многих временных масштабах или демонстрирующих изменения в мощности сигналов.

34 Приложение СВЧ колебания во вспышке BATSE Помимо гармоники 175 Гц спектральный анализ позволяет выявить в гамма вспышке СВЧ колебания в диапазоне сотен килогерц. Эти колебания детектируются только в районе максимума вспышки и наблюдаются в течение примерно двадцати миллисекунд. Природа столь высокочастотных колебаний остается полностью неясной. СВЧ колебания невозможно увидеть глазом, поскольку из-за дефицита квантов нельзя построить кривую блеска. Однако для спектрального анализа разреженные потоки квантов не являются помехой. Время регистрации квантов определяется с точностью около 1 микросекунды. Это позволяет детектировать переменность в диапазоне частот до 500 килогерц.

35 СВЧ колебания во вспышке BATSE (Пример микросекундной астрофизики) СВЧ колебания детектируются в каналах и keV. Фурье спектры, слаженные с применением прямоугольного окна и окна Тьюки, демонстрируют колебания на частотах около 190 и 310 Кгц. Данный пример показывает, что фактически достаточно всего несколько сотен квантов, чтобы обнаружить СВЧ переменность методами быстрой фотометрии.

36 Моделирование СВЧ переменности A=10; N=1e5; Lambda = *(1+A*cos(2*pi/10.52*1:N)).^2 +A*cos(2*pi/6.46*(1:N)).^2); for i =1:N; v5(i) = poissrnd (Lambda(i),1,1); end; % SIGNAL Ampl. = 10, Freq. = 190 & 310 kHz Моделирование позволяет сделать вывод, что разреженный модулированный поток квантов со средней интенсивностью, равной наблюдаемой (всего квантов за время выборки сигнала), демонстрирует колебания с частотой 190 и 310 Кгц (левый Рис.). Правый Рис. - шумовой, сигнал в выборке отсутствует.

37 Gamma-ray burst models proposed are compact mergers and massive stellar collapses. Numerical simulation of the coalescence of two equal mass black holes. About 1% of the total mass energy will emerge as gravitational waves and gamma-rays during the final stages of the collision (ring-down phase) lasting some milliseconds only.

38 Gravitational Radiation from Gamma-Ray Burst Progenitors Black hole – neutron star scenario: thin lines and curve: 3 and 1.4 M ; thick lines and curve: 12 and 1.4 M. In-spiral (solid line), merger (dashed dotted line) and ring-down (spike). Collapsar: blob merger (dashed dotted line) and ring- down (solid spike). S. Kobayashi and P. Mezaros, arXiv:astro-ph/ v2 11 Feb 2003 Мы можем заключить, что гармоника 175 Гц во время вспышки BATSE могла возникнуть в процессе коалесценции черной дыры и нейтронной звезды.