OCR

% ALGORITMO DE OCSR

 

% EN PRIMER LUGAR SE OBTIENEN LOS VALORES PROPIOS DE LA MATRIZ:

eig(H)

 

 

% CASO 1: H>0

x=-H\p

% es un mínimo global estricto, y por tanto único

 

 

% CASO 2: H>=0

% En este caso existen infinitos mínimos globales o bien no existe solución

U=null(H),U’*p

 

% Caso 2.1: U’*p es nulo –> Calculamos uno de entre los infinitos mínimos globales

x-H\p

 

% Caso 2.2 U’*p no es nulos –> No existe solución

 

% Caso 3  H tiene un valor propio negativo –> No existe solución

 

 

% SOLUCIÓN (cuando exista)

 

x,fx=0.5*x’*H*x+p’*x+q

 

 

 

 

 

 

% nulo viejo no hacer

b=-p;[Q,R,E]=qr(H);r=rank(R);[fil,col]=size(H);Z=Q(:,r+1:fil=;Z’*b

Y=Q(:,1:r);Rs=R(1:r,1:r);x=(Rs\(Y’*b);x=[x;zeros(col-r)];x=E*x

Anuncios

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s