Page 74 - 363_

75

31.0960
0.1284

c =
0.6299 0
d=
-0.0723
Необхідно обчислити значення полюсів та нулів відповідної функції
передачі. Ця задача зводиться до знаходження власних значень , які

задовольняють матричному рівнянню:

R  r     Q ; r 
 R b  Q 0
   r     . r 
 c d   0 0 

Розв'язок першого рівняння дозволяє обчислити полюси функції
передачі, а другого – нулі. Нижче наведено сукупність операторів, яка

приводить до розрахунку полюсів:

>> (АА, BB] = qz(R,-Q) % Приведення матриць до форми Шура
AA =
5.5039 + 2.7975i 24.8121 -25.364i
0.0000 - 0.0000i 5.5158 - 2.8036i
BB =
-0.6457 + 0.7622i -0.1337 + 0.1378i
0 -0.6471 - 0.7638i
>> diag(AA) ./dlag(BB) % Розрахунок полюсів
ans =
-1.4245 - 6.0143i
-1.4245 + 6.0143i

Розрахунок нулів здійснюється таким чином:

>> A = [-R b % Формування
c d] % першої матриці
A =
-1.1190 1.0000 0.1284
-36.4800 -1.5380 31.0960
0.6300 0 -0.0723
>> В = [ -Q zeros(size(b)) % Формування
zeros(size(c)) 0 ] % другої матриці
В =
-1.0000 0 0
-0.1920 -1.0000 0
0 0 0
>> [AA,BB) = qz(A,B) %Приведення матриць
AA = % до форми Шура
31.0963 -0.7169 -36.5109
0.0000 1.0647 0.9229
0 0.0000 0.5119