DCNLProjectHeartbeatModel/SlidingModeControl_script.m at main · AndreaMazzera/DCNLProjectHeartbeatModel · 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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
%% Sliding Mode Control
close all
% Parametri Sliding Mode Control
p1=0;
p2=1;
k=400;

%% Esecuzione

out=sim("SlidingModeControl.slx");

t=out.tout;
e=out.err_and_control.signals(1).values;
u=out.err_and_control.signals(2).values;
ref=out.statevar_and_ref.signals(1).values(:,1);
x1 =out.statevar_and_ref.signals(1).values(:,2);
x2 =out.statevar_and_ref.signals(2).values;
sigma=out.sigma.signals.values;

%% Plot Errore
figure(1)
subplot(2,1,1)
plot(t,e,'g','LineWidth',1.0);
title('Evolution of error e(t)','fontsize',14,'interpreter','latex')
xlabel('t [s]','fontsize',14,'interpreter','latex')
ylabel('e [mV]','fontsize',14,'interpreter','latex')
grid on
%% Plot Controllo
subplot(2,1,2)
plot(t,u,'LineWidth',1.0);
title('Evolution of control u(t)','fontsize',14,'interpreter','latex')
xlabel('t [s]','fontsize',14,'interpreter','latex')
ylabel('u [mV]','fontsize',14,'interpreter','latex')
grid on
%% Plot Variabili di Stato e Riferimento
figure(2)
subplot(2,1,1)
h=plot(t,ref,'--',t,x1,'k');
set(h,{'LineWidth'},{2.5;1.0});
set(h,{'Color'},{"#EDB120";"#000000"});
title('Evolution of reference and real ECG','fontsize',14,'interpreter','latex')
xlabel('t [s]','fontsize',14,'interpreter','latex')
ylabel('ref,$x_1$ [mV]','fontsize',14,'interpreter','latex')
legend('ECG ref','ECG real')
grid on

subplot(2,1,2)
plot(t,x2,'LineWidth',1.0)
title('Evolution of fiber length','fontsize',14,'interpreter','latex')
xlabel('t [s]','fontsize',14,'interpreter','latex')
ylabel('$x_2$ [$\mu$m]','fontsize',14,'interpreter','latex')
grid on

%% Plot Sigma
figure(3)
plot(t,sigma,'r')
title('Evolution of $\sigma$','fontsize',14,'interpreter','latex')
xlabel('t [s]','fontsize',14,'interpreter','latex')
ylabel('$\sigma$','fontsize',14,'interpreter','latex')
grid on;