1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
Nlevels_with_doppler_with_z_4wm
%% field propagation
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);
Npl=Nv/2+1;
Iv_out_v = reshape(I2_out_v_3(1:Npl,:,Nz),Npl,Nt);
for i=1:Npl
l_str{i}=num2str(v_3(i));
end
figure(1); hold off;
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(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)');
|
