поэтапный. WidebandTwoRayChannel

Широкополосный канал распространения 2D луча

Описание

Модели phased.WidebandTwoRayChannel широкополосный канал распространения 2D луча. Канал распространения 2D луча является самым простым типом многопутевого канала. Можно использовать канал 2D луча, чтобы моделировать распространение сигналов в гомогенном, изотропном носителе с одним контуром отражения. Этот тип носителя имеет два пути к распространению: угол обзора (прямой) путь к распространению от одной точки до другого и путь к лучу, отраженный от контура.

Можно использовать эту Систему object™ для ближнего радара и приложений мобильной связи, где сигналы распространяют вдоль прямых путей, и земля принята, чтобы быть плоской. Можно также использовать этот объект для приложений микрофона и гидролокатора. Для акустических приложений можно выбрать неполяризованные поля и настроить скорость распространения, чтобы быть скоростью звука в воздухе или воде. Можно использовать phased.WidebandTwoRayChannel для образцового распространения от нескольких точек одновременно.

Несмотря на то, что Системный объект работает на все частоты, модели затухания для атмосферных газов и дождя допустимы для электромагнитных сигналов в частотном диапазоне 1-1000 ГГц только. Модель затухания для вуали и облаков допустима для 10-1000 ГГц. Вне этих частотных диапазонов Системный объект использует самое близкое допустимое значение.

Системный объект phased.WidebandTwoRayChannel применяет зависимые областью значений задержки к сигналам, а также прибыли или убытки, сдвиги фазы и граничную отражательную потерю. Когда или источник или место назначения перемещаются, Системный объект применяет эффекты Доплера к сигналам.

Сигналы в канале вывод могут быть разделены или объединены. Если вы разделяете сигналы, оба сигнала прибывают к месту назначения отдельно и не объединены. Если вы принимаете решение объединить сигналы, два сигнала из источника распространяют отдельно, но когерентно суммированы в месте назначения в одно количество. Выберите эту опцию, когда различие между датчиком или усилениями массивов в направлениях этих двух путей будет незначительно.

В отличие от phased.WidebandFreeSpace и Системных объектов phased.WidebandLOSChannel, этот Системный объект не поддерживает двухстороннее распространение.

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

Задайте и настройте свой канал 2D луча. Смотрите Конструкцию.
Вызовите метод step, чтобы вычислить распространенный сигнал с помощью свойств Системного объекта phased.WidebandTwoRayChannel.

Примечание

Также вместо того, чтобы использовать метод step, чтобы выполнить операцию, заданную Системным объектом, можно вызвать объект с аргументами, как будто это была функция. Например, y = step(obj,x) и y = obj(x) выполняют эквивалентные операции.

Конструкция

channel = phased.WidebandTwoRayChannel создает Системный объект канала распространения 2D луча, channel.

channel = phased.WidebandTwoRayChannel(Name,Value) создает Системный объект, channel, с каждым заданным набором свойства Name к заданному Value. Можно задать дополнительное имя и аргументы пары значения в любом порядке как (Name1,Value1..., NameN,ValueN).

Свойства

развернуть все

`PropagationSpeed` — Скорость распространения сигнала
`physconst('LightSpeed')` (значение по умолчанию) | положительная скалярная величина

Скорость распространения сигнала, заданная как положительная скалярная величина. Модули исчисляются в метрах в секунду. Скорость распространения по умолчанию является значением, возвращенным physconst('LightSpeed'). Смотрите physconst.

Пример: 3e8

Типы данных: double

`OperatingFrequency` — Рабочая частота
`300e6` (значение по умолчанию) | положительная скалярная величина

Рабочая частота, заданная как положительная скалярная величина. Модули находятся в Гц.

Пример: 1e9

Типы данных: double

`SpecifyAtmosphere` — Включите атмосферную модель затухания
`false` (значение по умолчанию) | `true`

Опция, чтобы включить атмосферную модель затухания, заданную как false или true. Установите это свойство на true добавлять затухание сигнала, вызванное атмосферными газами, дождем, вуалью или облаками. Установите это свойство на false игнорировать атмосферные эффекты в распространении.

`NumSubbands` — Количество обработки поддиапазонов
`64` (значение по умолчанию) | положительное целое число

Количество обработки поддиапазонов, заданных как положительное целое число.

Пример: 128

Типы данных: double

`EnablePolarization` — Включите поляризованные поля
`false` (значение по умолчанию) | `true`

Опция, чтобы включить поляризованные поля, заданные как false или true. Установите это свойство на true включать поляризацию. Установите это свойство на false игнорировать поляризацию.

Типы данных: логический

`GroundReflectionCoefficient` — Оснуйте отражательный коэффициент
`-1` (значение по умолчанию) | скаляр с комплексным знаком | 1 с комплексным знаком N вектором - строкой

Оснуйте отражательный коэффициент для поля в отражательной точке, заданной как скаляр с комплексным знаком или 1 с комплексным знаком N вектором - строкой. Каждый коэффициент имеет абсолютное значение, меньше чем или равное одному. Количество N является количеством каналов 2D луча. Модули являются безразмерными. Используйте это свойство для неполяризованных сигналов модели. К поляризованным сигналам модели используйте свойство GroundRelativePermittivity.

Пример: -0.5

Зависимости

Чтобы включить это свойство, установите EnablePolarization на false.

Типы данных: double
Поддержка комплексного числа: Да

`GroundRelativePermittivity` — Оснуйте относительную проницаемость
`15` (значение по умолчанию) | положительный скаляр с действительным знаком | 1 с действительным знаком Nrow вектором положительных значений

