How do i setup coverage by longley rice model as transmitter to eight other receiver sites?

Hi,

Does anybody have experience in setting up a code to determine difference between transmitter antenna model? Or perhaps there's something wrong in the code?

1) Normal dipole antenna

% Define transmitter site at Ayer Rajah Telecoms

fq = 3.5e9; % 3.5 GHz

tx = txsite('Name','Ayer Rajah Telecoms', ...

'Latitude',1.308485, ...

'Longitude',103.776380, ...

'Antenna',design(dipole,fq), ...

'TransmitterFrequency',fq, ... % Units: Hz

'TransmitterPower',5); % Units: Watts

% Define receiver sites in several surrounding towns and cities

rxNames = {...

'SIT@SP', ...

'Tennis Courts', ...

'Dover International Court Bus Stop', ...

'Opposite ACS Independent', ...

'ACS Independent', ...

'Opposite ACS Boarding School', ...

'Fairfield Methodist Primary School',...

'Fusionopolis'};

rxLocations = [1.307522 103.777824; ... % SIT@SP

1.307118 103.778666; ... % Tennis Courts

1.306243 103.780141; ... % Dover International Court Bus Stop

1.305326 103.781079; ... % Opposite ACS Independent

1.303722 103.781820; ... % ACS Independent

1.302816 103.782823; ... % Opposite ACS Boarding School

1.300585 103.784496;... % Fairfield Methodist Primary School

1.299009 103.787652]; % Fusionopolis

% Define receiver sensitivity. Sensitivity is the minimum signal strength in

% power that is necessary for the receiver to accurately detect the signal.

rxSensitivity = -66.5; % Units: dBm

rxs = rxsite('Name',rxNames, ...

'Latitude',rxLocations(:,1), ...

'Longitude',rxLocations(:,2), ...

'Antenna',design(dipole,tx.TransmitterFrequency), ...

'ReceiverSensitivity',rxSensitivity,...

'AntennaHeight', 10); % Units: dBm

viewer = siteviewer;

show(tx)

show(rxs)

viewer.Basemap = "openstreetmap";

coverage(tx,'freespace', ...

'SignalStrengths',rxSensitivity)

link(rxs,tx,'freespace')

2) Longley Rice

% Define transmitter site at Ayer Rajah Telecoms

fq = 3.5e9; % 3.5 GHz

tx = txsite('Name','Ayer Rajah Telecoms', ...

'Latitude',1.308485, ...

'Longitude',103.776380, ...

'Antenna',design(dipole,fq), ...

'TransmitterFrequency',fq, ... % Units: Hz

'TransmitterPower',5); % Units: Watts

% Define receiver sites in several surrounding towns and cities

rxNames = {...

'SIT@SP', ...

'Tennis Courts', ...

'Dover International Court Bus Stop', ...

'Opposite ACS Independent', ...

'ACS Independent', ...

'Opposite ACS Boarding School', ...

'Fairfield Methodist Primary School',...

'Fusionopolis'};

rxLocations = [1.307522 103.777824; ... % SIT@SP

1.307118 103.778666; ... % Tennis Courts

1.306243 103.780141; ... % Dover International Court Bus Stop

1.305326 103.781079; ... % Opposite ACS Independent

1.303722 103.781820; ... % ACS Independent

1.302816 103.782823; ... % Opposite ACS Boarding School

1.300585 103.784496;... % Fairfield Methodist Primary School

1.299009 103.787652]; % Fusionopolis

% Define receiver sensitivity. Sensitivity is the minimum signal strength in

% power that is necessary for the receiver to accurately detect the signal.

rxSensitivity = -66.5; % Units: dBm

rxs = rxsite('Name',rxNames, ...

'Latitude',rxLocations(:,1), ...

'Longitude',rxLocations(:,2), ...

'Antenna',design(dipole,tx.TransmitterFrequency), ...

'ReceiverSensitivity',rxSensitivity,...

'AntennaHeight', 10); % Units: dBm

pm = propagationModel('longley-rice','TimeVariabilityTolerance',0.7)

viewer = siteviewer;

show(tx)

show(rxs)

viewer.Basemap = "openstreetmap";

ss = sigstrength(rxs,tx,pm)

%coverage(rx, 'PropagationModel', pm);

%coverage(tx,'freespace', ...

% 'SignalStrengths',rxSensitivity)

coverage(tx,rxs,pm,"SignalStrengths",-66.5:-5)

link(rxs,tx,'freespace')

Answer 1

[Curious]

Your code works well on longley rice model..

[Curious]

you just need to change freespace to longley-rice