step
Системный объект: twoRayChannel
Распространите сигнал от точки до точки с помощью модели канала 2D луча
Синтаксис
prop_sig = step(channel,sig,origin_pos,dest_pos,origin_vel,dest_vel)
Описание
Примечание
Запуск в R2016b, вместо того, чтобы использовать step метод, чтобы выполнить операцию, заданную Системой object™, можно вызвать объект с аргументами, как будто это была функция. Например, y = step(obj,x) и y = obj(x) выполните эквивалентные операции.
prop_sig = step(channel,sig,origin_pos,dest_pos,origin_vel,dest_vel)prop_sig, когда узкополосный сигнал, sig, распространяет через канал 2D луча от origin_pos положение к dest_pos положение. Любой origin_pos или dest_pos аргументы могут иметь несколько точек, но вы не можете задать обоих как имеющий несколько точек. Скорость источника сигнала задана в origin_vel и скорость места назначения сигнала задана в dest_vel. Размерности origin_vel и dest_vel должен согласиться с размерностями origin_pos и dest_pos, соответственно.
Электромагнитные поля, распространенные через канал 2D луча, могут быть поляризованы или не поляризованы. Поскольку, неполяризованные поля, такие как акустическое поле, распространение сигнализирует о поле, sig, вектор или матрица. Когда поля поляризованы, sig массив структур. Каждый элемент структуры представляет вектор электрического поля в Декартовой форме.
В среде 2D луча существует два пути прохождения сигнала, соединяющие каждый источник сигнала и целевую пару. Для источников сигнала N (или N сигнализируют местам назначения), существует количество 2N путей. Сигналы для каждой целевой источником пары не должны быть связаны. Сигналы вдоль этих двух путей для любой одной целевой источником пары могут также отличаться из-за фазы или амплитудных различий.
Можно сохранить два сигнала в месте назначения отдельными или объединенными — управляемый CombinedRaysOutput свойство. Объединенный означает, что сигналы в источнике распространяют отдельно вдоль этих двух путей, но когерентно суммированы в месте назначения в одно количество. Чтобы использовать отдельную опцию, установите CombinedRaysOutput к false. Чтобы использовать объединенную опцию, установите CombinedRaysOutput к true. Эта опция удобна, когда различие между датчиком или усилениями массивов в направлениях этих двух путей не является значительным и не должно быть учтено.
Примечание
Объект выполняет инициализацию в первый раз, когда объект выполняется. Эта инициализация блокирует ненастраиваемые свойства и входные технические требования, такие как размерности, сложность и тип данных входных данных. Если вы изменяете ненастраиваемое свойство или входную спецификацию, Системный объект выдает ошибку. Чтобы изменить ненастраиваемые свойства или входные параметры, необходимо сначала вызвать release метод, чтобы разблокировать объект.
Входные параметры
Выходные аргументы
Примеры
Сравните 2D луч с распространением свободного пространства
Распространите сигнал в среде канала 2D луча от радара в (0,0,10) метры к цели в (300,200,30) метры. Примите, что радар и цель являются стационарными и что передающая антенна имеет шаблон косинуса. Сравните объединенные сигналы от этих двух путей с одним сигналом, следующим из распространения свободного пространства. Установите CombinedRaysOutput к true произвести объединенный распространенный сигнал.
Создайте прямоугольный радиоимпульс
Установите частоту дискретизации на 2 МГц.
fs = 2e6;
waveform = phased.RectangularWaveform('SampleRate',fs);
wavfrm = waveform();Создайте передающую антенну и излучателя
Настройте phased.Radiator Система object™, чтобы передать от антенны косинуса
antenna = phased.CosineAntennaElement;
radiator = phased.Radiator('Sensor',antenna);Задайте передатчик и предназначайтесь для координат
posTx = [0;0;10]; posTgt = [300;200;30]; velTx = [0;0;0]; velTgt = [0;0;0];
Распространение свободного пространства
Вычислите направление передачи к цели для модели свободного пространства. Затем излучите сигнал.
[~,angFS] = rangeangle(posTgt,posTx); wavTx = radiator(wavfrm,angFS);
Распространите сигнал к цели.
fschannel = phased.FreeSpace('SampleRate',waveform.SampleRate);
yfs = fschannel(wavTx,posTx,posTgt,velTx,velTgt);
release(radiator);2D излучите распространение
Вычислите два угла передачи к цели для пути угла обзора (LOS) и отраженных путей. Вычислите направления передачи к цели для двух лучей. Затем излучите сигналы.
[~,angTwoRay] = rangeangle(posTgt,posTx,'two-ray');
wavTwoRay = radiator(wavfrm,angTwoRay);Распространите сигналы к цели.
channel = twoRayChannel('SampleRate',waveform.SampleRate,... 'CombinedRaysOutput',true); y2ray = channel(wavTwoRay,posTx,posTgt,velTx,velTgt);
Постройте распространенные сигналы
Постройте объединенный сигнал против сигнала свободного пространства
plot(abs([y2ray yfs])) legend('Two-ray','Free space') xlabel('Samples') ylabel('Signal Magnitude')
Поляризованное полевое распространение в канале 2D луча
Создайте поляризованное электромагнитное поле, состоящее из линейных импульсов формы волны FM. Распространите поле из стационарного источника с антенным элементом пересеченного диполя к стационарному приемнику на расстоянии приблизительно в 10 км. Передающая антенна составляет 100 метров над землей. Приемная антенна составляет 150 м над землей. Приемная антенна является также пересеченным диполем. Постройте полученный сигнал.
Установите радарные параметры формы волны
Примите, что ширина импульса и частота дискретизации составляет 10 МГц. Полоса пропускания импульса составляет 1 МГц. Примите 50%-й рабочий цикл, в котором ширина импульса является половиной импульсного интервала повторения. Создайте 2D импульсную волну, обучаются. Примите несущую частоту 100 МГц.
c = physconst('LightSpeed');
fs = 10e6;
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','Up','Envelope','Rectangular','SweepInterval',... 'Positive'); antenna = phased.CrossedDipoleAntennaElement(... 'FrequencyRange',[50,200]*1e6); radiator = phased.Radiator('Sensor',antenna,'OperatingFrequency',fc,... 'Polarization','Combined'); channel = twoRayChannel('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 Система object™ использует эти углы, чтобы применить отдельные усиления антенны к двум сигналам.
[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 Система object™ использует эти углы, чтобы применить отдельные усиления антенны к двум сигналам.
[~,angsRx] = rangeangle(posTx,posRx,laxRx,'two-ray');Соберите и объедините полученные лучи.
y = collector(prop_sig,angsRx,laxRx);
Постройте принятую форму волны
plot((0:(n-1))/fs*1000000,real(y)) xlabel('Time ({\mu}sec)') ylabel('Received Waveform')
2D излучите распространение формы волны LFM
Распространите линейный сигнал FM в канале 2D луча. Сигнал распространяет от передатчика, расположенного в (1000,10,10) метры в глобальной системе координат к приемнику в (10000,200,30) метры. Примите, что передатчик и приемник являются стационарными и что у них обоих есть шаблоны антенны косинуса. Постройте полученный сигнал.
Настройте радарный сценарий. Во-первых, создайте необходимые Системные объекты.
waveform = phased.LinearFMWaveform('SampleRate',1000000,... 'OutputFormat','Pulses','NumPulses',2); fs = waveform.SampleRate; antenna = phased.CosineAntennaElement; radiator = phased.Radiator('Sensor',antenna); collector = phased.Collector('Sensor',antenna); channel = twoRayChannel('SampleRate',fs,... 'CombinedRaysOutput',false,'GroundReflectionCoefficient',0.95);
Настройте геометрию сцены. Задайте положения передатчика и приемника и скорости. Передатчик и приемник являются стационарными.
posTx = [1000;10;10]; 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/fs; n = size(yR,1); plot((0:(n-1))*dt*1000000,real(yR)) xlabel('Time ({\mu}sec)') ylabel('Signal Magnitude')
2D излучите распространение формы волны LFM с атмосферными потерями
Распространите линейный сигнал FM на 100 МГц в канал 2D луча. Примите, что существует потеря сигнала, вызванная атмосферными газами и дождем. Сигнал распространяет от передатчика, расположенного в (0,0,0) метры в глобальной системе координат к приемнику в (10000,200,30) метры. Примите, что передатчик и приемник являются стационарными и что у них обоих есть шаблоны антенны косинуса. Постройте полученный сигнал. Установите сухое давление воздуха на 102,5 Па и уровень дождя к 5 мм/час.
Настройте радарный сценарий
waveform = phased.LinearFMWaveform('SampleRate',1e6,... 'OutputFormat','Pulses','NumPulses',2); antenna = phased.CosineAntennaElement; radiator = phased.Radiator('Sensor',antenna); collector = phased.Collector('Sensor',antenna); channel = twoRayChannel('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-direction и точки приемной антенны близко к - 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*1000000,real(yR)) xlabel('Time ({\mu}sec)') ylabel('Signal Magnitude')
Ссылки
[1] Proakis, J. Цифровая связь. Нью-Йорк: McGraw-Hill, 2001.
[2] Skolnik, M. Введение в радиолокационные системы, 3-го Эда. Нью-Йорк: McGraw-Hill
[3] Саакян, основные принципы распространения волны A.Radio. Норвуд, MA: дом Artech, 2011.
[4] Balanis, электромагнетизм разработки C.Advanced. Нью-Йорк: Wiley & Sons, 1989.
[5] Rappaport, коммуникации T.Wireless: принципы и практика, 2-й Эд Нью-Йорк: Prentice Hall, 2002.