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

 
 
> Увеличение разрешения по частоте, засчёт накопления
EKirshin
сообщение Jul 30 2008, 15:01
Сообщение #1


Участник
*

Группа: Свой
Сообщений: 41
Регистрация: 2-10-06
Из: Москва
Пользователь №: 20 904



Здравствуйте!

Такая задачка.

Есть стационарный сигнал: шум + синусоида.

Из этого сигнала периодически набирается N отсчётов.
Вопрос: можно ли, имея M кусков этого сигнала использовать их для получения спектра сигнала размерностью M*N?

Какой будет эффект, если на каждый из этих кусков наложить окно, "склеить" эти куски (получится сигнал длиной M*N), затем сделать fft на M*N точек?
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
GetSmart
сообщение Jul 31 2008, 10:40
Сообщение #2


.
******

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



Цитата(fontp)
А разрешение здесь не при чём. Маклауд использовал мощнейшую априорную информацию о сингулярности спектра сигнала. А разрешение определено только по отношению к сложным спектрам, в которых как минимум два горба
Хватит повторять друг за другом характеристику сигнала, не имеющую к теме отношения. Для сигнала, состоящего из одного тона она бессмыслена.

Цитата(fontp)
Если куски идут подряд то проблемы нет
(отказаться от когерентной обработки может только недостаток быстродействия или памяти, либо нестабильность частоты)

Если между кусками произвольные, но известные расстояния, то тоже проблемы нет - если ввести коррекцию по фазе для каждого куска и накоплять результаты ДПФ
S = Сумма(ДПФ(i)*exp(j*ф(i))) Принцип правильный, но формула неверная в том смысле, что фаза Ф ещё зависит и от частот в ДПФ, двойные суммы, это мне не написать нормально. В, общем, вы поняли


Если же расстояния не известны, то остаётся только некогерентное накопление. Собственно разница когерентного накопления от некогерентного собственно в том, что в одном случае накопляется результат ДПФ с правильной фазой, а в другом его квадрат модуля
Если куски идут подряд, то это уже не куски, а одно большое окно для FFT.

Если между кусками произвольные, но известные расстояния, то ничего у вас, дорогой мой, не получится. Т.к. для этой "операции" нужно знать ответ на вопрос: какая частота у тона, причём с идеальной точностью.

Так что остаётся метод один - некоггерентное накопление.


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post
fontp
сообщение Jul 31 2008, 10:51
Сообщение #3


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

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



Цитата(GetSmart @ Jul 31 2008, 14:40) *
Хватит повторять друг за другом характеристику сигнала, не имеющую к теме отношения. Для сигнала, состоящего из одного тона она бессмыслена.

Если куски идут подряд, то это уже не куски, а одно большое окно для FFT.

Если между кусками произвольные, но известные расстояния, то ничего у вас, дорогой мой, не получится. Т.к. для этой "операции" нужно знать ответ на вопрос: какая частота у тона, причём с идеальной точностью.

Так что остаётся метод один - некоггерентное накопление.


Я уже было ушёл, но вы опять... Отвечу коротко
1. Разрешение - это другое, да
2 Одно большое FFT - отлично
3.

А я уже зачеркнул ту формулу. Её просто было трудно написать - там двумерная сумма и всё такое.
Я не могу и не хочу даже учиться писать в форумах серьёзные двумерные суммы.
Но! Это не меняет ничего в принципе.
Принцип-то верный. Мне не нужно знать частоту синусоиды. Я беру много синусод - всю линейку ДПФ - в качестве эталонов. Для каждой синусоиды частоты F я знаю набег фазы между блоками - это
2*pi*dt*F. На каждом блоке я корректирую фазу. Что мне помешает считать "сплошной ДПФ"?
Можете сказать, что это не ДПФ. Может быть, назовите как хотите. Только это согласованая когерентная фильтрация по всем блокам

