SensorFDwithPITL/runSimulation_v1.m at main · fguc/SensorFDwithPITL · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
function [DMDys,DMDr,data] = runSimulation_v1(simData,numSim)

simName = simData.simName;
wnd = simData.DMDWndw;
simInputs = simData.simInputs;

% Simulation parameters
K = 1;      %gain
T = 0.1;    %sec
L = 1;  %delay

TsSystem = 0.01;
RefSignal = 0;

Ks = 0;
Ka = 0;
BL = 0;

da = 0;
ds = numSim;

% simOut{1} = sim(simName,'StopTime',num2str(RefSignal.Time(end)),'SrcWorkspace','current');
simOut{1} = sim(simName,'StopTime',num2str(20),'SrcWorkspace','current');

% Data Extraction
data.r{1} = simOut{1, 1}.logsout{1}.Values.signal1;
data.e{1} = simOut{1, 1}.logsout{1}.Values.signal2;
data.u{1} = simOut{1, 1}.logsout{1}.Values.signal3;
data.ua{1} = simOut{1, 1}.logsout{1}.Values.signal4;
data.y{1} = simOut{1, 1}.logsout{1}.Values.signal5;
data.ys{1} = simOut{1, 1}.logsout{1}.Values.signal6;
data.da{1} = simOut{1, 1}.logsout{1}.Values.signal7;
data.ds{1} = simOut{1, 1}.logsout{1}.Values.signal8;

% Data Prep
data.time = data.ys{1, 1}.Time;
dataDMDraw.ys = data.ys{1, 1}.Data;
dataDMDraw.r = data.r{1, 1}.Data;
dataDMD.ys = zeros(wnd,length(dataDMDraw.ys)-wnd);
dataDMD.r = zeros(wnd,length(dataDMDraw.r)-wnd);
for idx = 1:1:length(dataDMDraw.ys)-wnd
    dataDMD.ys(1:wnd,idx) = dataDMDraw.ys(idx:idx+wnd-1,1);
    dataDMD.r(1:wnd,idx) = dataDMDraw.r(idx:idx+wnd-1,1);
end

DMDys = dataDMD.ys;
DMDr = dataDMD.r;
end