NLCSmoothReg: Lplot(lambdavec,rho,eta) - 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.

Lplot(lambdavec,rho,eta)

function Lplot(lambdavec,rho,eta)

% make discrepancy plot
dp=figure;
set(dp, 'WindowButtonDownFcn', 'getposDP')
loglog(lambdavec,rho);hold on;loglog(lambdavec(1),rho(1),'LineStyle','none');
xlabel('\lambda'),ylabel('||S - K*g||');title('Discrepancy plot');

% make L-plot
lp=figure;
set(lp, 'WindowButtonDownFcn', 'getposLP')
loglog(rho,eta);hold on;loglog(rho(1),eta(1),'LineStyle','none');hold on;loglog(rho,lambdavec,'LineStyle','none')
set(gca,'XLim',[(0.99*min(rho)) max(rho)],'YLim',[min(eta) max(eta)]);xlabel('||S - K*g||'),ylabel('R[g]^0.5');title('L-curve')

Highlights from NLCSmoothReg

Highlights from
NLCSmoothReg