Считается скалярное произведение сигнальной синусоиды с синусоидой эталонной как при ДПФ. Но некоторые точки недоступны, их не считаем. Но положение всех сигнальных точек выборок известны - их считаем. Имеем когерентное накопление.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- EKirshin   Увеличение разрешения по частоте   Jul 30 2008, 15:01
- - GetSmart   А повысить разрешение методом интерполяции не кати...   Jul 30 2008, 15:56
- - Stanislav   Цитата(EKirshin @ Jul 30 2008, 19:01) Ест...   Jul 30 2008, 16:43
|- - GetSmart   Цитата(Stanislav @ Jul 30 2008, 22:43) Сп...   Aug 1 2008, 08:56
- - GetSmart   Цитата(Stanislav)Скажите, каким образом интерполяц...   Jul 30 2008, 17:17
|- - fontp   Цитата(GetSmart @ Jul 30 2008, 21:17) fon...   Jul 30 2008, 17:57
|- - Stanislav   Цитата(GetSmart @ Jul 30 2008, 21:17) fon...   Jul 30 2008, 18:18
|- - fontp   Цитата(Stanislav @ Jul 30 2008, 22:18) По...   Jul 30 2008, 18:27
|- - Stanislav   Цитата(fontp @ Jul 30 2008, 22:27) Отчаст...   Jul 30 2008, 18:49
|- - fontp   Цитата(Stanislav @ Jul 30 2008, 22:49) ЗЫ...   Jul 31 2008, 07:26
||- - GetSmart   Цитата(fontp @ Jul 31 2008, 13:26) Когере...   Jul 31 2008, 08:01
|||- - fontp   Цитата(GetSmart @ Jul 31 2008, 12:01) Как...   Jul 31 2008, 08:19
||- - Stanislav   Цитата(fontp @ Jul 31 2008, 11:26) Запрос...   Jul 31 2008, 09:05
||- - fontp   Цитата(Stanislav @ Jul 31 2008, 13:05) За...   Jul 31 2008, 09:15
||- - Stanislav   Цитата(fontp @ Jul 31 2008, 13:15) Если к...   Jul 31 2008, 09:31
||- - fontp   Цитата(Stanislav @ Jul 31 2008, 13:31) Те...   Jul 31 2008, 09:33
|- - EKirshin   Цитата(Stanislav @ Jul 30 2008, 22:49) От...   Jul 31 2008, 16:14
|- - shf_05   Цитата(EKirshin @ Jul 31 2008, 22:14) Про...   Aug 1 2008, 07:07
||- - fontp   Цитата(shf_05 @ Aug 1 2008, 11:07) я бы д...   Aug 1 2008, 08:56
||- - EKirshin   Цитата(shf_05 @ Aug 1 2008, 11:07) так ли...   Aug 1 2008, 09:44
||- - fontp   Цитата(EKirshin @ Aug 1 2008, 13:44) Нет,...   Aug 1 2008, 09:51
|- - fontp   Цитата(EKirshin @ Jul 31 2008, 20:14) Про...   Aug 1 2008, 07:07
|- - Stanislav   Цитата(EKirshin @ Jul 31 2008, 20:14) ......   Aug 1 2008, 08:31
|- - blackfin   Цитата(Stanislav @ Aug 1 2008, 12:31) Вы ...   Aug 1 2008, 08:56
|- - Stanislav   Цитата(blackfin @ Aug 1 2008, 12:56) А чт...   Aug 1 2008, 10:00
|- - fontp   Цитата(Stanislav @ Aug 1 2008, 14:00) Пот...   Aug 1 2008, 10:03
- - GetSmart   Цитата(Stanislav)Пожалуй, в 1001-й раз приводить е...   Jul 30 2008, 18:59
|- - Stanislav   Цитата(GetSmart @ Jul 30 2008, 22:59) Ну ...   Jul 30 2008, 19:18
- - GetSmart   fontp, согласитесь, что неограниченное увеличение ...   Jul 31 2008, 08:42
|- - fontp   Цитата(GetSmart @ Jul 31 2008, 12:42) fon...   Jul 31 2008, 08:53
|- - GetSmart   Цитата(fontp @ Jul 31 2008, 14:53) Теорет...   Jul 31 2008, 09:05
||- - fontp   Цитата(GetSmart @ Jul 31 2008, 13:05) Вы ...   Jul 31 2008, 09:25
|- - Stanislav   Цитата(fontp @ Jul 31 2008, 12:53) Теорет...   Jul 31 2008, 09:23
- - RadioJunior   а у меня более банальный вопрос: а как вообще оцен...   Jul 31 2008, 09:34
|- - fontp   Цитата(RadioJunior @ Jul 31 2008, 13:34) ...   Jul 31 2008, 10:02
|- - 729   Цитата(GetSmart @ Jul 31 2008, 14:40) Есл...   Jul 31 2008, 10:58
- - GetSmart   Цитата(fontp)Можете сказать, что это не ДПФ. Может...   Jul 31 2008, 11:44
|- - 729   Цитата(GetSmart @ Jul 31 2008, 15:44) Во-...   Jul 31 2008, 11:57
|- - fontp   Цитата(GetSmart @ Jul 31 2008, 15:44) При...   Jul 31 2008, 11:59
- - GetSmart   Цитата(fontp)Кстати так и делают во многих приложе...   Jul 31 2008, 12:09
|- - 729   Цитата(GetSmart @ Jul 31 2008, 16:09) Доп...   Jul 31 2008, 12:22
|- - fontp   Цитата(GetSmart @ Jul 31 2008, 16:09) Ког...   Jul 31 2008, 12:29
|- - AndrewN   Цитата(fontp @ Jul 31 2008, 15:29) Так ра...   Jul 31 2008, 12:38
|- - blackfin   Цитата(AndrewN @ Jul 31 2008, 16:38) Гм. ...   Jul 31 2008, 12:42
|- - fontp   Цитата(AndrewN @ Jul 31 2008, 16:38) Гм. ...   Jul 31 2008, 12:44
- - GetSmart   Цитата(729)Непонятно. Сами частоты (6 штук) можете...   Jul 31 2008, 12:25
|- - 729   Цитата(GetSmart @ Jul 31 2008, 16:25) 99....   Jul 31 2008, 13:13
|- - GetSmart   Цитата(729 @ Jul 31 2008, 19:13) Да, но п...   Jul 31 2008, 13:18
|- - 729   Цитата(GetSmart @ Jul 31 2008, 17:18) Сме...   Jul 31 2008, 13:22
|- - GetSmart   Цитата(729 @ Jul 31 2008, 19:22) А вот из...   Jul 31 2008, 13:34
|- - 729   Цитата(GetSmart @ Jul 31 2008, 17:34) Поч...   Jul 31 2008, 13:50
- - GetSmart   Цитата(fontp)Фома Вы не верующий. Так работают все...   Jul 31 2008, 12:40
- - Stanislav   Цитата(fontp @ Aug 1 2008, 12:56) Если ра...   Aug 1 2008, 10:16
|- - fontp   Цитата(Stanislav @ Aug 1 2008, 14:16) Пра...   Aug 1 2008, 10:20
|- - Stanislav   Цитата(fontp @ Aug 1 2008, 14:20) Вот Вы ...   Aug 1 2008, 10:24
|- - fontp   Цитата(Stanislav @ Aug 1 2008, 14:24) Пос...   Aug 1 2008, 10:30
- - EKirshin   Цитата(Stanislav @ Aug 1 2008, 12:31) Пон...   Aug 1 2008, 10:39
|- - 729   Цитата(EKirshin @ Aug 1 2008, 14:39) Част...   Aug 1 2008, 11:07
|- - GetSmart   Цитата(EKirshin @ Aug 1 2008, 16:39) Соот...   Aug 1 2008, 11:17
||- - 729   Цитата(GetSmart @ Aug 1 2008, 15:17) Чем ...   Aug 1 2008, 11:24
|- - Stanislav   Цитата(EKirshin @ Aug 1 2008, 14:39) Част...   Aug 1 2008, 12:09
|- - GetSmart   Цитата(Stanislav @ Aug 1 2008, 18:09) Кто...   Aug 1 2008, 13:19
||- - Stanislav   Цитата(GetSmart @ Aug 1 2008, 17:19) Stan...   Aug 1 2008, 14:27
||- - GetSmart   Цитата(Stanislav @ Aug 1 2008, 20:27) Это...   Aug 1 2008, 15:26
||- - 729   Цитата(GetSmart @ Aug 1 2008, 19:26) Вы и...   Aug 1 2008, 15:42
|- - blackfin   Цитата(Stanislav @ Aug 1 2008, 16:09) Сей...   Aug 1 2008, 14:57
|- - Stanislav   Цитата(blackfin @ Aug 1 2008, 18:57) Мы ж...   Aug 1 2008, 19:11
|- - shasik   Я, конечно, дико извиняюсь, что без стука Почитайт...   Aug 1 2008, 19:19
- - EKirshin   Цитата(GetSmart @ Aug 1 2008, 15:17) Чем ...   Aug 1 2008, 15:25
|- - alex_os   Цитата(EKirshin @ Aug 1 2008, 19:25) Испо...   Aug 1 2008, 17:42
||- - 729   Цитата(alex_os @ Aug 1 2008, 21:42) Дык е...   Aug 1 2008, 17:55
|||- - alex_os   [quote name='729' date='Aug 1 2008, 21...   Aug 1 2008, 18:36
|||- - Stanislav   Цитата(729 @ Aug 1 2008, 21:55) Фазировка...   Aug 3 2008, 07:28
||- - Stanislav   Цитата(alex_os @ Aug 1 2008, 21:42) Дык е...   Aug 1 2008, 19:23
||- - alex_os   Цитата(Stanislav @ Aug 1 2008, 23:23) А к...   Aug 1 2008, 19:32
||- - Stanislav   Цитата(alex_os @ Aug 1 2008, 23:32) Дык я...   Aug 1 2008, 19:37
|- - 729   Цитата(EKirshin @ Aug 1 2008, 19:25) Испо...   Aug 2 2008, 20:36
- - makc   За грубость и оффтопик трое участников уже получи...   Aug 1 2008, 15:55
- - Stanislav   Вот обещанная "рыба".Кодfunction phase_a...   Aug 2 2008, 00:00
|- - fontp   Цитата(Stanislav @ Aug 2 2008, 04:00) - ...   Aug 2 2008, 09:01
|- - GetSmart   Цитата(fontp @ Aug 2 2008, 15:01) " ...   Aug 2 2008, 10:01
|- - Stanislav   Цитата(fontp @ Aug 2 2008, 13:01) Маклауд...   Aug 2 2008, 15:15
|- - alex_os   Цитата- модуль спектра последовательности двойной ...   Aug 2 2008, 17:17
||- - Stanislav   Цитата(alex_os @ Aug 2 2008, 21:17) Фаза ...   Aug 2 2008, 17:31
||- - alex_os   Цитата(Stanislav @ Aug 2 2008, 21:31) Это...   Aug 2 2008, 17:41
||- - Stanislav   Цитата(alex_os @ Aug 2 2008, 21:41) А отк...   Aug 2 2008, 17:45
|- - fontp   Цитата(Stanislav @ Aug 2 2008, 19:15) Стр...   Aug 3 2008, 08:09
|- - Stanislav   Цитата(fontp @ Aug 3 2008, 12:09) Слишком...   Aug 3 2008, 17:29
|- - GetSmart   Цитата(Stanislav @ Aug 3 2008, 23:29) - ...   Aug 4 2008, 07:42
- - Stanislav   Вот, немного поправил программу, и подбавил аддити...   Aug 3 2008, 07:03
- - EKirshin   Цитата(Stanislav @ Aug 1 2008, 23:37) 2 E...   Aug 3 2008, 07:10
|- - alexkok   Цитата(EKirshin @ Aug 3 2008, 10:10) Да, ...   Aug 5 2008, 16:35
- - alex_os   Сделал модельку аналогичную модели Stanislav'а...   Aug 4 2008, 11:40
- - Stanislav   С комплексным сигналом работать проще - гильберта ...   Aug 4 2008, 11:48
- - alexkok   Цитата(EKirshin @ Jul 30 2008, 18:01) Воп...   Aug 4 2008, 20:06
- - Stanislav   Цитата(alexkok @ Aug 5 2008, 00:06) Я дум...   Aug 4 2008, 22:36
- - alexkok   Цитата(Stanislav @ Aug 5 2008, 01:36) Это...   Aug 5 2008, 01:05
- - Stanislav   Цитата(alexkok @ Aug 5 2008, 05:05) Повыш...   Aug 5 2008, 06:51
2 страниц V   1 2 >


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

 


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


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