|
Почему не работает теоремма Котельникова при F = 1/2Fs |
|
|
|
Oct 2 2010, 12:08
|

Местный
  
Группа: Участник
Сообщений: 228
Регистрация: 4-06-09
Пользователь №: 49 940

|
Всем привет. Я только приступаю к изучению ЦОС. Заметил такой факт- в теоремме написано, что частота дискретизации должны быть по крайней мере в 2 раза выше частоты спектра сигнала. Открываю матлаб, набираю там скриптик, где частота дискретизации ровно в 2 раза выше частоты сигнала Код t=0:1/8000:0.01;
f1=4000;
x = sin(2*pi*f1*t);
plot(t,x); Получаю следующий график. Вопрос- почему сигнала практически нет? (10^-14)
Эскизы прикрепленных изображений
|
|
|
|
|
 |
Ответов
|
Oct 3 2010, 19:00
|
Участник

Группа: Участник
Сообщений: 72
Регистрация: 7-06-08
Пользователь №: 38 128

|
http://ru.wikipedia.org/wiki/%D0%A2%D0%B5%...%BE%D0%B2%D0%B0 : Цитата …если аналоговый сигнал x(t)\; имеет ограниченный спектр, то он может быть восстановлен однозначно и без потерь по своим дискретным отсчётам, взятым с частотой строго большей удвоенной максимальной частоты спектра Всё строго и чётко. Восстановить можно при частоте строго большей. Синус, дискретизированной с частотой равной его собственной частоте является просто рядом нулей. Естессно, если не добавлять сдвиг фазы. Чем ближе подбирается полоса сигнала к частоте дискретизации, тем больше времени надо на переходные процессы при включении/выключении. Хорошее и простое объяснение: чем ближе полосу сигнала к частоте дискретизации вы хотите, тем более крутую переходную область у фильтра вам надо, тем длиннее его импульсная характеристика.
|
|
|
|
|
Oct 4 2010, 10:14
|
Местный
  
Группа: Участник
Сообщений: 312
Регистрация: 21-10-05
Пользователь №: 9 952

|
Хочу возразить всем, кто считает что в теореме Котельникова должно быть строгое неравенство. 1. Немного странно, но насколько мне известно, ни Котельников ни Шеннон свои первоначальные формулировки теоремы с равенством так и не скорректировали. 2. Теорема Котельникова для сигналов со спектром, содержащим ненулевые комплексные значения на F=1/2*Fs, работает. Возьмите отсчеты функции Котельникова с произвольным некратным 1/Fs смещением по времени (что дает отличный от нуля комплексный спректр в точках +-1/2*Fs) и подставьте в ряд Котельникова - после недолгих преобразований получите равномерно сходящийся ряд (сумма которого приведена в "Интегралы и ряды" Прудникова), суммой ряда будет исходная смещенная по времени функция Котельникова. 3. Очевидно, что дискретизация синуса частоты Fs/2 ведет к потере информации. Также очевидно, что такая потеря возникает из-за наложения спектров при дискретизации. К такому же эффекту приводит и дискретизация косинуса частоты Fs/2 - энергия отсчетов косинуса частоты Fs/2 ровно вдвое больше энергии отсчетов косинуса частоты, например, Fs/4. А если так, то ни синус частоты Fs/2 ни косинус частоты Fs/2 под условия теоремы не попадают. Почему - предлагаю подумать самостоятельно.
Сообщение отредактировал 729 - Oct 4 2010, 10:15
|
|
|
|
|
Oct 4 2010, 11:19
|
Местный
  
Группа: Участник
Сообщений: 312
Регистрация: 21-10-05
Пользователь №: 9 952

|
Цитата(тау @ Oct 4 2010, 14:51)  Комплексные отсчеты сигнала способны восстановить гармонику Fs/2, хотя бы по той простой причине что в спектре эта гармоника будет стоять только на одном краю (в отличие от вещественного сигнала) . А стало быть спектр получается чисто периодический и без алиасинга на краях. Математика просто обязана сработать, потому что главное в ней ( в теореме) - строгая периодичность спектров. С комплексными отсчетами она соблюдается. Я говорил про комплексность спектра дискретизируемой функции, но не про комплексность дискретизируемой функции. Про периодичность спектра - это отдельная тема. Но даже в рамках модели, которая периодизирует спектр при дискретизации, в точках +-Fs/2 ничего страшного не будет - значения ДВПФ от отсчетов косинуса в точках +-Fs/2 будут стремиться к бесконечности в 2 раза быстрее, чем ДВПФ от отсчетов косинуса в точках +-Fs/4 или +-Fs/3. Но ряд Котельникова эту двойку съест и чистый неудвоенный по амплитуде косинус выдаст.
|
|
|
|
|
Oct 4 2010, 11:36
|
Местный
  
