干扰源对信号源的干扰
Matlab实现
sinr的实现原理请看第二篇链接
%载入地图
dtedfile = "n39_w106_3arc_v2.dt1";
attribution = "SRTM 3 arc-second resolution. Data available " + ...
"from the U.S. Geological Survey.";
addCustomTerrain("southboulder",dtedfile,"Attribution",attribution);
viewer = siteviewer("Terrain","southboulder");
%发射机由射频电路和天线组成,射频电路用信号和功率激发天线。发射机的关键特性包括其输出功率、工作频率和天线辐射模式。
%fq=6e9
% %tx = txsite('Name','MathWorks Apple Hill',...
% 'Latitude',39.74356, ...
% 'Longitude',-105.04193,...
% 'TransmitterFrequency',fq);
%定义信号源
tx = txsite('Name','MathWorks Apple Hill',...
'Latitude',39.74356, ...
'Longitude',-105.04193);
show(tx);
% %接收器由射频电路和天线组成,天线拦截无线电波,并将其转换成电流,再由射频电路解码(例如解调)为信号。接收机的关键特性包括其灵敏度和天线辐射图。
% rx = rxsite('Name','MathWorks Lakeside', ...
% 'Latitude',39.54569, ...
% 'Longitude',-105.3764);
% show(rx);
% %使用距离函数计算发射点和接收点之间的距离。默认情况下,该函数在发送端和接收端之间绘制直线,以计算距离。这条路径被称为欧几里得路径,它忽略了包括地球在内的所有障碍物。
% dm = distance(tx,rx);
%
% dkm = dm / 1000;
% %使用角度函数来计算发射点和接收点之间的角度。
% azFromEast = angle(tx,rx); %从东方逆时针方向度
%
% azFromNorth = -azFromEast + 90;%从北转顺时针角度
% %接收功率和链接裕度
% %信号强度是接收机与发射机在特定距离和角度下可用的信号水平。它也被称为接收信号强度。信号强度可以是功率量或电压量。
% ss = sigstrength(rx,tx); %使用sigstrength函数计算接收端发射机的信号强度。可以计算两种信号强度:功率和电场(e-field)。
%
% %链路裕度是通信链路鲁棒性的量度,计算方法是从接收机可用的信号强度中减去所需的接收机灵敏度。
% margin = abs(rx.ReceiverSensitivity - ss);
% %使用link函数在地图上绘制两个站点之间的通信成功情况。链路的成功取决于接收机从发射机接收到的功率。
% %默认情况下,绿色链路表示接收功率满足或超过接收灵敏度。红色的链接表示发送端和接收端之间的通信失败。
% link(rx,tx);
% %发射场覆盖图
% % 覆盖图是发射机地理服务区域的可视化,即参考接收机接收到的信号强度满足其灵敏度的地方。覆盖图可以将信号强度描述为功率量(典型的dBm)或电压量(典型的dBμV/m)。
% % 使用覆盖功能计算发射机的覆盖范围。
% %使用覆盖功能计算发射机的覆盖范围。
% coverage(tx,'close-in', ...
% 'SignalStrengths',-100:5:-60);
% 新发射机位置
% 使用定位函数可以从现有的站点计算出新的接收站点或发射机站点的位置。
% 计算一个新的发射机位置在发射机基地以北1公里。
%[lat,lon] = location(tx,1000,90);
%定义干扰源
rxLocations = [...
39.14356 -105.05193; ... % Boston
39.14556 -105.04493; ... % Lexington
39.23356 -105.03173; ... % Concord
39.45356 -105.02153; ... % Marlborough
39.16456 -105.06133; ... % Hopkinton
39.25356 -105.03143; ... % Holliston
39.64356 -105.041338; ... % Foxborough
39.6756 -105.06183]; % Quincy
tx2 = txsite('Name','Transmitter2','Latitude',rxLocations(:,1),'Longitude',rxLocations(:,2),'AntennaHeight',30);
show(tx2);
%利用信噪比函数计算接收机的信噪比
sinr([tx,tx2]);
close(viewer);
removeCustomTerrain("southboulder")