|
Почему не работает теоремма Котельникова при 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, 10:36
|

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

|
QUOTE (729 @ Oct 4 2010, 14:14)  3. Очевидно, что дискретизация синуса частоты Fs/2 ведет к потере информации. Также очевидно, что такая потеря возникает из-за наложения спектров при дискретизации. К такому же эффекту приводит и дискретизация косинуса частоты Fs/2 - энергия отсчетов косинуса частоты Fs/2 ровно вдвое больше энергии отсчетов косинуса частоты, например, Fs/4. А если так, то ни синус частоты Fs/2 ни косинус частоты Fs/2 под условия теоремы не попадают. Почему - предлагаю подумать самостоятельно. 3. Авторы возражали в практическом смысле - что восстановить сигнал в действительности не получится, поскольку потребуется и бесконечное время и идеальные прямоугольные фильтры. И они не поспеют за real-time. Вы же ссылаетесь на математическую формулировку, что какой-то там бесконечный ряд не сходится (или сходится) равномерно А почему если дискретизация косинуса частоты Fs/2 - энергия отсчетов косинуса частоты Fs/2 ровно вдвое больше энергии отсчетов косинуса частоты, например, Fs/4 - то они под условие теоремы не попадают?  Или это просто так сказано, для красоты и ничего из этого не следует? Энергия косинуса частоты, например, Fs/4 или ещё лучше Fs/3 тоже зависит от фазы дискретизации и с физической энергией не совпадает. Да и в условиях теоремы о них не сказано ничего, как не напрягайся.
|
|
|
|
|
Oct 4 2010, 11:01
|
Местный
  
Группа: Участник
Сообщений: 312
Регистрация: 21-10-05
Пользователь №: 9 952

|
Цитата(fontp @ Oct 4 2010, 14:36)  3. Авторы возражали в практическом смысле - что восстановить сигнал в дествительности не получится, поскольку потребуется и бесконечное время и идеальные прямоугольные фильтры. Вы же ссылаетесь на математическую формулировку, что какой-то там бесконечный ряд не сходится равномерно Я ссылаюсь на математическую формулировку, ибо речь шла о ней (формулировке теоремы). Практический смысл только упоминался иногда. Согласитесь, что практического смысла и сам ряд Котельникова не имеет, ибо в общем случае он бесконечный. Цитата(fontp @ Oct 4 2010, 14:36)  Энергия косинуса частоты, например, Fs/4 или ещё лучше Fs/3 тоже зависит от фазы дискретизации и с физической энергией не совпадает. Имелась в виду энергия отсчетов косинуса частоты Fs/4 и косинуса частоты Fs/2 на интервале, где укладывается целое число периодов обеих частот. А привел только по тому, что удвоение энергии есть результат спектральных наложений при дискретизации косинуса частоты Fs/2. При дискретизации косинуса частоты Fs/4 или Fs/3 спектральных наложений нет. Физическая энергия здесь вообще не причем.
|
|
|
|
|
Oct 4 2010, 11:06
|

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

|
QUOTE (729 @ Oct 4 2010, 15:01)  Согласитесь, что практического смысла и сам ряд Котельникова не имеет, ибо в общем случае он бесконечный. Не соглашусь. Ряд Котельникова имеет практический физический смысл повсеместно. Там, где с нужной точностью можно приблизить сигнал конечной суммой членов ряда. Бесконечности вообще имеют практический смысл только в смысле предельного перехода, да и то хорошо бы всегда знать как любое епсилон зависит от того N,что когда-то найдётся)) Иначе бесконечности полностью бесполезны QUOTE (Oldring @ Oct 4 2010, 15:06)  Энергию для чистой гармоники сравнивать вообще нельзя, так как она бесконечна. А вот мощность неизбежно совпадает, так как базис Котельникова ортонормирован. Это зависит смотря на каком интервале, как объяснил автор. Практически зависит во всех трёх ваших предложениях)) Не отвлекайте
|
|
|
|
|
Oct 4 2010, 11:26
|
Местный
  
Группа: Участник
Сообщений: 312
Регистрация: 21-10-05
Пользователь №: 9 952

|
Цитата(fontp @ Oct 4 2010, 15:06)  Не соглашусь. Ряд Котельникова имеет практический физический смысл повсеместно. Там, где с нужной точностью можно приблизить сигнал конечной суммой членов ряда. Бесконечности вообще имеют практический смысл только в смысле предельного перехода, да и то хорошо бы всегда знать как любое епсилон зависит от того N,что когда-то найдётся)) Иначе бесконечности полностью бесполезны В такой трактовке и я с Вами соглашусь - имеет, конечно. Но сами бесконечности интересны только для теоретических выкладок.
|
|
|
|
Сообщений в этой теме
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      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    Oldring Цитата(fontp @ Oct 4 2010, 14:36) Энергия... Oct 4 2010, 11:06   тау Цитата(729 @ Oct 4 2010, 14:14) Хочу возр... Oct 4 2010, 10:51    729 Цитата(тау @ Oct 4 2010, 14:51) Комплексн... Oct 4 2010, 11:19     Oldring Цитата(729 @ Oct 4 2010, 15:19) Я говорил... Oct 4 2010, 11:26      729 Цитата(Oldring @ Oct 4 2010, 15:26) Мне к... Oct 4 2010, 11:36       Oldring Цитата(729 @ Oct 4 2010, 15:36) Под спект... Oct 4 2010, 11:47        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        729 Цитата(Oldring @ Oct 4 2010, 15:47) Я, на... Oct 4 2010, 12:14         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
|
|
|