This is machine translation

Translated by Microsoft
Mouseover text to see original. Click the button below to return to the English verison of the page.

Note: This page has been translated by MathWorks. Please click here
To view all translated materals including this page, select Japan from the country navigator on the bottom of this page.


State-space or transfer function plant augmentation for use in weighted mixed-sensitivity H and H2 loopshaping design


P = AUGW(G,W1,W2,W3)


P = AUGW(G,W1,W2,W3) computes a state-space model of an augmented LTI plant P(s) with weighting functions W1(s), W2(s), and W3(s) penalizing the error signal, control signal and output signal respectively (see block diagram) so that the closed-loop transfer function matrix is the weighted mixed sensitivity


where S, R and T are given by


The LTI systems S and T are called the sensitivity and complementary sensitivity, respectively.

Plant Augmentation

For dimensional compatibility, each of the three weights W1, W2 and W3 must be either empty, a scalar (SISO) or have respective input dimensions Ny, Nu, and Ny where G is Ny-by-Nu. If one of the weights is not needed, you may simply assign an empty matrix [ ]; e.g., P = AUGW(G,W1,[],W3) is P(s) as in the Algorithms section below, but without the second row (without the row containing W2).


collapse all

s = zpk('s'); 
G = (s-1)/(s+1);
W1 = 0.1*(s+100)/(100*s+1); 
W2 = 0.1; 
W3 = [];
P = augw(G,W1,W2,W3);

[K,CL,GAM] = hinfsyn(P); 
[K2,CL2,GAM2] = h2syn(P);

L = G*K; 
S = inv(1+L); 
T = 1-S; 

legend('S = 1/(1+L)','GAM/W1','T=L/(1+L)','GAM*G/W2','Location','Northwest')


The transfer functions G, W1, W2 and W3 must be proper, i.e., bounded as s or, in the discrete-time case, as z. Additionally, W1, W2 and W3 should be stable. The plant G should be stabilizable and detectable; else, P will not be stabilizable by any K.


The augmented plant P(s) produced by is


Partitioning is embedded via P=mktito(P,NY,NU), which sets the InputGroup and OutputGroup properties of P as follows

P.InputGroup  = struct('U1',1:c-NU,'U2',c-NU+1:c);
P.OutputGroup = struct('Y1',1:r-NY,'Y2',r-NY+1:r);

See Also

| | |

Introduced before R2006a

Was this topic helpful?