Относительная проницаемость земли в отражательной точке, заданной как положительный скаляр с действительным знаком или 1 N вектором - строкой с действительным знаком из положительных значений. Размерность N является количеством каналов 2D луча. Модули проницаемости являются безразмерными. Относительная проницаемость задана как отношение фактической наземной проницаемости к проницаемости свободного пространства. Это свойство применяется, когда вы устанавливаете свойство EnablePolarization на true. Используйте это свойство для поляризованных сигналов модели. К неполяризованным сигналам модели используйте свойство GroundReflectionCoefficient.

Пример 5

Зависимости

Чтобы включить это свойство, установите EnablePolarization на true.

Типы данных: double

`CombinedRaysOutput` — Опция, чтобы объединить два луча при выводе
`true` (значение по умолчанию) | `false`

Опция, чтобы объединить два луча в канале вывод, заданный как true или false. Когда этим свойством является true, объект когерентно добавляет угол обзора распространенный сигнал и отраженный сигнал пути при формировании выходного сигнала. Используйте этот режим, когда вы не должны будете включать направленное усиление антенны или массива в вашей симуляции.

Типы данных: логический

`MaximumDistanceSource` — Источник максимального одностороннего расстояния распространения
`'Auto'` (значение по умолчанию) | `'Property'`

Источник максимального одностороннего расстояния распространения, заданного как 'Auto' или 'Property'. Максимальное одностороннее расстояние распространения используется, чтобы выделить достаточную память для вычисления задержки сигнала. Когда вы устанавливаете это свойство на 'Auto', Системный объект автоматически выделяет память. Когда вы устанавливаете это свойство на 'Property', вы задаете максимальное одностороннее расстояние распространения с помощью значения свойства MaximumDistance.

Типы данных: char

`MaximumDistance` — Максимальное одностороннее расстояние распространения
`10000` (значение по умолчанию) | положительный скаляр с действительным знаком

Максимальное одностороннее расстояние распространения, заданное как положительный скаляр с действительным знаком. Модули исчисляются в метрах. Любой сигнал, который распространяет больше, чем максимальное одностороннее расстояние, проигнорирован. Максимальное расстояние должно быть больше, чем или равным самому большому расстоянию положения-к-позиционному.

Пример: 5000

Зависимости

Чтобы включить это свойство, установите свойство MaximumDistanceSource на 'Property'.

Типы данных: double

`MaximumNumInputSamplesSource` — Источник максимального количества выборок
`'Auto'` (значение по умолчанию) | `'Property'`

Источник максимального количества выборок входного сигнала, заданного как 'Auto' или 'Property'. Когда вы устанавливаете это свойство на 'Auto', модель распространения автоматически выделяет достаточно памяти, чтобы буферизовать входной сигнал. Когда вы устанавливаете это свойство на 'Property', вы задаете максимальное количество выборок во входном сигнале с помощью свойства MaximumNumInputSamples. Любой входной сигнал дольше, чем то значение является усеченным.

Когда вы используете модели распространения в MATLAB^® Function Block в Simulink^® с сигналами переменного размера, необходимо установить параметр MaximumNumInputSamples на 'Property' и установить параметр MaximumNumInputSamples соответственно.

Пример: 'Property'

Зависимости

Чтобы включить это свойство, установите MaximumDistanceSource на 'Property'.

Типы данных: char

`MaximumNumInputSamples` — Максимальное количество выборок входного сигнала
`100` (значение по умолчанию) | положительное целое число

Максимальное количество выборок входного сигнала, заданных как положительное целое число. Входной сигнал является первым аргументом метода step после самого Системного объекта. Размер входного сигнала является количеством строк во входной матрице. Любой входной сигнал дольше, чем этот номер является усеченным. К сигналам процесса полностью, гарантируйте, что это значение свойства больше, чем какая-либо максимальная длина входного сигнала.

Генерирующие форму волны Системные объекты определяют максимальный размер сигнала:

Для любой формы волны, если свойство OutputFormat формы волны установлено в 'Samples', максимальная длина сигнала является значением, заданным в свойстве NumSamples.
Для импульсных форм волны, если OutputFormat установлен в 'Pulses', длина сигнала является продуктом наименьшей импульсной частоты повторения, количеством импульсов и частотой дискретизации.
Для непрерывных форм волны, если OutputFormat установлен в 'Sweeps', длина сигнала является продуктом времени развертки, количеством разверток и частотой дискретизации.

Пример: 2048

Зависимости

Чтобы включить это свойство, установите MaximumNumInputSamplesSource на 'Property'.

Типы данных: double

Методы

сброс	Сбросьте состояния Системного объекта
шаг	Распространите широкополосный сигнал от точки до точки с помощью модели канала 2D луча

Характерный для всех системных объектов
`release`	Позвольте изменения значения свойства Системного объекта

Примеры

развернуть все

Скалярное широкополосное распространение сигнала в канале 2D луча

Скрипт Open Live Script

Этот пример иллюстрирует распространение 2D луча широкополосного сигнала, показывая, как сигналы от пути угла обзора и отраженного пути прибывают в получатель в разное время.

Примечание: можно заменить каждый вызов функции с эквивалентным синтаксисом step. Например, замените myObject(x) на step(myObject,x).

Создайте и постройте переданную форму волны

