47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
## Compute Gamma
Gamma = [];
for i=1:n_U
U = zeros(n_U,1); U(i,1) = 1; # Set up U_i
y_i = ppp_ystar (A,B,C,D,x_0,A_u,U,Tau);# Compute y* for ith input for each tau
y_i = y_i(:,i_y:n_y:n_y*n_tau); # Pluck out output i_y
Gamma = [Gamma [-y_i';y_i']]; # Put in parts for Min and max
endfor
## Compute gamma
gamma = [-Min';Max'];
endfunction
|
|
|
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
## Compute Gamma
Gamma = [];
for i=1:n_U
U = zeros(n_U,1); U(i,1) = 1; # Set up U_i
y_i = ppp_ystar (A,B,C,D,x_0,A_u,U,Tau);# Compute y* for ith input for each tau
y_i = y_i(i_y,:); # Pluck out output i_y
Gamma = [Gamma [-y_i';y_i']]; # Put in parts for Min and max
endfor
## Compute gamma
gamma = [-Min';Max'];
endfunction
|