Критический компонент в поэтапных системных приложениях массивов является способностью смоделировать движение на пробеле. Такое моделирование включает движение массивов, целей и источников интерференции. Для удобства можно проигнорировать различие между этими объектами и коллективно смоделировать движение платформы.
Расширенные тела могут подвергнуться и переводному и вращательному движению на пробеле. Моделирование поддержек программного обеспечения Phased Array System Toolbox™ переводного движения.
Моделирование переводного движения платформы требует спецификации радиус-вектора и вектора скорости. Спецификация радиус-вектора подразумевает систему координат. В Phased Array System Toolbox положение платформы и скорость заданы в Глобальной системе координат. Можно думать о позиции платформы вектора смещения от глобального источника или как координаты точки относительно глобального источника.
Позвольте r0 обозначить, что радиус-вектор во время 0 и v обозначает вектор скорости. Радиус-вектор платформы как функция времени, r(t):
Следующая фигура изображает векторную интерпретацию переводного движения.
Когда платформа представляет элемент датчика или массив, важно знать ориентацию элемента или осей локальной координаты массивов. Например, ориентация осей локальной координаты необходима, чтобы извлечь информацию об угле от инцидентных форм волны. Смотрите Системы Глобальной и Локальной координаты для описания систем глобальной и локальной координаты в программном обеспечении. Наконец, для платформ с переменной скоростью, необходимо смочь обновить вектор скорости в зависимости от времени.
Можно смоделировать положение платформы, скорость и локальную ориентацию осей с объектом phased.Platform.
При начале с простого примера смоделируйте движение платформы более чем десять временных шагов. Чтобы определить временной шаг, примите, что у вас есть импульсный передатчик с импульсной частотой повторения (PRF) 1 килогерца. Соответственно, временной интервал между каждым импульсом является 1 миллисекундой. Установите временной шаг, равный импульсному интервалу повторения.
PRF = 1e3; Tstep = 1/PRF; Nsteps = 10;
Затем, создайте объект платформы определение исходного положения платформы и скорости. Примите, что исходное положение платформы в 100 метрах от источника в (60,80,0). Примите, что скорость составляет приблизительно 30 метров в секунду (m/s) с постоянным вектором скорости, данным (15,25.98,0).
platform = phased.Platform('InitialPosition',[60;80;0],'Velocity',[15;25.98;0]);
Ориентация осей локальной координаты платформы является значением свойства InitialOrientationAxes. Можно просмотреть значение этого свойства путем ввода hplat.InitialOrientationAxes в командной строке MATLAB™. Поскольку свойство InitialOrientationAxes не задано в конструкции Системы phased.Platform object™, свойство присвоено его значение по умолчанию [1 0 0;0 1 0;0 0 1]. Используйте метод step, чтобы моделировать переводное движение платформы.
initialPos = platform.InitialPosition; for k = 1:Nsteps pos = platform(Tstep); end finalPos = pos + platform.Velocity*Tstep; distTravel = norm(finalPos - initialPos)
distTravel = 0.3000
Метод step возвращает текущее положение платформы и затем обновляет положение платформы на основе временного шага и скорости. Эквивалентно, в первый раз, когда вы вызываете метод step, вывод является положением платформы в t = 0.
Вспомните, что платформа перемещается с постоянной скоростью приблизительно 30 м/с. Общее истекшее время составляет 0,01 секунды. Вызов метода step возвращает текущее положение платформы и затем обновляет то положение. Соответственно, вы ожидаете, что конечное положение будет отличаться от исходного положения на 0,30 метра. Подтвердите это различие путем исследования значения distTravel.
Большинство платформ в поэтапных приложениях массивов не перемещается с постоянной скоростью. Если временной интервал, описанный количеством временных шагов, является маленьким относительно скорости платформы, можно часто аппроксимировать скорость как постоянную. Однако существуют ситуации, где необходимо обновлять скорость платформы в зависимости от времени. Можно сделать так с Системой phased.Platform object�, потому что свойство Velocity является настраиваемым.
Этот пример моделирует цель первоначально в покое. Начальный вектор скорости (0,0,0). Примите, что временной шаг является 1 миллисекундой. После 500 миллисекунд платформа начинает перемещаться со скоростью приблизительно 10 м/с. Вектор скорости (7.07,7.07,0). Платформа продолжается в этой скорости еще для 500 миллисекунд.
Tstep = 1e-3; Nsteps = 1/Tstep; platform = phased.Platform('InitialPosition',[100;100;0]); for k = 1:Nsteps/2 [tgtpos,tgtvel] = platform(Tstep); end platform.Velocity = [7.07; 7.07; 0]; for k = Nsteps/2+1:Nsteps [tgtpos,tgtvel] = platform(Tstep); end
Этот пример использует Систему phased.Platform object™, чтобы смоделировать изменение в области значений между стационарным радаром и движущейся целью. Радар расположен в (1000,1000,0) и имеет скорость (0,0,0). Цель имеет исходное положение (5000,8000,0) и перемещается с постоянной скоростью (? 30? 45,0). Импульсная частота повторения (PRF) составляет 1 кГц. Примите, что радар испускает десять импульсов.
PRF = 1e3; Tstep = 1/PRF; radar = phased.Platform('InitialPosition',[1000;1000;0]); target = phased.Platform('InitialPosition',[5000;8000;0],... 'Velocity',[-30;-45;0]);
Вычислите начальный целевой диапазон и угол.
[initRng, initAng] = rangeangle(target.InitialPosition,radar.InitialPosition);
Вычислите относительную радиальную скорость.
v = radialspeed(target.InitialPosition,target.Velocity,...
radar.InitialPosition);Моделируйте целевое движение.
Npulses = 10; for num = 1:Npulses tgtpos = target(Tstep); end
Вычислите последнее целевое положение.
tgtpos = tgtpos + target.Velocity*Tstep;
Вычислите итоговый целевой диапазон и угол.
[finalRng,finalAng] = rangeangle(tgtpos,radar.InitialPosition); deltaRng = finalRng - initRng
deltaRng = -0.5396
Постоянная скорость цели составляет приблизительно 54 м/с. Общее истекшее время составляет 0,01 секунды. Область значений между целью и радаром должна уменьшиться приблизительно на 54 сантиметра. Сравните начальную область значений цели, initRng, к итоговой области значений, finalRng, чтобы подтвердить, что это уменьшение происходит.