Создайте неполяризованное электромагнитное поле, состоящее из двух линейных импульсов формы волны FM в несущей частоте 100 МГц. Примите, что ширина импульса является 20 μs, и уровень выборки составляет 10 МГц. Пропускная способность импульса составляет 1 МГц. Примите 50%-й рабочий цикл так, чтобы ширина импульса была половиной импульсного интервала повторения. Создайте 2D импульсный train волны. Установите GroundReflectionCoefficient на –0.9 моделировать отражающую способность веского основания. Распространите поле со стационарного источника на стационарный получатель. Вертикальное разделение источника и получателя составляет приблизительно 10 км.

c = physconst('LightSpeed');
fs = 10e6;
pw = 20e-6;
pri = 2*pw;
PRF = 1/pri;
fc = 100e6;
lambda = c/fc;
bw = 1e6;
waveform = phased.LinearFMWaveform('SampleRate',fs,'PulseWidth',pw,...
    'PRF',PRF,'OutputFormat','Pulses','NumPulses',2,'SweepBandwidth',bw,...
    'SweepDirection','Down','Envelope','Rectangular','SweepInterval',...
    'Positive');
wav = waveform();
n = size(wav,1);
plot([0:(n-1)]/fs*1e6,real(wav),'b')
xlabel('Time (\mu s)')
ylabel('Waveform Magnitude')

Задайте местоположение источника и получателя

Поместите источник и получатель на расстоянии приблизительно в 1 км горизонтально и на расстоянии приблизительно в 5 км вертикально.

pos1 = [0;0;100];
pos2 = [1e3;0;5.0e3];
vel1 = [0;0;0];
vel2 = [0;0;0];

Создайте широкополосный системный объект канала 2D луча

Создайте Систему канала распространения 2D луча object™ и распространите сигнал вдоль обоих угол обзора и отраженные пути к лучу. Тот же сигнал распространен вдоль обоих путей.

channel = phased.WidebandTwoRayChannel('SampleRate',fs,...
    'GroundReflectionCoefficient',-0.9,'OperatingFrequency',fc,...
    'CombinedRaysOutput',false);
prop_signal = channel([wav,wav],pos1,pos2,vel1,vel2);

[rng2,angs] = rangeangle(pos2,pos1,'two-ray');

Вычислите задержки μs.

tm = rng2/c*1e6;
disp(tm)

   16.6815   17.3357

Отобразите расчетный азимут путей к распространению и углы повышения в градусах.

disp(angs)

         0         0
   78.4654  -78.9063

Постройте распространенные сигналы

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

n = size(prop_signal,1);
delay = [0:(n-1)]/fs*1e6;
subplot(3,1,1)
plot(delay,real([prop_signal(:,1)]),'b')
grid
xlabel('Time (\mu sec)')
ylabel('Real Part')
title('Direct Path')

subplot(3,1,2)
plot(delay,real([prop_signal(:,2)]),'b')
grid
xlabel('Time (\mu sec)')
ylabel('Real Part')
title('Reflected Path')

subplot(3,1,3)
plot(delay,real([prop_signal(:,1) + prop_signal(:,2)]),'b')
grid
xlabel('Time (\mu sec)')
ylabel('Real Part')
title('Combined Paths')

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

Сравните широкополосное распространение канала 2D луча со свободным пространством

Скрипт Open Live Script

Вычислите результат распространения широкополосного сигнала LFM в среде 2D луча от радара на 10 метров выше источника (0,0,10) к цели в (3000,2000,2000) метры. Примите, что радар и цель являются стационарными и что антенна передачи является изотропной. Объедините сигнал от этих двух путей и сравните сигнал с распространением сигнала в свободном пространстве. Система действует на уровне 300 МГц. Установите свойство CombinedRaysOutput на true комбинировать прямой путь и отраженные сигналы пути при формировании выходного сигнала.

Примечание: Этот пример запускается только в R2016b или позже. Если вы используете более ранний релиз, заменяете каждый вызов функции с эквивалентным синтаксисом step. Например, замените myObject(x) на step(myObject,x).

Создайте линейную форму волны FM.

fop = 300.0e6;
fs = 1.0e6;
waveform = phased.LinearFMWaveform();
x = waveform();

Задайте целевое положение и скорость.

posTx = [0; 0; 10];
posTgt = [3000; 2000; 2000];
velTx = [0;0;0];
velTgt = [0;0;0];

Смоделируйте распространение свободного пространства.

fschannel = phased.WidebandFreeSpace('SampleRate',waveform.SampleRate);
y_fs = fschannel(x,posTx,posTgt,velTx,velTgt);

Образцовое распространение 2D луча от положения радара к цели.

tworaychannel = phased.WidebandTwoRayChannel('SampleRate',waveform.SampleRate,...
    'CombinedRaysOutput',true);
y_tworay = tworaychannel(x,posTx,posTgt,velTx,velTgt);
plot(abs([y_tworay y_fs]))
legend('Wideband two-ray (Position 1)','Wideband free space (Position 1)',...
    'Location','best')
xlabel('Samples')
ylabel('Signal Magnitude')
hold on

Переместите радар на 10 метров горизонтально к второй позиции.

posTx = posTx + [10;0;0];
y_fs = fschannel(x,posTx,posTgt,velTx,velTgt);
y_tworay = tworaychannel(x,posTx,posTgt,velTx,velTgt);
plot(abs([y_tworay y_fs]))
legend('Wideband two-ray (Position 1)','Wideband free space (Position 1)',...
    'Wideband two-ray (Position 2)','Wideband free space (Position 2)',...
    'Location','best')
hold off

Потери распространения свободного пространства являются тем же самым и для первых и для вторых положений радара. Потери 2D луча отличаются из-за эффекта взаимодействия путей 2D луча.

