% === Load directional antenna pattern === vars = load('Directional12dBi.mat'); fns = fieldnames(vars); phasedAntenna = vars.(fns{1}); az = phasedAntenna.AzimuthAngles; el = phasedAntenna.ElevationAngles; magPattern_dBi = magPattern; antenna = phased.CustomAntennaElement( ... 'AzimuthAngles', az, ... 'ElevationAngles', el, ... 'MagnitudePattern', magPattern_dBi, ... 'PhasePattern', zeros(size(magPattern_dBi)) ); % === Create siteviewer === viewer = siteviewer('Name', 'NB-IoT 3-Sector Antenna View'); % === Loop through each tower and each sector === for i = 1:length(towerLat) for angle = [0, 120, 240 ] % 3-sector antennas tx = txsite("Name", sprintf("Tower %d - %d", i, angle), ... "Latitude", towerLat(i), ... "Longitude", towerLong(i), ... "AntennaHeight", 10, ... "TransmitterPower", 40, ... "TransmitterFrequency", 900e6, ... "Antenna", antenna, ... "AntennaAngle", angle); end end this is for the antenna right