Какова частота «Найквиста» для дискретизации производной сигнала?


12

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

Вопрос: Учитывая, что напряжение на конденсаторе ограничено полосой пропускания, и я делаю выборку производной от этого напряжения, какова минимальная частота дискретизации, необходимая для идеальной реконструкции сигнала напряжения из текущих выборок?

Если на этот вопрос нет постоянного ответа, было бы полезно все, что могло бы указать мне правильное направление. Спасибо заранее за любую помощь!!


1
Вы хотите «идеально восстановить» исходный сигнал из семплов? Что ты имеешь в виду?
Эллиот Олдерсон

1
Частота Найквиста в два раза выше самой высокой частоты в исходном сигнале.
Питер Карлсен

@Dweerberkitty, как упоминал Дэйв, сигнал - это просто сигнал :). Если серьезно, если вы используете системы реального измерения, то могут быть задержки, которые будут влиять на вашу производную работу. Таким образом, если вы учитываете их (если повезет, если система проста), вы можете аналитически вывести необходимый период выборки.
Раджа

«Напряжение на конденсаторе ограничено пропускной способностью». Почему?
Родриго де Азеведо

@RodrigodeAzevedo, это всего лишь предположение для упрощения постановки задачи. На самом деле, это не ограниченная полоса пропускания, но интересующий диапазон частот четко определен в этой проблеме. Благодарность!
ВЕНДЕРН

Ответы:


19

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

Таким образом, коэффициент Найквиста для производной такой же, как и для исходного сигнала.


3
Это правда в идеальном мире, в котором сигналы с абсолютно ограниченной полосой пропускания, идеальные фильтры нижних частот и вообще отсутствуют тепловые шумы.
Родриго де Азеведо

Весь баланс SNR изменяется. Небольшой высокочастотный компонент, который может иметь псевдоним, но мало что может сделать из-за своего размера, может стать значительным монстром, вызывающим уверенность в большой-низкой частоте-на-выборке.
Скотт Сейдман

-1

Взятие производной умножает преобразование на s, которое эффективно поворачивает график величины против часовой стрелки. Таким образом, вполне могут быть высокочастотные составляющие в производной. Более краткий способ выразить это заключается в том, что деривация усиливает высокочастотное содержимое.

1s+1

 bode(tf(1, [ 1 1 ])) 

введите описание изображения здесь

ss+1

bode(tf([1 0], [ 1 1 ])) 

введите описание изображения здесь

Производная в этом случае явно имеет более высокие частотные составляющие. Возможно, более правильно, он имеет гораздо большие высокочастотные компоненты, чем непроизводные. Можно с некоторой уверенностью выбрать выборку первого сигнала со скоростью 200 рад / с, поскольку при скорости Найквиста энергия очень мала, но сглаживание будет существенным, если вы производите выборку производной с той же скоростью.

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

РЕДАКТИРОВАТЬ: В ответ на downvote, я забью этот дом на конкретном примере. Позвольте мне взять синусоидальную волну и добавить к ней случайный нормальный шум (одну десятую величины синусоидальной волны).

введите описание изображения здесь

БПФ этого сигнала:

введите описание изображения здесь

Теперь позвольте мне взять производную сигнала: введите описание изображения здесь

и БПФ производной

введите описание изображения здесь

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


2
Я не уверен, что вы думаете, что вы планируете здесь, но это не ограниченные полосой сигналы.
Дэйв Твид

Преобразование Фурье сигнала и преобразование Фурье его производной.
Скотт Сейдман

Что это за язык?
Дэйв Твид

1
Ах. В этом случае tf()не представляет сигнал, он представляет передаточную функцию. Определенно не ограниченный диапазон.
Дэйв Твид

1
Вы все еще упускаете момент, когда сигнал ограничен. Вы добавляете в сигнал не ограниченный диапазон, чтобы выразить свою точку зрения, что выходит за рамки вопроса. Да, это практическое соображение, но вопрос (на мой взгляд) теоретический.
Дэйв Твид

-2

Ты не можешь

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

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

У вас также есть проблема в том, что, хотя ваши текущие выборочные измерения ограничены по Найквисту, фактический ток через конденсатор может не быть. Если вы не можете гарантировать, что ток через конденсатор имеет жесткий фильтр нижних частот где-то ниже предела Найквиста, вы никогда не сможете измерить ток достаточно точно, чтобы воспроизвести напряжение. Мне нужно прояснить, что это на самом деле математически невозможно, потому что это потребует частоты дискретизации бесконечности.

Но если вы знаете начальное напряжение и если фактический ток через конденсатор соответственно фильтруется по нижним частотам, то DaveTweed прав, что предел Найквиста для интеграла такой же, как и для выборочных данных.


Я не понимаю, почему вам нужно проводить различие между фактическим током через конденсатор и измеренным значением, ограниченным полосой. Что такого волшебного в этой ситуации, что хорошо известная линейность производных, фильтров и интеграции больше не применяется?
труба

@pipe Одним словом, выборка. Предположим, мы сэмплируем при 1 кГц. Теперь предположим, что у нас есть текущий всплеск 0,5 мс. Выборочная версия никогда не увидит всплеска, но фактическое напряжение на конденсаторе, безусловно, будет. Тогда у вас есть остаточные ошибки между каждой формой цифровой интеграции и фактическим значением. И я даже не начал ни с чего связанного с разрешением, которое является еще одной червей.
Грэм

Но энергия в этом импульсе будет распространяться на полосы, которые увидит сэмплер . Например: последовательность импульсов с очень короткими импульсами, после ограничения полосы, составит слегка повышенный уровень постоянного тока. Область вашего импульса останется прежней, и интеграция версии с ограниченной полосой частот приведет к тому же результату.
труба
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.