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

 
 
13 страниц V  « < 2 3 4 5 6 > »   
Reply to this topicStart new topic
> Измерение частоты основной гармоники (50 Гц) с точностью 0.01 Гц
blackfin
сообщение Sep 11 2015, 08:58
Сообщение #46


Гуру
******

Группа: Свой
Сообщений: 3 106
Регистрация: 18-04-05
Пользователь №: 4 261



Цитата(rudy_b @ Sep 11 2015, 03:30) *
Мгновенную частоту можно оценить только задним числом по достаточно большой выборке, но, при этом, само это понятие теряет смысл.

А ещё был такой чудак Гейзенберг, так тот вообще всем доказывал, что ошибка "∆F" при оценке измерении мгновенной частоты "F" равна бесконечности..

И даже формулу в свое оправдание привел:

∆E*∆t ≥ ћ/2,

и, с учетом того, что:

∆E == ћ*∆ω == ћ*2*pi*∆F,

такую же формулу для принципа неопределенности:

∆F*∆t ≥ 1/(4*pi).

Так что да, оценить мгновенную частоту сложно, а вот в какую-нить формулу запихнуть - легко!

biggrin.gif
Go to the top of the page
 
+Quote Post
rudy_b
сообщение Sep 11 2015, 14:26
Сообщение #47


Знающий
****

Группа: Свой
Сообщений: 888
Регистрация: 25-09-08
Из: Питер
Пользователь №: 40 458



Цитата(serjj @ Sep 11 2015, 10:52) *
f(t) = d/dt(arg(y(t))), простейший частотный дискриминатор, известный со времён старого доброго ЧМ, нет? )

Нет, вы не отличите изменение фазы от изменения амплитуды.

Цитата(serjj @ Sep 11 2015, 10:52) *
А ещё оказывается в современном ЦОС есть целый класс алгоритмов, т.н. DIFM, например
[attachment=95242:GE_47.pdf]
Вполне себе практические вещи делают.

Ага, но для этого нужно сделать выборку длиной не менее периода. И, только тогда, можно оценить "мгновенную" частоту.
Go to the top of the page
 
+Quote Post
serjj
сообщение Sep 11 2015, 14:32
Сообщение #48


Знающий
****

Группа: Участник
Сообщений: 527
Регистрация: 4-06-14
Из: Санкт-Петербург
Пользователь №: 81 866



Цитата
Нет, вы не отличите изменение фазы от изменения амплитуды.

Да вы что! biggrin.gif
Ну расскажите мне, где вы в формуле f(t) = d/dt(arg(y(t))) увидили изменение амплитуды? И как же это у нас FM радио тогда ловит?..
Цитата
Ага, но для этого нужно сделать выборку длиной не менее периода. И, только тогда, можно оценить "мгновенную" частоту.

Тогда она уже не будет мгновенной, вам не кажется? )
То о чём вы говорите, не мгновенная частота, а спектральная оценка. И для её определения со значительной точностью не обязательно иметь даже период колебания (в сущности всё определяется уровнем шумов и спектральным составом сигнала). Вот пруф, который я уже как-то приводил.

Сообщение отредактировал serjj - Sep 11 2015, 14:46
Go to the top of the page
 
+Quote Post
rudy_b
сообщение Sep 11 2015, 18:15
Сообщение #49


Знающий
****

Группа: Свой
Сообщений: 888
Регистрация: 25-09-08
Из: Питер
Пользователь №: 40 458



Цитата(serjj @ Sep 11 2015, 17:32) *
Да вы что! biggrin.gif
Ну расскажите мне, где вы в формуле f(t) = d/dt(arg(y(t))) увидили изменение амплитуды? И как же это у нас FM радио тогда ловит?..

Ну тогда уж сначала вы расскажите мне, как вы измеряете arg(y(t))?

Цитата(serjj @ Sep 11 2015, 17:32) *
Тогда она уже не будет мгновенной, вам не кажется? )
То о чём вы говорите, не мгновенная частота, а спектральная оценка. И для её определения со значительной точностью не обязательно иметь даже период колебания (в сущности всё определяется уровнем шумов и спектральным составом сигнала). Вот пруф, который я уже как-то приводил.

Как я уже говорил ранее, понятие "мгновенная частота" в метрологии отсутствует. Оно есть только в головах пишущих формулы теоретиков.