Широкополосное поляризованное полевое распространение в канале 2D луча

Скрипт Open Live Script

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

Установите радарные параметры формы волны

Примите, что ширина импульса $10 μ s$ и уровень выборки составляет 10 МГц. Пропускная способность импульса составляет 1 МГц. Примите 50%-й рабочий цикл, в котором ширина импульса является половиной импульсного интервала повторения. Создайте 2D импульсный train волны. Примите несущую частоту 100 МГц.

c = physconst('LightSpeed');
fs = 20e6;
pw = 10e-6;
pri = 2*pw;
PRF = 1/pri;
fc = 100e6;
bw = 1e6;
lambda = c/fc;

Настройте необходимые системные объекты

Используйте GroundRelativePermittivity 10.

waveform = phased.LinearFMWaveform('SampleRate',fs,'PulseWidth',pw,...
    'PRF',PRF,'OutputFormat','Pulses','NumPulses',2,'SweepBandwidth',bw,...
    'SweepDirection','Down','Envelope','Rectangular','SweepInterval',...
    'Positive');
antenna = phased.CrossedDipoleAntennaElement(...
    'FrequencyRange',[50,200]*1e6);
radiator = phased.Radiator('Sensor',antenna,'OperatingFrequency',fc,...
    'Polarization','Combined');
channel = phased.WidebandTwoRayChannel('SampleRate',fs,...
    'OperatingFrequency',fc,'CombinedRaysOutput',false,...
    'EnablePolarization',true,'GroundRelativePermittivity',10);
collector = phased.Collector('Sensor',antenna,'OperatingFrequency',fc,...
    'Polarization','Combined');

Настройте геометрию сцены

Задайте положения передатчика и получателя, скорости и ориентации. Поместите источник и получатель на расстоянии приблизительно в 1 000 м горизонтально и на расстоянии приблизительно в 50 м вертикально.

posTx = [0;100;100];
posRx = [1000;0;150];
velTx = [0;0;0];
velRx = [0;0;0];
laxRx = rotz(180);
laxTx = rotx(1)*eye(3);

Создайте и излучите сигналы от передатчика

Вычислите углы передачи для двух лучей, перемещающихся к получателю. Эти углы заданы относительно системы локальной координаты передатчика. phased.Radiator System object(TM) использует эти углы, чтобы применить отдельные усиления антенны к двум сигналам.

[rng,angsTx] = rangeangle(posRx,posTx,laxTx,'two-ray');
wav = waveform();

Постройте переданную форму волны.

n = size(wav,1);
plot([0:(n-1)]/fs*1000000,real(wav))
xlabel('Time ({\mu}sec)')
ylabel('Waveform')

sig = radiator(wav,angsTx,laxTx);

Распространите сигналы к получателю через канал 2D луча.

prop_sig = channel(sig,posTx,posRx,velTx,velRx);

Получите распространенный сигнал

Вычислите углы приема для двух лучей, прибывающих в получатель. Эти углы заданы относительно системы локальной координаты получателя. phased.Collector System object(TM) использует эти углы, чтобы применить отдельные усиления антенны к двум сигналам.

[rng1,angsRx] = rangeangle(posTx,posRx,laxRx,'two-ray');
delays = rng1/c*1e6

delays = 1×2

    3.3564    3.4544

Соберите и объедините полученные лучи.

y = collector(prop_sig,angsRx,laxRx);

Постройте полученную форму волны.

plot([0:(n-1)]/fs*1000000,real(y))
xlabel('Time ({\mu}sec)')
ylabel('Received Waveform')

2D излучите распространение широкополосной формы волны LFM с атмосферными потерями

Скрипт Open Live Script

Распространите широкополосный линейный сигнал FM в канале 2D луча. Пропускная способность сигнала составляет 15% несущей частоты. Примите, что существует потеря сигнала, вызванная атмосферными газами и дождем. Сигнал распространяет от передатчика, расположенного на уровне метров (0,0,0) в глобальной системе координат к получателю на уровне метров (10000,200,30). Примите, что передатчик и получатель являются стационарными и что у них обоих есть шаблоны антенны косинуса. Постройте полученный сигнал. Установите сухое давление воздуха на 102,0 Па и уровень дождя к 5 мм/час.

Установите радарные параметры формы волны

c = physconst('LightSpeed');
fs = 40e6;
pw = 10e-6;
pri = 2.5*pw;
PRF = 1/pri;
fc = 100e6;
bw = 15e6;
lambda = c/fc;

Настройте радарный сценарий

Создайте необходимые Системные объекты.

waveform = phased.LinearFMWaveform('SampleRate',fs,'PulseWidth',pw,...
    'PRF',PRF,'OutputFormat','Pulses','NumPulses',2,'SweepBandwidth',bw,...
    'SweepDirection','Down','Envelope','Rectangular','SweepInterval',...
    'Positive');
antenna = phased.CosineAntennaElement;
radiator = phased.Radiator('Sensor',antenna);
collector = phased.Collector('Sensor',antenna);
channel = phased.WidebandTwoRayChannel('SampleRate',waveform.SampleRate,...
    'CombinedRaysOutput',false,'GroundReflectionCoefficient',0.95,...
    'SpecifyAtmosphere',true,'Temperature',20,...
    'DryAirPressure',102.5,'RainRate',5.0);

Настройте геометрию сцены. Задайте положения передатчика и получателя и скорости. Передатчик и получатель являются стационарными.

posTx = [0;0;0];
posRx = [10000;200;30];
velTx = [0;0;0];
velRx = [0;0;0];

