BEISPIEL_2_4.M

Smith-McMillan-Form mit Matlab

Beispiel_2_4.m — Objective-C source code, 679 bytes

Dateiinhalt

%
% Beispiel 2.4 in Skript "Moderne Methoden der Regelungstechnik"
%
s = sym('s');
z11 = (s-1);
z12 = 4*(s+1);  
z21 = 4.5*(s+1);
z22 = 2*(s-1);
H_s = (s+1)*(s+2);
Z=[ z11 z12; z21 z22]
G=1/H_s*Z
%
%Bestimmung der Smith'schen Normalform mit der Maple-Funktion "Smith"
S=maple('smith',Z,s); 
S=factor(S)
M=1/H_s*S;
% simple(M);
M
%
k11=sym('k11');
k12=sym('k12');
k21=sym('k21');
k22=sym('k22');
K=[k11 k12; k21 k22]
%
% �bertragungsmatrix des geschlossenen Kreises zeigt, 
% da� die Nullstellen der Elemente der �bertragungsmatrix durch
% die R�ckf�hrung U(s) = -K Y(s) + Z(s) ver�ndert werden.
%
G_g=inv(eye(2)+G*K)*G;
G_g=simple(G_g);
pretty(G_g)