-
Notifications
You must be signed in to change notification settings - Fork 3
/
sim_model_amb_plot.m
53 lines (45 loc) · 1.73 KB
/
sim_model_amb_plot.m
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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
%
% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program. If not, see <http://www.gnu.org/licenses/>.
%
clear;
close all;
%profile on
% create an instance to generate the filter function
target = fbmc_chain();
% calculate the ambiguity inteferente with higher resolution
k_k_dot_diff_max = 3;
m_m_dot_diff_max = 5;
f_res = 50;
t_res = 64/2; % must be a devider of target.n_subc/2
[ambiguity_interference, ...
k_k_dot_diffs, ...
m_m_dot_diffs] = fbmc_lib.ambiguity_interference(target, k_k_dot_diff_max, m_m_dot_diff_max, f_res, t_res);
% take the logarithm of the absolute
ambiguity_interference_abs = abs(ambiguity_interference);
ambiguity_interference_abs_log = 20*log10(ambiguity_interference_abs);
% limit to a minimum
minimum = -100;
ambiguity_interference_abs_log(ambiguity_interference_abs_log < minimum) = minimum;
% plot 2d picture, we only need one vector
h = pcolor(m_m_dot_diffs(1,:), k_k_dot_diffs(:,1), ambiguity_interference_abs_log);
set(h, 'EdgeColor', 'none');
colorbar()
colormap(jet)
% labels
title('Ambiguity Function')
xlabel('\tau / T/2')
ylabel('\nu / carrier spacing 1/K')
%profile off
%profile viewer
disp('Simulation finished')