реклама на сайте
подробности

 
 
4 страниц V  < 1 2 3 4 >  
Reply to this topicStart new topic
> Уточнение момента прихода сигнала
alex_os
сообщение Sep 7 2009, 10:30
Сообщение #31


Знающий
****

Группа: Свой
Сообщений: 521
Регистрация: 12-05-06
Пользователь №: 17 030



Цитата(GetSmart @ Sep 7 2009, 13:53) *
Это по-вашему, гипотетически. Ну а практически - реально меньше.

Практическое подтверждение smile.gif.

Девиация положения максимума
std_max = 0.1613
Девиация положения точки пересечения порога 0.5
std_edge = 2.5813
Так что природу не обманешь...


Код
x = zeros(1000,1);
x = zeros(1000,1);

for i=1:10:1000
    if(randn(1,1)>0)     x(i:i+9) = 1;
    else        x(i:i+9) = -1;    end
end

pos_max  = zeros(1,1000);
pos_edge = zeros(1,1000);
plot(y);

for j=1:1000
        y=xcorr(x,x+3*randn(1000,1));        
        [v, m] = max(y);
        I=find(y>500);
        pos_max(j)  = m;
        pos_edge(j) = I(1);
end        

std_max  = std(pos_max)
std_edge = std(pos_edge)


--------------------
ну не художники мы...
Go to the top of the page
 
+Quote Post
fontp
сообщение Sep 7 2009, 10:42
Сообщение #32


Эксперт
*****

Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183



Цитата(alex_os @ Sep 7 2009, 14:30) *
Практическое подтверждение smile.gif.

Девиация положения максимума
std_max = 0.1613
Девиация положения точки пересечения порога 0.5
std_edge = 2.5813
Так что природу не обманешь...


Зависит от того как сформированы импульсы. Есть приложения, где корреляционная функция треугольная или даже может и трапеция и применяют геометрический первый момент (центр масс) функции корреляции импульса или даже медиану (вес справа = весу слева), как самое устойчивое решение - для выявления максимума. Почему обязательно уровень пересечения 0.5? Есть разные способы максимум искать и все они правильные, и не всегда подогнать параболу лучший ))
Go to the top of the page
 
+Quote Post
GetSmart
сообщение Sep 7 2009, 11:48
Сообщение #33


.
******

Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753



Цитата(petrov @ Sep 7 2009, 15:45) *
максимум на выходе согласованного фильтра это и есть самый луший обнаружитель

А шо за фильтр, который умеет фильтровать треугольники? Все фильтры предназначены для синусов. А треугольники - это уже из области геометрии.


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post
fontp
сообщение Sep 7 2009, 11:52
Сообщение #34


Эксперт
*****

Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183



Цитата(GetSmart @ Sep 7 2009, 15:48) *
А шо за фильтр, который умеет фильтровать треугольники? Все фильтры предназначены для синусов. А треугольники - это уже из области геометрии.


Так всякий сигнал обычно формируется в передатчике посредством фильтрации.
И в модемах и в радиолокаторах. Вот о чём речь. В функции корреляции выходит скругление.

Треугольник - это результат корреляции двух идеально прямоугольных импульсов. Идеально прямоугольные импульсы бывают только в идеальных устройствах, например в мозге )) Ultra Wide Band
Go to the top of the page
 
+Quote Post
GetSmart
сообщение Sep 7 2009, 12:01
Сообщение #35


.
******

Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753



Цитата(fontp @ Sep 7 2009, 17:52) *
Так всякий сигнал обычно формируется в передатчике посредством фильтрации.
И в модемах и в радиолокаторах. Вот о чём речь. В функции корреляции выходит скругление.

Я так понял, что petrov предлагал фильтровать результат корреляции чтобы на отфильтрованном результате проще/качественней найти пик. Но синусоидальный фильтр не предназначен для этого. Кроме того, фильтр должен иметь частоту среза соизмеримую с длительностью опорного сигнала. При этом треугольник реально исказится, т.к. к нему применили "недопустимую операцию" biggrin.gif

