216849 (588966), страница 4

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

Текст из файла (страница 4)

-частота модуляции,

-глубина модуляции,

-число тактов обработки,

-частотный диапазон прослушивания,

-угол падения фронта волны на элементы АР,

- частотный сдвиг для выбранного ЧД,

- степень расширения полосы частот (количество отсчетов),

- спектральное окно.

1. Формирование процессов на выходе приемных элементов АР. Имитационное моделирование входного процесса включает в себя задание моделей полей сигнала, формирование временной реализации сигнала, фильтрацию ее в заданном диапазоне частот, формирование плоского фронта сигнала на элементах АР.

2. Перечень процедур, проводимых в этом функциональном блоке:

• расчет координат приемных элементов в АР,

• формирование сигнала заданной частоты в поле,

• расчет временных задержек фронта волны сигнала исходя из направления его прихода,

• формирование дискретизированных с частотой fd=24000Гц сигналов на приемных элементах АР.

1. Функциональные программы расчета выходного эффекта тракта прослушивания:

   • Последовательный набор по n временных отсчетов сигналов на элементах АР с перекрытием 25% - так называемый такт обработки;

Дальнейшие процедуры выполняются на каждом такте обработки:

• Переход в частотную область с использованием процедуры БПФ;

• Формирование ПК на частотах рабочей полосы;

• Вырезание полосы частот выбранного частотного диапазона (ЧД);

• Сдвиг полосы частот на заранее рассчитанное для каждого ЧД число спектральных отсчетов;

• Умножение отсчетов спектра сформированного канала на спектральное окно;

• Восстановление сигнала во временной области (процедура обратного ПФ);

• Отбрасывание некорректных отчетов из реализации – по n/8 отсчетов в начале и в конце реализации;

• набор реализации большой длительности.

6. Выдача сигнала на ЦАП, вычисление спектра сигнала с высоким разрешением по частоте, графическое представление результатов обработки.

4.3 Расчет выходных эффектов тракта

Сформированный указанным выше способом сигнал на элементах антенны «нарезается» на временные фрагменты длиной n=512 точек с учетом перекрытия, которое вводится в модель тракта как параметр. Каждый фрагмент подвергается процедуре БПФ, в результате чего получаются спектры процессов на выходе приемных элементов. Из получившегося набора частот в дальнейшую обработку берутся только номера частот f_k, соответствующие заданной полосе обработки (f_н÷f_в).

Как было показано выше, вырезание полосы фильтром с прямоугольной частотной характеристикой вносит дополнительные искажения в восстановленный сигнал. Поэтому в модели предусмотрено расширение полосы обработки от значения (0,3-4,5) кГц до (0,05-5) кГц.

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

Рис. 19 Структурная схема программного макета тракта прослушивания

Рис.20 Обработка в канале прослушивания

После ФХН выполняется вырезание полосы частот выбранного частотного диапазона и сдвиг ее на заранее рассчитанное для каждого ЧД число спектральных отсчетов. Затем проводится процедура ОБПФ, на выходе которой получаем временную реализацию. Как было показано ранее, для устранения эффектов, вызванных вырезанием полосы частот, необходимо перед восстановлением сигнала во временную область ввести частотное окно, вид которого в модели можно выбирать из заранее созданных файлов либо из имеющегося в системе MatLab банка частотных фильтров, а также задавать их параметры. В модели тракта прослушивания параллельно формируются и обрабатываются сигналы, восстановленные с применением частотного окна и без него.

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

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

5 Результаты расчетов

Чтобы сформировать сигнал на элементах АР нужно рассчитать задержки на каждом элементе и сымитировать их.

Исходя из заданных нами параметров, система формирует приходящий на АР сигнал под определенным углом.

Рис.22 Приход сигнала под углом на элементы АР

Обработка в канале прослушивания

Блок основных вычислительных операций подробно будет рассмотрен далее.

Вывод графических и аудио результатов

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

