NLCSmoothReg: getposLP.m - File Exchange

Code covered by the BSD License

Highlights from
NLCSmoothReg

Lplot(lambdavec,rho,eta) make discrepancy plot
Lt=trans2d(dimr,dimc);
[Lgfit,rho,eta]=LCurve(s,K,la... determines a set of solutions $g$ according to a user-defined set of regularization parameters.
[f,d]=chi2spectrum(gcol,guess... calculate deviation vector in spectral domain
[gfit,rho,eta]=NLCSmoothReg(s... determines the solution $g$ according to the user-defined regularization parameter, which has to be estimated
ddL2=L2_2d(n,dimr,dimc); calculates the nth derivative matrix with respect to each columnvector of a
get_l(n,d) GET_L Compute discrete derivative operators.
plotODF(ODF,nocontours) generates a contourplot of a two dimensionally function with
getposDP.m
getposLP.m
View all files

from NLCSmoothReg by Michael Wendlandt
Computes a smooth regularized solution of an ill-posed discrete linear inverse problem.

getposLP.m

pt=get(gca,'Currentpoint');
xdat_pt=pt(1,1);
h=findobj(gcf,'Type','line');

xdat=get(h(1),'XData');
ydat0=get(h(3),'YData');
ydat=get(h(1),'YData');

ydat_pt=interp1(xdat,ydat,xdat_pt);
ydat0_pt=interp1(xdat,ydat0,xdat_pt);

cla
loglog(xdat,ydat0);hold on;loglog(xdat_pt,ydat0_pt,'x');hold on;loglog(xdat,ydat,'LineStyle','none');
set(gca,'XLim',[(0.99*min(xdat)) max(xdat)]);
xlabel('||S - K*g||'),ylabel('R[g]^0.5');title(['L-curve, \lambda = ',num2str(ydat_pt)])

Highlights from NLCSmoothReg

Highlights from
NLCSmoothReg