Модель Guidance для БПЛА фиксированного крыла
fixedwing объект представляет модель руководства уменьшаемого порядка для беспилотного воздушного транспортного средства (UAV). Модель аппроксимирует поведение системы с обратной связью, состоящей из контроллера автопилота и фиксированного крыла кинематическая модель для 3-D движения.
Для БПЛА мультиротора смотрите multirotor.
model = fixedwing создает модель движения фиксированного крыла с double значения точности для входных параметров, выходных параметров и параметров конфигурации модели руководства.
model = fixedwing(DataType) задает точность типа данных (DataType свойство) для входных параметров, выходных параметров и параметров настроек модели руководства.
control | Управляйте командами для UAV |
derivative | Производная времени состояний UAV |
environment | Экологические входные параметры для UAV |
state | Вектор состояния UAV |
В этом примере показано, как использовать fixedwing модель руководства, чтобы симулировать изменение в состоянии UAV из-за ввода команд.
Создайте модель руководства фиксированного крыла.
model = fixedwing;
Установите воздушную скорость транспортного средства путем изменения структуры от state функция.
s = state(model);
s(4) = 5; % 5 m/sЗадайте команду управления, u, это обеспечивает воздушную скорость и дает крен pi/12.
u = control(model); u.RollAngle = pi/12; u.AirSpeed = 5;
Создайте стандартную среду без ветра.
e = environment(model);
Вычислите производную времени состояния, учитывая текущее состояние, управляйте командой и средой.
sdot = derivative(model,s,u,e);
Симулируйте состояние UAV использование ode45 интегрирование. y поле выводит фиксированное крыло состояния UAV на основе этой симуляции.
simOut = ode45(@(~,x)derivative(model,x,u,e), [0 50], s); size(simOut.y)
ans = 1×2
8 904
Постройте изменение в крене на основе симуляции выход. Крен является 7-й строкой simOut.y вывод .
plot(simOut.y(7,:))
Можно также построить траекторию фиксированного крыла с помощью plotTransforms. Создайте векторы перевода и вращения из симулированного состояния. Downsample (каждый 30-й элемент) и транспонирует simOut элементы, и преобразуют Углы Эйлера в кватернионы. Задайте mesh как fixedwing.stl файл и положительное Z-направление как "down". Отображенное представление показывает UAV, делающий постоянный поворот на основе постоянного крена.
downsample = 1:30:size(simOut.y,2); translations = simOut.y(1:3,downsample)'; % xyz-position rotations = eul2quat([simOut.y(5,downsample)',simOut.y(6,downsample)',simOut.y(7,downsample)']); % ZYX Euler plotTransforms(translations,rotations,... 'MeshFilePath','fixedwing.stl','InertialZDirection',"down") hold on plot3(simOut.y(1,:),-simOut.y(2,:),simOut.y(3,:),'--b') % full path xlim([-10.0 10.0]) ylim([-20.0 5.0]) zlim([-0.5 4.00]) view([-45 90]) hold off
Для БПЛА фиксированного крыла следующие уравнения используются, чтобы задать модель руководства UAV. Используйте derivative функция, чтобы вычислить производную времени состояния UAV, использующего эти управляющие уравнения. Задайте входные параметры с помощью state, control, и environment функции.
Положением UAV в наземной системе координат является [xe, ye, h] с ориентацией как угол рыскания, угол угла тангажа и крен, [χ, γ, ϕ] в радианах.
В модели принимается, что UAV летит при условии скоординированного поворота с нулевым заносом. Автопилот управляет скоростью полета, высотой и креном. Соответствующие уравнения движения:
Va и Vg обозначают воздух UAV и скорости относительно земли.
Скорость ветра задана как [Vwn, Vwe, Vwd] для севера, востока, и вниз направлений. Чтобы сгенерировать структуру для этих входных параметров, используйте environment функция.
k* является усилениями контроллера. Чтобы задать эти усиления, используйте Configuration свойство fixedwing объект.
От этих управляющих уравнений модель дает следующие переменные: 
Эти переменные совпадают с выходом state функция.
[1] Рэндал В. Бирд и Тимоти В. Маклэйн. "Глава 9". Маленькая беспилотная теория самолета и практика, NJ: Издательство Принстонского университета, 2012.
ode45 | control | derivative | environment | state | plotTransforms