Сообщение отредактировал GetSmart - Sep 7 2009, 12:06


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post
thermit
сообщение Sep 7 2009, 12:03
Сообщение #36


Знающий
****

Группа: Участник
Сообщений: 781
Регистрация: 3-08-09
Пользователь №: 51 730



Цитата
GetSmart:
Все фильтры предназначены для синусов


Отакот век живи - век учись...
Go to the top of the page
 
+Quote Post
petrov
сообщение Sep 7 2009, 12:06
Сообщение #37


Гуру
******

Группа: Свой
Сообщений: 2 220
Регистрация: 21-10-04
Из: Balakhna
Пользователь №: 937



Цитата(GetSmart @ Sep 7 2009, 16:01) *
Я так понял, что petrov предлагал фильтровать результат кореляции чтобы на отфильтрованном результате проще/качественней найти пик. Но синусоидальный фильтр не предназначен для этого.


Ничего такого я не предлагал. Имеете понятие о согласованной фильтрации? На выходе согласованного фильтра для вашей последовательности псевдослучайной ищете максимум, всё.
Go to the top of the page
 
+Quote Post
DRUID3
сообщение Sep 7 2009, 13:39
Сообщение #38


山伏
*****

Группа: Свой
Сообщений: 1 827
Регистрация: 3-08-06
Из: Kyyiv
Пользователь №: 19 294



Цитата(DMax @ Sep 6 2009, 11:47) *
Ну хорошо, расскажите нам тогда в каком месте по-вашему я наступил на горло Котельникову или же Найквисту, чтобы они вертелись в гробах?

ОК. Если призадумаетесь, то теорема Котельникова/Найквиста устанавливает ограничение на скорость изменения функции относительно количества отсчетов на определенном интервале... Частотный критерий - это дух той эпохи - гармонический подход ко всякого рода явлениям. Надеюсь Вы понимаете, что никаких частот(частостей, вейвлетов, чирплетов, разложений по Тейлору и т.д) "внутри" сигнала нет. Вот именно поэтому то и работает стробоскопический осциллограф и оцифровка радиосигнала 1.0..1.1 GHz с частотой 200 MHz.

Далее Вы сами же приходите к противоречию. К стационарному сигналу Вы подмешиваете априори неизвестный сигнал (в виде шума) получая при этом непредсказуемую скорость изменения и задаете вопрос следующего плана:
"ребята сколько будет 2+Y если раньше было 18, а еще раньше -7".

Так вот это Y Вы можете только прогнозировать но никак не знать заранее. У шума есть 2-е характеристики - мощность на интервале и функция распределения(количество тех или иных значений отсчетов) на интервале - и никакая статистика Вам не скажет чему равен такой-то по-порядку отсчет - иначе это уже не шум будет.

Далее вообще супер... Поздравляю, Вы - на грани открытия(чур я соаФФтАр wink.gif ):

Установка которой Вы пытаетесь "угадать" момент прихода радиоимпульса является отличнейшим прибором для изучения статистических свойств тех или иных шумов/помех. Эдакий "идеологический туннельный микроскоп" biggrin.gif - задаете правило аппроксимации/предсказания - и на довольно длительном интервале собираете статистику имея однозначный параметр - коэффициент корреляции. Так практически можно "снять"(как ВАХ для лампы, транзистора или полевика) любую статистистическую функцию распределения. Кажется, а зачем можно и непосредственно подсчитать? Но Вашим способом удобно делать быстрые шаблонные тесты - выдвинули гипотезу, ага пошло разбегание, новая гипотеза т т.д. smile.gif


Цитата(fontp @ Sep 7 2009, 13:13) *
По первой части вопроса, есть теоретически достижимая оценка максимального правдоподобия
Крамера-Рао
http://www.cplire.ru/joined/mac/lection5/text.html
http://www.ndt.net/article/ultragarsas/63-...%20Dumbrava.pdf

