Newton

% Método de Newton para OGSR   % Se comprueba la derivada (el gradiente y el hessiano) ERROR ABSOLUTO   x=rand(n,1);dejer(x)-jacnum(‘ejer’,x)’,A=jacnum(‘dejer’,x);A=(A+A’)/2;hejer(x)-A   % Opcion del ERROR RELATIVO   x=rand(n,1);gfx=dejer(x);v=jacnum(‘ejer’,x)’;norm(gfx-v)/max(1,norm(gfx)) hfx=hejer(x);A=jacnum(‘dejer’,x)’;A=(A+A’)/2;norm(hfx-A,1)/max(1,norm(hfx,1))       % Se ejecuta la siguiente línea una sola vez em=eps/2;e0=sqrt(em);x=[0 0]’;fx=ejer(x);gfx=dejer(x);A=hejer(x);[R,p]=chol(A);if p==0,r=rcond(A);if r<e0,p=1;end,end,if p>0,V=eig(A);M=max(V);m=min(V);if M==m,A=eye(size(A));else lambda=(e0*M-m)/(1-e0); A=A+lambda*eye(size(A));end,end;d=-A\gfx;[xn,fxn,ro,info]=armijo(‘ejer’,x,d,fx,gfx);if info==0,gfxn=dejer(xn);[fx-fxn,norm(xn-x),norm(gfxn)],else info,end;   % Actualización […]

BFGS

% Método de BFGS para OGSR   % Se comprueba el gradiente ERROR ABSOLUTO   x=rand(n,1);dejer(x)-jacnum(‘ejer’,x)’   % Se comprueba el gradiente ERROR RELATIVO   x=rand(n,1);gfx=dejer(x);v=jacnum(‘ejer’,x)’;norm(gfx-v)/max(1,norm(gfx))   % Se ejecuta la siguiente línea una sola vez   em=eps/2;e0=sqrt(em);x=[0 0]’;fx=ejer(x);gfx=dejer(x);B=eye(length(x));d=-B\gfx;ro=-2*max(1,abs(fx))/(gfx’*d);[xn,ro,fxn,gfxn,info]=wolfe(x,d,fx,gfx,ro,’ejer’,’dejer’); if info==0, s=xn-x;y=gfxn-gfx;[fx-fxn,norm(s),norm(gfxn)], else,info,end; sy=s’*y;bs=B*s;sbs=s’*bs; if sy<0.2*sbs, theta=0.8*sbs/(sbs-sy);y=theta*y+(1-theta)*bs;sy=0.2*sbs; end, alpha=(y’*y)/sy;B=alpha*B;bs=alpha*bs;sbs=alpha*sbs;B=B+(y*y’)/sy-(bs*bs’)/sbs; x=xn;fx=fxn;gfx=gfxn;d=-B\gfx;ro=1;[xn,ro,fxn,gfxn,info]=wolfe(x,d,fx,gfx,ro,’ejer’,’dejer’); if info==0, s=xn-x;y=gfxn-gfx;[fx-fxn,norm(s),norm(gfxn)], else,info,end;   […]

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 […]

Yo #soy15m

Como parte del #15m me declaro una persona pacífica y condeno radicalmente todo tipo de violencia: la de los violentos infiltrados en nuestras manifestaciones, y la del Estado, que ha causado más dolor y heridos. Además, condeno la manipulación mediática que enfatiza la información sesgada, parcial o errónea con el propósito de demonizar a los […]