Задайте передачу и получение радарных ориентаций антенны относительно глобальных координат. Антенна передачи указывает вдоль положительного направления X и точек антенны получения близко к отрицательному направлению X.

laxTx = eye(3);
laxRx = rotx(5)*rotz(170);

Вычислите углы передачи, которые являются углами, под которыми два луча, перемещающиеся к получателю, оставляют передатчик. Система phased.Radiator object™ использует эти углы, чтобы применить отдельные усиления антенны к двум сигналам. Поскольку усиления антенны зависят от направления контура, необходимо передать и получить два луча отдельно.

[~,angTx] = rangeangle(posRx,posTx,laxTx,'two-ray');

Создайте и излучите сигналы от передатчика

Излучите сигналы вдоль направлений передачи.

wavfrm = waveform();
wavtrans = radiator(wavfrm,angTx);

Распространите сигналы к получателю через канал 2D луча.

wavrcv = channel(wavtrans,posTx,posRx,velTx,velRx);

Соберите сигнал в получателе

Вычислите угол, под которым два луча, перемещающиеся из передатчика, прибывают в получатель. Система phased.Collector object™ использует эти углы, чтобы применить отдельные усиления антенны к двум сигналам.

[~,angRcv] = rangeangle(posTx,posRx,laxRx,'two-ray');

Соберите и объедините два полученных луча.

yR = collector(wavrcv,angRcv);

Постройте полученный сигнал

dt = 1/waveform.SampleRate;
n = size(yR,1);
plot([0:(n-1)]*dt*1e6,real(yR))
xlabel('Time ({\mu}sec)')
ylabel('Signal Magnitude')

Больше о

развернуть все

2D излучите пути к распространению

Канал распространения 2D луча является следующим, подходят в сложности от канала свободного пространства, и самый простой случай многопутевой среды распространения. Канал свободного пространства моделирует прямолинейный путь угла обзора от точки 1 до точки 2. В канале 2D луча носитель задан как гомогенный, изотропный носитель с отражающимся плоским контуром. Контур всегда устанавливается в z = 0. Существует самое большее два распространения лучей от точки 1 до точки 2. Первый путь к лучу распространяет вдоль того же пути угла обзора как в канале свободного пространства (см. Системный объект phased.FreeSpace). Путь угла обзора часто называется прямым путем. Второй луч отражается от контура прежде, чем распространить к точке 2. Согласно Закону Отражения, угол отражения равняется углу падения. В ближних симуляциях, таких как системы сотовой связи и автомобильные радары, можно принять, что отражающаяся поверхность, земля или океанская поверхность, является плоской.

phased.TwoRayChannel и задержка распространения модели Системных объектов phased.WidebandTwoRayChannel, сдвиг фазы, эффект Доплера и эффекты потерь для обоих путей. Для отраженного пути эффекты потерь включают отражательную потерю в контур.

Фигура иллюстрирует два пути к распространению. От исходного положения, _ss, и положения получателя, _sr, можно вычислить углы падения обоих путей, _θ′los и _θ′rp. Углы падения являются повышением и углами азимута прибывающего излучения относительно системы локальной координаты. В этом случае система локальной координаты совпадает с глобальной системой координат. Можно также вычислить углы передачи, _θlos и _θrp. В глобальных координатах угол отражения на контуре совпадает с углами _θrp и _θ′rp. Отражательный угол важен, чтобы знать, когда вы используете зависимые углом данные отражательной потери. Можно определить отражательный угол при помощи функции rangeangle и установки ссылочных осей к глобальной системе координат. Общую длину пути для пути угла обзора показывает в фигуре _Rlos, который равен геометрическому расстоянию между источником и получателем. Общей длиной пути для отраженного пути является _{Rrp= R1 + R2}. Количество L является наземной областью значений между источником и получателем.

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

$\begin{array}{l} \vec{R} = {\vec{x}}_{s} - {\vec{x}}_{r} \\ R_{l o s} = | \vec{R} | = \sqrt{{(z_{r} - z_{s})}^{2} + L^{2}} \\ R_{1} = \frac{z_{r}}{z_{r} + z_{z}} \sqrt{{(z_{r} + z_{s})}^{2} + L^{2}} \\ R_{2} = \frac{z_{s}}{z_{s} + z_{r}} \sqrt{{(z_{r} + z_{s})}^{2} + L^{2}} \\ R_{r p} = R_{1} + R_{2} = \sqrt{{(z_{r} + z_{s})}^{2} + L^{2}} \\ загар θ_{l o s} = \frac{(z_{s} - z_{r})}{L} \\ загар θ_{r p} = - \frac{(z_{s} + z_{r})}{L} \\ {θ^{'}}_{l o s} = - θ_{l o s} \\ {θ^{'}}_{r p} = θ_{r p} \end{array}$

2D излучите затухание

Потеря затухания или пути в канале 2D луча является продуктом пяти компонентов, _{L = Ltworay LG Lg Lc Lr}, где

_Ltworay является 2D лучом геометрическое затухание пути
_LG является наземным затуханием отражения
_Lg является атмосферным затуханием пути
_Lc является вуалью и затуханием пути к облаку
_Lr является затуханием пути к дождю

Каждый компонент находится в единицах значения, не в дБ.

Оснуйте потерю отражения и распространения

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

$E_{l o s} = E_{0} (\frac{R_{0}}{R_{l o s}}) e^{i ω (t - R_{l o s} / c)}$

для пути угла обзора. Можно выразить отраженный о земле E - поле как