Попробуйте почитать книжку, указанную в посте
Цитата(rudy_b @ Sep 11 2015, 03:30) *
...Финк Л. М. Сигналы, помехи, ошибки…

Go to the top of the page
 
+Quote Post
serjj
сообщение Sep 11 2015, 20:08
Сообщение #50


Знающий
****

Группа: Участник
Сообщений: 527
Регистрация: 4-06-14
Из: Санкт-Петербург
Пользователь №: 81 866



Цитата
Ну тогда уж сначала вы расскажите мне, как вы измеряете arg(y(t))?

Арктангенсом, а вы?

Цитата
И, только тогда, можно оценить "мгновенную" частоту.
...
Как я уже говорил ранее, понятие "мгновенная частота" в метрологии отсутствует. Оно есть только в головах пишущих формулы теоретиков.

То у вас её можно оценить, то она только в "головах у теоретиков". Вам приводят конкретные статьи и примеры, а вы говорите какие-то общие фразы и не хотите их количественно подтвердить. Речь же о точном предмете )
Go to the top of the page
 
+Quote Post
rudy_b
сообщение Sep 12 2015, 16:40
Сообщение #51


Знающий
****

Группа: Свой
Сообщений: 888
Регистрация: 25-09-08
Из: Питер
Пользователь №: 40 458



Цитата(serjj @ Sep 11 2015, 23:08) *
Арктангенсом, а вы?

Что-то я не знаю приборов измеряющих арктангенс произвольного аналогового сигнала, может просветите?

Кроме того, арктангенс - разрывная функция и брать производную от него...

Несерьезно как-то, лучше книжку почитайте.
Go to the top of the page
 
+Quote Post
анатолий
сообщение Sep 13 2015, 19:44
Сообщение #52


Местный
***

Группа: Свой
Сообщений: 221
Регистрация: 10-12-05
Из: Украина
Пользователь №: 12 052



Цитата
Хочется выделить основную гармонику (цифровым фильтром) и измерить программно её частоту с точностью 0,01 Гц.

Что я могу:
С помощью ДПФ посчитать амплитуду основной гармоники.
С помощью БПФ посчитать амплитуды всех гармоник.
С помощью цифрового фильтра выделить полосу около 50 Гц и получить отфильтрованную выборку из 64-х точек.


Когда-то в 1993г. такую задачу решал для измерения токов в рельсах метро на контроллере 8051 на частотах 50, 75, 125 ... Гц.
Решалась просто: сигнал умножался на косинус, синус 50 Гц и накапливался как скользящее среднее. Результат - как фильтрация узким фильтром. Если сигнал строго =50 Гц. то выходной сигнал - постоянная, если 49 Гц - вых. сигнал 1 Гц, если 51 Гц - сигнал 1 Гц, но крутящийся в противоположную сторону. И т.п. Соответственно, если на выходе синусоида с периодом 100 сек. то сигнал 50,01 Гц.
Go to the top of the page
 
+Quote Post
Krys
сообщение Sep 14 2015, 04:10
Сообщение #53


Гуру
******

Группа: Свой
Сообщений: 2 002
Регистрация: 17-01-06
Из: Томск, Россия
Пользователь №: 13 271



Цитата(анатолий @ Sep 14 2015, 02:44) *
И т.п. Соответственно, если на выходе синусоида с периодом 100 сек. то сигнал 50,01 Гц.
дак это надо тогда 100 секунд ждать, пока хоть 1 период наловишь? А чтобы принять надёжное решение, нужен ещё и не один период? Видимо хотелось бы побыстрее...


--------------------
Зная себе цену, нужно ещё и пользоваться спросом...
Go to the top of the page
 
+Quote Post
serjj
сообщение Sep 14 2015, 07:04
Сообщение #54


Знающий
****

Группа: Участник
Сообщений: 527
Регистрация: 4-06-14
Из: Санкт-Петербург
Пользователь №: 81 866



Цитата
Что-то я не знаю приборов измеряющих арктангенс произвольного аналогового сигнала, может просветите?

Кроме того, арктангенс - разрывная функция и брать производную от него...


Просвещайтесь на здоровье
Go to the top of the page
 
+Quote Post
blackfin
сообщение Sep 14 2015, 07:13
Сообщение #55


Гуру
******

Группа: Свой
Сообщений: 3 106
Регистрация: 18-04-05
Пользователь №: 4 261



