Цифровой согласованный обнаружитель сигналов

Основным назначением любого приемника является выделение (обнаружение) полезного сигнала или его параметров из действующей на входе приемника аддитивной смеси сигнала и помехи (шума). Среди прочих задач можно выделить задачу обнаружения полезного сигнала, она состоит в том, чтобы определить, имеется ли в действующем на входе приемника колебании полезный сигнал или оно образовано только помехой (шумом). Приемник в результате решения этой задачи должен дать ответ типа тАЬдатАЭ или тАЬнеттАЭ, т.е. имеется ли полезный сигнал или нет.

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

Для решения задач обнаружения строят оптимальные приемники, в состав которых входят линейные фильтры, а именно согласованные. Выбор критерия оптимальности определяется решаемой задачей.

В данной курсовой работе входными сигналами является пачка, состоящая из линейных частотно тАУ манипулированных (ЛЧМ) сигналов с заданными значением частоты девиации.

1 Анализ технического задания

Согласно заданию, необходимо спроектировать цифровой согласованный обнаружитель сигналов для пачки, состоящей из пяти ЛЧМ тАУ импульсов с заданным значением частоты девиации, а именно девиация равна 1 МГц, а частота входного сигнала 4 МГц, тогда по теореме Котельникова (частота дискретизации должна быть как минимум в 2 раза больше максимальной частоты сигнала) возьмем частоту дискретизации равной 10 МГц.

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

В· передаваемые сигналы полностью известны

В· канал связи имеет известные параметры

В· помеха носит аддитивный характер

В· синхронизация сигналов является идеальной

Согласно этим условиям и осуществляется когерентный прием.

Так как параметры ЛЧМ - сигнала должны быть известны, то примем изначально, что фаза равна нулю.

Структуру обнаружителя можно определить следующим образом:

Рисунок 1- Структурная схема цифрового согласованного обнаружителя сигналов

2 Линейные частотно тАУ манипулированные сигналы

Подобный сигнал изображен на рисунке 2,а, а закон изменения частоты заполнения импульса тАУ на рисунке 2,б.

Рисунок 2 тАУ ЛЧМ тАУ импульс (а) и

изменение частоты его заполнения(б).

Мгновенную частоту заполнения Ваможно определить выражением:

ВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВа (1)

где ВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВа (2)

есть скорость линейного изменения частоты внутри импульса. Тогда мгновенное значение колебания, представленного на рисунке 2, а, можно записать в виде:

Ва,ВаВаВа ВаВаВаВаВаВаВаВаВаВаВаВаВаВа Ва(3)

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

является основным параметром ЛЧМ тАУ сигнала.[1]

3 Согласованный фильтр

Оптимальный по критерию максимума отношения:

Ва, ВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВа (5)

где Ва- возможный максимум (пика) сигнал y(nT);

T тАУ период дискретизации;

Ва- среднеквадратичное (эффективное) значение шума линейный фильтр называется согласованным (i) с сигналом.

Согласованный фильтр, являясь линейным, полностью описывается импульсной hi(nT) и частотной Hi(ejwt) характеристиками, которые связаны между собой преобразованием Фурье:

ВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВаВа (6)

Импульсная характеристика i является Влзеркальным отражениемВ» согласованного с ним сигнала:

Соответственно амлитудо тАУ частотная характеристика (АЧХ) определится как модуль частотной:

а фазочастотная характеристика(ФЧХ) как аргумент частотной. [2]

4 Моделирование

Моделирование будем выполнять в программной среде MATLAB 7.0. Для начала необходимо создать один ЛЧМ - сигнал, а затем сформировать из пяти таких импульсов пачку, пропустить через i, подать на накопитель для формирования сжатого сигнала, после чего детектировать с помощью детектора. Затеи отправить на пороговое устройство, чтобы определить наличие полезного сигнала и дать решение о том, сигнал был или нет. Рассмотрим каждый из этих этапов более подробно.

4.1 Создание ЛЧМ импульса

Исходными данными являются частота входного сигнала f0= 4МГц, девиация равна 1 МГц, это значит, что частота линейно изменяется от 3 МГц до 5 МГц. Длительность импульса равна 50 мкс.

Для того чтобы построить ЛЧМ тАУ импульс, необходимо определить и построить закон изменения частоты. В результате получаем , изобразим его:

Рисунок 3 тАУ Закон изменения частоты заполнения

