0001 function plotdata=plBeamSize(lindata,r,dpp,varargin)
0002
0003
0004
0005
0006 [r2,radindex,~]=atradon(r);
0007 [envelope,~,~]=ohmienvelope(r2,radindex,1:length(r2)+1);
0008
0009 s=real(cat(2,envelope.Sigma))*1e3;
0010
0011
0012
0013
0014 beta=cat(1,lindata.beta);
0015 plotdata(1).values=beta;
0016 plotdata(1).labels={'\beta_x','\beta_z'};
0017 plotdata(1).axislabel='\beta [m]';
0018
0019 dispersion=cat(2,lindata.Dispersion)';
0020
0021 plotdata(2).values=[dispersion(:,1) s'];
0022 plotdata(2).labels={'\eta_x ','\sigma_x [mm]','\sigma_y [mm]'};
0023 plotdata(2).axislabel='dispersion [m] and \sigma [mm]';
0024
0025
0026
0027
0028
0029
0030 return