Commit 7f911b98 authored by beamdyn's avatar beamdyn
Browse files

more figures in qemsvd_mod

parent b958b1c4
function dq=qemsvd_mod(a,b,neig,plot)
function dq=qemsvd_mod(a,b,neig,plot_flag)
%function dq=qemsvd_mod(a,b,neig,plot)
% given response matrix a and vector b to be corrected the function finds
% the vector dq so that b-a*dq=0
......@@ -22,10 +22,10 @@ dq=v(:,1:neig)*eigscor;
% plot correction effect to set appropriate number of eigenvectors.
if nargin<4
plot=0;
plot_flag=0;
end
if plot
if plot_flag
numeigen0=neig;
neig=sort([1:5:nmax,numeigen0]);
dqstd=zeros(size(neig));
......@@ -44,10 +44,29 @@ if plot
set(get(AX(1),'Ylabel'),'String','rms vector to correct')
set(get(AX(2),'Ylabel'),'String','rms correctors')
%set(AX(1),'YScale','log');
set(AX(2),'YScale','log');
% set(AX(2),'YScale','log');
title('Correctors strength and residual vs number of eigenvectors')
saveas(gca,[datestr(now,'yyyymmddTHHMMSS') '_eigplot.fig']);
export_fig([datestr(now,'yyyymmddTHHMMSS') '_eigplot.jpg']);
close(f);
% close(f);
figure;
plot(lambda)
ax=gca;
set(ax,'YScale','log');
figure;
surf(v);
shading flat;view(2); colormap('jet');
xlabel('eigenvector #');
ylabel('corrector #');
figure;
surf(abs(fft(v(1:end,:))));
shading flat;view(2); colormap('jet');
xlabel('eigenvector #');
ylabel('corrector #');
end
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment