Leader (and sub Leader) Election per uniformare e ... - Automatica

E.4
Codice MATLAB CAP 7.
E.4.1 Implementazione del metodo di controllo
ValoreFinaleMinimo.
function soluzione = Msistema(Mg,tempo,alfa)
dim = size (Mg);
% Trovo la matrice da risolvere
S = [];
for i = 1:dim(1,1)
rad = zeros(dim(1,1),1);
rad(i,1) = 1;
O = MaFigliEquispaziati(Mg,rad);
O = O * ones(dim(1,1),1);
S = [S O];
end
Q = inv(S) * ones(dim(1,1),1) - quadprog (2*S’*S ,
zeros(dim(1,1),1) , eye(dim(1,1)) , (inv(S) * ones(dim(1,1),1)));
% Normalizzo la soluzione in funzione del tempo di simulazione
lambda = (alfa+1) * (tempo+1) + alfa * (tempo+1) / (ones(1,dim(1,1)) * Q);
Q = Q * (lambda - (1+alfa) * (tempo+1)) / alfa;
% Approssimo la soluzione a valori naturali
vet = zeros(dim(1,1),1);
for i = 1: dim(1,1)
vet(i,1) = fix (Q(i,1));
end
n = tempo - vet’ * ones(dim(1,1),1); % Quanti valori perdo con la
% approssimazione. Devo rimetterli
for i = 1: n
vet (i,1) = vet (i,1) + 1;
end
soluzione = vet;
63