По сформированным входным воздействиям реализуется обработка в канале прослушивания .

Сформированный сигнал на элементах АР

Исходя из заданных нами параметров, система формирует приходящий на АР сигнал под определенным углом.

БПФ с перекрытием

Выполняется БПФ с перекрытием, перекрытие задается произвольным, либо кратно числу точек БПФ (как и выяснилось в результате моделирования), но задать перекрытие мы можем произвольно. Выполняется перекрытие на каждом приемном элементе АР.

Умножение на фазирующие коэффициенты

Этап процедуры формирования. внесение задержек в частотной области

Суммирование по элементам АР

Вырезание полосы

Поскольку у нас частотное окно имеет спады на краях диапазона, мы расширяем полосу сигнала, чтобы не сильно зарезать энергетику сигнала.

Сдвиг

В первом частотном диапазоне полоса частот 1-2,5 кГц это соответствует спектральным отсчетам: К нижнее =21, К верхнее =53. В соответствии с такой закономерностью сдвигать при выбранной величине перекрытия мы можем только порциями по 4Δf значит первый спектральный отсчет ,для того чтобы попасть в полосу 3-3.5 кГц должен быть сдвинут на 4 отсчетов . Соответственно для второго диапазона 2-5 кГц : К нижнее =42, К верхнее =107 , 32 отсчетов. Для третьего 4-8 Кгц: К нижнее =85, К верхнее =171, 76 отсчетов.

Рассчитывать коэффициент сдвига будем проводить по формуле

Где нижнее значение частотного диапазона,

= коэффициент фильтра, - число сдвига спектральных отсчетов.

7 – коэффициент минимальной частоты переноса спектра = 300 ГЦ

Умножение на спектральное окно

Для каждого из трех диапазонов мы выбираем свои параметры фильтра . В данной случае будем менять лишь параметры фильтра Ханна (53 81 99). Полученные по формуле:

K=kB-kH+2*kl+1

ОБПФ

Обратное преобразование Фурье с количеством точек равным в БПФ

Отбрасывание некорректных отсчетов

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

Стыковка реализаций

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

MATLAB - это интерактивная система, в которой основным элементом данных является массив. Это позволяет решать различные задачи, связанные с техническими вычислениями, особенно в которых используются матрицы и вектора, в несколько раз быстрее, чем при написании программ с использованием "скалярных" языков программирования, таких как Си или Фортран.

В результате работы программы мы получаем следующие результаты :

Для первого частотного диапазона (1-2.5) Khz

Каждый частотный диапазон соответствует спектральным отсчетам k1=[21] k2=[53]

Угол прихода сигнала берем равный Pi/4

Сдвиг полосы в область от 0.3 , 4Δf = 187 Гц

Умножение на спектральное окно K=53

Частота сигнала 2000 Гц

Спектр исходного сигнала Спектр восстановленного сигнала

со сдвигом

Исходный сигнал во временной области Восстановленный сигнал во

временной области

Для второго частотного диапазона (2-5) Khz

Каждый частотный диапазон соответствует спектральным отсчетам k2=[42] k2=[107]

Угол прихода сигнала берем равный Pi/4

Сдвиг полосы в область от 0.3 , 32Δf = 1500 Гц

Умножение на спектральное окно K=81

Частота сигнала 3500 Гц

Спектр исходного сигнала Спектр восстановленного

сигнала со сдвигом

Исходный сигнал во Восстановленный сигнал

временной области во временной области

Для третьего частотного диапазона (4-8) Khz

Каждый частотный диапазон соответствует спектральным отсчетам k2=[85] k2=[171]

Угол прихода сигнала берем равный Pi/4

Сдвиг полосы в область от 0.3 , 76Δf = 3570 Гц

Умножение на спектральное окно K=99

Частота сигнала 5000 Гц

Спектр исходного сигнала Спектр восстановленного

Характеристики

Список файлов ВКР