Построим сам импульс в соответствии с данным законом изменения частоты заполнения.

Рисунок 4 тАУ ЛЧМ тАУ импульс

Теперь сформируем пачку из пяти таких импульсов периодом следования Т=500 мкс и добавим в канал аддитивный шум, так как он необходим для когерентного приема. Получаем:

Рисунок 5 тАУ Изображение шума

Рисунок 6 тАУ ЛЧМ тАУ импульс с шумом

Рисунок 7 тАУ Пачка из пяти ЛЧМ тАУ импульсов с шумом

Как видно из рисунка 5, амплитуда шума больше амплитуды сигнала, поэтому отношение С/Ш в данном случае равно 0.0683.

4.2 Создание согласованного фильтра и его импульсной характеристики

При синтезе исходим из того, что импульсная характеристика согласованного фильтра (i) должна представлять собой тАЬзеркальнуютАЭ копию выделяемого сигнала с обращенным во времени порядком следования отдельных позиций. Получаем:

Рисунок 8 тАУ Изображение импульсной характеристики i без шума

Рисунок 9 тАУ Изображение импульсной характеристики i с шумом

Все дальнейшие операции, то есть накопление и детектирование, будем производить при наличии аддитивного шума, изображенного на рисунке 5.

4.3 Прохождение через согласованный фильтр

Как было сказано выше, в канале присутствует шум.

Рисунок 10 тАУ Реакция i на один ЛЧМ тАУ импульс

4.4 Создание накопителя и прохождение через него

обнаружитель сигнал частотный фильтр

Накопитель создадим с помощью блоков повторения через период Т=500мкс, таким образом, их будет четыре, так как необходимо получить пачку из пяти импульсов, и с помощью сумматора, который суммирует их и получит сжатый сигнал . После всех этих операций получаем:

Рисунок 11 тАУ Сигнал на выходе накопителя

4.5 Создание детектора и прохождение через него

Как было показано в структурной схеме, детектор состоит из устройства взятия модуля и ФНЧ. Для подавления лепестков сжатого сигнала используем метод взвешивания импульсной характеристики. Новая импульсная характеристика формируется по правилу:

h' (nT) = W (nt) * h (nT)

где W (nT) - весовая функция или "сглаживающее окно".

Находят применение различные типы окон, например "окно ХэммингатАЭ:

W=0.42+0.5*cos(2*pi*n/N)+0.08*cos(4*pi*n/N); с помощью этого окна ВлвзвесимВ» импульсную характеристику ФНЧ. На выходе детектора получаем следующий сигнал:

Рисунок 12 тАУ Сигнал на выходе детектора

4.6 Создание порогового устройства и анализ полученного результата

Значение порога выбираем экспериментально, возьмем значение порога равным пяти, сравним возможный максимум (пик) полученного сигнала с этим порогом, при этом примем, что если полезный сигнал присутствует, то на выходе порогового устройства выдается прямоугольный импульс с амплитудой 1 на заданном интервале времени. Получаем следующий результат.

Рисунок 13 тАУ Сигнал на выходе порогового устройства

По рисунку 13, определяем, что в канале в интервале от 48мкс до 52 мкс присутствует полезный сигнал.

Здесь продемонстрирован прием лишь при одном значении С/Ш, использовали это значение, так как оно показывает наиболее реальную ситуацию, когда шум больше сигнала. Таким образом, можно сделать вывод, сто чем больше отношения С/Ш, тем лучше будет прием сигналов.

5 Функциональная схема цифрового согласованного обнаружителя сигналов

ВаВл1В»

Рисунок 14 - Функциональная схема цифрового согласованного обнаружителя сигналов

Заключение

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

Была сформирована пачка из пяти ЛЧМ тАУ импульсов , для этих импульсов были заданны длительность и периодом повторения, а также частота входного сигнала и девиация, была получена импульсная характеристика используемого согласованного фильтра.

Согласно заданию было необходимо организовать согласованный прием каждого импульса, что и было сделано с помощью среды моделирования MATLAB.

Приложение 1 : Текст написанной программы

clc;

clear all;

Fd=10; % Частота дискретизации, МГц

f0=4; % Частота сигнала, МГц

t=0:1/Fd:2500; % время в микросекундах

T=1/Fd; % период дискретизации

f=0.04*t+3;%построим закон изменения частоты

figure (1);

