MPC control of aircraft wingtip acceleration and control surface

1、 Content abstract

A little
372- Can communicate 、 consulting 、 Answering question

2、 Content description

clc
close all
clear
% Open file 
path = pwd;
pathFull = strcat(path,'\ System state space ');
addpath(pathFull);
fid = fopen('180.DAT'); 
disp(['Read file ' '180.DAT']) 
% Read and display the header 
head=fgetl(fid); disp(head) 
V=str2num(head(36:46)); ro=str2num(head(56:66)); 
% Read dimensions 
tmp=fgetl(fid); tmp=sscanf(tmp,'%i'); 
if length(tmp)==1 
n=tmp(1); m=0; l=0; 
elseif length(tmp)==3 
n=tmp(1); m=tmp(2); l=tmp(3); 
else 
disp(['FATAL ERROR: File ' filename ' corupt.']); return 
end 
% Read state-space matrices 
ABCD = zeros(n+l,n+m); 
ABCD = fscanf(fid,'%g',[n+l,n+m]);
A=ABCD(1:n,1:n); 
B=ABCD(1:n,n+1:n+m); 
C=ABCD(n+1:n+l,1:n); 
D=ABCD(n+1:n+l,n+1:n+m); 
% Read gain matrix 
if head(74:80)=='VEHICLE' 
fgetl(fid);line = fgetl(fid); 
G = fscanf(fid,'%g',[m,l]); 
else 
G = []; 
end 
% Read gust state-space matrices 
% 
fgetl(fid);line = fgetl(fid); 
if ~isstr(line) 
Bw=[]; CG=[]; return 
else 
nG2= fscanf(fid,'%g',1); 
Bw = fscanf(fid,'%g',[n,nG2]); 
end 
fgetl(fid);line = fgetl(fid); 
if ~isstr(line), CG=[]; return, 
else CG = fscanf(fid,'%g',[l,nG2]); 
end 

% %% establish MPC Model
% sys = ss(A,B,C,D);
% Ts=0.05;     % Sampling time
% p=15;       % Predicted time domain length
% m=2;        % Control the time domain length
% MPC1=mpc(sys,Ts,p,m);    % Definition M1 and M2 The separation of MPC Model
% % Set limits ,F The size of cannot be 0,F The absolute value of slope is the maximum 1000
% MPC1.MV=struct('Min',-1,'Max',1,'RateMin',-1e2,'RateMax',1e2);
% MPC1.OV=struct('Min',-1000,'Max',1000);
% MPC1.OV.MaxECR=1;  
%% establish MPC Model
close all
clc
% sys = ss(A,B,C,D);
% Ts=0.05;     % Sampling time
% p=15;       % Predicted time domain length
% m=5;        % Control the time domain length
% MPC1=mpc(sys,Ts,p,m);    % Definition M1 and M2 The separation of MPC Model
% % Set limits ,F The size of cannot be 0,F The absolute value of slope is the maximum 1000
% MPC1.MV=struct('Min',-1,'Max',1,'RateMin',-1e1,'RateMax',1e1);
% MPC1.OV=struct('Min',-1000,'Max',1000);
% MPC1.OV.MaxECR=1;  
sys = ss(A,B,C,D);
Ts=0.05;     % Sampling time
p=15;       % Predicted time domain length
m=5;        % Control the time domain length
MPC1=mpc(sys,Ts,p,m);    % Definition M1 and M2 The separation of MPC Model
% Set limits ,F The size of cannot be 0,F The absolute value of slope is the maximum 1000
MPC1.MV=struct('Min',-1,'Max',1,'RateMin',-1e1,'RateMax',1e1);
MPC1.OV=struct('Min',-10,'Max',10);
MPC1.OV.MaxECR=1;  
in = 10;
sim('model_mpc1.slx')
figure
plot(tout,yout(:,1))
xlabel Time /s
ylabel Wingtip acceleration sensor

figure
plot(tout,yout(:,2),tout,yout(:,3))
xlabel Time /s
ylabel Control the deflection angle of the face

3、 Simulation analysis

 

4、 Reference paper

A little