$E_{r p} = L_{G} E_{0} (\frac{R_{0}}{R_{r p}}) e^{i ω (t - R_{r p} / c)}$

где _Rrp является отраженным расстоянием пути. Количество _LG представляет потерю из-за отражения в наземной плоскости. Чтобы задать _LG, используйте свойство GroundReflectionCoefficient. В целом _LG зависит от угла установки поля. Если у вас есть эмпирическая информация об угловой зависимости _LG, можно использовать rangeangle, чтобы вычислить угол установки отраженного пути. Общее поле в месте назначения является суммой полей line-of-sight и reflected-path.

Для электромагнитных волн более сложная, но более реалистическая модель использует векторное представление поляризованного поля. Можно разложить инцидентное электрическое поле на два компонента. Один компонент, _Ep, параллелен плоскости падения. Другой компонент, _Es, перпендикулярен плоскости падения. Наземные коэффициенты отражения для этих компонентов отличаются и могут быть записаны с точки зрения наземной проницаемости и угла установки.

$\begin{array}{l} G_{p} = \frac{Z_{1} потому что θ_{1} - Z_{2} потому что θ_{2}}{Z_{1} потому что θ_{1} + Z_{2} потому что θ_{2}} = \frac{потому что θ_{1} - \frac{Z_{2}}{Z_{1}} потому что θ_{2}}{потому что θ_{1} + \frac{Z_{2}}{Z_{1}} потому что θ_{2}} \\ G_{s} = \frac{Z_{2} потому что θ_{1} - Z_{1} потому что θ_{2}}{Z_{2} потому что θ_{1} + Z_{1} потому что θ_{2}} = \frac{потому что θ_{2} - \frac{Z_{2}}{Z_{1}} потому что θ_{1}}{потому что θ_{2} + \frac{Z_{2}}{Z_{1}} потому что θ_{1}} \\ Z_{1} = \sqrt{\frac{μ_{1}}{ε_{1}}} \\ Z_{2} = \sqrt{\frac{μ_{2}}{ε_{2}}} \end{array}$

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

$\begin{array}{l} G_{p} = \frac{\sqrt{ρ} потому что θ_{1} - потому что θ_{2}}{\sqrt{ρ} потому что θ_{1} + потому что θ_{2}} \\ G_{s} = \frac{\sqrt{ρ} потому что θ_{2} - потому что θ_{1}}{\sqrt{ρ} потому что θ_{2} + потому что θ_{1}} \end{array}$

где количество _{ρ = ε2/ε1} является наземной проницаемостью родственника, установленной свойством GroundRelativePermittivity. Угол _θ1 является углом установки и углом _θ2, является углом преломления на контуре. Можно определить _θ2 с помощью закона Поводка преломления.

После отражения все поле восстановлено от параллельных и перпендикулярных компонентов. Общее наземное затухание плоскости, _LG, является комбинацией _Gs и _Gp.

Когда источник и место назначения являются стационарными друг относительно друга, можно записать вывод Y step как Y(t) = F(t-τ)/L. Количество τ является задержкой сигнала и L, является потерей пути свободного пространства. Задержка τ дана R/c. R является или расстоянием пути к распространению угла обзора или отраженным расстоянием пути, и c является скоростью распространения. Потеря пути

$L_{t w o r a y} = \frac{{(4 π R)}^{2}}{λ^{2}},$

где λ является длиной волны сигнала.

Атмосферная газовая модель затухания

Эта модель вычисляет затухание сигналов, которые распространяют через атмосферные газы.

Электромагнитные сигналы затухают, когда они распространяют через атмосферу. Этот эффект должен, в основном, к строкам резонанса поглощения кислорода и водяного пара с меньшими вкладами, поступающими из газа азота. Модель также включает непрерывный спектр поглощения ниже 10 ГГц. Рекомендация модели ITU ITU-R P.676-10: Затухание атмосферными газами используется. Модель вычисляет определенное затухание (затухание на километр) как функция температуры, давления, плотности водяного пара и частоты сигнала. Атмосферная газовая модель допустима для частот от 1-1000 ГГц и применяется к поляризованным и неполяризованным полям.

Формула для определенного затухания на каждой частоте

$γ = γ_{o} (f) + γ_{w} (f) = 0.1820 f N^{″} (f) .$

Количество N"() является мнимой частью комплексного атмосферного явления преломления и состоит из спектрального компонента строки и непрерывного компонента:

$N^{″} (f) = \sum_{i} S_{i} F_{i} + {N^{″}}_{D}^{} (f)$

Спектральный компонент состоит из суммы дискретных условий спектра, состоявших из локализованной функции пропускной способности частоты, F(f) _i, умноженный на спектральную силу строки, S _i. Для атмосферного кислорода каждая спектральная сила строки

