70829

Децимация и интерполяция

Лабораторная работа

Коммуникация, связь, радиоэлектроника и цифровые приборы

Выполнение процедуры децимации (уменьшения частоты дискретизации в заданное целое число раз) приводит к уменьшению частоты дискретизации исходной последовательности. В процессе децимации исходная последовательность обрабатывается НЧ фильтром, после чего производится выборка с необходимой частотой.

Русский

2014-10-27

104 KB

19 чел.

МИНИСТЕРСТВО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

Марийский государственный технический университет

Факультет информатики и вычислительной техники

Кафедра ИВС

Отчет по лабораторной работе №4

по дисциплине "ТДЛС"

"Децимация и интерполяция"

Выполнили: ст. гр. ВМ-41

Дудоладов В.А.

Алдушкин А.А.

Венедиктов Р.В.

Проверила: Малашкевич И.А.

Йошкар-Ола

2003


Текст программы, которая строит график исходного сигнала, находит его спектр, выполняет децимацию и интерполяцию:

t=0:0.0002:1

figure

x=sin(2*pi*30*t)+sin(2*pi*60*t)

stem(x(1:120))

figure

y=decimate(x,4)

stem(y(1:30))

figure

z=interp(y,4)

stem(z(1:120))

figure

fft_out=mot563_fft(x,'56301')

mag=abs(fft_out)

freq=(0:127)*1000/256

plot(freq,mag(1:128))

Частота дискретизации сигнала = 5кГц

Исходный сигнал

x=sin(2*pi*30*t)+sin(2*pi*60*t)

Спектр сигнала

Так как исходный сигнал сложный, то есть является суммой 2-х синусоид, то в спектре сигнала хорошо просматриваются две частотные составляющие:

Понятие децимации сигнала

Выполнение процедуры децимации (уменьшения частоты дискретизации в заданное целое число раз) приводит к уменьшению частоты дискретизации исходной последовательности. В процессе децимации исходная последовательность обрабатывается НЧ фильтром, после чего производится выборка с необходимой частотой.

Функция у = decimate(x, r) уменьшает в г раз частоту дискретизации выборки, находящейся в векторе х. Результирующая выборка возвращается в векторе y, длина которого в г раз меньше длины вектора х. По умолчанию используется НЧ фильтр Чебышева типа 18-го порядка. Для исключения фазовых искажений фильтрация выполняется в прямом и обратном направлениях.

Понятие интерполяции сигнала

Выполнение процедуры интерполяции (увеличения частоты дискретизации в заданное целое число раз) приводит к увеличению частоты дискретизации исходной последовательности. Функция interp производит интерполяцию за два шага. Сначала в исходную последовательность вносятся нулевые элементы, а затем полученная последовательность обрабатывается специальным НЧ фильтром.

Функция z = interp(y,r) увеличивает в г раз частоту дискретизации выборки, находящейся в векторе y. Результирующая выборка возвращается в векторе z, длина которого в г раз больше длины вектора y.

 Практическая реализация децимации и интерполяции

Для изучения децимации и интерполяции в каждом случае использовался один и тот же исходный сигнал. На каждом шаге для этого сигнала производилось уменьшение частоты дискретизации в n раз (децимация), а затем увеличение частоты децимированного сигнала также в n раз (интерполяция). После этого производилось сравнение графиков сигналов, полученных при помощи программной среды MatLab: конечного сигнала (сигнала после интерполяции) и исходного сигнала.

Исходный сигнал:

  (Fд=5000 Гц)

 Уменьшение частоты в 4 раза (Fд=1250 Гц):         Увеличение частоты в 4 раза (Fд=5000 Гц):

   Уменьшение частоты в 8 раз (Fд=625 Гц):          Увеличение частоты в 8 раз (Fд=5000 Гц):

Уменьшение частоты в 16 раз (Fд=312,5 Гц):        Увеличение частоты в 16 раз (Fд=5000 Гц):

Исходный сигнал:

   (Fд=5000 Гц)

Уменьшение частоты в 24 раза (Fд≈208,3 Гц):      Увеличение частоты в 24 раза (Fд=5000 Гц):

  Уменьшение частоты в 30 раз (Fд≈166,7 Гц):        Увеличение частоты в 30 раз (Fд=5000 Гц):

Выводы

 Согласно теореме Котельникова сигнал подлежит корректному восстановлению тогда, когда частота дискретизации больше максимальной частоты исходного сигнала в 2 и более раз. Так как исходный сигнал имеет максимальную частоту  ≈ 380 Гц (это видно из его спектра), то минимальная частота дискретизации, при которой сигнал еще может быть восстановлен корректно, должна быть больше 380*2=760 Гц. Имея частоту исходного сигнала 5000 Гц, определим во сколько раз можно максимально уменьшить частоту дискретизации сигнала: 5000/760≈6,6≈6 раз (округление должно выполняться в меньшую сторону до целого числа).


 

А также другие работы, которые могут Вас заинтересовать

22599. Контроль клавіатурного вводу 32 KB
  Скенкод це однобайтне число молодші 7 бітів якого представляють ідентифікаційний номер призначений кожній клавіші. На всіх машинах крім AT старший біт коду говорить про те чи була клавіша натиснута біт = 1 код натискання або відпущена біт = 0 код звільнення. Наприклад 7бітный скенкод клавіші В 48 або 110000 в двійковій системі.
22600. Управління відеоадаптером IBM PC 35.5 KB
  Однак вона також встановлює режим екрана управляє курсором і для кольорового графічного адаптора управляє кольором. Розмір і розташування цих буферів міняється з системою режимом екрана а також кількістю заздалегідь відведеної пам'яті. Коли в буфері зберігається декілька образів екрана то кожний окремий образ називають екранною сторінкою. Цього досить для відображення одного графічного екрана без сторінок або від чотирьох до восьми екранів тексту в залежності від числа символів в рядку 40 або 80.
22602. SQL. Знайти всі відомості про всіх постачальників 56 KB
  SELECT FROM П; Знайти номери постачальників з статусом більш 20 що живуть у місті N. SELECT КП FROM П WHERE місто = ‘N’AND статус 20; Знайти прізвища постачальників які постачають деталь Д1. SELECT Прізвище FROM П WHERE КП IN SELECT КП FROM ОПД WHERE КД = ‘Д1’; Знайти прізвища постачальників які постачають принаймні одну червону деталь. SELECT Прізвище FROM П WHERE КП IN SELECT КП FROM ОПД WHERE КД IN SELECT КД FROM ОПД WHERE колір = ‘червоний’; Дужки обов’язково ставляться якщо є неоднозначність.
22603. Political System of Federal Republic Germany 2.03 MB
  0 – December 6 2004 CONTENTS Introduction I. December 2004 06:20:44 Bundesland. December 2004 06:20:45 Candidate. December 2004 06:20:45 City.
22605. Психологические особенности внедрения систем электронного документооборота 267 KB
  Рассмотреть внедрение системы электронного документооборота как проект, рассмотреть технические особенности и сложности внедрения системы электронного документооборота, рассмотреть психологические особенности внедрения системы электронного документооборота, рассмотреть существующие психологические типы сотрудников
22606. Реляційне числення. Мова “Альфа” 58.5 KB
  RANGE ОПД X GET WП.КД = Д3 RANGE ОПД X – оператор декларації ОПД – тип X – змінна. Перший варіант: RANGE Д X GET WОПД.колір = червоний RANGE ОПД Y GET W2П.