Грубо говоря стандартное отклонение позиционирования во времени не может быть для одиночного импульса лучше чем
T / sqrt(SNR)
, T - продолжительность импульса, SNR - сигнал/шум
Для псевдослучайной последовательности длиной N должно получиться что-то ~ T/sqrt(N*SNR)

Поищите Крамера-Рао или CRLB в Гугле. Это в экзаменах есть ))

Цитата(fontp @ Sep 7 2009, 13:42) *
Зависит от того как сформированы импульсы. Есть приложения, где корреляционная функция треугольная или даже может и трапеция и применяют геометрический первый момент (центр масс) функции корреляции импульса или даже медиану (вес справа = весу слева), как самое устойчивое решение - для выявления максимума. Почему обязательно уровень пересечения 0.5? Есть разные способы максимум искать и все они правильные, и не всегда подогнать параболу лучший ))

Позволил себе добавить эти цитаты как под другим углом раскрывающие мою позицию.

Цитата(GetSmart @ Sep 6 2009, 11:54) *
А про Котельникова не переживайте. Его ТК здесь никаким боком не требуется.

biggrin.gif Ай-яй-яй. Такой большой, а обманывает младших!


Цитата(GetSmart @ Sep 6 2009, 11:54) *
Не, я простой гений smile.gif

+
Цитата(GetSmart @ Sep 7 2009, 14:48) *
А шо за фильтр, который умеет фильтровать треугольники? Все фильтры предназначены для синусов. А треугольники - это уже из области геометрии.

А продижи местного разлива часом не выступал научным консультантом в проекте "Булава"? smile.gif

P.S.: Еще в гробу крутится и Шенон. Полоса+соотношение шума и сигнала - вероятность достоверного приема - это все к нему wink.gif ...


--------------------
Нас помнят пока мы мешаем другим...
//--------------------------------------------------------
Хороший блатной - мертвый...
//--------------------------------------------------------
Нет старик, это те дроиды которых я ищу...
Go to the top of the page
 
+Quote Post
GetSmart
сообщение Sep 7 2009, 14:17
Сообщение #39


.
******

Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753



Цитата(DRUID3 @ Sep 7 2009, 19:39) *
Частотный критерий - это дух той эпохи - гармонический подход ко всякого рода явлениям. Надеюсь Вы понимаете, что никаких частот(частостей, вейвлетов, чирплетов, разложений по Тейлору и т.д) "внутри" сигнала нет.
Во-первых, корреляция не раскладывает сигнал на спектры, поэтому как я и сказал - проблем с ТК нет.

Цитата
Так вот это Y Вы можете только прогнозировать но никак не знать заранее. У шума есть 2-е характеристики - мощность на интервале и функция распределения(количество тех или иных значений отсчетов) на интервале - и никакая статистика Вам не скажет чему равен такой-то по-порядку отсчет - иначе это уже не шум будет.

По-вашему существование "шумоподавителей" можно считать мистикой.
---
Во-вторых, шумы имеют свойства, отличные от полезных сигналов и их можно отделять от оных. Кроме того, есстественные шумы имеют реально ограниченную амплитуду.


По-поводу корреляции. Жаль, что ширина её завязана на ширину опорной последовательности, т.к при увеличении ширины интеграла уровень шумов в результате всегда пропорционально падает. Допустим, что шумы распределяются одинаково на треугольнике кор.функции (я сам не проверял). чтобы сгладить эти шумы вполне естественно взять как можно больше точек в окресностях пика треугольника. На двух или трёх точках разброс будет куда больший.


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post
thermit
сообщение Sep 7 2009, 14:38
Сообщение #40


Знающий
****

Группа: Участник
Сообщений: 781
Регистрация: 3-08-09
Пользователь №: 51 730