Группа: Участник
Сообщений: 312
Регистрация: 21-10-05
Пользователь №: 9 952

|
Цитата(Oldring @ Oct 4 2010, 15:26)  Мне кажется, вы заговариваетесь. Начнем с того, что если говорить про "спектр" как про энергетический спектр, то он действителен неизбежно. А если говорить про преобразование Фурье, то оно действительно только для четных действительных функций. Можете сформулировать свои утверждения четко? Под спектром всегда понимал спектральную функцию, но не модули и прочие производные от неё. Утверждение простое. Спектр sin(pi*Fs*t)/(pi*Fs*t) - функции Котельникова - действителен. Спектр смещенной на тау функции Котельникова - sin(pi*Fs*(t-тау))/(pi*Fs*(t-тау)) уже комплексный в точках +-Fs/2, если тау не равно k*1/Fs. Вы это утверждение имели в виду? Цитата(Oldring @ Oct 4 2010, 15:29)  Так а в чём же тогда конструктивный вклад ваших постов в обсуждаемую тему? В том, что: 1 Равенство в теореме есть и оно работает. 2 Что синус и косинус Fs/2 условиям теоремы не удовлетворяют. Цитата(Oldring @ Oct 4 2010, 15:29)  Про Fs/2 речь не шла. Ну вот... Извиняюсь. Я встрял в Ваш с fontp диалог.
|
|
|
|
|
Oct 4 2010, 11:47
|

Гуру
     
Группа: Свой
Сообщений: 3 041
Регистрация: 10-01-05
Из: Москва
Пользователь №: 1 874

|
Цитата(729 @ Oct 4 2010, 15:36)  Под спектром всегда понимал спектральную функцию Эту? http://www.answers.com/topic/spectral-functionЯ, например, привык под спектром понимать спектр стационарного случайного процесса. Называть "спектром" просто результат преобразование Фурье, наверное, не вполне корректно. Цитата(729 @ Oct 4 2010, 15:36)  В том, что: 1 Равенство в теореме есть и оно работает. 2 Что синус и косинус Fs/2 условиям теоремы не удовлетворяют. Так речь в теме шла про строгое неравенство в части границ частотного диапазона. Ваше второе утверждение и означает, что это неравенство строгое. А первое утверждение про какое-то другое равенство-неравенство, наверное. Ссылку на текст работы Котельникова я привел выше. Давайте отталкиваться от него.
--------------------
Пишите в личку.
|
|
|
|
|
Oct 4 2010, 12:14
|
Местный
  
Группа: Участник
Сообщений: 312
Регистрация: 21-10-05
Пользователь №: 9 952

