diff options
| -rw-r--r-- | xmds2/Nlevels_with_doppler_with_z_4wm/pp_doppler.m | 19 |
1 files changed, 14 insertions, 5 deletions
diff --git a/xmds2/Nlevels_with_doppler_with_z_4wm/pp_doppler.m b/xmds2/Nlevels_with_doppler_with_z_4wm/pp_doppler.m index 8d1a219..d98cc0e 100644 --- a/xmds2/Nlevels_with_doppler_with_z_4wm/pp_doppler.m +++ b/xmds2/Nlevels_with_doppler_with_z_4wm/pp_doppler.m @@ -1,8 +1,10 @@ Nlevels_with_doppler_with_z_4wm %% field propagation -z_3=z_1*100; % z in cm -t_3=t_1*1e6; % time now measured in uS +z_3=z_3*100; % z in cm +z_1=z_1*100; % z in cm +t_3=t_3*1e6; % time now measured in uS +t_1=t_1*1e6; % time now measured in uS [Nv, Nt, Nz]=size(I2_out_v_3); @@ -15,12 +17,19 @@ end figure(1); hold off; -plot(Iv_out_v'); hold on -plot(I2_out_1(:,end), '-b', 'LineWidth',2); +plot(t_3, Iv_out_v'); hold on + +plot(t_1, I2_out_1(:,end), '-b', 'LineWidth',2); l_str{Npl+1}=num2str('Doppler averaged out'); -plot(I2_out_1(:,1), '-r', 'LineWidth',2); + +plot(t_1, I2_out_1(:,1), '-r', 'LineWidth',2); l_str{Npl+2}=num2str('Doppler averaged in'); + legend(l_str); +xlabel('time (uS)'); + figure(2); hold off; imagesc(z_1, t_1, I2_out_1) +ylabel('time (uS)'); +xlabel('z (cm)'); |