xmds2/Nlevels_with_doppler_with_z_4wm/pp_doppler.m


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

Nlevels_with_doppler_with_z_4wm

%% field propagation
z_1=z_1*100; % z in cm
t_1=t_1*1e6; % time now measured in uS

[Nv, Nt, Nz]=size(I2_out_1);

Iv_out = reshape(I2_out_1(1:Nv/2,:,Nz),Nv/2,Nt);

for i=1:Nv/2
	l_str{i}=num2str(v_1(i));
end
l_str{Nv/2+1}=num2str('Doppler averaged');


figure(1); hold off;
plot(Iv_out'); hold on
plot(I2_out_avgd_3(:,end), 'LineWidth',2);
legend(l_str);