|
Цитата(Oldring @ Oct 4 2010, 15:47)  Я, например, привык под спектром понимать спектр стационарного случайного процесса. Называть "спектром" просто результат преобразование Фурье, наверное, не вполне корректно. Давайте для определённости считать спектром просто непрерывное преобразование Фурье. Цитата(Oldring @ Oct 4 2010, 15:47)  Так речь в теме шла про строгое неравенство в части границ частотного диапазона. Ваше второе утверждение и означает, что это неравенство строгое. А первое утверждение про какое-то другое равенство-неравенство, наверное. Ссылку на текст работы Котельникова я привел выше. Давайте отталкиваться от него. Моё утверждение просто до безобразия. 1. Если спектр функции отличен от нуля вне отрезка [-Fh,Fh] (точки +-Fh в этот отрезок входят), то частота дискретизации Fs>=1/(2Fh). 2. Функции синус и косинус Fs/2 дискретизируются со спектральными наложениями, следовательно спектры этих функций за отрезок [-Fh,Fh] "вылазят".
Сообщение отредактировал 729 - Oct 4 2010, 12:16
|
|
|
|
Сообщений в этой теме
Skaf Почему не работает теоремма Котельникова при F = 1/2Fs Oct 2 2010, 12:08 MrYuran Ну и где удвоенная частота? Oct 2 2010, 12:13 Andron_ потому что функцию неправильно задаете.
y = sin (... Oct 2 2010, 12:15 Skaf Странно, нам на лекциях давали именно так... Oct 2 2010, 12:18 Taradov Alexander QUOTE (Skaf @ Oct 2 2010, 16:18) Странно,... Oct 2 2010, 13:26 V_G Вообще-то в теореме Котельникова неравенство строг... Oct 2 2010, 14:04 Microwatt Цитата(V_G @ Oct 2 2010, 17:04) Вообще-то... Oct 2 2010, 16:04 Skaf ЦитатаА нули получаются из-за того, что вы всегда ... Oct 2 2010, 15:50 Aner Теорема Котельникова работает, точка! Как то д... Oct 2 2010, 16:22 Dr.Alex Цитата(Aner @ Oct 2 2010, 20:22) Теорема ... Oct 2 2010, 17:37 rezident Шо, опять???!!! Была обширная тема в... Oct 2 2010, 19:13 тау Цитата(rezident @ Oct 2 2010, 23:13) То б... Oct 2 2010, 20:11  rezident Цитата(тау @ Oct 3 2010, 02:11) для абсол... Oct 2 2010, 20:33   тау Цитата(rezident @ Oct 3 2010, 00:33) Угу.... Oct 2 2010, 22:37    Oldring Теорема работает, но автор, очевидно, использует е... Oct 3 2010, 11:43     Lisitsin Цитата(Oldring @ Oct 3 2010, 15:43) Теоре... Oct 3 2010, 15:29      Tanya Цитата(Lisitsin @ Oct 3 2010, 19:29) Можн... Oct 3 2010, 16:12   Tanya Цитата(729 @ Oct 4 2010, 14:14) 1. Немног... Oct 4 2010, 10:31   fontp QUOTE (729 @ Oct 4 2010, 14:14) 3. Очевид... Oct 4 2010, 10:36    729 Цитата(fontp @ Oct 4 2010, 14:36) 3. Авто... Oct 4 2010, 11:01     fontp QUOTE (729 @ Oct 4 2010, 15:01) Согласите... Oct 4 2010, 11:06      Oldring Цитата(fontp @ Oct 4 2010, 15:06) Это зав... Oct 4 2010, 11:21       fontp QUOTE (Oldring @ Oct 4 2010, 15:21) Сигна... Oct 4 2010, 11:21      729 Цитата(fontp @ Oct 4 2010, 15:06) Не согл... Oct 4 2010, 11:26    Oldring Цитата(fontp @ Oct 4 2010, 14:36) Энергия... Oct 4 2010, 11:06        fontp QUOTE (Oldring @ Oct 4 2010, 15:47) Эту? ... Oct 4 2010, 11:59         Oldring Цитата(fontp @ Oct 4 2010, 15:59) Нет, вп... Oct 4 2010, 12:16          729 Цитата(Oldring @ Oct 4 2010, 16:16) Прежд... Oct 4 2010, 12:20           Oldring Цитата(729 @ Oct 4 2010, 16:20) В чем же?... Oct 4 2010, 12:25            729 Цитата(Oldring @ Oct 4 2010, 16:25) В том... Oct 4 2010, 13:17             fontp QUOTE (729 @ Oct 4 2010, 17:17) 2. Спектр... Oct 4 2010, 13:28              729 Цитата(fontp @ Oct 4 2010, 17:28) К налож... Oct 4 2010, 13:37               fontp QUOTE (729 @ Oct 4 2010, 17:37) Вот и пол... Oct 4 2010, 13:40              Oldring Цитата(fontp @ Oct 4 2010, 17:28) К налож... Oct 4 2010, 13:45               fontp QUOTE (Oldring @ Oct 4 2010, 17:45) Это к... Oct 4 2010, 13:47                Oldring Цитата(fontp @ Oct 4 2010, 17:47) Хвосты ... Oct 4 2010, 13:51                 fontp QUOTE (Oldring @ Oct 4 2010, 17:51) Если ... Oct 4 2010, 13:52                  Oldring Цитата(fontp @ Oct 4 2010, 17:52) Это дру... Oct 4 2010, 13:53             Oldring Цитата(729 @ Oct 4 2010, 17:17) 2. Спектр... Oct 4 2010, 13:38         fontp QUOTE (729 @ Oct 4 2010, 16:14) Давайте д... Oct 4 2010, 12:56          729 Цитата(fontp @ Oct 4 2010, 16:56) 1. Если... Oct 4 2010, 13:30           fontp QUOTE (729 @ Oct 4 2010, 17:30) Действите... Oct 4 2010, 13:33         GetSmart Цитата(729 @ Oct 4 2010, 17:14) 2. Функци... Oct 4 2010, 16:51          bahurin я для себя случай представления гармоники сигнала ... Oct 4 2010, 18:42           GetSmart Цитата(bahurin @ Oct 4 2010, 23:42) Это Н... Oct 4 2010, 19:03            fontp QUOTE (GetSmart @ Oct 4 2010, 23:03) Спор... Oct 5 2010, 03:40             GetSmart Цитата(fontp @ Oct 5 2010, 08:40) Нет ник... Oct 5 2010, 06:13   Oldring Цитата(729 @ Oct 4 2010, 14:14) 2. Теорем... Oct 4 2010, 10:55    729 Цитата(Oldring @ Oct 4 2010, 14:55) По по... Oct 4 2010, 11:24     Oldring Цитата(729 @ Oct 4 2010, 15:24) Вы невним... Oct 4 2010, 11:29 Tanya Цитата(Designer56 @ Oct 3 2010, 20:20) Не... Oct 4 2010, 03:23  DRUID3 Цитата(Tanya @ Oct 4 2010, 06:23) В этом ... Oct 4 2010, 04:03  Designer56 Цитата(Tanya @ Oct 4 2010, 09:23) В этом ... Oct 4 2010, 08:24   DRUID3 Цитата(Designer56 @ Oct 4 2010, 11:24) Не... Oct 4 2010, 08:30 thermit Цитата729:
Спектры функций (в смысле ПФ обобщенных... Oct 4 2010, 13:33 729 Цитата(thermit @ Oct 4 2010, 17:33) С это... Oct 4 2010, 14:19  Oldring Цитата(729 @ Oct 4 2010, 17:55) посчитаем... Oct 4 2010, 14:27   729 Цитата(Oldring @ Oct 4 2010, 18:27) ... Oct 4 2010, 14:41    Oldring Цитата(729 @ Oct 4 2010, 18:41) Это там, ... Oct 4 2010, 14:45     729 Цитата(Oldring @ Oct 4 2010, 18:45) У мен... Oct 4 2010, 14:49      Oldring Цитата(729 @ Oct 4 2010, 18:49) Да и не в... Oct 4 2010, 14:58       729 Цитата(Oldring @ Oct 4 2010, 18:58) Всё р... Oct 4 2010, 15:05        Oldring Цитата(729 @ Oct 4 2010, 19:05) Именно та... Oct 4 2010, 15:26         729 Цитата(Oldring @ Oct 4 2010, 19:26) То ес... Oct 4 2010, 20:14          Oldring Цитата(729 @ Oct 5 2010, 00:14) Это как в... Oct 4 2010, 21:44           729 Цитата(Oldring @ Oct 5 2010, 01:44) Давай... Oct 5 2010, 10:26            Oldring Цитата(729 @ Oct 5 2010, 14:04) Дельта им... Oct 5 2010, 10:35 thermit Блин. Ну дык кто-нить может ответить на 1 частный ... Oct 5 2010, 04:41 Tanya Цитата(thermit @ Oct 5 2010, 08:41) Блин.... Oct 5 2010, 05:26  Oldring Цитата(Tanya @ Oct 5 2010, 09:26) Не меша... Oct 5 2010, 09:12   тау Цитата(Oldring @ Oct 5 2010, 13:12) Разже... Oct 5 2010, 09:36    Oldring Цитата(тау @ Oct 5 2010, 13:36) Oldring, ... Oct 5 2010, 09:58     729 Цитата(Oldring @ Oct 5 2010, 13:58) Поэто... Oct 5 2010, 10:34      Oldring Цитата(729 @ Oct 5 2010, 14:34) Интересно... Oct 5 2010, 10:44       729 Цитата(Oldring @ Oct 5 2010, 14:44) Отвра... Oct 5 2010, 10:50        Oldring Цитата(729 @ Oct 5 2010, 14:50) Имелись в... Oct 5 2010, 11:07 thermit ЦитатаЦитата(thermit @ Oct 5 2010, 09:41) *
Почему... Oct 5 2010, 06:50 thermit ЦитатаOldring:
Разжевывать, что синусоиду нельзя в... Oct 5 2010, 09:32 Oldring Цитата(thermit @ Oct 5 2010, 13:32) Гы. А... Oct 5 2010, 09:39 thermit Цитататау:
Неужели в допущении при переходе от инт... Oct 5 2010, 09:57 thermit Цитата729:
Не удовлетворяет этот cos условиям теор... Oct 5 2010, 10:32 729 Цитата(thermit @ Oct 5 2010, 14:32) Каком... Oct 5 2010, 10:42  Oldring Цитата(729 @ Oct 5 2010, 14:42) То есть, ... Oct 5 2010, 10:56 thermit Цитата729:
Условию равенства нулю спектра за преде... Oct 5 2010, 10:43 729 Цитата(thermit @ Oct 5 2010, 14:43) Т е д... Oct 5 2010, 10:56 thermit ЦитатаOldring:
Потому что вы зафиксировали фазу, в... Oct 5 2010, 11:44 Oldring Цитата(thermit @ Oct 5 2010, 15:44) Ну и ... Oct 5 2010, 11:49  GetSmart Цитата(Oldring @ Oct 5 2010, 16:49) Добав... Oct 5 2010, 14:30  bahurin Цитата(Oldring @ Oct 5 2010, 15:49) В слу... Oct 5 2010, 18:17   rezident Цитата(bahurin @ Oct 6 2010, 00:17) Ну эт... Oct 5 2010, 18:20    bahurin Цитата(rezident @ Oct 5 2010, 22:20) Пурк... Oct 5 2010, 18:28
2 страниц
1 2 >
2 чел. читают эту тему (гостей: 2, скрытых пользователей: 0)
Пользователей: 0
|
|
|