plot (t(1:500),f(1:500));

grid on;

xlabel('time,microsec');

ylabel ('frequency,MegaHz');

title('Zakon izmenenyia');

% создание шума

Noise=randn(1,25001);

figure (2); % вывод на экран изображение шума

plot (t(1:25000),Noise (1:25000));

xlabel ('Time, microsec');

ylabel ('Amplitude, V');

grid on;

title ('Noise');

% % создание ЛЧМ-импульса

ВаS=cos(2*pi*(0.04*t+3).*t)+ Noise;

% S=cos(2*pi*(0.04*t+3).*t)

figure (3);

plot(t(1:500),S(1:500));

grid on;

xlabel('time, microsec');

ylabel ('amplitude, V');

title('Signal')

% создаем соглаcованный фильтр

for q=1:1:500

h_SF(q)=S(501-q);

end

figure (4);

plot (t(1:500),h_SF(1:500));

grid on;

xlabel('frequency,MegaHz');

ylabel ('amplitude, V');

title('Impulse haracneristic')

% реакция согласованного фильтра на один ЛЧМ-импульс

ReactionSF1=filter (h_SF,1,S)./(0.5*50*10);

figure (5); % вывод на экран реакции

plot (t(1:1000),ReactionSF1(1:1000));

grid on;

xlabel('time, microsec');

ylabel ('amplitude, V');

title('Reaction for 1 imp');

% создание пачки из 5 ЛЧМ-импульсов

P=zeros(1,25000);

for i=1:5000:25000

P(i+1:i+500) = S(1:500);

end

figure(6);

plot(t(1:25000),P(1:25000));

grid on;

xlabel('time, microsec');

ylabel ('amplitude, V');

% Вычисляем уровень сигнал\шум

Noise_v_kv = power (Noise,2);

z=trapz(t,Noise_v_kv); % вычисление среднеквадратичного значения шума

P_N=max(abs(P))./sqrt(z);

disp (P_N);

% формируем накопитель

for q=1:1:500;

for i=1:5000:25000;

Ваh_1(q) = P((501-q)+i);

Ваend

end

ReactionSF = filter(h_1,1,P)./(0.5*50*10);

figure (7);

plot(t(1:25000),ReactionSF(1:25000));

grid on;

xlabel ('time, microsec');

ylabel ('amplitude, V');

title('Reaction for 5 imp');

for i = 5000:5000:25000

ВаReaction=ReactionSF(1:1000)+ReactionSF(i+1:i+1000);

figure (8);

plot(t(1:1000),Reaction(1:1000));

grid on;

xlabel ('time, microsec');

ylabel ('amplitude');

title('On exit of nakopitel')

%создание детектора

reaction=abs(Reaction);

figure (9);

plot(t(1:1000),reaction(1:1000));

%создание ФНЧ, входящего в состав детектора

N =10;

f1=0.2; f2=0.3;

f0=(f1+f2)/2;

d1=3; d2=50;

n = -(N-1)/2:1:(N-1)/2; %отсчеты

N0=(f2-f1)/2;

Hn =sin(f0*2*pi*n)./(n*pi); %задание ИХ

[H,w]=freqz(Hn,1,1024);% АЧХ фильтра

W1=HAMMING(N)'; % применяем окно в качестве взвешивающей функции

[Hw1,ww1]=freqz(W1,1,1024); % часть характеристики окна

h1=Hn.*W1;

[H1,w1]=freqz(h1,1,1024);

ВаF_h=filter(h1,1,reaction);

figure(10);

Ваplot(t(1:1000),F_h(1:1000));

grid on;

xlabel ('Time, microsec')

ylabel ('Amplitude');

title('On exit of detector');

Ва% пороговое устройство

priem=zeros (1, 1000);

Ваif max (reaction)> 2; % порог

Ваpriem(480:520) = 1;% max (vuhod_modul);

Ваdisp('сигнал был');% вывод на экран

Ваend

Ваfigure (11);

Ваplot(t(1:1000),priem(1:1000));

Ваgrid on;

xlabel ('Time, microsec');

ylabel ('Amplitude, V');

title ('Priem');

Вместе с этим смотрят:

IP-телефония. Особенности цифровой офисной связи

РЖсторiя звтАЩязку та його розвиток

Анализ режимов автоматического управления

Аргоновый лазер

Архитектуры реализации корпоративных информационных систем