Цитата
GetSmart:
Во-первых, корреляция не раскладывает сигнал на спектры, поэтому как я и сказал - проблем с ТК нет.


Ы-ы-ы-ы-ы...
А как же широко известная в узких кругах теорема винера нашего, хинчина?
Go to the top of the page
 
+Quote Post
GetSmart
сообщение Sep 7 2009, 14:57
Сообщение #41


.
******

Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753



Это как в том анекдоте - "во, валит!..." :-)))

Сообщение отредактировал GetSmart - Sep 7 2009, 14:58


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post
DRUID3
сообщение Sep 7 2009, 15:22
Сообщение #42


山伏
*****

Группа: Свой
Сообщений: 1 827
Регистрация: 3-08-06
Из: Kyyiv
Пользователь №: 19 294



Цитата(GetSmart @ Sep 7 2009, 17:17) *
Во-первых, корреляция не раскладывает сигнал на спектры, поэтому как я и сказал - проблем с ТК нет.

Повторяю специально для местных гениев. Спектры действительно нипричем. А вот скорость изменения относительно количества взятых отсчетов - очень даже причем. И эту скорость можно измерять частотным критеием biggrin.gif о чем и писали Найквист/Котельников ...гы... Так вот подмешивая случайным образом дополнительные веса для каждого отсчета мы как бы размазываем детерминированную картинку нашего стробоскопического осциллографа намного быстрее изменяющейся функцией шума. И все, уже ничего точно между отсчетов мы никогда не зафиксируем. Остается применить некий статистический фильтр, который из определенного количества измерений будет давать тот или иной процент адекватный в некотором диапазоне значений... Преднамеренно получить единичное точное значение - исключено.

Цитата(GetSmart @ Sep 7 2009, 17:17) *
По-вашему существование "шумоподавителей" можно считать мистикой.

По-моему мистикой можно считать существование универсальных и 100%-ных шумоподавителей, предсказателей, прогнозов погоды и т.д smile.gif Тем лучше будет изучена/предсказана природа помехи, тем лучше ее подавят. Все эти методы предсказателей/интерполяций предлагаемые автору так или иначе хорошо себя зарекомендовали при определенных условиях...

Цитата(GetSmart @ Sep 7 2009, 17:17) *
Во-вторых, шумы имеют свойства, отличные от полезных сигналов и их можно отделять от оных. Кроме того, есстественные шумы имеют реально ограниченную амплитуду.

Ага... вот и я про то...


--------------------
Нас помнят пока мы мешаем другим...
//--------------------------------------------------------
Хороший блатной - мертвый...
//--------------------------------------------------------
Нет старик, это те дроиды которых я ищу...
Go to the top of the page
 
+Quote Post
GetSmart
сообщение Sep 7 2009, 15:42
Сообщение #43


.
******

Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753



Цитата(DRUID3 @ Sep 7 2009, 21:22) *
А вот скорость изменения относительно количества взятых отсчетов - очень даже причем. И эту скорость можно измерять частотным критеием biggrin.gif о чем и писали Найквист/Котельников ...гы... Так вот подмешивая случайным образом дополнительные веса для каждого отсчета мы как бы размазываем детерминированную картинку нашего стробоскопического осциллографа намного быстрее изменяющейся функцией шума. И все, уже ничего точно между отсчетов мы никогда не зафиксируем. Остается применить некий статистический фильтр, который из определенного количества измерений будет давать тот или иной процент адекватный в некотором диапазоне значений... Преднамеренно получить единичное точное значение - исключено.

Котельников/Найквист писали, что в дискретном сигнале есть ортогональные частоты. И о том, что сигнал можно полностью восстановить если он сверху ограничен по частоте равной половине дискретизации (пусть чуть меньше). Со вторым я никогда не спорил. А первое тут нипричём. Шум мы не подмешиваем, он уже подмешан "природой". Перед оцифровкой сигнал ессно отфильтрован по критерию Котельникова. Остальные промежуточные точки сигнала вместе с шумом находятся фильтрацией, причём идеально. Какие проблемы?!
Тут проблема только в том, что сам шум даёт корреляцию с сигналом, тем самым искажая треугольник.

