6 votos

Convexidad de matriz rastro funciones ATr(Bf(A))ATr(Bf(A))

Cualquier función de f:RR se puede extender a una función de auto-adjoint n×n de las matrices f(A)=ni=1f(ai)vivi donde A=ni=1aivivi es la descomposición espectral de A.

Si f:RR es convexa, entonces la función de ATr(f(A)) también es convexa como una función de las matrices. (Ver, por ejemplo, el Teorema 3.27 en la Introducción al Análisis Matricial y Aplicaciones, o ver la prueba a continuación).

Mi pregunta: Si f:RR es convexa, es la función \etiqueta1UnTr(Bf(A)) convexo para cualquier auto-adjunto positivo de la matriz B0?

Me pregunto si las siguientes pruebas pueden ser adaptados para mostrar la convexidad de (1). Aquí os muestro que ATrf(A) es una función convexa mientras f es convexa. No puedo sacar una otra prueba similar a la obra para mostrar la convexidad de ATr(Bf(A)), ni puedo encontrar un contraejemplo.

Lema. Supongamos f:RR es convexo y deje A ser un auto-adjoint n×n matriz. Para cualquier base ortonormales {u1,,un}Cn, sostiene que Trf(A)nj=1f(ujAuj). Prueba. Deje A=ni=1aivivi ser la descomposición espectral de A. A continuación, nj=1|ujvi|2=1 por cada j y

Trf(A)=nj=1ni=1f(ai)|ujvi|2j=1f(ni=1ai|ujvi|2)=j=1f(ujAuj) donde la desigualdad es debido a la convexidad de f.

Teorema. Deje f:RR ser convexo. A continuación, ATrf(A) es convexa como una función de las matrices.

Prueba. Deje A B ser uno mismo-adjoint n×n matrices y deje t(0,1). Vamos a mostrar que tTrf(A)+(1t)Trf(B)Tr(f(tA+(1t)B)). Let {u1,,un} be the eigenbasis of tA+(1t)B. Por el Lema, tTrf(A)+(1t)Trf(B)tnj=1f(ujAuj)+(1t)nj=1f(ujBuj)=nj=1(tf(ujAuj)+(1t)f(ujBuj))nj=1f(tujAuj+(1t)ujBuj)=nj=1f(uj(tA+(1t))Buj)=Tr(f(tA+(1t)B)), como deseado, donde la segunda desigualdad es debido a la convexidad de f.

5voto

Faramir66103 Puntos 6

Después de leer user1551 la respuesta, tengo posteriormente se dio cuenta de que mi pregunta puede ser contestada en sentido negativo de la siguiente manera.

Para n×n matrices, podemos escribir la A0 si A es positivo semidefinite, y AB si AB0. Decimos que una función de f:RR es la matriz convexo si, para todos los enteros n, sostiene que f(tA+(1t)B)(tf)(A)+(1t)f(B) para todos los n×n hermitian matrices A,B y todos los t(0,1), donde la desigualdad es el respeto a la ordenación de positivo semidefinite matrices.

Teorema Deje f:RR ser una función. Los siguientes son equivalentes:

  1. f matrix es convexo
  2. Para todos los enteros n y todos los n×n matrices B0, la función de ATr(Bf(A)) es convexa.

Prueba. Si f matrix es convexo, entonces el resultado claramente de la siguiente manera. Así que supongo que f no es la matriz convexo. Entonces existe un n n×n hermitian matrices A B t(0,1) tal que f(tA+(1t)B)\no(tf)(A)+(1t)f(B). Por lo tanto no es un vector vCn tal que v(f(tA+(1t)B))v>tvf(A)v+(1t)vf(B)v que es equivalente a Tr(vvf(tA+(1t)B))>tTr(vvf(A))+(1t)Tr(vvf(B)), por lo tanto la función de ATr(vvf(A)) no es convexa.

2voto

Chris Ballance Puntos 17329

La respuesta es no. Ya que cada positivos semidefinite matriz es un nonnegatively suma ponderada de las proyecciones ortogonales, podemos suponer que la B es una proyección ortogonal. Con esta simplificación, no es difícil ver que la desigualdad en cuestión deben mantener si f es el operador convexo.

Así que, para encontrar un contraejemplo a la desigualdad, se debe buscar una función de f que es convexa, pero no operador convexo. Tome f(x)=|x|. Deje n=2B=diag(1,0), por lo que el tr(Bf(A)) es sólo el (1,1)-ésima de a f(A). Utilizamos los siguientes Octave/Matlab script para ver si la primera entrada de f(A) es el punto medio convexo:

[U,S,V]=svd(rand(2,2)); % this generates two orthogonal matrices U,V
D1=diag(2*rand(2,1)-1); % random diagonal matrix
D2=diag(2*rand(2,1)-1); % random diagonal matrix
A1=U*D1*U';      % random self-adjoint matrix
A2=V*D2*V';      % random self-adjoint matrix
F1=U*abs(D1)*U'; % f(A1)
F2=V*abs(D2)*V'; % f(A2)
A3=(A1+A2)./2;   % mean of A1 and A2
[Q,D]=eig(A3);   % A3 = Q*D*Q' is a spectral decomposition
F3=Q*abs(D)*Q';  % f(A3)
% If the (1,1)-th entry of f(A) is convex, we always have d>=0:
d = (F1(1,1)+F2(1,1))/2 - F3(1,1)

La respuesta resulta ser "no":

d = -0.013076

A1 =
   0.24718   0.63317
   0.63317   0.25794

A2 =
  -0.18503   0.21871
   0.21871  -0.61285

A3 =
   0.031074   0.425940
   0.425940  -0.177452

F1 =
   0.63104   0.25255
   0.25255   0.63534

F2 =
   0.18503  -0.21871
  -0.21871   0.61285

F3 =
   0.421113  -0.071090
  -0.071090   0.455917

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X