Filter
Top Products
296 Chapter 6. Function Reference
# that mu is not equal to its upper bound for
# more than three full blocks.
gamma = 3 + sqrt(3); beta = sqrt(3) -1
a = sqrt(2/gamma); b = 1/sqrt(gamma)
c = 1/sqrt(gamma); d = -sqrt(beta/gamma)
f = (1+jay)*sqrt(1/(gamma*beta))
psi1 = -pi/2; psi2 = pi
U = [a,0; b,b; c,jay*c; d,f]
V = [0,a; b,-b; c,-jay*c; f*exp(jay*psi1), d*exp(jay*psi2)]
scl = diagonal(random(4,1)+0.1*ones(4,1))
M = scl*U*V*’*inv(scl)
# Consider 4 1x1 blocks
blk1 = [1,1; 1,1; 1,1; 1,1]
[mubnds1,D1,Dinv1,Delta1] = mu(M,blk1)
max(svd(M))?
ans (a scalar) = 2.81264
max(svd(D1*M*Dinv1))?
ans (a scalar) = 1
mubnds1?
mubnds1 (a column vector) =
1
0.860682
# Consider 1 4x4 block (equivalent to max sing. val.)
blk2 = [4,4]
[mubnds2,D2,Dinv2,Delta2] = mu(M,blk2)
# Note that the perturbation is such that
# det(I-M Delta) = 0.
max(svd(D2*M*Dinv2))?