1. Проектирование фнч методом весовых окон: (Требуемое ослабление 65дБ)

Рисунок 3.1 – Импульсная, амплитудно-частотная и фазовая характеристики ФНЧ (методом весовых окон)

fd=16000; % частота дискретизации

fc=4800; % частота среза Гц

fp=3200; % полоса пропускания

wcp=2*pi*fc/fd;% норм. Частоты среза

wper=2*fp*pi/fd; %норм. переходная полоса

M=round((65-7.95)/(2.285*wper)); %порядок фильтра для окна Хэмминга

Wn=2*fc/fd;

Wpod=2*(fc+fp)/fd;

h=fir1(M,Wn,'low',window);

freqz(h);

2. Проектирование фнч методом частотной выборки:

Рисунок 3.2 – Импульсная, амплитудно-частотная и фазовая характеристики ФНЧ (методом частотной выборки)

M=23

f=0:fd/(2*M):fd/2;

f1=2*f./fd;

b=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0];

h1=fir2(M,f1,b);

freqz(h1);

Добились требуемого ослабления на частоте fpod=8000Гц.

Чтобы добиться ослабления на частоте подавления (8000 Гц) при использовании метода частотной выборки нам понадобился более высокий порядок фильтра. (М=23)

3. Проверка фильтра и подача на его вход смеси сигналов

Подадим на вход фильтра, например: sin(2*pi*f1*t)+cos(2*pi*f2*t)+0.1*randn(size(t)), где f1=500 Гц, f2=8000 Гц.

Используем наш ФНЧ, чтобы подавить частоту 8000 Гц.

Используя команду sptool и и пропускаем смесь через наш ФНЧ, спроектированный методом весовых окон:

Рисунок 3.3 – Частотный спектр отфильтрованного сигнала (ФНЧ с весовым окном)

Добились подавления нижней частоты, следовательно, спроектированный нами ФНЧ методом весового окна Хэмминга правильный.

Проверим фильтр, спроектированный методом частотной выборки.

Рисунок 3.4 – Частотный спектр отфильтрованного сигнала (ФНЧ методом частотной выборки)

Фильтр также удовлетворяет нашим требованиям.

Вывод: как видно из работы , чтобы добиться ослабления на частоте подавления (500 Гц) при использовании метода частотной выборки нам понадобился более высокий порядок фильтра. (М=23) в отличии от метода весовых окон, который в этом случае является более практичным.

Часть 4.

_{1) Провести операцию интерполяции смеси гармонических сигналов и шума (увеличить частоту дискретизации)}

_{2) Спроектировать ФНЧ, для последующего его использования при децимации}

_{3) Провести операцию децимации (уменьшение частоты дискретизации)}

Децима́ция (от лат. decimatio, от decem — «десять») — уменьшение частоты дискретизации дискретного во времени сигнала путем удаления его отсчетов.

При децимации из исходной последовательности отсчетов

a0, a1, a2, …

берется каждый N-й отсчет (N — целое число):

a0, aN, a2N, … ; N > 1

остальные отсчеты отбрасываются. Преобразование спектра при децимации существенно зависит от спектра исходного сигнала:

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

Если исходный сигнал содержит частоты, превышающие частоту Найквиста децимированного сигнала, то при децимации будет иметь место алиасинг (наложение спектров).

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

Интерполяция это процесс обратный децимации (увеличение частоты дискретизации)

Дано:

№ по списку в журнале 17.

Fd=8400 Гц, f1=2000 Гц, f2=3800, l=4, m=5

F₁, F₂– несущие частоты сигнла,

l – коэффициентинтероляции,

m– коэффициент децимации,

F_d – частота дискритизации.

Рисунок 4.1 – Структурная схема проектируемого фильтра.

4.1. Зададим сигнал в виде смеси гаромнических и шума:

clear all

Fs=8400

f1=2000

f2=3800

t1=1

td=1/Fs

t=0:td:t1

x=sin(2*pi*f1*t)+cos(2*pi*f2*t)+0.1*randn(size(t))

Рисунок 4.2 - Временная диаграмма нашей смеси

Рисунок 4.3 - Частотный спектр смеси, при помощи программы sptool

2. Проводим интерполяцию: увеличиваем частоту дискретизации в 3 раза командой:

w=upfirdn(x,1,4,1)

Рисунок 4.4 - Частотный спектр, после интерполяции

2. Проектируем фильтр и подготавливаем сигнал для децимации

Выбираем частоту среза исходя из соотношения: fd*3/16=>2662 Гц < fc. Выберем fc =2100 Гц

На основе этого проектируем фильтр с заданной частотой дискретизации и частотой среза при помощи команды fdatool

Рисунок 4.5 - Амплитудно-частотная характеристика, требуемого ФНЧ

Рисунок 4.6 - Импульсная характеристика нашего ФНЧ

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

Уменьшение частоты дискретизации в 8 раз), при помощи команды: w1=upfirdn(w,Num,1,5), где Num – коэффициенты испульсной характеристики нашего ФНЧ.

Рисунок 4.7 - Временная диаграмма сигнала и его спектр, после децимации и фильтрации

На спектре видно, что теперь наша частота f1=2000 Гц, передается уже на новой частоте дискретизации.

Вывод.

Сравнивая спектры, полученные на последних рисунках, можно сделать вывод о некорректной работе функции upfirdn(x,h,l,m). Как мы можем видеть одна из зеркальных составляющих спектра последнего рисунка не была подавлена. Мы же смогли устранить эффект Гиббса с помощью промежуточного фильтра. В итоге, хоть и часть информации была утеряна, но искажения проявились в меньшей степени (см. временные диаграммы).

Заключение

Все понятия и методы, использованные в данной курсовой работе, образуют взаимосвязанное единство по анализу и синтезу любых ЦУ.

В данной работе был произведён ряд вычислений, наиболее полно отражающий содержание курса. В частности научились дискретизировать различные сигналы и получали их спектры (ДПФ). Так же научились по спектру дискретного сигнала получать сигнал во временной области (ОДПФ). Научились проектировать цифровые фильтры в среде МatLAB с использованием пакетов sptool и fdatool. Выполнили преобразование частоты дискретизации.

В данной работе я выполнил анализ большого объема информации по ЦОС, радиоэлектронике. В процессе работы, были использованы следующие математические пакеты: Mathcad 14, MatLab2006, AdvancedGrapher.Все они мне помогли сделать работу точнее, быстрее и нагляднее.

Все расчеты имеют практический характер в реально жизни. Примером может стать цифровое телевидение, которое в момент написания курсовой работы бурно вводится в России (DVB-T). Все эти преобразования используются там, но конечно сигналы и спектры там намного сложнее.

Список использованной литературы

1. Конспект лекций по МоЦОС: СибГУТИ.

2. ЦОС. Часть 1. Методические указания к лабораторным работам/Рязань. Гос РТ университет./Сост: В.В Витязев,2009, 36с.

3. Радиотехника: Энциклопедия / Под ред. Ю. Л. Мазова, Е. А. Мачусского, В.И. Правды.- 2-е изд., стер.- М.: Издательский дом «Додэка-XXI», 2009.-944 с.

4. Зеленцов Б. П. Математика в формулах и таблицах – Новосибирск: СибГУТИ, 2005. – 117 с.

5. Макаров Е.Г. Mathcad 14: учебный курс.-Спб.: Питер, 2009.-384 с.: ил.