-
Notifications
You must be signed in to change notification settings - Fork 0
/
par_v_damp.m
139 lines (115 loc) · 20.9 KB
/
par_v_damp.m
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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
% conclusão do teste
% não convém paralelizar este laço para mzk_mobile
global ncontrol tcromossomo opoints alfa
w1= pop(1,:);
ncontrol= 16;
tcromossomo= 3*ncontrol;
opoints= 24;
alfa= 0.1;
run('D:\Sérgio\Pesquisas\PSO revista\alg\op\load_NENY.m')
%define as matrizes A, B e C, de acordo com os Alfas, Betas e Gamas
%Ac(i)= [ -alfa(i) 0 0
% gama(i)-alfaBeta(i) -gama(i) 0
% beta(i)Gama(i)-alfaBeta(i)^2 gama(i)-Beta(i)Gama(i) -gama(i)]
%Bc(i)= [ 1; beta(i); beta(i)^2]
%Cc(i)= [0 0 K]
k=1; %reinicia controladores
j=1; %reinicia j
while j <= tcromossomo
a(1,k).c= [ -alfa 0 0;
(w1(1,j+2)-alfa*w1(1,j+1)) -w1(1,j+2) 0;
(w1(1,j+1)*w1(1,j+2)-alfa*(w1(1,j+1))^2) (w1(1,j+2)*(1-w1(1,j+1))) -w1(1,j+2)];
b(1,k).c= [ 1; w1(1,j+1); w1(1,j+1)*w1(1,j+1)];
cc(1,k).c= [ 0 0 w1(1,j)];
k= k+1;
j= j+3;
end
% define as matrizes Ac, Bc e Cc
% Ac(1,1).c =[ a(1,1).c zeros zeros zeros ...
% zeros a(1,2).c zeros zeros ...
% ... ... ... ... ... ]
Ac(1,1).c= [ a(1,1).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) a(1,2).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,3).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,4).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,5).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,6).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,7).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,8).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,9).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,10).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,11).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,12).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,13).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,14).c zeros(size(a(1,1).c)) zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,15).c zeros(size(a(1,1).c));
zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) zeros(size(a(1,1).c)) a(1,16).c];
Bc(1,1).c= [ b(1,1).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) b(1,2).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,3).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,4).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,5).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,6).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,7).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,8).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,9).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,10).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,11).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,12).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,13).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,14).c zeros(size(b(1,1).c)) zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,15).c zeros(size(b(1,1).c));
zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) zeros(size(b(1,1).c)) b(1,16).c];
Cc(1,1).c= [ cc(1,1).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) cc(1,2).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,3).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,4).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,5).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,6).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,7).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,8).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,9).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,10).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,11).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,12).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,13).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,14).c zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,15).c zeros(size(cc(1,1).c));
zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) zeros(size(cc(1,1).c)) cc(1,16).c];
% define matriz x~' e salva na estrutura xtl.c (x-tiu-linha).c
% par 30x+
for j=1:opoints
xtl(1,j).c= [ A(1,j).po (B(1,j).po*Cc(1,1).c);
(Bc(1,1).c*C(1,j).po*A(1,j).po) (Ac(1,1).c+Bc(1,1).c*C(1,j).po*B(1,j).po*Cc(1,1).c)];
end
clear Ac Bc Cc a b cc
x=0; w=0; z=0; %cria variáveis temporárias
% verifica amortecimentos da matriz x~' e salva na estrutura ax.c
% par um pouco mais rápido (as vezes) -> não compensa
for j=1:opoints
[w,x,z]= damp(xtl(1,j).c);
ax(1,j).c= x;
end
%%% fim da verificação
clear xtl.c;
x=0; w=0; z=0; %zera armazenadores temporários
[size_ax w]= size(ax(1,1).c);
%numera os amortecimentos
%este laço serve apenas para conferir o funcionamento do programa
%na etapa da programação
% par 200x+
for k=1:opoints
for j=1:size_ax
ax(1,k).c(j,2)= j;
end %j
end %k
% verifica o menor amortecimento entre os pontos de operação
% os menores amortecimentos são guardados na matriz amt
w=0;
% par 1000x+
for k=1:opoints
w(k)= min(ax(1,k).c(:,1));
amt(1,k)= w(k);
end %k
md=0;
md= min(amt(1,:));