$S_{i} = a_{1} \times 10^{- 7} {(\frac{300}{T})}^{3} \exp [a_{2} (1 - (\frac{300}{T})] P .$

Для атмосферного водяного пара каждая спектральная сила строки

$S_{i} = b_{1} \times 10^{- 1} {(\frac{300}{T})}^{3.5} \exp [b_{2} (1 - (\frac{300}{T})] W .$

P является сухим давлением воздуха, W является парциальным давлением водяного пара, и T является температурой окружающей среды. Единицы давления находятся в гектопаскалях (гПа), и температурой является в градусах Келвин. Парциальное давление водяного пара, W, связано с плотностью водяного пара, ρ,

$W = \frac{ρ T}{216.7} .$

Общим атмосферным давлением является P + W.

Для каждой кислородной строки _Si зависит от двух параметров, _a1 и _a2. Точно так же каждая строка водяного пара зависит от двух параметров, _b1 и _b2. Документация ITU, процитированная в конце этого раздела, содержит табулирование этих параметров как функции частоты.

Локализованные функции пропускной способности частоты _Fi(f) являются сложными функциями частоты, описанной в ссылках ITU, процитированных ниже. Функции зависят от эмпирических параметров модели, которые также сведены в таблицу в ссылке.

Чтобы вычислить общее затухание для узкополосных сигналов вдоль пути, функция умножает определенное затухание на длину пути, R. Затем общим затуханием является _{Lg= R(γo + γw)}.

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

Вуаль и модель затухания облака

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

Вуаль и затухание облака являются тем же атмосферным явлением. Модель ITU, Рекомендация ITU-R P.840-6: Затухание из-за облаков и вуали используется. Модель вычисляет определенное затухание (затухание на километр), сигнала как функция жидкой водной плотности, частоты сигнала и температуры. Модель применяется к поляризованным и неполяризованным полям. Формула для определенного затухания на каждой частоте

$γ_{c} = K_{l} (f) M,$

где M является жидкой водной плотностью в ^gm/m3. Количество _Kl(f) является определенным коэффициентом затухания и зависит от частоты. Модель затухания облака и вуали допустима для частот 10-1000 ГГц. Модули для определенного коэффициента затухания (дБ/км) / ^(g/m3).

Чтобы вычислить общее затухание для узкополосных сигналов вдоль пути, функция умножает определенное затухание на длину пути R. Общим затуханием является _{Lc = Rγc}.

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

Модель затухания ливня

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

Электромагнитные сигналы ослаблены при распространении через область ливня. Затухание ливня вычисляется согласно модели Recommendation ITU-R P.838-3 ливня ITU: Определенная модель затухания для дождя для использования в методах прогноза. Модель вычисляет определенное затухание (затухание на километр) сигнала как функция уровня ливня, частоты сигнала, поляризации и угла повышения пути. Вычислить затухание, эта модель использование

$γ_{r} = k r^{α},$

где r является уровнем дождя в мм/час. Параметр k и экспонента α зависит от частоты, вида поляризации и угла повышения пути прохождения сигнала. Определенная модель затухания допустима для частот от 1-1000 ГГц.

Чтобы вычислить общее затухание для узкополосных сигналов вдоль пути, функция умножает определенное затухание на расстояние распространения, R. Затем общим затуханием является _{Lr = Rγr}. Вместо того, чтобы использовать геометрическую область значений в качестве расстояния распространения, тулбокс использует измененную область значений. Измененная область значений является геометрической областью значений, умноженной на фактор области значений

$\frac{1}{1 + \frac{R}{R_{0}}}$

где

$R_{0} = 35 e^{- 0.015 r}$

эффективная длина пути в километрах (см. Seybold, J. Введение в Распространение РФ.), Когда нет никакого дождя, эффективная длина пути составляет 35 км. Когда уровень дождя - например, 10 мм/час, эффективная длина пути составляет 30,1 км. С близкого расстояния расстояние распространения является приблизительно геометрической областью значений. Для более длинных областей значений расстояние распространения асимптотически приближается к эффективной длине пути.

Обработка частоты поддиапазона

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

При использовании широкополосных Системных объектов частоты или блоков, вы задаете количество поддиапазонов, N _B, в котором можно анализировать широкополосный сигнал. Частоты центра поддиапазона и ширины автоматически вычисляются из общей пропускной способности и количества поддиапазонов. Общий диапазон частот сосредоточен на поставщике услуг или рабочей частоте, _fc. Полная пропускная способность дана частотой дискретизации, _fs. Ширинами поддиапазона частоты является Δf = f _s/NB. Центральные частоты поддиапазонов

$f_{m} = {\begin{matrix} f_{c} - \frac{f_{s}}{2} + (m - 1) Δ f, N_{B} даже \\ f_{c} - \frac{(N_{B} - 1) f_{s}}{2 N_{B}} + (m - 1) Δ f, N_{B} нечетный \end{matrix}, m = 1, \dots, N_{B}$

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

Ссылки

[1] Proakis, J. Цифровая связь. Нью-Йорк: McGraw-Hill, 2001.

[2] Skolnik, M. Введение в радиолокационные системы, 3-го Эда. Нью-Йорк: McGraw-Hill.

[3] Саакян, A. Основные принципы распространения радиоволны. Норвуд, MA: дом Artech, 2011.

[4] Balanis, C. Усовершенствованный технический электромагнетизм. Нью-Йорк: Wiley & Sons, 1989.

[5] Rappaport, T. Радиосвязи: принципы и практика, 2-й Эд Нью-Йорк: Prentice Hall, 2002.

[6] Сектор радиосвязи Международного союза электросвязи. Рекомендация ITU-R P.676-10: Затухание атмосферными газами. 2013.

[7] Сектор радиосвязи Международного союза электросвязи. Рекомендация ITU-R P.840-6: Затухание из-за облаков и вуали. 2013.

[8] Сектор радиосвязи Международного союза электросвязи. Рекомендация ITU-R P.838-3: Определенная модель затухания для дождя для использования в методах прогноза. 2005.

Расширенные возможности

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.

Указания и ограничения по применению:

Смотрите системные объекты в Генерации кода MATLAB (MATLAB Coder).

Смотрите также

Функции

fogpl | fspl | gaspl | rainpl | rangeangle

Системные объекты

phased.FreeSpace | phased.LOSChannel | phased.TwoRayChannel | phased.WidebandBackscatterRadarTarget | phased.WidebandFreeSpace | phased.WidebandLOSChannel

Документация

поэтапный. WidebandTwoRayChannel

Описание

Примечание

Конструкция

Свойства

PropagationSpeed — Скорость распространения сигнала physconst('LightSpeed') (значение по умолчанию) | положительная скалярная величина

OperatingFrequency — Рабочая частота 300e6 (значение по умолчанию) | положительная скалярная величина

SpecifyAtmosphere — Включите атмосферную модель затухания false (значение по умолчанию) | true

Temperature — Температура окружающей среды 15 (значение по умолчанию) | скаляр с действительным знаком

Зависимости

DryAirPressure — Атмосферное сухое давление воздуха 101.325e3 (значение по умолчанию) | положительный скаляр с действительным знаком

Зависимости

WaterVapourDensity — Атмосферная плотность водяного пара 7.5 (значение по умолчанию) | положительный скаляр с действительным знаком

Зависимости

LiquidWaterDensity — Жидкая водная плотность 0.0 (значение по умолчанию) | неотрицательный скаляр с действительным знаком

Зависимости

RainRate — Уровень ливня 0.0 (значение по умолчанию) | неотрицательный скаляр

Зависимости

SampleRate — Частота дискретизации сигнала 1e6 (значение по умолчанию) | положительная скалярная величина

NumSubbands — Количество обработки поддиапазонов 64 (значение по умолчанию) | положительное целое число

EnablePolarization — Включите поляризованные поля false (значение по умолчанию) | true

GroundReflectionCoefficient — Оснуйте отражательный коэффициент -1 (значение по умолчанию) | скаляр с комплексным знаком | 1 с комплексным знаком N вектором - строкой

Зависимости

Зависимости

CombinedRaysOutput — Опция, чтобы объединить два луча при выводе true (значение по умолчанию) | false

MaximumDistanceSource — Источник максимального одностороннего расстояния распространения 'Auto' (значение по умолчанию) | 'Property'

MaximumDistance — Максимальное одностороннее расстояние распространения 10000 (значение по умолчанию) | положительный скаляр с действительным знаком

Зависимости

MaximumNumInputSamplesSource — Источник максимального количества выборок 'Auto' (значение по умолчанию) | 'Property'

Зависимости

MaximumNumInputSamples — Максимальное количество выборок входного сигнала 100 (значение по умолчанию) | положительное целое число

Зависимости

Методы

Примеры

Скалярное широкополосное распространение сигнала в канале 2D луча

Сравните широкополосное распространение канала 2D луча со свободным пространством

Широкополосное поляризованное полевое распространение в канале 2D луча

2D излучите распространение широкополосной формы волны LFM с атмосферными потерями

Больше о

2D излучите пути к распространению

2D излучите затухание

Оснуйте потерю отражения и распространения

Атмосферная газовая модель затухания

Вуаль и модель затухания облака

Модель затухания ливня

Обработка частоты поддиапазона

Ссылки

Расширенные возможности

Генерация кода C/C++ Генерация кода C и C++ с помощью MATLAB® Coder™.

Смотрите также

Функции

Системные объекты

Введенный в R2017b

Документация Phased Array System Toolbox

Поддержка

`PropagationSpeed` — Скорость распространения сигнала
`physconst('LightSpeed')` (значение по умолчанию) | положительная скалярная величина

`OperatingFrequency` — Рабочая частота
`300e6` (значение по умолчанию) | положительная скалярная величина

`SpecifyAtmosphere` — Включите атмосферную модель затухания
`false` (значение по умолчанию) | `true`

`Temperature` — Температура окружающей среды
`15` (значение по умолчанию) | скаляр с действительным знаком

`DryAirPressure` — Атмосферное сухое давление воздуха
`101.325e3` (значение по умолчанию) | положительный скаляр с действительным знаком

`WaterVapourDensity` — Атмосферная плотность водяного пара
`7.5` (значение по умолчанию) | положительный скаляр с действительным знаком

`LiquidWaterDensity` — Жидкая водная плотность
`0.0` (значение по умолчанию) | неотрицательный скаляр с действительным знаком

`RainRate` — Уровень ливня
`0.0` (значение по умолчанию) | неотрицательный скаляр

`SampleRate` — Частота дискретизации сигнала
`1e6` (значение по умолчанию) | положительная скалярная величина

`NumSubbands` — Количество обработки поддиапазонов
`64` (значение по умолчанию) | положительное целое число

`EnablePolarization` — Включите поляризованные поля
`false` (значение по умолчанию) | `true`

`GroundReflectionCoefficient` — Оснуйте отражательный коэффициент
`-1` (значение по умолчанию) | скаляр с комплексным знаком | 1 с комплексным знаком N вектором - строкой

`CombinedRaysOutput` — Опция, чтобы объединить два луча при выводе
`true` (значение по умолчанию) | `false`

`MaximumDistanceSource` — Источник максимального одностороннего расстояния распространения
`'Auto'` (значение по умолчанию) | `'Property'`

`MaximumDistance` — Максимальное одностороннее расстояние распространения
`10000` (значение по умолчанию) | положительный скаляр с действительным знаком

`MaximumNumInputSamplesSource` — Источник максимального количества выборок
`'Auto'` (значение по умолчанию) | `'Property'`

`MaximumNumInputSamples` — Максимальное количество выборок входного сигнала
`100` (значение по умолчанию) | положительное целое число

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.