Задание 2. Фильтрация сигналов для демодуляции амплитудно-манипулированных сигналов в гауссовских каналах связи
КАЗАНСКИЙ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ
ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ им. А.Н. Туполева - КАИ
Институт радиоэлектроники и телекоммуникаций
Практическая работа по дисциплине
«Статистические методы в теории связи»
Вариант ___
Выполнили:
Студенты группы 5179
Проверил:
Д.т.н., проф., Р.Р. Файзуллин
Казань 2018
Задание 1. Демодуляция АМ-сигнала
Порядок выполнения работы:
1) Зададим временной массив t: 0≤t≤1;
2) Зададим линейные частоты модулирующего 
и несущего 
колебаний.
3) Зададим амплитуду несущегосигнала 
. Для модулирующего сигнала зададим амплитуду 
, постоянную составляющую 
, и начальную фазу колебаний 
.
4) Зададим коэффициент модуляции m: 0<m≤1 (при m>1 наступает перемодуляция(избыточная модуляция)), простые демодуляторы (типа квадратичного детектора) демодулируют такой сигнал с сильными искажениями.
5) Сгенерируем несущий сигнал и модулирующий сигнал, используя входные параметры:
;
;
6) Сгенерируем амплитудно-модулированный сигнал 
:
;
Т.к. , то:
;
7) Сгенерируем массив случайных величин (соразмерный смассивом t),распределенных по равномерному закону в интервале (0, 1)[белый шум 
]. Зададим его амплитуду 
. В канале связи будет создаваться аддитивная помеха, т.е. сигнал на входе приемника будет иметь вид:
;
8) Проведем фильтрацию принимаемого сигнала .
Операция линейной дискретной фильтрации в общем случае описывается следующим образом:
|
|
|
;
— отсчеты входного сигнала;
— отсчеты выходного сигнала;
—постоянные коэффициенты;
Максимальное из чисел 
и 
называется порядком фильтра.
Рис 1. Цифровой рекурсивный фильтр
Передаточная функция фильтра имеет вид:
;
В среде MatLab используется функция 
.
9) Зададим параметры фильтра:
;
;
И проведем операцию демодуляции принимаемого сигнала.
10) Вычислим отношение сигнал/шум по формуле:
;
Листинг программы
clear
t=0:0.001:1; %Задаем массив отсчетов времени
fc=100; %Задаем частоту несущего сигнала
fm=10; %Задаем частоту модулирующего сигнала
Uc0=3; %Задаем амплитуду несущего колебания
Um0=2; %Задаем амплитуду модулирующего сигнала
U0=3; %Задаем постоянную составляющую модулирующего сигнала
Un0=0.1; %Задаем амплитуду белого шума
fi0=pi/3; %Задаем начальную фазу модулирующего колебания
m=1; %Задаем коэффициент модуляции
Uc=Uc0*cos(2*pi*fc*t); %Задаем несущее гармоническое колебание
Um=U0+Um0*cos(2*pi*fm*t+fi0); %Задаем модулирующий гармонический сигнал
Uam=Uc0.*(1+m*Um/max(abs(Um))).*cos(2*pi*fc*t); %Выражение для АМ-сигнала
Un=Un0*randn(size(t)); %Генерирует белый шум с МО=0 и СКО=1
Ud=Uam+Un; %Сигнал на входе приемника (аддитивная помеха)
b=[1 1 1]; %Коэффиенты b цифрового фильтра
|
|
|
a=1; %Коэффициенты a цифрового фильтра
Udmd=filter(b,a,abs(Ud));
subplot(3,2,1)
plot(t,Um)
gridon
xlabel('Время, с')
ylabel('Амплитуда, B')
title ('Модулирующийсигнал')
subplot(3,2,2)
plot(t,Uc)
gridon
xlabel('Время, с')
ylabel('Амплитуда, B')
title ('Несущееколебание')
subplot(3,2,3)
plot(t,Uam)
gridon
xlabel('Время, с')
ylabel('Амплитуда, B')
title ('АМсигнал')
subplot(3,2,4)
plot(t,Un)
gridon
xlabel('Время, с')
ylabel('Амплитуда, B')
title ('Белыйшум')
subplot(3,2,5)
plot(t,Ud)
gridon
xlabel('Время, с')
ylabel('Амплитуда, B')
title ('Сигнал на входе приемника')
subplot(3,2,6)
plot(t,Udmd)
gridon
xlabel('Время, с')
ylabel('Амплитуда, B')
title ('Отфильтрованный сигнал')
Задание 2. Фильтрация сигналов для демодуляции амплитудно-манипулированных сигналов в гауссовских каналах связи.
Проделаем задание 2 аналогично заданию 1 с той лишь разницей, что в качестве модулирующего сигнала будет выступать битовая последовательность.
В задании использовался модифицированный код из предыдущего задания.
Рис 2. Временное представление модулирующей битовой последовательности, детектируемого сигнала и демодулированного сигналов 
;
Задание 3. Обнаружение детерминированного импульсного сигнала 
на фоне АБГШ.
Порядок выполнения работы:
1) Введем известные данные и посчитаем порог Байеса для принятия решения 
:
|
|
|
;
;
;
;
;
– гипотеза о том, что в сигнале 
присутствует 
;
– гипотеза о том,в сигнале отсутствует ;
2) Из теории проверки статистических гипотез имеем:
Необходимо проверить соотношение:
Так как помеха, сгенерированная в задании 2, имеет математическое ожидание МО=0 и среднеквадратическое отклонение СКО=1, то:
Листинг программы
t=0:0.001:1; %Задаем массив отсчетов времени
fc=100; %Задаем частоту несущего сигнала
fm=10; %Задаем частоту модулирующего сигнала
Uc0=3; %Задаем амплитуду несущего колебания
Um0=2; %Задаем амплитуду модулирующего сигнала
U0=3; %Задаем постоянную составляющую модулирующего сигнала
fi0=pi/3; %Задаем начальную фазу модулирующего колебания
m=1; %Задаем коэффициент модуляции
Un0=0.3; %Задаем амплитуду белого шума
Uc=Uc0*cos(2*pi*fc*t); %Задаем несущее гармоническое колебание
Um=U0+Um0*cos(2*pi*fm*t+fi0); %Задаем модулирующий гармонический сигнал
Uam=Uc0.*(1+m*Um/max(abs(Um))).*cos(2*pi*fc*t); %Выражение для АМ-сигнала
Un=Un0*randn(size(t)); %Генерирует белый шум с МО=0 и СКО=1
Ud=Uam+Un; %Сигнал на входе приемника (аддитивная помеха)
U1=Un.^2;
U2=Um.^2;
U3=U1+U2;
U4=U3';
L=2.71828^(sum(U4));
disp('Отношение правдоподобий')
disp(L)
, следовательно, на основании критерия Байеса принимается гипотеза (в сигнале присутствует ).
|
|
|
Дата добавления: 2018-02-28; просмотров: 341; Мы поможем в написании вашей работы! |
Мы поможем в написании ваших работ!