308 {-0.21, 0.977701, 0, 0, 0, 0},
309 {-0.48, 0.459445, 0.747335, 0, 0, 0},
310 {-0.05, 0.377927, 0.28416, 0.879729, 0, 0},
311 {-0.02, 0.691213, 0.258017, 0.073265, 0.670734, 0},
312 {-0.31, -0.00521632, -0.115615, 0.0788023, 0.00218104, 0.940368},
317 {-0.25, 0.968246, 0, 0, 0, 0},
318 {-0.52, 0.35115, 0.778648, 0, 0, 0},
319 {-0.04, 0.371806, 0.345008, 0.860889, 0, 0},
320 {0, 0.743613, 0.281102, 0.0424415, 0.605161, 0},
321 {-0.32, 0.0206559, -0.0689057, 0.154832, 0.061865, 0.929852},
326 {-0.21, 0.977701, 0, 0, 0, 0},
327 {-0.52, 0.450853, 0.725487, 0, 0, 0},
328 {-0.04, 0.288023, 0.260989, 0.920504, 0, 0},
329 {0.01, 0.697657, 0.386856, 0.0418183, 0.601472, 0},
330 {-0.33, 0.0416283, -0.0694268, 0.166137, 0.139937, 0.915075},
335 {-0.26, 0.965609, 0, 0, 0, 0},
336 {-0.53, 0.395813, 0.749955, 0, 0, 0},
337 {-0.04, 0.299914, 0.320139, 0.897754, 0, 0},
338 {0.01, 0.696556, 0.372815, 0.0580784, 0.610202, 0},
339 {-0.33, 0.0457742, -0.0173584, 0.154417, 0.129332, 0.920941},
344 {-0.25, 0.968246, 0, 0, 0, 0},
345 {-0.57, 0.420864, 0.705672, 0, 0, 0},
346 {-0.03, 0.229797, 0.235501, 0.943839, 0, 0},
347 {0.03, 0.679063, 0.384466, 0.0681379, 0.6209, 0},
348 {-0.41, -0.147173, -0.229228, 0.270707, 0.293002, 0.773668},
353 {-0.2, 0.979796, 0, 0, 0, 0},
354 {-0.53, 0.473568, 0.703444, 0, 0, 0},
355 {-0.05, 0.204124, 0.109225, 0.971547, 0, 0},
356 {0.03, 0.68994, 0.411073, 0.0676935, 0.591202, 0},
357 {-0.4, -0.224537, -0.292371, 0.275609, 0.301835, 0.732828},
362 {-0.19, 0.981784, 0, 0, 0, 0},
363 {-0.5, 0.524555, 0.689088, 0, 0, 0},
364 {-0.03, 0.228462, 0.18163, 0.955989, 0, 0},
365 {-0.02, 0.637818, 0.428725, 0.00608114, 0.639489, 0},
366 {-0.36, -0.18171, -0.282523, 0.106726, 0.123808, 0.854894},
371 {-0.2, 0.979796, 0, 0, 0, 0},
372 {-0.49, 0.502145, 0.712566, 0, 0, 0},
373 {-0.01, 0.232702, 0.151916, 0.960558, 0, 0},
374 {-0.05, 0.612372, 0.376106, 0.0206792, 0.693265, 0},
375 {-0.37, -0.320475, -0.365405, -0.00376264, 0.0364343, 0.790907},
380 {-0.19, 0.981784, 0, 0, 0, 0},
381 {-0.38, 0.58852, 0.713613, 0, 0, 0},
382 {-0.03, 0.360874, 0.12082, 0.924269, 0, 0},
383 {-0.12, 0.526796, 0.34244, 0.0594196, 0.766348, 0},
384 {-0.33, -0.257389, -0.24372, -0.257035, -0.176521, 0.817451},
455 {-0.36, 0.932952, 0, 0, 0, 0},
456 {0.45, 0.516639, 0.728412, 0, 0, 0},
457 {0.02, 0.329277, 0.371881, 0.867687, 0, 0},
458 {-0.06, 0.59853, 0.436258, -0.0324062, 0.668424, 0},
459 {-0.07, 0.0373009, 0.305087, -0.0280496, -0.225204, 0.921481},
464 {-0.39, 0.920815, 0, 0, 0, 0},
465 {0.52, 0.426579, 0.740021, 0, 0, 0},
466 {0, 0.347518, -0.0381664, 0.936896, 0, 0},
467 {-0.04, 0.710675, 0.172483, 0.116993, 0.670748, 0},
468 {-0.17, -0.0394216, 0.115154, 0.243458, -0.0702635, 0.944498},
473 {-0.41, 0.912086, 0, 0, 0, 0},
474 {0.54, 0.49491, 0.680782, 0, 0, 0},
475 {0, 0.350844, -0.152231, 0.923977, 0, 0},
476 {-0.04, 0.694672, 0.0702137, 0.0832998, 0.709903, 0},
477 {-0.19, -0.0854087, 0.0805978, 0.283811, -0.137441, 0.922318},
482 {-0.37, 0.929032, 0, 0, 0, 0},
483 {0.53, 0.480177, 0.698949, 0, 0, 0},
484 {0.01, 0.434538, 0.00864797, 0.900556, 0, 0},
485 {-0.05, 0.765851, -0.0303947, 0.0421641, 0.63896, 0},
486 {-0.17, -0.16458, 0.0989022, 0.158081, -0.150425, 0.941602},
491 {-0.4, 0.916515, 0, 0, 0, 0},
492 {0.55, 0.403703, 0.731111, 0, 0, 0},
493 {0.02, 0.499719, -0.0721341, 0.862947, 0, 0},
494 {-0.06, 0.835775, -0.156481, 0.0373835, 0.521534, 0},
495 {-0.19, -0.301141, 0.145082, 0.144564, -0.0238067, 0.911427},
500 {-0.41, 0.912086, 0, 0, 0, 0},
501 {0.56, 0.339442, 0.755764, 0, 0, 0},
502 {0.02, 0.436582, -0.0256617, 0.899076, 0, 0},
503 {-0.07, 0.856608, -0.12116, 0.0715303, 0.491453, 0},
504 {-0.2, -0.331109, 0.15136, 0.036082, 0.031313, 0.908391},
509 {-0.4, 0.916515, 0, 0, 0, 0},
510 {0.56, 0.386246, 0.732949, 0, 0, 0},
511 {0.04, 0.573913, -0.0601289, 0.815726, 0, 0},
512 {-0.11, 0.813953, -0.0720183, 0.0281118, 0.565158, 0},
513 {-0.19, -0.432071, 0.236423, -0.0247788, -0.0557206, 0.847113},
518 {-0.46, 0.887919, 0, 0, 0, 0},
519 {0.58, 0.469412, 0.665772, 0, 0, 0},
520 {0.01, 0.309262, -0.286842, 0.90663, 0, 0},
521 {-0.05, 0.762457, -0.268721, -0.0467443, 0.584605, 0},
522 {-0.23, -0.580909, 0.399665, 0.0403629, 0.326208, 0.584698},
527 {-0.3, 0.953939, 0, 0, 0, 0},
528 {0.47, 0.81871, 0.329868, 0, 0, 0},
529 {0.06, 0.0712834, -0.595875, 0.797654, 0, 0},
530 {-0.1, 0.408831, -0.0233859, 0.0412736, 0.905873, 0},
531 {-0.13, -0.407783, 0.439436, -0.0768289, -0.212875, 0.756631},
550 {-0.36, 0.932952, 0, 0, 0, 0},
551 {0.45, 0.527358, 0.72069, 0, 0, 0},
552 {0.02, 0.350715, 0.355282, 0.866241, 0, 0},
553 {-0.07, 0.562515, 0.478504, 0.0162932, 0.670406, 0},
554 {-0.06, 0.0411597, 0.270982, 0.0121094, -0.159927, 0.946336},
559 {-0.38, 0.924986, 0, 0, 0, 0},
560 {0.52, 0.473088, 0.711188, 0, 0, 0},
561 {0, 0.367573, -0.0617198, 0.927944, 0, 0},
562 {-0.04, 0.68628, 0.149228, 0.115257, 0.701332, 0},
563 {-0.16, -0.0441088, 0.118207, 0.251641, -0.0752458, 0.943131},
568 {-0.42, 0.907524, 0, 0, 0, 0},
569 {0.54, 0.48131, 0.690464, 0, 0, 0},
570 {0, 0.363627, -0.137613, 0.921324, 0, 0},
571 {-0.04, 0.686704, 0.117433, 0.104693, 0.708581, 0},
572 {-0.19, -0.0438556, 0.0922685, 0.269877, -0.136292, 0.928469},
577 {-0.36, 0.932952, 0, 0, 0, 0},
578 {0.53, 0.483197, 0.696865, 0, 0, 0},
579 {0.01, 0.464761, -0.0285153, 0.88492, 0, 0},
580 {-0.05, 0.763169, 0.140255, 0.0562856, 0.626286, 0},
581 {-0.16, -0.126051, 0.0942905, 0.195354, -0.217188, 0.92967},
586 {-0.39, 0.920815, 0, 0, 0, 0},
587 {0.55, 0.406705, 0.729446, 0, 0, 0},
588 {0.01, 0.503793, -0.123923, 0.854831, 0, 0},
589 {-0.06, 0.821664, -0.207246, 0.0245302, 0.526988, 0},
590 {-0.19, -0.254231, 0.10679, 0.190931, -0.0665276, 0.920316},
595 {-0.42, 0.907524, 0, 0, 0, 0},
596 {0.56, 0.391395, 0.730213, 0, 0, 0},
597 {0.02, 0.427978, -0.0393147, 0.902712, 0, 0},
598 {-0.06, 0.820694, -0.119986, 0.105509, 0.545281, 0},
599 {-0.2, -0.279882, 0.180145, 0.0563477, -0.0121631, 0.919723},
604 {-0.36, 0.932952, 0, 0, 0, 0},
605 {0.54, 0.519212, 0.662434, 0, 0, 0},
606 {0.04, 0.412025, -0.0234416, 0.909992, 0, 0},
607 {-0.09, 0.758452, -0.0682296, 0.0214276, 0.64151, 0},
608 {-0.17, -0.387158, 0.306169, -0.0291255, -0.109344, 0.845378},
613 {-0.44, 0.897998, 0, 0, 0, 0},
614 {0.57, 0.43519, 0.696928, 0, 0, 0},
615 {0.01, 0.316705, -0.248988, 0.915207, 0, 0},
616 {-0.06, 0.805793, -0.296262, -0.0419182, 0.507514, 0},
617 {-0.22, -0.497551, 0.289742, 0.0785823, 0.328773, 0.711214},
622 {-0.27, 0.96286, 0, 0, 0, 0},
623 {0.46, 0.741748, 0.488067, 0, 0, 0},
624 {0.04, 0.0735309, -0.374828, 0.923308, 0, 0},
625 {-0.08, 0.517624, 0.128779, 0.0795063, 0.838308, 0},
626 {-0.11, -0.321646, 0.0802763, -0.131981, -0.193429, 0.907285},
668 {10, {-7.12, 0.8, -3.06, 0.48, 0.94, 0.7, 0.82, 0.03, -2.52, 0.5, 4.4,
669 3.3, 2.5, 24.4, 3.8, 3.0, 20.0, 3.9, 0.0, 11.0, 7.0, 3.0}},
670 {20, {-7.28, 0.67, -2.68, 0.36, 0.87, 0.66, 0.5, 0.09, -2.29, 0.53, 9.0,
671 6.6, 2.5, 23.6, 4.7, 3.0, 20.0, 3.9, 0.0, 11.0, 7.0, 3.0}},
672 {30, {-7.45, 0.68, -2.51, 0.38, 0.92, 0.68, 0.82, 0.05, -2.19, 0.58, 9.3,
673 6.1, 2.5, 23.2, 4.6, 3.0, 20.0, 3.9, 0.0, 11.0, 7.0, 3.0}},
674 {40, {-7.73, 0.66, -2.4, 0.32, 0.79, 0.64, 1.23, 0.03, -2.24, 0.51, 7.9,
675 4.0, 2.5, 22.6, 4.9, 3.0, 20.0, 3.9, 0.0, 11.0, 7.0, 3.0}},
676 {50, {-7.91, 0.62, -2.31, 0.33, 0.72, 0.63, 1.43, 0.06, -2.3, 0.46, 7.4,
677 3.0, 2.5, 21.8, 5.7, 3.0, 20.0, 3.9, 0.0, 11.0, 7.0, 3.0}},
678 {60, {-8.14, 0.51, -2.2, 0.39, 0.6, 0.54, 1.56, 0.05, -2.48, 0.35, 7.0,
679 2.6, 2.5, 20.5, 6.9, 3.0, 20.0, 3.9, 0.0, 11.0, 7.0, 3.0}},
680 {70, {-8.23, 0.45, -2.0, 0.4, 0.55, 0.52, 1.66, 0.05, -2.64, 0.31, 6.9,
681 2.2, 2.5, 19.3, 8.1, 3.0, 20.0, 3.9, 0.0, 11.0, 7.0, 3.0}},
682 {80, {-8.28, 0.31, -1.64, 0.32, 0.71, 0.53, 1.73, 0.02, -2.68, 0.39, 6.5,
683 2.1, 2.5, 17.4, 10.3, 3.0, 20.0, 3.9, 0.0, 11.0, 7.0, 3.0}},
684 {90, {-8.36, 0.08, -0.63, 0.53, 0.81, 0.62, 1.79, 0.01, -2.61, 0.28, 6.8,
685 1.9, 2.5, 12.3, 15.2, 3.0, 20.0, 3.9, 0.0, 11.0, 7.0, 3.0}},
689 {10, {-7.43, 0.9, -3.43, 0.54, 0.65, 0.82, 0.82, 0.05, -2.75, 0.55, 6.1,
690 2.6, 2.5, 24.7, 2.1, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
691 {20, {-7.62, 0.78, -3.06, 0.41, 0.53, 0.78, 0.47, 0.11, -2.64, 0.64, 13.7,
692 6.8, 2.5, 24.4, 2.8, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
693 {30, {-7.76, 0.8, -2.91, 0.42, 0.6, 0.83, 0.8, 0.05, -2.49, 0.69, 12.9,
694 6.0, 2.5, 24.4, 2.7, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
695 {40, {-8.02, 0.72, -2.81, 0.34, 0.43, 0.78, 1.23, 0.04, -2.51, 0.57, 10.3,
696 3.3, 2.5, 24.2, 2.7, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
697 {50, {-8.13, 0.61, -2.74, 0.34, 0.36, 0.77, 1.42, 0.1, -2.54, 0.5, 9.2,
698 2.2, 2.5, 23.9, 3.1, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
699 {60, {-8.3, 0.47, -2.72, 0.7, 0.16, 0.84, 1.56, 0.06, -2.71, 0.37, 8.4,
700 1.9, 2.5, 23.3, 3.9, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
701 {70, {-8.34, 0.39, -2.46, 0.4, 0.18, 0.64, 1.65, 0.07, -2.85, 0.31, 8.0,
702 1.5, 2.5, 22.6, 4.8, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
703 {80, {-8.39, 0.26, -2.3, 0.78, 0.24, 0.81, 1.73, 0.02, -3.01, 0.45, 7.4,
704 1.6, 2.5, 21.2, 6.8, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
705 {90, {-8.45, 0.01, -1.11, 0.51, 0.36, 0.65, 1.79, 0.01, -3.08, 0.27, 7.6,
706 1.3, 2.5, 17.6, 12.7, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
719 {10, {-6.84, 0.82, -2.08, 0.87, 1.0, 1.6, 1.0, 0.63, -2.08, 0.58,
720 2.3, 23.8, 4.4, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
721 {20, {-6.81, 0.61, -1.68, 0.73, 1.44, 0.87, 0.94, 0.65, -1.66, 0.5,
722 2.3, 21.9, 6.3, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
723 {30, {-6.94, 0.49, -1.46, 0.53, 1.54, 0.64, 1.15, 0.42, -1.48, 0.4,
724 2.3, 19.7, 8.1, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
725 {40, {-7.14, 0.49, -1.43, 0.5, 1.53, 0.56, 1.35, 0.28, -1.46, 0.37,
726 2.3, 18.1, 9.3, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
727 {50, {-7.34, 0.51, -1.44, 0.58, 1.48, 0.54, 1.44, 0.25, -1.53, 0.47,
728 2.3, 16.3, 11.5, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
729 {60, {-7.53, 0.47, -1.33, 0.49, 1.39, 0.68, 1.56, 0.16, -1.61, 0.43,
730 2.3, 14.0, 13.3, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
731 {70, {-7.67, 0.44, -1.31, 0.65, 1.42, 0.55, 1.64, 0.18, -1.77, 0.5,
732 2.3, 12.1, 14.9, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
733 {80, {-7.82, 0.42, -1.11, 0.69, 1.38, 0.6, 1.7, 0.09, -1.9, 0.42,
734 2.3, 8.7, 17.0, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
735 {90, {-7.84, 0.55, -0.11, 0.53, 1.23, 0.6, 1.7, 0.17, -1.99, 0.5,
736 2.3, 6.4, 12.3, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
740 {10, {-6.86, 0.81, -2.12, 0.94, 1.02, 1.44, 1.01, 0.56, -2.11, 0.59,
741 2.3, 23.7, 4.5, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
742 {20, {-6.84, 0.61, -1.74, 0.79, 1.44, 0.77, 0.96, 0.55, -1.69, 0.51,
743 2.3, 21.8, 6.3, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
744 {30, {-7.0, 0.56, -1.56, 0.66, 1.48, 0.7, 1.13, 0.43, -1.52, 0.46,
745 2.3, 19.6, 8.2, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
746 {40, {-7.21, 0.56, -1.54, 0.63, 1.46, 0.6, 1.3, 0.37, -1.51, 0.43,
747 2.3, 18.0, 9.4, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
748 {50, {-7.42, 0.57, -1.45, 0.56, 1.4, 0.59, 1.4, 0.32, -1.54, 0.45,
749 2.3, 16.3, 11.5, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
750 {60, {-7.86, 0.55, -1.64, 0.78, 0.97, 1.27, 1.41, 0.45, -1.84, 0.63,
751 2.3, 15.9, 12.4, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
752 {70, {-7.76, 0.47, -1.37, 0.56, 1.33, 0.56, 1.63, 0.17, -1.86, 0.51,
753 2.3, 12.3, 15.0, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
754 {80, {-8.07, 0.42, -1.29, 0.76, 1.12, 1.04, 1.68, 0.14, -2.16, 0.74,
755 2.3, 10.5, 15.7, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
756 {90, {-7.95, 0.59, -0.41, 0.59, 1.04, 0.63, 1.7, 0.17, -2.21, 0.61,
757 2.3, 10.5, 15.7, 4.0, 20.0, 3.9, 0.0, 15.0, 7.0, 3.0}},
767static const std::map<std::string, std::map<int, std::array<float, 22>>>
NTNUrbanLOS{
770 {10, {-7.97, 1.0, -2.6, 0.79, 0.18, 0.74, -0.63, 2.6, -2.54, 2.62, 31.83,
771 13.84, 2.5, 8.0, 4.0, 4.0, 20.0, 3.9, 0.09, 12.55, 1.25, 3.0}},
772 {20, {-8.12, 0.83, -2.48, 0.8, 0.42, 0.9, -0.15, 3.31, -2.67, 2.96, 18.78,
773 13.78, 2.5, 8.0, 4.0, 3.0, 20.0, 3.9, 0.09, 12.76, 3.23, 3.0}},
774 {30, {-8.21, 0.68, -2.44, 0.91, 0.41, 1.3, 0.54, 1.1, -2.03, 0.86, 10.49,
775 10.42, 2.5, 8.0, 4.0, 3.0, 20.0, 3.9, 0.12, 14.36, 4.39, 3.0}},
776 {40, {-8.31, 0.48, -2.6, 1.02, 0.18, 1.69, 0.35, 1.59, -2.28, 1.19, 7.46,
777 8.01, 2.5, 8.0, 4.0, 3.0, 20.0, 3.9, 0.16, 16.42, 5.72, 3.0}},
778 {50, {-8.37, 0.38, -2.71, 1.17, -0.07, 2.04, 0.27, 1.62, -2.48, 1.4, 6.52,
779 8.27, 2.5, 8.0, 4.0, 3.0, 20.0, 3.9, 0.2, 17.13, 6.17, 3.0}},
780 {60, {-8.39, 0.24, -2.76, 1.17, -0.43, 2.54, 0.26, 0.97, -2.56, 0.85, 5.47,
781 7.26, 2.5, 8.0, 4.0, 3.0, 20.0, 3.9, 0.28, 19.01, 7.36, 3.0}},
782 {70, {-8.38, 0.18, -2.78, 1.2, -0.64, 2.47, -0.12, 1.99, -2.96, 1.61, 4.54,
783 5.53, 2.5, 8.0, 4.0, 3.0, 20.0, 3.9, 0.44, 19.31, 7.3, 3.0}},
784 {80, {-8.35, 0.13, -2.65, 1.45, -0.91, 2.69, -0.21, 1.82, -3.08, 1.49, 4.03,
785 4.49, 2.5, 8.0, 4.0, 3.0, 20.0, 3.9, 0.9, 22.39, 7.7, 3.0}},
786 {90, {-8.34, 0.09, -2.27, 1.85, -0.54, 1.66, -0.07, 1.43, -3.0, 1.09, 3.68,
787 3.14, 2.5, 8.0, 4.0, 3.0, 20.0, 3.9, 2.87, 27.8, 9.25, 3.0}},
791 {10, {-8.52, 0.92, -3.18, 0.79, -0.4, 0.77, -0.67, 2.22, -2.61, 2.41, 40.18,
792 16.99, 2.5, 8.0, 4.0, 4.0, 20.0, 1.6, 0.09, 11.8, 1.14, 3.0}},
793 {20, {-8.59, 0.79, -3.05, 0.87, -0.15, 0.97, -0.34, 3.04, -2.82, 2.59, 23.62,
794 18.96, 2.5, 8.0, 4.0, 3.0, 20.0, 1.6, 0.09, 11.6, 2.78, 3.0}},
795 {30, {-8.51, 0.65, -2.98, 1.04, -0.18, 1.58, 0.07, 1.33, -2.48, 1.02, 12.48,
796 14.23, 2.5, 8.0, 4.0, 3.0, 20.0, 1.6, 0.11, 13.05, 3.87, 3.0}},
797 {40, {-8.49, 0.48, -3.11, 1.06, -0.31, 1.69, -0.08, 1.45, -2.76, 1.27, 8.56,
798 11.06, 2.5, 8.0, 4.0, 3.0, 20.0, 1.6, 0.15, 14.56, 4.94, 3.0}},
799 {50, {-8.48, 0.46, -3.19, 1.12, -0.58, 2.13, -0.21, 1.62, -2.93, 1.38, 7.42,
800 11.21, 2.5, 8.0, 4.0, 3.0, 20.0, 1.6, 0.18, 15.35, 5.41, 3.0}},
801 {60, {-8.44, 0.34, -3.25, 1.14, -0.9, 2.51, -0.25, 1.06, -3.05, 0.96, 5.97,
802 9.47, 2.5, 8.0, 4.0, 3.0, 20.0, 1.6, 0.27, 16.97, 6.31, 3.0}},
803 {70, {-8.4, 0.27, -3.33, 1.25, -1.16, 2.47, -0.61, 1.88, -3.45, 1.51, 4.88,
804 7.24, 2.5, 8.0, 4.0, 3.0, 20.0, 1.6, 0.42, 17.96, 6.66, 3.0}},
805 {80, {-8.37, 0.19, -3.22, 1.35, -1.48, 2.61, -0.79, 1.87, -3.66, 1.49, 4.22,
806 5.79, 2.5, 8.0, 4.0, 3.0, 20.0, 1.6, 0.86, 20.68, 7.31, 3.0}},
807 {90, {-8.35, 0.14, -2.83, 1.62, -1.14, 1.7, -0.58, 1.19, -3.56, 0.89, 3.81,
808 4.25, 2.5, 8.0, 4.0, 3.0, 20.0, 1.6, 2.55, 25.08, 9.23, 3.0}},
818static const std::map<std::string, std::map<int, std::array<float, 22>>>
NTNUrbanNLOS{
821 {10, {-7.24, 1.26, -1.58, 0.89, 0.13, 2.99, -1.13, 2.66, -2.87, 2.76, 0.0,
822 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.08, 14.72, 1.57, 3.0}},
823 {20, {-7.7, 0.99, -1.67, 0.89, 0.19, 3.12, 0.49, 2.03, -2.68, 2.76, 0.0,
824 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.1, 14.62, 4.3, 3.0}},
825 {30, {-7.82, 0.86, -1.84, 1.3, 0.44, 2.69, 0.95, 1.54, -2.12, 1.54, 0.0,
826 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.14, 16.4, 6.64, 3.0}},
827 {40, {-8.04, 0.75, -2.02, 1.15, 0.48, 2.45, 1.15, 1.02, -2.27, 1.77, 0.0,
828 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.22, 17.86, 9.21, 3.0}},
829 {50, {-8.08, 0.77, -2.06, 1.23, 0.56, 2.17, 1.14, 1.61, -2.5, 2.36, 0.0,
830 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.31, 19.74, 10.32, 3.0}},
831 {60, {-8.1, 0.76, -1.99, 1.02, 0.55, 1.93, 1.13, 1.84, -2.47, 2.33, 0.0,
832 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.49, 19.73, 10.3, 3.0}},
833 {70, {-8.16, 0.73, -2.19, 1.78, 0.48, 1.72, 1.16, 1.81, -2.83, 2.84, 0.0,
834 0.0, 2.3, 7.0, 3.0, 2.0, 20.0, 1.6, 0.97, 20.5, 10.2, 3.0}},
835 {80, {-8.03, 0.79, -1.88, 1.55, 0.53, 1.51, 1.28, 1.35, -2.82, 2.87, 0.0,
836 0.0, 2.3, 7.0, 3.0, 2.0, 20.0, 1.6, 1.52, 26.16, 12.27, 3.0}},
837 {90, {-8.33, 0.7, -2.0, 1.4, 0.32, 1.2, 1.42, 0.6, -4.55, 4.27, 0.0,
838 0.0, 2.3, 7.0, 3.0, 2.0, 20.0, 1.6, 5.36, 25.83, 12.75, 3.0}},
842 {10, {-7.24, 1.26, -1.58, 0.89, 0.13, 2.99, -1.13, 2.66, -2.87, 2.76, 0.0,
843 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.08, 14.72, 1.57, 3.0}},
844 {20, {-7.7, 0.99, -1.67, 0.89, 0.19, 3.12, 0.49, 2.03, -2.68, 2.76, 0.0,
845 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.1, 14.62, 4.3, 3.0}},
846 {30, {-7.82, 0.86, -1.84, 1.3, 0.44, 2.69, 0.95, 1.54, -2.12, 1.54, 0.0,
847 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.14, 16.4, 6.64, 3.0}},
848 {40, {-8.04, 0.75, -2.02, 1.15, 0.48, 2.45, 1.15, 1.02, -2.27, 1.77, 0.0,
849 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.22, 17.86, 9.21, 3.0}},
850 {50, {-8.08, 0.77, -2.06, 1.23, 0.56, 2.17, 1.14, 1.61, -2.5, 2.36, 0.0,
851 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.31, 19.74, 10.32, 3.0}},
852 {60, {-8.1, 0.76, -1.99, 1.02, 0.55, 1.93, 1.13, 1.84, -2.47, 2.33, 0.0,
853 0.0, 2.3, 7.0, 3.0, 3.0, 20.0, 1.6, 0.49, 19.73, 10.3, 3.0}},
854 {70, {-8.16, 0.73, -2.19, 1.78, 0.48, 1.72, 1.16, 1.81, -2.83, 2.84, 0.0,
855 0.0, 2.3, 7.0, 3.0, 2.0, 20.0, 1.6, 0.97, 20.5, 10.2, 3.0}},
856 {80, {-8.03, 0.79, -1.88, 1.55, 0.53, 1.51, 1.28, 1.35, -2.82, 2.87, 0.0,
857 0.0, 2.3, 7.0, 3.0, 2.0, 20.0, 1.6, 1.52, 26.16, 12.27, 3.0}},
858 {90, {-8.33, 0.7, -2.0, 1.4, 0.32, 1.2, 1.42, 0.6, -4.55, 4.27, 0.0,
859 0.0, 2.3, 7.0, 3.0, 2.0, 20.0, 1.6, 5.36, 25.83, 12.75, 3.0}},
869static const std::map<std::string, std::map<int, std::array<float, 22>>>
NTNSuburbanLOS{
872 {10, {-8.16, 0.99, -3.57, 1.62, 0.05, 1.84, -1.78, 0.62, -1.06, 0.96, 11.4,
873 6.26, 2.2, 21.3, 7.6, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
874 {20, {-8.56, 0.96, -3.8, 1.74, -0.38, 1.94, -1.84, 0.81, -1.21, 0.95, 19.45,
875 10.32, 3.36, 21.0, 8.9, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
876 {30, {-8.72, 0.79, -3.77, 1.72, -0.56, 1.75, -1.67, 0.57, -1.28, 0.49, 20.8,
877 16.34, 3.5, 21.2, 8.5, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
878 {40, {-8.71, 0.81, -3.57, 1.6, -0.59, 1.82, -1.59, 0.86, -1.32, 0.79, 21.2,
879 15.63, 2.81, 21.1, 8.4, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
880 {50, {-8.72, 1.12, -3.42, 1.49, -0.58, 1.87, -1.55, 1.05, -1.39, 0.97, 21.6,
881 14.22, 2.39, 20.7, 9.2, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
882 {60, {-8.66, 1.23, -3.27, 1.43, -0.55, 1.92, -1.51, 1.23, -1.36, 1.17, 19.75,
883 14.19, 2.73, 20.6, 9.8, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
884 {70, {-8.38, 0.55, -3.08, 1.36, -0.28, 1.16, -1.27, 0.54, -1.08, 0.62, 12.0,
885 5.7, 2.07, 20.3, 10.8, 2.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
886 {80, {-8.34, 0.63, -2.75, 1.26, -0.17, 1.09, -1.28, 0.67, -1.31, 0.76, 12.85,
887 9.91, 2.04, 19.8, 12.2, 2.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
888 {90, {-8.34, 0.63, -2.75, 1.26, -0.17, 1.09, -1.28, 0.67, -1.31, 0.76, 12.85,
889 9.91, 2.04, 19.1, 13.0, 2.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
893 {10, {-8.07, 0.46, -3.55, 0.48, 0.89, 0.67, 0.63, 0.35, -3.37, 0.28, 8.9,
894 4.4, 2.5, 23.2, 5.0, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
895 {20, {-8.61, 0.45, -3.69, 0.41, 0.31, 0.78, 0.76, 0.3, -3.28, 0.27, 14.0,
896 4.6, 2.5, 23.6, 4.5, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
897 {30, {-8.72, 0.28, -3.59, 0.41, 0.02, 0.75, 1.11, 0.28, -3.04, 0.26, 11.3,
898 3.7, 2.5, 23.5, 4.7, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
899 {40, {-8.63, 0.17, -3.38, 0.35, -0.1, 0.65, 1.37, 0.23, -2.88, 0.21, 9.0,
900 3.5, 2.5, 23.4, 5.2, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
901 {50, {-8.54, 0.14, -3.23, 0.35, -0.19, 0.55, 1.53, 0.23, -2.83, 0.18, 7.5,
902 3.0, 2.5, 23.2, 5.7, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
903 {60, {-8.48, 0.15, -3.19, 0.43, -0.54, 0.96, 1.65, 0.17, -2.86, 0.17, 6.6,
904 2.6, 2.5, 23.3, 5.9, 3.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
905 {70, {-8.42, 0.09, -2.83, 0.33, -0.24, 0.43, 1.74, 0.11, -2.95, 0.1, 5.9,
906 1.7, 2.5, 23.4, 6.2, 2.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
907 {80, {-8.39, 0.05, -2.66, 0.44, -0.52, 0.93, 1.82, 0.05, -3.21, 0.07, 5.5,
908 0.7, 2.5, 23.2, 7.0, 2.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
909 {90, {-8.37, 0.02, -1.22, 0.31, -0.15, 0.44, 1.87, 0.02, -3.49, 0.24, 5.4,
910 0.3, 2.5, 23.1, 7.6, 2.0, 20.0, 1.6, 0.0, 11.0, 7.0, 3.0}},
920static const std::map<std::string, std::map<int, std::array<float, 22>>>
NTNSuburbanNLOS{
923 {10, {-7.43, 0.5, -2.89, 0.41, 1.49, 0.4, 0.81, 0.36, -3.09, 0.32, 0.0,
924 0.0, 2.3, 22.5, 5.0, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
925 {20, {-7.63, 0.61, -2.76, 0.41, 1.24, 0.82, 1.06, 0.41, -2.93, 0.47, 0.0,
926 0.0, 2.3, 19.4, 8.5, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
927 {30, {-7.86, 0.56, -2.64, 0.41, 1.06, 0.71, 1.12, 0.4, -2.91, 0.46, 0.0,
928 0.0, 2.3, 15.5, 10.0, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
929 {40, {-7.96, 0.58, -2.41, 0.52, 0.91, 0.55, 1.14, 0.39, -2.78, 0.54, 0.0,
930 0.0, 2.3, 13.9, 10.6, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
931 {50, {-7.98, 0.59, -2.42, 0.7, 0.98, 0.58, 1.29, 0.35, -2.7, 0.45, 0.0,
932 0.0, 2.3, 11.7, 10.0, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
933 {60, {-8.45, 0.47, -2.53, 0.5, 0.49, 1.37, 1.38, 0.36, -3.03, 0.36, 0.0,
934 0.0, 2.3, 9.8, 9.1, 3.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
935 {70, {-8.21, 0.36, -2.35, 0.58, 0.73, 0.49, 1.36, 0.29, -2.9, 0.42, 0.0,
936 0.0, 2.3, 10.3, 9.1, 3.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
937 {80, {-8.69, 0.29, -2.31, 0.73, -0.04, 1.48, 1.38, 0.2, -3.2, 0.3, 0.0,
938 0.0, 2.3, 15.6, 9.1, 3.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
939 {90, {-8.69, 0.29, -2.31, 0.73, -0.04, 1.48, 1.38, 0.2, -3.2, 0.3, 0.0,
940 0.0, 2.3, 15.6, 9.1, 3.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
944 {10, {-7.43, 0.5, -2.89, 0.41, 1.49, 0.4, 0.81, 0.36, -3.09, 0.32, 0.0,
945 0.0, 2.3, 22.5, 5.0, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
946 {20, {-7.63, 0.61, -2.76, 0.41, 1.24, 0.82, 1.06, 0.41, -2.93, 0.47, 0.0,
947 0.0, 2.3, 19.4, 8.5, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
948 {30, {-7.86, 0.56, -2.64, 0.41, 1.06, 0.71, 1.12, 0.4, -2.91, 0.46, 0.0,
949 0.0, 2.3, 15.5, 10.0, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
950 {40, {-7.96, 0.58, -2.41, 0.52, 0.91, 0.55, 1.14, 0.39, -2.78, 0.54, 0.0,
951 0.0, 2.3, 13.9, 10.6, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
952 {50, {-7.98, 0.59, -2.42, 0.7, 0.98, 0.58, 1.29, 0.35, -2.7, 0.45, 0.0,
953 0.0, 2.3, 11.7, 10.0, 4.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
954 {60, {-8.45, 0.47, -2.53, 0.5, 0.49, 1.37, 1.38, 0.36, -3.03, 0.36, 0.0,
955 0.0, 2.3, 9.8, 9.1, 3.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
956 {70, {-8.21, 0.36, -2.35, 0.58, 0.73, 0.49, 1.36, 0.29, -2.9, 0.42, 0.0,
957 0.0, 2.3, 10.3, 9.1, 3.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
958 {80, {-8.69, 0.29, -2.31, 0.73, -0.04, 1.48, 1.38, 0.2, -3.2, 0.3, 0.0,
959 0.0, 2.3, 15.6, 9.1, 3.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
960 {90, {-8.69, 0.29, -2.31, 0.73, -0.04, 1.48, 1.38, 0.2, -3.2, 0.3, 0.0,
961 0.0, 2.3, 15.6, 9.1, 3.0, 20.0, 1.6, 0.0, 15.0, 7.0, 3.0}},
971static const std::map<std::string, std::map<int, std::array<float, 22>>>
NTNRuralLOS{
974 {10, {-9.55, 0.66, -3.42, 0.89, -9.45, 7.83, -4.2, 6.3, -6.03, 5.19, 24.72,
975 5.07, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.39, 10.81, 1.94, 3.0}},
976 {20, {-8.68, 0.44, -3.0, 0.63, -4.45, 6.86, -2.31, 5.04, -4.31, 4.18, 12.31,
977 5.75, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.31, 8.09, 1.83, 3.0}},
978 {30, {-8.46, 0.28, -2.86, 0.52, -2.39, 5.14, -0.28, 0.81, -2.57, 0.61, 8.05,
979 5.46, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.29, 13.7, 2.28, 3.0}},
980 {40, {-8.36, 0.19, -2.78, 0.45, -1.28, 3.44, -0.38, 1.16, -2.59, 0.79, 6.21,
981 5.23, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.37, 20.05, 2.93, 3.0}},
982 {50, {-8.29, 0.14, -2.7, 0.42, -0.99, 2.59, -0.38, 0.82, -2.59, 0.65, 5.04,
983 3.95, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.61, 24.51, 2.84, 3.0}},
984 {60, {-8.26, 0.1, -2.66, 0.41, -1.05, 2.42, -0.46, 0.67, -2.65, 0.52, 4.42,
985 3.75, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.9, 26.35, 3.17, 3.0}},
986 {70, {-8.22, 0.1, -2.53, 0.42, -0.9, 1.78, -0.49, 1.0, -2.69, 0.78, 3.92,
987 2.56, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 1.43, 31.84, 3.88, 3.0}},
988 {80, {-8.2, 0.05, -2.21, 0.5, -0.89, 1.65, -0.53, 1.18, -2.65, 1.01, 3.65,
989 1.77, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 2.87, 36.62, 4.17, 3.0}},
990 {90, {-8.19, 0.06, -1.78, 0.91, -0.81, 1.26, -0.46, 0.91, -2.65, 0.71, 3.59,
991 1.77, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 5.48, 36.77, 4.29, 3.0}},
995 {10, {-9.68, 0.46, -4.03, 0.91, -9.74, 7.52, -5.85, 6.51, -7.45, 5.3, 25.43,
996 7.04, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.36, 4.63, 0.75, 3.0}},
997 {20, {-8.86, 0.29, -3.55, 0.7, -4.88, 6.67, -3.27, 5.36, -5.25, 4.42, 12.72,
998 7.47, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.3, 6.83, 1.25, 3.0}},
999 {30, {-8.59, 0.18, -3.45, 0.55, -2.6, 4.63, -0.88, 0.93, -3.16, 0.68, 8.4,
1000 7.18, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.25, 12.91, 1.93, 3.0}},
1001 {40, {-8.46, 0.19, -3.38, 0.52, -1.92, 3.45, -0.93, 0.96, -3.15, 0.73, 6.52,
1002 6.88, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.35, 18.9, 2.37, 3.0}},
1003 {50, {-8.36, 0.14, -3.33, 0.46, -1.56, 2.44, -0.99, 0.97, -3.2, 0.77, 5.24,
1004 5.28, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.53, 22.44, 2.66, 3.0}},
1005 {60, {-8.3, 0.15, -3.29, 0.43, -1.66, 2.38, -1.04, 0.83, -3.27, 0.61, 4.57,
1006 4.92, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 0.88, 25.69, 3.23, 3.0}},
1007 {70, {-8.26, 0.13, -3.24, 0.46, -1.59, 1.67, -1.17, 1.01, -3.42, 0.74, 4.02,
1008 3.4, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 1.39, 27.95, 3.71, 3.0}},
1009 {80, {-8.22, 0.03, -2.9, 0.44, -1.58, 1.44, -1.19, 1.01, -3.36, 0.79, 3.7,
1010 2.22, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 2.7, 31.45, 4.17, 3.0}},
1011 {90, {-8.21, 0.07, -2.5, 0.82, -1.51, 1.13, -1.13, 0.85, -3.35, 0.65, 3.62,
1012 2.28, 3.8, 12.0, 4.0, 2.0, 20.0, 0.0, 4.97, 28.01, 4.14, 3.0}},
1022static const std::map<std::string, std::map<int, std::array<float, 22>>>
NTNRuralNLOS{
1025 {10, {-9.01, 1.59, -2.9, 1.34, -3.33, 6.22, -0.88, 3.26, -4.92, 3.96, 0.0,
1026 0.0, 1.7, 7.0, 3.0, 3.0, 20.0, 0.0, 0.03, 18.16, 2.32, 3.0}},
1027 {20, {-8.37, 0.95, -2.5, 1.18, -0.74, 4.22, -0.07, 3.29, -4.06, 4.07, 0.0,
1028 0.0, 1.7, 7.0, 3.0, 3.0, 20.0, 0.0, 0.05, 26.82, 7.34, 3.0}},
1029 {30, {-8.05, 0.92, -2.12, 1.08, 0.08, 3.02, 0.75, 1.92, -2.33, 1.7, 0.0,
1030 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.07, 21.99, 8.28, 3.0}},
1031 {40, {-7.92, 0.92, -1.99, 1.06, 0.32, 2.45, 0.72, 1.92, -2.24, 2.01, 0.0,
1032 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.1, 22.86, 8.76, 3.0}},
1033 {50, {-7.92, 0.87, -1.9, 1.05, 0.53, 1.63, 0.95, 1.45, -2.24, 2.0, 0.0,
1034 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.15, 25.93, 9.68, 3.0}},
1035 {60, {-7.96, 0.87, -1.85, 1.06, 0.33, 2.08, 0.97, 1.62, -2.22, 1.82, 0.0,
1036 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.22, 27.79, 9.94, 3.0}},
1037 {70, {-7.91, 0.82, -1.69, 1.14, 0.55, 1.58, 1.1, 1.43, -2.19, 1.66, 0.0,
1038 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.5, 28.5, 8.9, 3.0}},
1039 {80, {-7.79, 0.86, -1.46, 1.16, 0.45, 2.01, 0.97, 1.88, -2.41, 2.58, 0.0,
1040 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 1.04, 37.53, 13.74, 3.0}},
1041 {90, {-7.74, 0.81, -1.32, 1.3, 0.4, 2.19, 1.35, 0.62, -2.45, 2.52, 0.0,
1042 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 2.11, 29.23, 12.16, 3.0}},
1046 {10, {-9.13, 1.91, -2.9, 1.32, -3.4, 6.28, -1.19, 3.81, -5.47, 4.39, 0.0,
1047 0.0, 1.7, 7.0, 3.0, 3.0, 20.0, 0.0, 0.03, 18.21, 2.13, 3.0}},
1048 {20, {-8.39, 0.94, -2.53, 1.18, -0.51, 3.75, -0.11, 3.33, -4.06, 4.04, 0.0,
1049 0.0, 1.7, 7.0, 3.0, 3.0, 20.0, 0.0, 0.05, 24.08, 6.52, 3.0}},
1050 {30, {-8.1, 0.92, -2.16, 1.08, 0.06, 2.95, 0.72, 1.93, -2.32, 1.54, 0.0,
1051 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.07, 22.06, 7.72, 3.0}},
1052 {40, {-7.96, 0.94, -2.04, 1.09, 0.2, 2.65, 0.69, 1.91, -2.19, 1.73, 0.0,
1053 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.09, 21.4, 8.45, 3.0}},
1054 {50, {-7.99, 0.89, -1.99, 1.08, 0.4, 1.85, 0.84, 1.7, -2.16, 1.5, 0.0,
1055 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.16, 24.26, 8.92, 3.0}},
1056 {60, {-8.05, 0.87, -1.95, 1.06, 0.32, 1.83, 0.99, 1.27, -2.24, 1.64, 0.0,
1057 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.22, 24.15, 8.76, 3.0}},
1058 {70, {-8.01, 0.82, -1.81, 1.17, 0.46, 1.57, 0.95, 1.86, -2.29, 1.66, 0.0,
1059 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.51, 25.99, 9.0, 3.0}},
1060 {80, {-8.05, 1.65, -1.56, 1.2, 0.33, 1.99, 0.92, 1.84, -2.65, 2.86, 0.0,
1061 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 0.89, 36.07, 13.6, 3.0}},
1062 {90, {-7.91, 0.76, -1.53, 1.27, 0.24, 2.18, 1.29, 0.59, -2.23, 1.12, 0.0,
1063 0.0, 1.7, 7.0, 3.0, 2.0, 20.0, 0.0, 1.68, 24.51, 10.56, 3.0}},
1225 double hUT = std::min(aMob->GetPosition().z, bMob->GetPosition().z);
1226 double hBS = std::max(aMob->GetPosition().z, bMob->GetPosition().z);
1228 double distance2D = sqrt(pow(aMob->GetPosition().x - bMob->GetPosition().x, 2) +
1229 pow(aMob->GetPosition().y - bMob->GetPosition().y, 2));
1238 bool los = channelCondition->IsLos();
1239 bool o2i = channelCondition->IsO2i();
1248 table3gpp->m_numOfCluster = 11;
1249 table3gpp->m_raysPerCluster = 20;
1250 table3gpp->m_uLgDS = -7.49;
1251 table3gpp->m_sigLgDS = 0.55;
1252 table3gpp->m_uLgASD = 0.90;
1253 table3gpp->m_sigLgASD = 0.38;
1254 table3gpp->m_uLgASA = 1.52;
1255 table3gpp->m_sigLgASA = 0.24;
1256 table3gpp->m_uLgZSA = 0.47;
1257 table3gpp->m_sigLgZSA = 0.40;
1258 table3gpp->m_uLgZSD = 0.34;
1259 table3gpp->m_sigLgZSD =
1260 std::max(-1.0, -0.17 * (distance2D / 1000.0) - 0.01 * (hUT - 1.5) + 0.22);
1261 table3gpp->m_offsetZOD = 0;
1262 table3gpp->m_cDS = 3.91e-9;
1263 table3gpp->m_cASD = 2;
1264 table3gpp->m_cASA = 3;
1265 table3gpp->m_cZSA = 3;
1266 table3gpp->m_uK = 7;
1267 table3gpp->m_sigK = 4;
1268 table3gpp->m_rTau = 3.8;
1269 table3gpp->m_uXpr = 12;
1270 table3gpp->m_sigXpr = 4;
1271 table3gpp->m_perClusterShadowingStd = 3;
1273 for (uint8_t row = 0; row < 7; row++)
1275 for (uint8_t column = 0; column < 7; column++)
1277 table3gpp->m_sqrtC[row][column] =
sqrtC_RMa_LOS[row][column];
1281 else if (!los && !o2i)
1283 table3gpp->m_numOfCluster = 10;
1284 table3gpp->m_raysPerCluster = 20;
1285 table3gpp->m_uLgDS = -7.43;
1286 table3gpp->m_sigLgDS = 0.48;
1287 table3gpp->m_uLgASD = 0.95;
1288 table3gpp->m_sigLgASD = 0.45;
1289 table3gpp->m_uLgASA = 1.52;
1290 table3gpp->m_sigLgASA = 0.13;
1291 table3gpp->m_uLgZSA = 0.58;
1292 table3gpp->m_sigLgZSA = 0.37;
1293 table3gpp->m_uLgZSD =
1294 std::max(-1.0, -0.19 * (distance2D / 1000.0) - 0.01 * (hUT - 1.5) + 0.28);
1295 table3gpp->m_sigLgZSD = 0.30;
1296 table3gpp->m_offsetZOD = atan((35 - 3.5) / distance2D) - atan((35 - 1.5) / distance2D);
1297 table3gpp->m_cDS = 3.91e-9;
1298 table3gpp->m_cASD = 2;
1299 table3gpp->m_cASA = 3;
1300 table3gpp->m_cZSA = 3;
1301 table3gpp->m_uK = 0;
1302 table3gpp->m_sigK = 0;
1303 table3gpp->m_rTau = 1.7;
1304 table3gpp->m_uXpr = 7;
1305 table3gpp->m_sigXpr = 3;
1306 table3gpp->m_perClusterShadowingStd = 3;
1308 for (uint8_t row = 0; row < 6; row++)
1310 for (uint8_t column = 0; column < 6; column++)
1318 table3gpp->m_numOfCluster = 10;
1319 table3gpp->m_raysPerCluster = 20;
1320 table3gpp->m_uLgDS = -7.47;
1321 table3gpp->m_sigLgDS = 0.24;
1322 table3gpp->m_uLgASD = 0.67;
1323 table3gpp->m_sigLgASD = 0.18;
1324 table3gpp->m_uLgASA = 1.66;
1325 table3gpp->m_sigLgASA = 0.21;
1326 table3gpp->m_uLgZSA = 0.93;
1327 table3gpp->m_sigLgZSA = 0.22;
1328 table3gpp->m_uLgZSD =
1329 std::max(-1.0, -0.19 * (distance2D / 1000.0) - 0.01 * (hUT - 1.5) + 0.28);
1330 table3gpp->m_sigLgZSD = 0.30;
1331 table3gpp->m_offsetZOD = atan((35 - 3.5) / distance2D) - atan((35 - 1.5) / distance2D);
1332 table3gpp->m_cDS = 3.91e-9;
1333 table3gpp->m_cASD = 2;
1334 table3gpp->m_cASA = 3;
1335 table3gpp->m_cZSA = 3;
1336 table3gpp->m_uK = 0;
1337 table3gpp->m_sigK = 0;
1338 table3gpp->m_rTau = 1.7;
1339 table3gpp->m_uXpr = 7;
1340 table3gpp->m_sigXpr = 3;
1341 table3gpp->m_perClusterShadowingStd = 3;
1343 for (uint8_t row = 0; row < 6; row++)
1345 for (uint8_t column = 0; column < 6; column++)
1347 table3gpp->m_sqrtC[row][column] =
sqrtC_RMa_O2I[row][column];
1356 table3gpp->m_numOfCluster = 12;
1357 table3gpp->m_raysPerCluster = 20;
1358 table3gpp->m_uLgDS = -6.955 - 0.0963 * log10(fcGHz);
1359 table3gpp->m_sigLgDS = 0.66;
1360 table3gpp->m_uLgASD = 1.06 + 0.1114 * log10(fcGHz);
1361 table3gpp->m_sigLgASD = 0.28;
1362 table3gpp->m_uLgASA = 1.81;
1363 table3gpp->m_sigLgASA = 0.20;
1364 table3gpp->m_uLgZSA = 0.95;
1365 table3gpp->m_sigLgZSA = 0.16;
1366 table3gpp->m_uLgZSD =
1367 std::max(-0.5, -2.1 * distance2D / 1000.0 - 0.01 * (hUT - 1.5) + 0.75);
1368 table3gpp->m_sigLgZSD = 0.40;
1369 table3gpp->m_offsetZOD = 0;
1370 table3gpp->m_cDS = std::max(0.25, -3.4084 * log10(fcGHz) + 6.5622) * 1e-9;
1371 table3gpp->m_cASD = 5;
1372 table3gpp->m_cASA = 11;
1373 table3gpp->m_cZSA = 7;
1374 table3gpp->m_uK = 9;
1375 table3gpp->m_sigK = 3.5;
1376 table3gpp->m_rTau = 2.5;
1377 table3gpp->m_uXpr = 8;
1378 table3gpp->m_sigXpr = 4;
1379 table3gpp->m_perClusterShadowingStd = 3;
1381 for (uint8_t row = 0; row < 7; row++)
1383 for (uint8_t column = 0; column < 7; column++)
1385 table3gpp->m_sqrtC[row][column] =
sqrtC_UMa_LOS[row][column];
1391 double uLgZSD = std::max(-0.5, -2.1 * distance2D / 1000.0 - 0.01 * (hUT - 1.5) + 0.9);
1393 double afc = 0.208 * log10(fcGHz) - 0.782;
1395 double cfc = -0.13 * log10(fcGHz) + 2.03;
1396 double efc = 7.66 * log10(fcGHz) - 5.96;
1398 double offsetZOD = efc - std::pow(10, afc * log10(std::max(bfc, distance2D)) + cfc);
1402 table3gpp->m_numOfCluster = 20;
1403 table3gpp->m_raysPerCluster = 20;
1404 table3gpp->m_uLgDS = -6.28 - 0.204 * log10(fcGHz);
1405 table3gpp->m_sigLgDS = 0.39;
1406 table3gpp->m_uLgASD = 1.5 - 0.1144 * log10(fcGHz);
1407 table3gpp->m_sigLgASD = 0.28;
1408 table3gpp->m_uLgASA = 2.08 - 0.27 * log10(fcGHz);
1409 table3gpp->m_sigLgASA = 0.11;
1410 table3gpp->m_uLgZSA = -0.3236 * log10(fcGHz) + 1.512;
1411 table3gpp->m_sigLgZSA = 0.16;
1412 table3gpp->m_uLgZSD =
uLgZSD;
1413 table3gpp->m_sigLgZSD = 0.49;
1414 table3gpp->m_offsetZOD = offsetZOD;
1415 table3gpp->m_cDS = std::max(0.25, -3.4084 * log10(fcGHz) + 6.5622) * 1e-9;
1416 table3gpp->m_cASD = 2;
1417 table3gpp->m_cASA = 15;
1418 table3gpp->m_cZSA = 7;
1419 table3gpp->m_uK = 0;
1420 table3gpp->m_sigK = 0;
1421 table3gpp->m_rTau = 2.3;
1422 table3gpp->m_uXpr = 7;
1423 table3gpp->m_sigXpr = 3;
1424 table3gpp->m_perClusterShadowingStd = 3;
1426 for (uint8_t row = 0; row < 6; row++)
1428 for (uint8_t column = 0; column < 6; column++)
1436 table3gpp->m_numOfCluster = 12;
1437 table3gpp->m_raysPerCluster = 20;
1438 table3gpp->m_uLgDS = -6.62;
1439 table3gpp->m_sigLgDS = 0.32;
1440 table3gpp->m_uLgASD = 1.25;
1441 table3gpp->m_sigLgASD = 0.42;
1442 table3gpp->m_uLgASA = 1.76;
1443 table3gpp->m_sigLgASA = 0.16;
1444 table3gpp->m_uLgZSA = 1.01;
1445 table3gpp->m_sigLgZSA = 0.43;
1446 table3gpp->m_uLgZSD =
uLgZSD;
1447 table3gpp->m_sigLgZSD = 0.49;
1448 table3gpp->m_offsetZOD = offsetZOD;
1449 table3gpp->m_cDS = 11e-9;
1450 table3gpp->m_cASD = 5;
1451 table3gpp->m_cASA = 8;
1452 table3gpp->m_cZSA = 3;
1453 table3gpp->m_uK = 0;
1454 table3gpp->m_sigK = 0;
1455 table3gpp->m_rTau = 2.2;
1456 table3gpp->m_uXpr = 9;
1457 table3gpp->m_sigXpr = 5;
1458 table3gpp->m_perClusterShadowingStd = 4;
1460 for (uint8_t row = 0; row < 6; row++)
1462 for (uint8_t column = 0; column < 6; column++)
1464 table3gpp->m_sqrtC[row][column] =
sqrtC_UMa_O2I[row][column];
1474 table3gpp->m_numOfCluster = 12;
1475 table3gpp->m_raysPerCluster = 20;
1476 table3gpp->m_uLgDS = -0.24 * log10(1 + fcGHz) - 7.14;
1477 table3gpp->m_sigLgDS = 0.38;
1478 table3gpp->m_uLgASD = -0.05 * log10(1 + fcGHz) + 1.21;
1479 table3gpp->m_sigLgASD = 0.41;
1480 table3gpp->m_uLgASA = -0.08 * log10(1 + fcGHz) + 1.73;
1481 table3gpp->m_sigLgASA = 0.014 * log10(1 + fcGHz) + 0.28;
1482 table3gpp->m_uLgZSA = -0.1 * log10(1 + fcGHz) + 0.73;
1483 table3gpp->m_sigLgZSA = -0.04 * log10(1 + fcGHz) + 0.34;
1484 table3gpp->m_uLgZSD =
1485 std::max(-0.21, -14.8 * distance2D / 1000.0 + 0.01 * std::abs(hUT - hBS) + 0.83);
1486 table3gpp->m_sigLgZSD = 0.35;
1487 table3gpp->m_offsetZOD = 0;
1488 table3gpp->m_cDS = 5e-9;
1489 table3gpp->m_cASD = 3;
1490 table3gpp->m_cASA = 17;
1491 table3gpp->m_cZSA = 7;
1492 table3gpp->m_uK = 9;
1493 table3gpp->m_sigK = 5;
1494 table3gpp->m_rTau = 3;
1495 table3gpp->m_uXpr = 9;
1496 table3gpp->m_sigXpr = 3;
1497 table3gpp->m_perClusterShadowingStd = 3;
1499 for (uint8_t row = 0; row < 7; row++)
1501 for (uint8_t column = 0; column < 7; column++)
1503 table3gpp->m_sqrtC[row][column] =
sqrtC_UMi_LOS[row][column];
1510 std::max(-0.5, -3.1 * distance2D / 1000.0 + 0.01 * std::max(hUT - hBS, 0.0) + 0.2);
1511 double offsetZOD = -1 * std::pow(10, -1.5 * log10(std::max(10.0, distance2D)) + 3.3);
1514 table3gpp->m_numOfCluster = 19;
1515 table3gpp->m_raysPerCluster = 20;
1516 table3gpp->m_uLgDS = -0.24 * log10(1 + fcGHz) - 6.83;
1517 table3gpp->m_sigLgDS = 0.16 * log10(1 + fcGHz) + 0.28;
1518 table3gpp->m_uLgASD = -0.23 * log10(1 + fcGHz) + 1.53;
1519 table3gpp->m_sigLgASD = 0.11 * log10(1 + fcGHz) + 0.33;
1520 table3gpp->m_uLgASA = -0.08 * log10(1 + fcGHz) + 1.81;
1521 table3gpp->m_sigLgASA = 0.05 * log10(1 + fcGHz) + 0.3;
1522 table3gpp->m_uLgZSA = -0.04 * log10(1 + fcGHz) + 0.92;
1523 table3gpp->m_sigLgZSA = -0.07 * log10(1 + fcGHz) + 0.41;
1524 table3gpp->m_uLgZSD =
uLgZSD;
1525 table3gpp->m_sigLgZSD = 0.35;
1526 table3gpp->m_offsetZOD = offsetZOD;
1527 table3gpp->m_cDS = 11e-9;
1528 table3gpp->m_cASD = 10;
1529 table3gpp->m_cASA = 22;
1530 table3gpp->m_cZSA = 7;
1531 table3gpp->m_uK = 0;
1532 table3gpp->m_sigK = 0;
1533 table3gpp->m_rTau = 2.1;
1534 table3gpp->m_uXpr = 8;
1535 table3gpp->m_sigXpr = 3;
1536 table3gpp->m_perClusterShadowingStd = 3;
1538 for (uint8_t row = 0; row < 6; row++)
1540 for (uint8_t column = 0; column < 6; column++)
1548 table3gpp->m_numOfCluster = 12;
1549 table3gpp->m_raysPerCluster = 20;
1550 table3gpp->m_uLgDS = -6.62;
1551 table3gpp->m_sigLgDS = 0.32;
1552 table3gpp->m_uLgASD = 1.25;
1553 table3gpp->m_sigLgASD = 0.42;
1554 table3gpp->m_uLgASA = 1.76;
1555 table3gpp->m_sigLgASA = 0.16;
1556 table3gpp->m_uLgZSA = 1.01;
1557 table3gpp->m_sigLgZSA = 0.43;
1558 table3gpp->m_uLgZSD =
uLgZSD;
1559 table3gpp->m_sigLgZSD = 0.35;
1560 table3gpp->m_offsetZOD = offsetZOD;
1561 table3gpp->m_cDS = 11e-9;
1562 table3gpp->m_cASD = 5;
1563 table3gpp->m_cASA = 8;
1564 table3gpp->m_cZSA = 3;
1565 table3gpp->m_uK = 0;
1566 table3gpp->m_sigK = 0;
1567 table3gpp->m_rTau = 2.2;
1568 table3gpp->m_uXpr = 9;
1569 table3gpp->m_sigXpr = 5;
1570 table3gpp->m_perClusterShadowingStd = 4;
1572 for (uint8_t row = 0; row < 6; row++)
1574 for (uint8_t column = 0; column < 6; column++)
1576 table3gpp->m_sqrtC[row][column] =
sqrtC_UMi_O2I[row][column];
1584 NS_ASSERT_MSG(!o2i,
"The indoor scenario does out support outdoor to indoor");
1587 table3gpp->m_numOfCluster = 15;
1588 table3gpp->m_raysPerCluster = 20;
1589 table3gpp->m_uLgDS = -0.01 * log10(1 + fcGHz) - 7.692;
1590 table3gpp->m_sigLgDS = 0.18;
1591 table3gpp->m_uLgASD = 1.60;
1592 table3gpp->m_sigLgASD = 0.18;
1593 table3gpp->m_uLgASA = -0.19 * log10(1 + fcGHz) + 1.781;
1594 table3gpp->m_sigLgASA = 0.12 * log10(1 + fcGHz) + 0.119;
1595 table3gpp->m_uLgZSA = -0.26 * log10(1 + fcGHz) + 1.44;
1596 table3gpp->m_sigLgZSA = -0.04 * log10(1 + fcGHz) + 0.264;
1597 table3gpp->m_uLgZSD = -1.43 * log10(1 + fcGHz) + 2.228;
1598 table3gpp->m_sigLgZSD = 0.13 * log10(1 + fcGHz) + 0.30;
1599 table3gpp->m_offsetZOD = 0;
1600 table3gpp->m_cDS = 3.91e-9;
1601 table3gpp->m_cASD = 5;
1602 table3gpp->m_cASA = 8;
1603 table3gpp->m_cZSA = 9;
1604 table3gpp->m_uK = 7;
1605 table3gpp->m_sigK = 4;
1606 table3gpp->m_rTau = 3.6;
1607 table3gpp->m_uXpr = 11;
1608 table3gpp->m_sigXpr = 4;
1609 table3gpp->m_perClusterShadowingStd = 6;
1611 for (uint8_t row = 0; row < 7; row++)
1613 for (uint8_t column = 0; column < 7; column++)
1621 table3gpp->m_numOfCluster = 19;
1622 table3gpp->m_raysPerCluster = 20;
1623 table3gpp->m_uLgDS = -0.28 * log10(1 + fcGHz) - 7.173;
1624 table3gpp->m_sigLgDS = 0.1 * log10(1 + fcGHz) + 0.055;
1625 table3gpp->m_uLgASD = 1.62;
1626 table3gpp->m_sigLgASD = 0.25;
1627 table3gpp->m_uLgASA = -0.11 * log10(1 + fcGHz) + 1.863;
1628 table3gpp->m_sigLgASA = 0.12 * log10(1 + fcGHz) + 0.059;
1629 table3gpp->m_uLgZSA = -0.15 * log10(1 + fcGHz) + 1.387;
1630 table3gpp->m_sigLgZSA = -0.09 * log10(1 + fcGHz) + 0.746;
1631 table3gpp->m_uLgZSD = 1.08;
1632 table3gpp->m_sigLgZSD = 0.36;
1633 table3gpp->m_offsetZOD = 0;
1634 table3gpp->m_cDS = 3.91e-9;
1635 table3gpp->m_cASD = 5;
1636 table3gpp->m_cASA = 11;
1637 table3gpp->m_cZSA = 9;
1638 table3gpp->m_uK = 0;
1639 table3gpp->m_sigK = 0;
1640 table3gpp->m_rTau = 3;
1641 table3gpp->m_uXpr = 10;
1642 table3gpp->m_sigXpr = 4;
1643 table3gpp->m_perClusterShadowingStd = 3;
1645 for (uint8_t row = 0; row < 6; row++)
1647 for (uint8_t column = 0; column < 6; column++)
1656 if (channelCondition->IsLos())
1660 table3gpp->m_numOfCluster = 12;
1661 table3gpp->m_raysPerCluster = 20;
1662 table3gpp->m_uLgDS = -0.2 * log10(1 + fcGHz) - 7.5;
1663 table3gpp->m_sigLgDS = 0.1;
1664 table3gpp->m_uLgASD = -0.1 * log10(1 + fcGHz) + 1.6;
1665 table3gpp->m_sigLgASD = 0.1;
1666 table3gpp->m_uLgASA = -0.1 * log10(1 + fcGHz) + 1.6;
1667 table3gpp->m_sigLgASA = 0.1;
1668 table3gpp->m_uLgZSA = -0.1 * log10(1 + fcGHz) + 0.73;
1669 table3gpp->m_sigLgZSA = -0.04 * log10(1 + fcGHz) + 0.34;
1670 table3gpp->m_uLgZSD = -0.1 * log10(1 + fcGHz) + 0.73;
1671 table3gpp->m_sigLgZSD = -0.04 * log10(1 + fcGHz) + 0.34;
1672 table3gpp->m_offsetZOD = 0;
1673 table3gpp->m_cDS = 5;
1674 table3gpp->m_cASD = 17;
1675 table3gpp->m_cASA = 17;
1676 table3gpp->m_cZSA = 7;
1677 table3gpp->m_uK = 3.48;
1678 table3gpp->m_sigK = 2;
1679 table3gpp->m_rTau = 3;
1680 table3gpp->m_uXpr = 9;
1681 table3gpp->m_sigXpr = 3;
1682 table3gpp->m_perClusterShadowingStd = 4;
1684 for (uint8_t row = 0; row < 7; row++)
1686 for (uint8_t column = 0; column < 7; column++)
1688 table3gpp->m_sqrtC[row][column] =
sqrtC_UMi_LOS[row][column];
1692 else if (channelCondition->IsNlos())
1694 table3gpp->m_numOfCluster = 19;
1695 table3gpp->m_raysPerCluster = 20;
1696 table3gpp->m_uLgDS = -0.3 * log10(1 + fcGHz) - 7;
1697 table3gpp->m_sigLgDS = 0.28;
1698 table3gpp->m_uLgASD = -0.08 * log10(1 + fcGHz) + 1.81;
1699 table3gpp->m_sigLgASD = 0.05 * log10(1 + fcGHz) + 0.3;
1700 table3gpp->m_uLgASA = -0.08 * log10(1 + fcGHz) + 1.81;
1701 table3gpp->m_sigLgASA = 0.05 * log10(1 + fcGHz) + 0.3;
1702 table3gpp->m_uLgZSA = -0.04 * log10(1 + fcGHz) + 0.92;
1703 table3gpp->m_sigLgZSA = -0.07 * log10(1 + fcGHz) + 0.41;
1704 table3gpp->m_uLgZSD = -0.04 * log10(1 + fcGHz) + 0.92;
1705 table3gpp->m_sigLgZSD = -0.07 * log10(1 + fcGHz) + 0.41;
1706 table3gpp->m_offsetZOD = 0;
1707 table3gpp->m_cDS = 11;
1708 table3gpp->m_cASD = 22;
1709 table3gpp->m_cASA = 22;
1710 table3gpp->m_cZSA = 7;
1711 table3gpp->m_uK = 0;
1712 table3gpp->m_sigK = 0;
1713 table3gpp->m_rTau = 2.1;
1714 table3gpp->m_uXpr = 8;
1715 table3gpp->m_sigXpr = 3;
1716 table3gpp->m_perClusterShadowingStd = 4;
1718 for (uint8_t row = 0; row < 6; row++)
1720 for (uint8_t column = 0; column < 6; column++)
1726 else if (channelCondition->IsNlosv())
1728 table3gpp->m_numOfCluster = 19;
1729 table3gpp->m_raysPerCluster = 20;
1730 table3gpp->m_uLgDS = -0.4 * log10(1 + fcGHz) - 7;
1731 table3gpp->m_sigLgDS = 0.1;
1732 table3gpp->m_uLgASD = -0.1 * log10(1 + fcGHz) + 1.7;
1733 table3gpp->m_sigLgASD = 0.1;
1734 table3gpp->m_uLgASA = -0.1 * log10(1 + fcGHz) + 1.7;
1735 table3gpp->m_sigLgASA = 0.1;
1736 table3gpp->m_uLgZSA = -0.04 * log10(1 + fcGHz) + 0.92;
1737 table3gpp->m_sigLgZSA = -0.07 * log10(1 + fcGHz) + 0.41;
1738 table3gpp->m_uLgZSD = -0.04 * log10(1 + fcGHz) + 0.92;
1739 table3gpp->m_sigLgZSD = -0.07 * log10(1 + fcGHz) + 0.41;
1740 table3gpp->m_offsetZOD = 0;
1741 table3gpp->m_cDS = 11;
1742 table3gpp->m_cASD = 22;
1743 table3gpp->m_cASA = 22;
1744 table3gpp->m_cZSA = 7;
1745 table3gpp->m_uK = 0;
1746 table3gpp->m_sigK = 4.5;
1747 table3gpp->m_rTau = 2.1;
1748 table3gpp->m_uXpr = 8;
1749 table3gpp->m_sigXpr = 3;
1750 table3gpp->m_perClusterShadowingStd = 4;
1752 for (uint8_t row = 0; row < 6; row++)
1754 for (uint8_t column = 0; column < 6; column++)
1756 table3gpp->m_sqrtC[row][column] =
sqrtC_UMi_LOS[row][column];
1767 if (channelCondition->IsLos())
1769 table3gpp->m_numOfCluster = 12;
1770 table3gpp->m_raysPerCluster = 20;
1771 table3gpp->m_uLgDS = -8.3;
1772 table3gpp->m_sigLgDS = 0.2;
1773 table3gpp->m_uLgASD = 1.4;
1774 table3gpp->m_sigLgASD = 0.1;
1775 table3gpp->m_uLgASA = 1.4;
1776 table3gpp->m_sigLgASA = 0.1;
1777 table3gpp->m_uLgZSA = -0.1 * log10(1 + fcGHz) + 0.73;
1778 table3gpp->m_sigLgZSA = -0.04 * log10(1 + fcGHz) + 0.34;
1779 table3gpp->m_uLgZSD = -0.1 * log10(1 + fcGHz) + 0.73;
1780 table3gpp->m_sigLgZSD = -0.04 * log10(1 + fcGHz) + 0.34;
1781 table3gpp->m_offsetZOD = 0;
1782 table3gpp->m_cDS = 5;
1783 table3gpp->m_cASD = 17;
1784 table3gpp->m_cASA = 17;
1785 table3gpp->m_cZSA = 7;
1786 table3gpp->m_uK = 9;
1787 table3gpp->m_sigK = 3.5;
1788 table3gpp->m_rTau = 3;
1789 table3gpp->m_uXpr = 9;
1790 table3gpp->m_sigXpr = 3;
1791 table3gpp->m_perClusterShadowingStd = 4;
1793 for (uint8_t row = 0; row < 7; row++)
1795 for (uint8_t column = 0; column < 7; column++)
1797 table3gpp->m_sqrtC[row][column] =
sqrtC_UMi_LOS[row][column];
1801 else if (channelCondition->IsNlosv())
1803 table3gpp->m_numOfCluster = 19;
1804 table3gpp->m_raysPerCluster = 20;
1805 table3gpp->m_uLgDS = -8.3;
1806 table3gpp->m_sigLgDS = 0.3;
1807 table3gpp->m_uLgASD = 1.5;
1808 table3gpp->m_sigLgASD = 0.1;
1809 table3gpp->m_uLgASA = 1.5;
1810 table3gpp->m_sigLgASA = 0.1;
1811 table3gpp->m_uLgZSA = -0.04 * log10(1 + fcGHz) + 0.92;
1812 table3gpp->m_sigLgZSA = -0.07 * log10(1 + fcGHz) + 0.41;
1813 table3gpp->m_uLgZSD = -0.04 * log10(1 + fcGHz) + 0.92;
1814 table3gpp->m_sigLgZSD = -0.07 * log10(1 + fcGHz) + 0.41;
1815 table3gpp->m_offsetZOD = 0;
1816 table3gpp->m_cDS = 11;
1817 table3gpp->m_cASD = 22;
1818 table3gpp->m_cASA = 22;
1819 table3gpp->m_cZSA = 7;
1820 table3gpp->m_uK = 0;
1821 table3gpp->m_sigK = 4.5;
1822 table3gpp->m_rTau = 2.1;
1823 table3gpp->m_uXpr = 8.0;
1824 table3gpp->m_sigXpr = 3;
1825 table3gpp->m_perClusterShadowingStd = 4;
1827 for (uint8_t row = 0; row < 6; row++)
1829 for (uint8_t column = 0; column < 6; column++)
1831 table3gpp->m_sqrtC[row][column] =
sqrtC_UMi_LOS[row][column];
1835 else if (channelCondition->IsNlos())
1838 "The fast fading parameters for the NLOS condition in the Highway scenario are not "
1839 "defined in TR 37.885, use the ones defined in TDoc R1-1803671 instead");
1841 table3gpp->m_numOfCluster = 19;
1842 table3gpp->m_raysPerCluster = 20;
1843 table3gpp->m_uLgDS = -0.3 * log10(1 + fcGHz) - 7;
1844 table3gpp->m_sigLgDS = 0.28;
1845 table3gpp->m_uLgASD = -0.08 * log10(1 + fcGHz) + 1.81;
1846 table3gpp->m_sigLgASD = 0.05 * log10(1 + fcGHz) + 0.3;
1847 table3gpp->m_uLgASA = -0.08 * log10(1 + fcGHz) + 1.81;
1848 table3gpp->m_sigLgASA = 0.05 * log10(1 + fcGHz) + 0.3;
1849 table3gpp->m_uLgZSA = -0.04 * log10(1 + fcGHz) + 0.92;
1850 table3gpp->m_sigLgZSA = -0.07 * log10(1 + fcGHz) + 0.41;
1851 table3gpp->m_uLgZSD = -0.04 * log10(1 + fcGHz) + 0.92;
1852 table3gpp->m_sigLgZSD = -0.07 * log10(1 + fcGHz) + 0.41;
1853 table3gpp->m_offsetZOD = 0;
1854 table3gpp->m_cDS = 11;
1855 table3gpp->m_cASD = 22;
1856 table3gpp->m_cASA = 22;
1857 table3gpp->m_cZSA = 7;
1858 table3gpp->m_uK = 0;
1859 table3gpp->m_sigK = 0;
1860 table3gpp->m_rTau = 2.1;
1861 table3gpp->m_uXpr = 8;
1862 table3gpp->m_sigXpr = 3;
1863 table3gpp->m_perClusterShadowingStd = 4;
1865 for (uint8_t row = 0; row < 6; row++)
1867 for (uint8_t column = 0; column < 6; column++)
1880 std::string freqBand = (fcGHz < 13) ?
"S" :
"Ka";
1882 double elevAngle = 0;
1883 bool isSatellite =
false;
1901 if (aNTNMob->GetGeographicPosition().z <
1902 bNTNMob->GetGeographicPosition().z)
1904 elevAngle = aNTNMob->GetElevationAngle(bNTNMob);
1905 if (bNTNMob->GetGeographicPosition().z > 50000)
1912 elevAngle = bNTNMob->GetElevationAngle(aNTNMob);
1913 if (aNTNMob->GetGeographicPosition().z > 50000)
1921 NS_FATAL_ERROR(
"Mobility Models needs to be of type Geocentric for NTN scenarios");
1925 int elevAngleQuantized = (elevAngle < 10) ? 10 : round(elevAngle / 10) * 10;
1929 if (channelCondition->IsLos())
1931 table3gpp->m_uLgDS =
1933 table3gpp->m_sigLgDS =
1937 table3gpp->m_uLgASD =
1940 table3gpp->m_sigLgASD =
1943 table3gpp->m_uLgASA =
1945 table3gpp->m_sigLgASA =
1947 table3gpp->m_uLgZSA =
1949 table3gpp->m_sigLgZSA =
1953 table3gpp->m_uLgZSD =
1956 table3gpp->m_sigLgZSD =
1967 table3gpp->m_sigXpr =
1984 for (uint8_t row = 0; row < 7; row++)
1986 for (uint8_t column = 0; column < 7; column++)
1992 else if (channelCondition->IsNlos())
1995 table3gpp->m_uLgDS =
1997 table3gpp->m_sigLgDS =
1999 table3gpp->m_uLgASD =
2003 table3gpp->m_uLgASA =
2007 table3gpp->m_uLgZSA =
2011 table3gpp->m_uLgZSD =
2019 table3gpp->m_sigXpr =
2036 for (uint8_t row = 0; row < 6; row++)
2038 for (uint8_t column = 0; column < 6; column++)
2047 if (channelCondition->IsLos())
2049 table3gpp->m_uLgDS =
2051 table3gpp->m_sigLgDS =
2053 table3gpp->m_uLgASD =
2055 table3gpp->m_sigLgASD =
2057 table3gpp->m_uLgASA =
2059 table3gpp->m_sigLgASA =
2061 table3gpp->m_uLgZSA =
2063 table3gpp->m_sigLgZSA =
2065 table3gpp->m_uLgZSD =
2067 table3gpp->m_sigLgZSD =
2077 table3gpp->m_sigXpr =
2079 table3gpp->m_numOfCluster =
2081 table3gpp->m_raysPerCluster =
NTNUrbanLOS.at(freqBand).at(
2091 table3gpp->m_perClusterShadowingStd =
NTNUrbanLOS.at(freqBand).at(
2094 for (uint8_t row = 0; row < 7; row++)
2096 for (uint8_t column = 0; column < 7; column++)
2102 else if (channelCondition->IsNlos())
2104 table3gpp->m_uLgDS =
2106 table3gpp->m_sigLgDS =
2108 table3gpp->m_uLgASD =
2110 table3gpp->m_sigLgASD =
2112 table3gpp->m_uLgASA =
2114 table3gpp->m_sigLgASA =
2116 table3gpp->m_uLgZSA =
2118 table3gpp->m_sigLgZSA =
2120 table3gpp->m_uLgZSD =
2122 table3gpp->m_sigLgZSD =
2132 table3gpp->m_sigXpr =
2134 table3gpp->m_numOfCluster =
2136 table3gpp->m_raysPerCluster =
NTNUrbanNLOS.at(freqBand).at(
2146 table3gpp->m_perClusterShadowingStd =
NTNUrbanNLOS.at(freqBand).at(
2149 for (uint8_t row = 0; row < 6; row++)
2151 for (uint8_t column = 0; column < 6; column++)
2153 table3gpp->m_sqrtC[row][column] =
2161 if (channelCondition->IsLos())
2163 table3gpp->m_uLgDS =
2165 table3gpp->m_sigLgDS =
2167 table3gpp->m_uLgASD =
2169 table3gpp->m_sigLgASD =
2171 table3gpp->m_uLgASA =
2173 table3gpp->m_sigLgASA =
2175 table3gpp->m_uLgZSA =
2177 table3gpp->m_sigLgZSA =
2179 table3gpp->m_uLgZSD =
2181 table3gpp->m_sigLgZSD =
2191 table3gpp->m_sigXpr =
2205 table3gpp->m_perClusterShadowingStd =
NTNSuburbanLOS.at(freqBand).at(
2208 for (uint8_t row = 0; row < 7; row++)
2210 for (uint8_t column = 0; column < 7; column++)
2216 else if (channelCondition->IsNlos())
2218 table3gpp->m_uLgDS =
2220 table3gpp->m_sigLgDS =
2222 table3gpp->m_uLgASD =
2224 table3gpp->m_sigLgASD =
2226 table3gpp->m_uLgASA =
2228 table3gpp->m_sigLgASA =
2230 table3gpp->m_uLgZSA =
2232 table3gpp->m_sigLgZSA =
2234 table3gpp->m_uLgZSD =
2236 table3gpp->m_sigLgZSD =
2246 table3gpp->m_sigXpr =
2260 table3gpp->m_perClusterShadowingStd =
NTNSuburbanNLOS.at(freqBand).at(
2263 for (uint8_t row = 0; row < 6; row++)
2265 for (uint8_t column = 0; column < 6; column++)
2274 if (channelCondition->IsLos())
2276 table3gpp->m_uLgDS =
2278 table3gpp->m_sigLgDS =
2280 table3gpp->m_uLgASD =
2282 table3gpp->m_sigLgASD =
2284 table3gpp->m_uLgASA =
2286 table3gpp->m_sigLgASA =
2288 table3gpp->m_uLgZSA =
2290 table3gpp->m_sigLgZSA =
2292 table3gpp->m_uLgZSD =
2294 table3gpp->m_sigLgZSD =
2304 table3gpp->m_sigXpr =
2306 table3gpp->m_numOfCluster =
2308 table3gpp->m_raysPerCluster =
NTNRuralLOS.at(freqBand).at(
2318 table3gpp->m_perClusterShadowingStd =
NTNRuralLOS.at(freqBand).at(
2321 for (uint8_t row = 0; row < 7; row++)
2323 for (uint8_t column = 0; column < 7; column++)
2329 else if (channelCondition->IsNlos())
2331 table3gpp->m_uLgDS =
2333 table3gpp->m_sigLgDS =
2335 table3gpp->m_uLgASD =
2337 table3gpp->m_sigLgASD =
2339 table3gpp->m_uLgASA =
2341 table3gpp->m_sigLgASA =
2343 table3gpp->m_uLgZSA =
2345 table3gpp->m_sigLgZSA =
2347 table3gpp->m_uLgZSD =
2349 table3gpp->m_sigLgZSD =
2359 table3gpp->m_sigXpr =
2361 table3gpp->m_numOfCluster =
2363 table3gpp->m_raysPerCluster =
NTNRuralNLOS.at(freqBand).at(
2373 table3gpp->m_perClusterShadowingStd =
NTNRuralNLOS.at(freqBand).at(
2376 if (freqBand ==
"S")
2378 for (uint8_t row = 0; row < 6; row++)
2380 for (uint8_t column = 0; column < 6; column++)
2382 table3gpp->m_sqrtC[row][column] =
2387 else if (freqBand ==
"Ka")
2389 for (uint8_t row = 0; row < 6; row++)
2391 for (uint8_t column = 0; column < 6; column++)
2393 table3gpp->m_sqrtC[row][column] =
2404 table3gpp->m_uLgASD = std::numeric_limits<double>::min();
2405 table3gpp->m_sigLgASD = 0;
2406 table3gpp->m_uLgZSD = std::numeric_limits<double>::min();
2407 table3gpp->m_sigLgZSD = 0;
2586 channelParams->m_nodeIds =
2588 channelParams->m_losCondition = channelCondition->GetLosCondition();
2589 channelParams->m_o2iCondition = channelCondition->GetO2iCondition();
2594 uint8_t paramNum = 6;
2601 for (uint8_t iter = 0; iter < paramNum; iter++)
2605 for (uint8_t row = 0; row < paramNum; row++)
2608 for (uint8_t column = 0; column < paramNum; column++)
2610 temp += table3gpp->m_sqrtC[row][column] * LSPsIndep[column];
2612 LSPs.push_back(temp);
2624 kFactor = LSPs[1] * table3gpp->m_sigK + table3gpp->m_uK;
2625 DS = pow(10, LSPs[2] * table3gpp->m_sigLgDS + table3gpp->m_uLgDS);
2626 ASD = pow(10, LSPs[3] * table3gpp->m_sigLgASD + table3gpp->m_uLgASD);
2627 ASA = pow(10, LSPs[4] * table3gpp->m_sigLgASA + table3gpp->m_uLgASA);
2628 ZSD = pow(10, LSPs[5] * table3gpp->m_sigLgZSD + table3gpp->m_uLgZSD);
2629 ZSA = pow(10, LSPs[6] * table3gpp->m_sigLgZSA + table3gpp->m_uLgZSA);
2633 DS = pow(10, LSPs[1] * table3gpp->m_sigLgDS + table3gpp->m_uLgDS);
2634 ASD = pow(10, LSPs[2] * table3gpp->m_sigLgASD + table3gpp->m_uLgASD);
2635 ASA = pow(10, LSPs[3] * table3gpp->m_sigLgASA + table3gpp->m_uLgASA);
2636 ZSD = pow(10, LSPs[4] * table3gpp->m_sigLgZSD + table3gpp->m_uLgZSD);
2637 ZSA = pow(10, LSPs[5] * table3gpp->m_sigLgZSA + table3gpp->m_uLgZSA);
2639 ASD = std::min(ASD, 104.0);
2640 ASA = std::min(ASA, 104.0);
2641 ZSD = std::min(ZSD, 52.0);
2642 ZSA = std::min(ZSA, 52.0);
2645 channelParams->m_DS = DS;
2646 channelParams->m_K_factor = kFactor;
2648 NS_LOG_INFO(
"K-factor=" << kFactor <<
", DS=" << DS <<
", ASD=" << ASD <<
", ASA=" << ASA
2649 <<
", ZSD=" << ZSD <<
", ZSA=" << ZSA);
2653 double minTau = 100.0;
2654 for (uint8_t cIndex = 0; cIndex < table3gpp->m_numOfCluster; cIndex++)
2656 double tau = -1 * table3gpp->m_rTau * DS * log(
m_uniformRv->GetValue(0, 1));
2661 clusterDelay.push_back(tau);
2664 for (uint8_t cIndex = 0; cIndex < table3gpp->m_numOfCluster; cIndex++)
2666 clusterDelay[cIndex] -= minTau;
2668 std::sort(clusterDelay.begin(), clusterDelay.end());
2675 double powerSum = 0;
2676 for (uint8_t cIndex = 0; cIndex < table3gpp->m_numOfCluster; cIndex++)
2679 exp(-1 * clusterDelay[cIndex] * (table3gpp->m_rTau - 1) / table3gpp->m_rTau / DS) *
2681 -1 *
m_normalRv->GetValue() * table3gpp->m_perClusterShadowingStd / 10.0);
2683 clusterPower.push_back(power);
2685 channelParams->m_clusterPower = clusterPower;
2687 double powerMax = 0;
2689 for (uint8_t cIndex = 0; cIndex < table3gpp->m_numOfCluster; cIndex++)
2691 channelParams->m_clusterPower[cIndex] =
2692 channelParams->m_clusterPower[cIndex] / powerSum;
2699 double kLinear = pow(10, kFactor / 10.0);
2701 for (uint8_t cIndex = 0; cIndex < table3gpp->m_numOfCluster; cIndex++)
2705 clusterPowerForAngles.push_back(channelParams->m_clusterPower[cIndex] /
2707 kLinear / (1 + kLinear));
2711 clusterPowerForAngles.push_back(channelParams->m_clusterPower[cIndex] /
2714 if (powerMax < clusterPowerForAngles[cIndex])
2716 powerMax = clusterPowerForAngles[cIndex];
2722 for (uint8_t cIndex = 0; cIndex < table3gpp->m_numOfCluster; cIndex++)
2724 clusterPowerForAngles.push_back(channelParams->m_clusterPower[cIndex]);
2725 if (powerMax < clusterPowerForAngles[cIndex])
2727 powerMax = clusterPowerForAngles[cIndex];
2734 double thresh = 0.0032;
2735 for (uint8_t cIndex = table3gpp->m_numOfCluster; cIndex > 0; cIndex--)
2737 if (clusterPowerForAngles[cIndex - 1] < thresh * powerMax)
2739 clusterPowerForAngles.erase(clusterPowerForAngles.begin() + cIndex - 1);
2740 channelParams->m_clusterPower.erase(channelParams->m_clusterPower.begin() + cIndex - 1);
2741 clusterDelay.erase(clusterDelay.begin() + cIndex - 1);
2745 NS_ASSERT(channelParams->m_clusterPower.size() < UINT8_MAX);
2746 channelParams->m_reducedClusterNumber = channelParams->m_clusterPower.size();
2751 0.7705 - 0.0433 * kFactor + 2e-4 * pow(kFactor, 2) + 17e-6 * pow(kFactor, 3);
2752 for (uint8_t cIndex = 0; cIndex < channelParams->m_reducedClusterNumber; cIndex++)
2754 clusterDelay[cIndex] = clusterDelay[cIndex] / cTau;
2764 switch (table3gpp->m_numOfCluster)
2809 double cPhi = cNlos;
2813 cPhi *= (1.1035 - 0.028 * kFactor - 2e-3 * pow(kFactor, 2) +
2814 1e-4 * pow(kFactor, 3));
2818 switch (table3gpp->m_numOfCluster)
2854 double cTheta = cNlos;
2855 if (channelCondition->IsLos())
2857 cTheta *= (1.3086 + 0.0339 * kFactor - 0.0077 * pow(kFactor, 2) +
2858 2e-4 * pow(kFactor, 3));
2865 for (uint8_t cIndex = 0; cIndex < channelParams->m_reducedClusterNumber; cIndex++)
2867 double logCalc = -1 * log(clusterPowerForAngles[cIndex] / powerMax);
2868 double angle = 2 * sqrt(logCalc) / 1.4 / cPhi;
2869 clusterAoa.push_back(ASA * angle);
2870 clusterAod.push_back(ASD * angle);
2871 angle = logCalc / cTheta;
2872 clusterZoa.push_back(ZSA * angle);
2873 clusterZod.push_back(ZSD * angle);
2876 Angles sAngle(bMob->GetPosition(), aMob->GetPosition());
2877 Angles uAngle(aMob->GetPosition(), bMob->GetPosition());
2879 for (uint8_t cIndex = 0; cIndex < channelParams->m_reducedClusterNumber; cIndex++)
2886 clusterAoa[cIndex] = clusterAoa[cIndex] * Xn + (
m_normalRv->GetValue() * ASA / 7.0) +
2888 clusterAod[cIndex] = clusterAod[cIndex] * Xn + (
m_normalRv->GetValue() * ASD / 7.0) +
2890 if (channelCondition->IsO2i())
2892 clusterZoa[cIndex] =
2893 clusterZoa[cIndex] * Xn + (
m_normalRv->GetValue() * ZSA / 7.0) + 90;
2897 clusterZoa[cIndex] = clusterZoa[cIndex] * Xn + (
m_normalRv->GetValue() * ZSA / 7.0) +
2900 clusterZod[cIndex] = clusterZod[cIndex] * Xn + (
m_normalRv->GetValue() * ZSD / 7.0) +
2902 table3gpp->m_offsetZOD;
2914 for (uint8_t cIndex = 0; cIndex < channelParams->m_reducedClusterNumber; cIndex++)
2916 clusterAoa[cIndex] -= diffAoa;
2917 clusterAod[cIndex] -= diffAod;
2918 clusterZoa[cIndex] -= diffZsa;
2919 clusterZod[cIndex] -= diffZsd;
2923 double sizeTemp = clusterZoa.size();
2924 for (uint8_t ind = 0; ind < 4; ind++)
2930 angleDegree = clusterAoa;
2933 angleDegree = clusterZoa;
2936 angleDegree = clusterAod;
2939 angleDegree = clusterZod;
2944 for (uint8_t nIndex = 0; nIndex < sizeTemp; nIndex++)
2946 while (angleDegree[nIndex] > 360)
2948 angleDegree[nIndex] -= 360;
2951 while (angleDegree[nIndex] < 0)
2953 angleDegree[nIndex] += 360;
2956 if (ind == 1 || ind == 3)
2958 if (angleDegree[nIndex] > 180)
2960 angleDegree[nIndex] = 360 - angleDegree[nIndex];
2967 clusterAoa = angleDegree;
2970 clusterZoa = angleDegree;
2973 clusterAod = angleDegree;
2976 clusterZod = angleDegree;
2987 for (uint8_t cInd = 0; cInd < channelParams->m_reducedClusterNumber; cInd++)
2989 channelParams->m_clusterPower[cInd] =
2990 channelParams->m_clusterPower[cInd] / pow(10, attenuationDb[cInd] / 10.0);
2995 attenuationDb.push_back(0);
2999 channelParams->m_attenuation_dB = attenuationDb;
3004 channelParams->m_reducedClusterNumber,
3008 channelParams->m_reducedClusterNumber,
3012 channelParams->m_reducedClusterNumber,
3016 channelParams->m_reducedClusterNumber,
3020 const double pow10_uLgZSD = pow(10, table3gpp->m_uLgZSD);
3021 for (uint8_t nInd = 0; nInd < channelParams->m_reducedClusterNumber; nInd++)
3023 for (uint8_t mInd = 0; mInd < table3gpp->m_raysPerCluster; mInd++)
3025 double tempAoa = clusterAoa[nInd] + table3gpp->m_cASA *
offSetAlpha[mInd];
3026 double tempZoa = clusterZoa[nInd] + table3gpp->m_cZSA *
offSetAlpha[mInd];
3027 std::tie(rayAoaRadian[nInd][mInd], rayZoaRadian[nInd][mInd]) =
3030 double tempAod = clusterAod[nInd] + table3gpp->m_cASD *
offSetAlpha[mInd];
3031 double tempZod = clusterZod[nInd] + 0.375 * pow10_uLgZSD *
offSetAlpha[mInd];
3032 std::tie(rayAodRadian[nInd][mInd], rayZodRadian[nInd][mInd]) =
3037 for (uint8_t cIndex = 0; cIndex < channelParams->m_reducedClusterNumber; cIndex++)
3039 Shuffle(&rayAodRadian[cIndex][0], &rayAodRadian[cIndex][table3gpp->m_raysPerCluster]);
3040 Shuffle(&rayAoaRadian[cIndex][0], &rayAoaRadian[cIndex][table3gpp->m_raysPerCluster]);
3041 Shuffle(&rayZodRadian[cIndex][0], &rayZodRadian[cIndex][table3gpp->m_raysPerCluster]);
3042 Shuffle(&rayZoaRadian[cIndex][0], &rayZoaRadian[cIndex][table3gpp->m_raysPerCluster]);
3046 channelParams->m_rayAodRadian = rayAodRadian;
3047 channelParams->m_rayAoaRadian = rayAoaRadian;
3048 channelParams->m_rayZodRadian = rayZodRadian;
3049 channelParams->m_rayZoaRadian = rayZoaRadian;
3055 auto& crossPolarizationPowerRatios = channelParams->m_crossPolarizationPowerRatios;
3057 auto& clusterPhase = channelParams->m_clusterPhase;
3059 const double uXprLinear = pow(10, table3gpp->m_uXpr / 10.0);
3060 const double sigXprLinear = pow(10, table3gpp->m_sigXpr / 10.0);
3063 clusterPhase.resize(channelParams->m_reducedClusterNumber);
3064 crossPolarizationPowerRatios.resize(channelParams->m_reducedClusterNumber);
3065 for (uint8_t nInd = 0; nInd < channelParams->m_reducedClusterNumber; nInd++)
3067 clusterPhase[nInd].resize(table3gpp->m_raysPerCluster);
3068 crossPolarizationPowerRatios[nInd].resize(table3gpp->m_raysPerCluster);
3069 for (uint8_t mInd = 0; mInd < table3gpp->m_raysPerCluster; mInd++)
3071 clusterPhase[nInd][mInd].resize(4);
3073 crossPolarizationPowerRatios[nInd][mInd] =
3074 std::pow(10, (
m_normalRv->GetValue() * sigXprLinear + uXprLinear) / 10.0);
3075 for (uint8_t pInd = 0; pInd < 4; pInd++)
3078 clusterPhase[nInd][mInd][pInd] =
m_uniformRv->GetValue(-1 * M_PI, M_PI);
3083 uint8_t cluster1st = 0;
3084 uint8_t cluster2nd = 0;
3085 double maxPower = 0;
3086 for (uint8_t cIndex = 0; cIndex < channelParams->m_reducedClusterNumber; cIndex++)
3088 if (maxPower < channelParams->m_clusterPower[cIndex])
3090 maxPower = channelParams->m_clusterPower[cIndex];
3091 cluster1st = cIndex;
3094 channelParams->m_cluster1st = cluster1st;
3096 for (uint8_t cIndex = 0; cIndex < channelParams->m_reducedClusterNumber; cIndex++)
3098 if (maxPower < channelParams->m_clusterPower[cIndex] && cluster1st != cIndex)
3100 maxPower = channelParams->m_clusterPower[cIndex];
3101 cluster2nd = cIndex;
3104 channelParams->m_cluster2nd = cluster2nd;
3106 NS_LOG_INFO(
"1st strongest cluster:" << +cluster1st
3107 <<
", 2nd strongest cluster:" << +cluster2nd);
3110 if (cluster1st == cluster2nd)
3112 clusterDelay.push_back(clusterDelay[cluster1st] + 1.28 * table3gpp->m_cDS);
3113 clusterDelay.push_back(clusterDelay[cluster1st] + 2.56 * table3gpp->m_cDS);
3115 clusterAoa.push_back(clusterAoa[cluster1st]);
3116 clusterAoa.push_back(clusterAoa[cluster1st]);
3118 clusterZoa.push_back(clusterZoa[cluster1st]);
3119 clusterZoa.push_back(clusterZoa[cluster1st]);
3121 clusterAod.push_back(clusterAod[cluster1st]);
3122 clusterAod.push_back(clusterAod[cluster1st]);
3124 clusterZod.push_back(clusterZod[cluster1st]);
3125 clusterZod.push_back(clusterZod[cluster1st]);
3131 if (cluster1st < cluster2nd)
3141 clusterDelay.push_back(clusterDelay[min] + 1.28 * table3gpp->m_cDS);
3142 clusterDelay.push_back(clusterDelay[min] + 2.56 * table3gpp->m_cDS);
3143 clusterDelay.push_back(clusterDelay[max] + 1.28 * table3gpp->m_cDS);
3144 clusterDelay.push_back(clusterDelay[max] + 2.56 * table3gpp->m_cDS);
3146 clusterAoa.push_back(clusterAoa[min]);
3147 clusterAoa.push_back(clusterAoa[min]);
3148 clusterAoa.push_back(clusterAoa[max]);
3149 clusterAoa.push_back(clusterAoa[max]);
3151 clusterZoa.push_back(clusterZoa[min]);
3152 clusterZoa.push_back(clusterZoa[min]);
3153 clusterZoa.push_back(clusterZoa[max]);
3154 clusterZoa.push_back(clusterZoa[max]);
3156 clusterAod.push_back(clusterAod[min]);
3157 clusterAod.push_back(clusterAod[min]);
3158 clusterAod.push_back(clusterAod[max]);
3159 clusterAod.push_back(clusterAod[max]);
3161 clusterZod.push_back(clusterZod[min]);
3162 clusterZod.push_back(clusterZod[min]);
3163 clusterZod.push_back(clusterZod[max]);
3164 clusterZod.push_back(clusterZod[max]);
3167 channelParams->m_delay = clusterDelay;
3168 channelParams->m_angle.clear();
3169 channelParams->m_angle.push_back(clusterAoa);
3170 channelParams->m_angle.push_back(clusterZoa);
3171 channelParams->m_angle.push_back(clusterAod);
3172 channelParams->m_angle.push_back(clusterZod);
3175 channelParams->m_cachedAngleSincos.resize(channelParams->m_angle.size());
3176 for (
size_t direction = 0; direction < channelParams->m_angle.size(); direction++)
3178 channelParams->m_cachedAngleSincos[direction].resize(
3179 channelParams->m_angle[direction].size());
3180 for (
size_t cluster = 0; cluster < channelParams->m_angle[direction].size(); cluster++)
3182 channelParams->m_cachedAngleSincos[direction][cluster] = {
3183 sin(channelParams->m_angle[direction][cluster] *
DEG2RAD),
3184 cos(channelParams->m_angle[direction][cluster] *
DEG2RAD)};
3204 uint8_t updatedClusterNumber = (channelParams->m_reducedClusterNumber == 1)
3205 ? channelParams->m_reducedClusterNumber + 2
3206 : channelParams->m_reducedClusterNumber + 4;
3208 for (uint8_t cIndex = 0; cIndex < updatedClusterNumber; cIndex++)
3217 dopplerTermAlpha.push_back(alpha);
3218 dopplerTermD.push_back(D);
3220 channelParams->m_alpha = dopplerTermAlpha;
3221 channelParams->m_D = dopplerTermD;
3223 return channelParams;
3242 channelMatrix->m_nodeIds =
3245 bool isSameDirection = (channelParams->m_nodeIds == channelMatrix->m_nodeIds);
3256 if (isSameDirection)
3258 rayAodRadian = channelParams->m_rayAodRadian;
3259 rayAoaRadian = channelParams->m_rayAoaRadian;
3260 rayZodRadian = channelParams->m_rayZodRadian;
3261 rayZoaRadian = channelParams->m_rayZoaRadian;
3265 rayAodRadian = channelParams->m_rayAoaRadian;
3266 rayAoaRadian = channelParams->m_rayAodRadian;
3267 rayZodRadian = channelParams->m_rayZoaRadian;
3268 rayZoaRadian = channelParams->m_rayZodRadian;
3274 size_t uSize = uAntenna->GetNumElems();
3275 size_t sSize = sAntenna->GetNumElems();
3281 uint16_t numOverallCluster = (channelParams->m_cluster1st != channelParams->m_cluster2nd)
3282 ? channelParams->m_reducedClusterNumber + 4
3283 : channelParams->m_reducedClusterNumber + 2;
3285 NS_ASSERT(channelParams->m_reducedClusterNumber <= channelParams->m_clusterPhase.size());
3286 NS_ASSERT(channelParams->m_reducedClusterNumber <= channelParams->m_clusterPower.size());
3287 NS_ASSERT(channelParams->m_reducedClusterNumber <=
3288 channelParams->m_crossPolarizationPowerRatios.size());
3289 NS_ASSERT(channelParams->m_reducedClusterNumber <= rayZoaRadian.size());
3290 NS_ASSERT(channelParams->m_reducedClusterNumber <= rayZodRadian.size());
3291 NS_ASSERT(channelParams->m_reducedClusterNumber <= rayAoaRadian.size());
3292 NS_ASSERT(channelParams->m_reducedClusterNumber <= rayAodRadian.size());
3293 NS_ASSERT(table3gpp->m_raysPerCluster <= channelParams->m_clusterPhase[0].size());
3294 NS_ASSERT(table3gpp->m_raysPerCluster <=
3295 channelParams->m_crossPolarizationPowerRatios[0].size());
3296 NS_ASSERT(table3gpp->m_raysPerCluster <= rayZoaRadian[0].size());
3297 NS_ASSERT(table3gpp->m_raysPerCluster <= rayZodRadian[0].size());
3298 NS_ASSERT(table3gpp->m_raysPerCluster <= rayAoaRadian[0].size());
3299 NS_ASSERT(table3gpp->m_raysPerCluster <= rayAodRadian[0].size());
3301 double x = sMob->GetPosition().x - uMob->GetPosition().x;
3302 double y = sMob->GetPosition().y - uMob->GetPosition().y;
3303 double distance2D = sqrt(x * x + y * y);
3306 double hUt = std::min(sMob->GetPosition().z, uMob->GetPosition().z);
3307 double hBs = std::max(sMob->GetPosition().z, uMob->GetPosition().z);
3309 double distance3D = std::sqrt(distance2D * distance2D + (hBs - hUt) * (hBs - hUt));
3311 Angles sAngle(uMob->GetPosition(), sMob->GetPosition());
3312 Angles uAngle(sMob->GetPosition(), uMob->GetPosition());
3324 for (
size_t polSa = 0; polSa < sAntenna->GetNumPols(); ++polSa)
3326 for (
size_t polUa = 0; polUa < uAntenna->GetNumPols(); ++polUa)
3328 raysPreComp[std::make_pair(polSa, polUa)] =
3329 Complex2DVector(channelParams->m_reducedClusterNumber, table3gpp->m_raysPerCluster);
3334 sinCosA.resize(channelParams->m_reducedClusterNumber);
3335 sinSinA.resize(channelParams->m_reducedClusterNumber);
3336 cosZoA.resize(channelParams->m_reducedClusterNumber);
3337 sinCosD.resize(channelParams->m_reducedClusterNumber);
3338 sinSinD.resize(channelParams->m_reducedClusterNumber);
3339 cosZoD.resize(channelParams->m_reducedClusterNumber);
3340 for (uint8_t nIndex = 0; nIndex < channelParams->m_reducedClusterNumber; nIndex++)
3342 sinCosA[nIndex].resize(table3gpp->m_raysPerCluster);
3343 sinSinA[nIndex].resize(table3gpp->m_raysPerCluster);
3344 cosZoA[nIndex].resize(table3gpp->m_raysPerCluster);
3345 sinCosD[nIndex].resize(table3gpp->m_raysPerCluster);
3346 sinSinD[nIndex].resize(table3gpp->m_raysPerCluster);
3347 cosZoD[nIndex].resize(table3gpp->m_raysPerCluster);
3350 for (uint8_t nIndex = 0; nIndex < channelParams->m_reducedClusterNumber; nIndex++)
3352 for (uint8_t mIndex = 0; mIndex < table3gpp->m_raysPerCluster; mIndex++)
3354 DoubleVector initialPhase = channelParams->m_clusterPhase[nIndex][mIndex];
3356 double k = channelParams->m_crossPolarizationPowerRatios[nIndex][mIndex];
3360 for (uint8_t polUa = 0; polUa < uAntenna->GetNumPols(); ++polUa)
3362 auto [rxFieldPatternPhi, rxFieldPatternTheta] = uAntenna->GetElementFieldPattern(
3363 Angles(channelParams->m_rayAoaRadian[nIndex][mIndex],
3364 channelParams->m_rayZoaRadian[nIndex][mIndex]),
3366 for (uint8_t polSa = 0; polSa < sAntenna->GetNumPols(); ++polSa)
3368 auto [txFieldPatternPhi, txFieldPatternTheta] =
3369 sAntenna->GetElementFieldPattern(
3370 Angles(channelParams->m_rayAodRadian[nIndex][mIndex],
3371 channelParams->m_rayZodRadian[nIndex][mIndex]),
3373 raysPreComp[std::make_pair(polSa, polUa)](nIndex, mIndex) =
3374 std::complex<double>(cos(initialPhase[0]), sin(initialPhase[0])) *
3375 rxFieldPatternTheta * txFieldPatternTheta +
3376 std::complex<double>(cos(initialPhase[1]), sin(initialPhase[1])) *
3377 std::sqrt(1.0 / k) * rxFieldPatternTheta * txFieldPatternPhi +
3378 std::complex<double>(cos(initialPhase[2]), sin(initialPhase[2])) *
3379 std::sqrt(1.0 / k) * rxFieldPatternPhi * txFieldPatternTheta +
3380 std::complex<double>(cos(initialPhase[3]), sin(initialPhase[3])) *
3381 rxFieldPatternPhi * txFieldPatternPhi;
3387 double sinRayZoa = sin(rayZoaRadian[nIndex][mIndex]);
3388 double sinRayAoa = sin(rayAoaRadian[nIndex][mIndex]);
3389 double cosRayAoa = cos(rayAoaRadian[nIndex][mIndex]);
3390 sinCosA[nIndex][mIndex] = sinRayZoa * cosRayAoa;
3391 sinSinA[nIndex][mIndex] = sinRayZoa * sinRayAoa;
3392 cosZoA[nIndex][mIndex] = cos(rayZoaRadian[nIndex][mIndex]);
3396 double sinRayZod = sin(rayZodRadian[nIndex][mIndex]);
3397 double sinRayAod = sin(rayAodRadian[nIndex][mIndex]);
3398 double cosRayAod = cos(rayAodRadian[nIndex][mIndex]);
3399 sinCosD[nIndex][mIndex] = sinRayZod * cosRayAod;
3400 sinSinD[nIndex][mIndex] = sinRayZod * sinRayAod;
3401 cosZoD[nIndex][mIndex] = cos(rayZodRadian[nIndex][mIndex]);
3407 uint8_t numSubClustersAdded = 0;
3408 for (uint8_t nIndex = 0; nIndex < channelParams->m_reducedClusterNumber; nIndex++)
3410 for (
size_t uIndex = 0; uIndex < uSize; uIndex++)
3412 Vector uLoc = uAntenna->GetElementLocation(uIndex);
3414 for (
size_t sIndex = 0; sIndex < sSize; sIndex++)
3416 Vector sLoc = sAntenna->GetElementLocation(sIndex);
3419 if (nIndex != channelParams->m_cluster1st && nIndex != channelParams->m_cluster2nd)
3421 std::complex<double> rays(0, 0);
3422 for (uint8_t mIndex = 0; mIndex < table3gpp->m_raysPerCluster; mIndex++)
3425 double rxPhaseDiff =
3427 (sinCosA[nIndex][mIndex] * uLoc.x + sinSinA[nIndex][mIndex] * uLoc.y +
3428 cosZoA[nIndex][mIndex] * uLoc.z);
3430 double txPhaseDiff =
3432 (sinCosD[nIndex][mIndex] * sLoc.x + sinSinD[nIndex][mIndex] * sLoc.y +
3433 cosZoD[nIndex][mIndex] * sLoc.z);
3436 rays += raysPreComp[std::make_pair(sAntenna->GetElemPol(sIndex),
3437 uAntenna->GetElemPol(uIndex))](nIndex,
3439 std::complex<double>(cos(rxPhaseDiff), sin(rxPhaseDiff)) *
3440 std::complex<double>(cos(txPhaseDiff), sin(txPhaseDiff));
3443 sqrt(channelParams->m_clusterPower[nIndex] / table3gpp->m_raysPerCluster);
3444 hUsn(uIndex, sIndex, nIndex) = rays;
3448 std::complex<double> raysSub1(0, 0);
3449 std::complex<double> raysSub2(0, 0);
3450 std::complex<double> raysSub3(0, 0);
3452 for (uint8_t mIndex = 0; mIndex < table3gpp->m_raysPerCluster; mIndex++)
3456 double rxPhaseDiff =
3458 (sinCosA[nIndex][mIndex] * uLoc.x + sinSinA[nIndex][mIndex] * uLoc.y +
3459 cosZoA[nIndex][mIndex] * uLoc.z);
3461 double txPhaseDiff =
3463 (sinCosD[nIndex][mIndex] * sLoc.x + sinSinD[nIndex][mIndex] * sLoc.y +
3464 cosZoD[nIndex][mIndex] * sLoc.z);
3466 std::complex<double> raySub =
3467 raysPreComp[std::make_pair(sAntenna->GetElemPol(sIndex),
3468 uAntenna->GetElemPol(uIndex))](nIndex,
3470 std::complex<double>(cos(rxPhaseDiff), sin(rxPhaseDiff)) *
3471 std::complex<double>(cos(txPhaseDiff), sin(txPhaseDiff));
3495 sqrt(channelParams->m_clusterPower[nIndex] / table3gpp->m_raysPerCluster);
3497 sqrt(channelParams->m_clusterPower[nIndex] / table3gpp->m_raysPerCluster);
3499 sqrt(channelParams->m_clusterPower[nIndex] / table3gpp->m_raysPerCluster);
3500 hUsn(uIndex, sIndex, nIndex) = raysSub1;
3503 channelParams->m_reducedClusterNumber + numSubClustersAdded) = raysSub2;
3506 channelParams->m_reducedClusterNumber + numSubClustersAdded + 1) =
3511 if (nIndex == channelParams->m_cluster1st || nIndex == channelParams->m_cluster2nd)
3513 numSubClustersAdded += 2;
3520 std::complex<double> phaseDiffDueToDistance(cos(-2 * M_PI * distance3D / lambda),
3521 sin(-2 * M_PI * distance3D / lambda));
3525 const double sinUAngleAz = sin(uAngle.
GetAzimuth());
3526 const double cosUAngleAz = cos(uAngle.
GetAzimuth());
3529 const double sinSAngleAz = sin(sAngle.
GetAzimuth());
3530 const double cosSAngleAz = cos(sAngle.
GetAzimuth());
3532 for (
size_t uIndex = 0; uIndex < uSize; uIndex++)
3534 Vector uLoc = uAntenna->GetElementLocation(uIndex);
3535 double rxPhaseDiff = 2 * M_PI *
3536 (sinUAngleIncl * cosUAngleAz * uLoc.x +
3537 sinUAngleIncl * sinUAngleAz * uLoc.y + cosUAngleIncl * uLoc.z);
3539 for (
size_t sIndex = 0; sIndex < sSize; sIndex++)
3541 Vector sLoc = sAntenna->GetElementLocation(sIndex);
3542 std::complex<double> ray(0, 0);
3543 double txPhaseDiff =
3545 (sinSAngleIncl * cosSAngleAz * sLoc.x + sinSAngleIncl * sinSAngleAz * sLoc.y +
3546 cosSAngleIncl * sLoc.z);
3548 auto [rxFieldPatternPhi, rxFieldPatternTheta] = uAntenna->GetElementFieldPattern(
3550 uAntenna->GetElemPol(uIndex));
3551 auto [txFieldPatternPhi, txFieldPatternTheta] = sAntenna->GetElementFieldPattern(
3553 sAntenna->GetElemPol(sIndex));
3555 ray = (rxFieldPatternTheta * txFieldPatternTheta -
3556 rxFieldPatternPhi * txFieldPatternPhi) *
3557 phaseDiffDueToDistance *
3558 std::complex<double>(cos(rxPhaseDiff), sin(rxPhaseDiff)) *
3559 std::complex<double>(cos(txPhaseDiff), sin(txPhaseDiff));
3561 double kLinear = pow(10, channelParams->m_K_factor / 10.0);
3563 hUsn(uIndex, sIndex, 0) =
3564 sqrt(1.0 / (kLinear + 1)) * hUsn(uIndex, sIndex, 0) +
3565 sqrt(kLinear / (1 + kLinear)) * ray /
3567 channelParams->m_attenuation_dB[0] / 10.0);
3568 for (
size_t nIndex = 1; nIndex < hUsn.
GetNumPages(); nIndex++)
3570 hUsn(uIndex, sIndex, nIndex) *=
3571 sqrt(1.0 / (kLinear + 1));
3577 NS_LOG_DEBUG(
"Husn (sAntenna, uAntenna):" << sAntenna->GetId() <<
", " << uAntenna->GetId());
3578 for (
size_t cIndex = 0; cIndex < hUsn.
GetNumPages(); cIndex++)
3580 for (
size_t rowIdx = 0; rowIdx < hUsn.
GetNumRows(); rowIdx++)
3582 for (
size_t colIdx = 0; colIdx < hUsn.
GetNumCols(); colIdx++)
3584 NS_LOG_DEBUG(
" " << hUsn(rowIdx, colIdx, cIndex) <<
",");
3589 NS_LOG_INFO(
"size of coefficient matrix (rows, columns, clusters) = ("
3592 channelMatrix->m_channel = hUsn;
3593 return channelMatrix;