Цитата(serjj @ Sep 14 2015, 10:04) *

А чего это они в формулах (4) и (7) после второго знака равенства "arctan" выкинули?
Go to the top of the page
 
+Quote Post
serjj
сообщение Sep 14 2015, 07:23
Сообщение #56


Знающий
****

Группа: Участник
Сообщений: 527
Регистрация: 4-06-14
Из: Санкт-Петербург
Пользователь №: 81 866



Цитата
А чего это они в формулах (4) и (7) после второго знака равенства "arctan" выкинули?

Арктангенс равен аргументу комплексного числа, записанного в (3) как две отдельные квадратуры, обычное приравнивание. Если кому-то интересно, то вот о процессе угловой модуляции. Тогда всё сразу становится понятнее )

Сообщение отредактировал serjj - Sep 14 2015, 07:24
Go to the top of the page
 
+Quote Post
blackfin
сообщение Sep 14 2015, 07:40
Сообщение #57


Гуру
******

Группа: Свой
Сообщений: 3 106
Регистрация: 18-04-05
Пользователь №: 4 261



Цитата(serjj @ Sep 14 2015, 10:23) *
Если кому-то интересно, то вот о процессе угловой модуляции. Тогда всё сразу становится понятнее )

Да читали мы все это. И даже экзамены сдавали по радиолокации.. Просто очень давно.. wacko.gif
Go to the top of the page
 
+Quote Post
rudy_b
сообщение Sep 14 2015, 12:34
Сообщение #58


Знающий
****

Группа: Свой
Сообщений: 888
Регистрация: 25-09-08
Из: Питер
Пользователь №: 40 458



И где же тут прямое измерение "мгновенной" частоты? После фильтров ФНЧ?
Вы просто подтвердили мои слова:
Цитата
Мгновенную частоту можно оценить только задним числом по достаточно большой выборке, но, при этом, само это понятие теряет смысл.

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

И, как обычно, вы перепутали реально существующую и оцениваемую по достаточно длинной выборке разность фаз двух сигналов с мифической фазой произвольного ОДИНОЧНОГО аналогового сигнала.

Измеряется не фаза, а разность фаз, и за несколько периодов.

И, если вы прочитали всю приведенную ссылку, в ней, опять же, подтверждают мои слова о разрывности дифференцируемого вами арктангенса
Цитата
Кроме того, арктангенс - разрывная функция и брать производную от него...

Поэтому, в реальной обработке, арктангенс старательно исключают.

Расскажите мне, как определить фазу сигнала по одному отсчету АЦП? Ну ладно, учитывая нарисованные вами производные - пусть будет целых три последовательных отсчета с интервалом 1 мксек: -1В; 0В; +1В.
Ну, хотя бы, любимую "мгновенную" частоту определите?

Не получается? Тогда прекратите увиливать, лучше читайте книжки и думайте.
Go to the top of the page
 
+Quote Post
Fat Robot
сообщение Sep 14 2015, 13:34
Сообщение #59


ʕʘ̅͜ʘ̅ʔ
*****

Группа: Свой
Сообщений: 1 008
Регистрация: 3-05-05
Пользователь №: 4 691



Нет ничего проще: На основании именно этих полученных от АЦП отсчетов можно смело сказать, что на его входе присутствует синусоидальный сигнал с частотой 250 kHz, c амплитудой 1В. Центральному отсчету соответствует фаза 0. Мы считаем, что амплитуда сигнала нормирована под полный "раскрыв" АЦП +/-1В.

Цитата(rudy_b @ Sep 14 2015, 13:34) *
Расскажите мне, как определить фазу сигнала по одному отсчету АЦП? Ну ладно, учитывая нарисованные вами производные - пусть будет целых три последовательных отсчета с интервалом 1 мксек: -1В; 0В; +1В.
Ну, хотя бы, любимую "мгновенную" частоту определите?
Go to the top of the page
 
+Quote Post
thermit
сообщение Sep 14 2015, 13:45
Сообщение #60


Знающий
****

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



Цитата
rudy_b:
подтверждают мои слова о разрывности дифференцируемого вами арктангенса


Прошу прощения за нескромный вопрос. А в каком месте atan(x) разрывается?

Сообщение отредактировал thermit - Sep 14 2015, 13:52
Go to the top of the page
 
+Quote Post

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

 


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


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