Сообщение отредактировал GetSmart - Sep 7 2009, 15:49


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post
DMax
сообщение Sep 7 2009, 15:57
Сообщение #44


Частый гость
**

Группа: Свой
Сообщений: 161
Регистрация: 19-01-08
Из: Питер
Пользователь №: 34 234



Цитата(DRUID3 @ Sep 7 2009, 17:39) *
ОК. Если призадумаетесь, то теорема Котельникова/Найквиста устанавливает ограничение на скорость изменения функции относительно количества отсчетов на определенном интервале... Частотный критерий - это дух той эпохи - гармонический подход ко всякого рода явлениям. Надеюсь Вы понимаете, что никаких частот(частостей, вейвлетов, чирплетов, разложений по Тейлору и т.д) "внутри" сигнала нет. Вот именно поэтому то и работает стробоскопический осциллограф и оцифровка радиосигнала 1.0..1.1 GHz с частотой 200 MHz.


Ёп!!! А мужики-то не знают! Лепят уже по второму десятку лет гетеродины и прочую аналоговую лабуду... побегу им ща рассказывать.

Цитата(DRUID3 @ Sep 7 2009, 17:39) *
Далее Вы сами же приходите к противоречию. К стационарному сигналу Вы подмешиваете априори неизвестный сигнал (в виде шума) получая при этом непредсказуемую скорость изменения и задаете вопрос следующего плана:
"ребята сколько будет 2+Y если раньше было 18, а еще раньше -7".

Так вот это Y Вы можете только прогнозировать но никак не знать заранее. У шума есть 2-е характеристики - мощность на интервале и функция распределения(количество тех или иных значений отсчетов) на интервале - и никакая статистика Вам не скажет чему равен такой-то по-порядку отсчет - иначе это уже не шум будет.

Далее вообще супер... Поздравляю, Вы - на грани открытия(чур я соаФФтАр wink.gif ):

Установка которой Вы пытаетесь "угадать" момент прихода радиоимпульса является отличнейшим прибором для изучения статистических свойств тех или иных шумов/помех. Эдакий "идеологический туннельный микроскоп" biggrin.gif - задаете правило аппроксимации/предсказания - и на довольно длительном интервале собираете статистику имея однозначный параметр - коэффициент корреляции. Так практически можно "снять"(как ВАХ для лампы, транзистора или полевика) любую статистистическую функцию распределения. Кажется, а зачем можно и непосредственно подсчитать? Но Вашим способом удобно делать быстрые шаблонные тесты - выдвинули гипотезу, ага пошло разбегание, новая гипотеза т т.д. smile.gif


Сколько слов, но нет банального ответа на первоначально заданный вопрос: "В каком месте, по-вашему, я наступил на горло Котельникову и/или же Найквисту, чтобы они вертелись в гробах?". Пока пахнет демагогией.
Go to the top of the page
 
+Quote Post
bahurin
сообщение Sep 7 2009, 16:16
Сообщение #45


Местный
***

Группа: Участник
Сообщений: 240
Регистрация: 20-09-08
Пользователь №: 40 347



Думаю, что если шум аддитивный с нулевым среднем, то можно получить любую точность оценки временного положения корреляционного пика если частота дискретизации была выбрана по Найквисту. Все упирается в возможную длительность наблюдения. При бесконечном наблюдении можно получить истинную оценку времени задержки.
Go to the top of the page
 
+Quote Post

4 страниц V  < 1 2 3 4 >
Reply to this topicStart new topic
2 чел. читают эту тему (гостей: 2, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 15th June 2025 - 18:12
Рейтинг@Mail.ru


Страница сгенерированна за 0.01528 секунд с 7
ELECTRONIX ©2004-2016