卫星星座对通信链路的干扰

这个示例使用:

卫星通信工具箱卫星通信工具箱

这个例子展示了如何分析从中地球轨道(MEO)卫星星座到位于太平洋的地面站的下行链路上的干扰。干扰星座由40颗近地轨道卫星(LEO)组成。这个例子确定了下行链路关闭的时间、载波噪声加干扰比和链路裕度。

本示例需要“卫星通信工具箱”。如果您也有天线工具箱™，您可以使用本示例学习如何从天线工具箱导入天线到卫星场景。如果您也有相控阵系统工具箱™，您可以使用本示例学习如何从相控阵系统工具箱导入天线到卫星场景，并使用波束形成来提高载波噪声干扰比(CNIR)。

创建卫星场景

创建一个卫星场景。定义场景的开始时间和结束时间。将采样时间设置为60秒。

开始时间= datetime(2021 3, 17日,22日,52岁,0);% 17 2021年3月10:52 PM UTCstopTime = startTime + minutes(10);% 17 2021年3月11:02 PM UTCsampleTime = 60;% ssc = satelliteScenario(开始时间、stopTime sampleTime);

添加中地球轨道卫星

通过指定其开普勒轨道元素，在MEO中添加一颗卫星。该卫星是数据下行链接的卫星。

semiMajorAxis = 12000000;在m %离心率= 0;倾向= 8;%的度raan = 0;上升节点的右经，以度数为单位argOfPeriapsis = 0;%的度trueAnomaly = 343.9391;%的度meoSat =卫星(sc，半maijoraxis，.．.怪癖,.．.倾向,.．.raan,.．.argOfPeriapsis,.．.trueAnomaly,.．.Name =“MEO卫星”，.．.OrbitPropagator =“two-body-keplerian”）;

增加干扰卫星星座

从双行元素(TLE)文件中添加干扰卫星星座。这些卫星被放置在低地球轨道上。

interferingSat =卫星(sc,“leoSatelliteConstellation.tle”）;

向MEO卫星添加发射机

给MEO卫星加一个发射机。这个发射机用于下行链路。定义天线规格，设置工作载波频率为3ghz。

txMEOFreq = 3 e9;%在赫兹txMEOSat =发射机(meoSat,.．.频率= txMEOFreq,.．.%在赫兹功率= 11);%在瓦分贝gaussianAntenna (txMEOSat.．.DishDiameter = 1);在m %

在近地轨道卫星中添加发射机

为LEO星座中的每颗卫星添加发射机，然后定义天线规格。这些发射器会干扰来自MEO卫星的下行链路。设置干扰卫星工作载波频率为2.99 GHz。该示例为每个干扰卫星分配了10到20 dBW范围内的随机功率。

interferenceFreq = 2.99 e9;%在赫兹rng (“默认”）;txInterferingSat =发射机(interferingSat,.．.频率= interferenceFreq,.．.%在赫兹功率= 10 + 10 *兰特(1,元素个数(interferingSat)));%在瓦分贝gaussianAntenna (txInterferingSat.．.DishDiameter = 0.2);在m %

添加地面站

通过指定地面站的纬度和经度，向卫星场景添加地面站。

gs = groundStation (sc,.．.0，.．.%纬度(以度为单位)180年,.．.%经度(以度表示)Name =“地面站”）;

指定地面站天线类型

对于本例，您可以从以下天线中选择一种:

高斯天线
抛物面反射器从天线工具箱
相控阵系统工具箱中的均匀矩形阵列

选择所需的地面站天线。groundStationAntennaType =“高斯天线”；

将接收器添加到地面站

给地面站增加一个接收器。如果你选择高斯天线或抛物线反射器，将接收器连接到一个万向节上，该万向节依次连接到地面站。配置万向节跟踪MEO卫星，使天线也跟踪MEO卫星。如果您选择了均匀矩形阵列，请直接将接收器连接到地面站。确定万向节和接收机的安装位置和安装角度，天线的规格。

开关groundStationAntennaType情况下｛“高斯天线”，“抛物面反射器”｝当选择高斯天线或抛物面反射器时，附加一个到地面站的万向节。gim =万向节(gs,.．.MountingLocation = (0, 0; 5),.．.在m %MountingAngles = [0; 180; 0]);%的度设置万向节跟踪MEO卫星。。pointAt (gim meoSat);如果groundStationAntennaType = =“高斯天线”%当选择高斯天线时创建接收器对象并将其添加到万向节。rxGs =接收机(gim时,.．.MountingLocation = [0, 0, 1]);在m %提供高斯天线的规格gaussianAntenna (rxGs.．.DishDiameter = 0.8);在m %其他的%当抛物面反射器被选中根据频率调整天线的大小%需要天线工具箱(TM)蚂蚁=设计(reflectorParabolic txMEOFreq);创建接收器对象并将其添加到万向节。rxGs =接收机(gim时,.．.天线=蚂蚁,.．.MountingLocation = [0, 0, 1]);在m %结束情况下“统一矩形数组”%当选择“均匀矩形阵列”时确定下行信号的波长c = physconst (“光速”）;λ= c / txMEOFreq;定义数组大小nrow = 8;ncol = 8;定义元素间距卓尔精灵=λ/ 2;dcol =λ/ 2;创建一个背折的6 × 6天线阵列%需要相控阵系统工具箱(TM)蚂蚁=分阶段。市建局(Size = [nrow ncol]，.．.ElementSpacing = [drow dcol]);ant.Element.BackBaffled = true;创建接收对象并将其添加到地面站。rxGs =接收机(gs,.．.天线=蚂蚁,.．.MountingAngles = [0; 90; 0]);%的度结束

建立干扰卫星星座与地面站之间的访问分析

增加干扰星座中每个卫星与地面站之间的访问分析。这种分析使得在稍后将发射的卫星场景查看器中可视化干扰成为可能。任何时候，只要地面站看到星座中的卫星，就会受到来自该可见卫星的某种程度的干扰。

ac =访问(interferingSat, gs);ac.LineColor = [1 1 0];%的黄色

为卫星设定跟踪目标

让卫星追踪地面站。这确保了每颗卫星上的发射机天线都能跟踪地面站。将干扰卫星发射机设置为跟踪地面站，会导致下行链路出现最坏情况的干扰。

pointAt ([meoSat interferingSat], gs);

均匀矩形阵列权值计算

如果选择“均匀矩形阵列”作为地面站天线，计算将主瓣指向MEO卫星所需的权值，并将空值指向干扰卫星所需的权值，从而消除干扰。使用指定计算的权重pointAt．

如果groundStationAntennaType = =“统一矩形数组”找到在地面视线内的近地轨道卫星。%。这些卫星是潜在的干扰者。currentInterferingSat =干扰ingsat (accessStatus(ac,sc.StartTime) == true);计算MEO卫星相对于的方向%的数组。这是瞭望方向。[azd,古人]= aer (rxGs meoSat, sc。开始时间,CoordinateFrame =“身体”）;计算可能干扰卫星的方向。%相对于数组。这些是零方向。[下跌,eln] = aer (rxGs currentInterferingSat, sc。开始时间,CoordinateFrame =“身体”）;计算瞭望方向的转向矢量。。需要相控阵系统工具箱。wd = steervec (getElementPosition (ant) /λ,[wrapTo180 (azd);高龄]);为零方向计算转向矢量。需要相控阵系统工具箱。wn = steervec (getElementPosition (ant) /λ,[wrapTo180(下跌);民族解放军的]);在零方向上计算所需转向的响应。。rn = (wn * wn) \ (wn * wd);旁瓣消去器-消除零方向的响应。w = rn wd-wn *;为相控阵分配权重。pointAt (rxGs、重量= w);结束

创建所需的下行

从MEO卫星上的发射机到地面站上的接收器建立一条下行链路。这条链路是遇到LEO星座干扰的下行链路。

下行=链接(txMEOSat rxGs);

创建干扰链接

在LEO星座中每个卫星上的发射机和地面站上的接收器之间建立一个连接。这些链路是与所需下行链路的干扰链路。

lnkInterference =链接(txInterferingSat rxGs);

发射卫星场景查看器

启动卫星场景查看器ShowDetails设置为false。当ShowDetails属性设置为假，只显示卫星、地面站、入口和链接。标签和轨道将被隐藏。将鼠标移到卫星和地面站上以显示它们的标签。单击MEO卫星，使其轨道投影到场景中StopTime它的标签是可见的，不需要鼠标移动。单击地面站，使其标签可见，而无需鼠标悬停。在MEO卫星上的发射机和地面站上的接收机之间出现绿线，意味着假设不存在来自卫星星座的干扰，下行链路可以成功关闭。星座中给定卫星与地面站之间的黄线表示它们可以相互访问，因此，该卫星的干扰存在。

v = satelliteScenarioViewer (sc ShowDetails = false);

下行链路中天线辐射图的可视化

可视化MEO卫星上的发射天线和地面站上的接收器的辐射模式。

模式(txMEOSat,.．.大小= 1000000);在m %模式(rxGs txMEOFreq,.．.大小= 1000000);在m %

MEO卫星天线的辐射谱图

设置摄像头为“查看地面站天线辐射图”

设置摄像头位置和方向以查看地面站天线。%的辐射模式。坎波斯(v, -8172, 2500000);camheading (v, 40);campitch (v, -60);

高斯天线

抛物面反射器

均匀的矩形阵列

模拟场景和可视化

带高斯天线或抛物面反射器

如果您选择了高斯天线或抛物面反射器(需要天线工具箱)，请使用玩使情景形象化开始时间来StopTime．这将在回放可视化之前自动模拟场景。注意当万向节跟踪MEO卫星时天线的指向是如何变化的。

如果groundStationAntennaType = =“高斯天线”| | groundStationAntennaType = =“抛物面反射器”玩(sc);坎波斯(v, -8172, 2500000);camheading (v, 40);campitch (v, -60);结束

与高斯天线

用抛物面反射器

均匀矩形阵列

如果选择了均匀矩形阵列(需要相控阵系统工具箱)，则必须手动分步进行模拟，以便能够基于MEO卫星和干扰LEO卫星的新位置重新计算每个时间步的权重。要手动步进模拟，首先设置AutoSimulate为假。接下来，您可以调用推进将模拟移动一个时间步。第一个电话推进将计算模拟状态在开始时间．随后的调用将使时间提前一步SampleTime然后计算相应的状态。

如果groundStationAntennaType = =“统一矩形数组”将自动模拟设置为false。sc.AutoSimulate = false;手动步进模拟。而推进(sc)确定每个LEO卫星的访问状态历史%，对应当前的SimulationTime。acStatusHistory = accessStatus (ac);acStatus = acStatusHistory(:,结束);确定地面可见的近地轨道卫星。%。这些卫星可能会干扰地面站目前的模拟。%的时间。currentInterferingSat =干扰ingsat (acStatus == true);的体帧内确定MEO卫星的方向均匀矩形阵列。这是瞭望方向%的数组。[azdHistory, eldHistory] = aer (rxGs、meoSat CoordinateFrame =“身体”）;azd = azdHistory(:,结束);古人= eldHistory(:,结束);确定这些干扰卫星的方向均匀矩形阵列的主体框架。这些都是%数组必须指向null的方向。[aznHistory, elnHistory] = aer (rxGs、currentInterferingSat CoordinateFrame =“身体”）;有= aznHistory(:,结束);eln = elnHistory(:,结束);计算瞭望方向的转向矢量。。需要相控阵系统工具箱。wd = steervec (getElementPosition (ant) /λ,[wrapTo180 (azd);高龄]);为零方向计算转向矢量。需要相控阵系统工具箱。wn = steervec (getElementPosition (ant) /λ,[wrapTo180(下跌);民族解放军的]);计算所需转向在零方向的响应。。rn = (wn * wn) \ (wn * wd);旁瓣消去器-删除零方向的响应。w = rn wd-wn *;为相控阵分配权重。pointAt (rxGs、重量= w);结束结束

绘制忽略干扰的下行链路关闭状态

从MEO卫星确定所需下行链路的关闭状态。的linkStatus功能忽略来自其他发射器的干扰。当下行链路关闭时，状态为true。否则，状态为false。在图中，状态分别用1和0表示。

[downlinkStatus t] = linkStatus(下行);情节(t, downlinkStatus“g”、线宽= 2);包含(“时间”）;ylabel (“下行关闭状态”）;标题(“连接状态与时间的函数”）;网格在；

图中包含一个axes对象。标题为“Link Status as a Function of Time”的axes对象包含一个类型为line的对象。

计算有干扰的下行链路关闭状态

通过计算地面站接收机输入的MEO下行链路和干扰信号功率等级，计算有干扰的下行链路关闭状态sigstrength。接收到的功率测量和损耗的位置在下面的地面站接收图中说明。

计算与下行链路对应的接收端输入功率MEO卫星。[~, downlinkPowerRxInput] = sigstrength(下行);%在瓦分贝计算每个对应的接收机输入处的干扰功率% LEO卫星。[~, interferencePowerRxInput] = sigstrength (lnkInterference);%在瓦分贝

计算接收机输入处的总干扰信号功率。通过将干扰LEO卫星的功率水平以瓦特为单位相加得到这个量。

interferencePowerRxInputW = 10 ^ (interferencePowerRxInput / 10);% WinterferencePowerRxInputSumW =总和(interferencePowerRxInputW);% W

按照以下步骤计算在信号带宽中造成干扰的总干扰信号功率。

1)计算信号带宽与干扰带宽的重叠部分。本例认为干扰卫星和MEO卫星的传输功率在各自的MEO卫星和干扰卫星的整个带宽上是恒定的。

2)计算对信号带宽起干扰作用的干扰功率。

此图为功率谱密度(PSD)图，图中显示了传输带宽和干扰带宽重叠时的实际干扰功率和模拟干扰功率。实际的干扰功率是重叠带宽区域的干扰功率密度所占的面积。这个实际的干扰功率被分散到整个传输带宽中，并假定为类噪声。

本例假设MEO卫星的传输(或信号)带宽为30 MHz，干扰信号的带宽为20 MHz。

txBandwidth = 30 e6;%在赫兹interferenceBandWidth = 20 e6;%在赫兹的带宽与干扰带宽的重叠部分%的兴趣。假设它们之间有相同的干扰功率%整个带宽。overlapFactor = getOverlapFactor (txMEOFreq txBandwidth,.．.interferenceFreq interferenceBandWidth);获取对信号产生干扰的干扰功率总干扰功率的%interferencePowerRxInputActual = interferencePowerRxInputSumW * overlapFactor;在W %

将干扰信号功率在重叠带宽中的贡献视为噪声来建模干扰。相应地，将此量加到地面站接收端输入的热噪声中。注意，干扰和噪声功率级必须以瓦为单位。

%计算地面站接收端输入的热噪声。T = HelperGetNoiseTemperature (txMEOFreq rxGs);%在Kkb = physconst (“波尔兹曼”）;thermalNoise = kb * T * txBandwidth;在W %%计算在接收机输入端的噪声加上干扰功率。noisePlusInterferencePowerW = thermalNoise + interferencePowerRxInputActual;在W %noisePlusInterferencePower = 10 * log10 (noisePlusInterferencePowerW);%在瓦分贝

计算解调器输入端载波噪声加干扰功率谱密度比如下:

$C / （ N_{0} + 我_{0} ）＝ P_{R x 我 n p u t} - （ N + 我） + 10 {日志}_{10} T x B 一个 n d w 我 d t h - l O 年代 {年代}_{R x}$ ，

地点:

$C / （ N_{0} + 我_{0} ）$ 为解调器输入端载波噪声加干扰功率密度比(单位为dB)。
$P_{R x 我 n p u t}$ 为在地面站接收器输入处测量到的MEO卫星接收到的下行功率(dBW)。
$（ N + 我）$ 为接收机系统热噪声与干扰信号功率在接收机输入处测量的重叠带宽中的贡献之和(单位为dBW)。
$T x B 一个 n d w 我 d t h$ 为MEO卫星下行传输带宽，单位为Hz。
$l O 年代 {年代}_{R x}$ 为发生在接收机输入和解调器输入之间的损耗(单位为dB)。

计算发生在接收机输入和解调器之间的损耗%的输入。rxGsLoss = rxGs。年代ystemLoss - rxGs.PreReceiverLoss;在解调器输入端计算C/(N0+I0)。CNoPlusInterference = downlinkPowerRxInput -.．.噪声加干扰功率+ 10*log10(txBandwidth) - rxGsLoss;

计算在解调器输入处的比特能量与噪声加干扰功率谱密度比如下:

$E_{b} / （ N_{0} + 我_{0} ）＝ C / （ N_{0} + 我_{0} ） - 10 {日志}_{10} B 我 T R 一个 T E - 60$ ，

地点:

$E_{b} / （ N_{0} + 我_{0} ）$ 为每比特能量对噪声加上在解调器输入处的干扰功率谱密度之比(单位为dB)。
$B 我 T R 一个 T E$ 为MEO卫星下行链路的比特率，单位为Mbps。

比特率= txMEOSat.BitRate;ebNoPlusInterference = CNoPlusInterference - 10*log10(bitRate) - 60;

链路裕量的计算方法如下:

$米一个 R G 我 N ＝ E_{b} / （ N_{0} + 我_{0} ） - {（ E_{b} / N_{0} ）}_{R e 问 u 我 r e d}$

地点:

$米一个 R G 我 N$ 为链接裕量，单位为dB。
${（ E_{b} / N_{0} ）}_{R e 问 u 我 r e d}$ 为解调器输入端闭合链路所需的每比特最小接收能量与噪声功率谱密度之比(单位为dB)。

marginWithInterference = ebNoPlusInterference - rxGs.RequiredEbNo;

计算有干扰的下行链路关闭状态。当链路裕度大于等于0db时，该状态为true。

downlinkStatusWithInterference = marginWithInterference >= 0;

计算每比特能量与噪声功率谱密度比

计算解调器输入端下行链路和干扰链路的比特能量与噪声功率谱密度比(Eb/N0)，供后续分析。

ebnoDownlink = ebno(下行);在dB %ebnoInterference = ebno (lnkInterference);在dB %

绘制有干扰的下行链路关闭状态

绘制考虑干扰的新的下行链路关闭状态。将新的链路状态与之前忽略干扰的情况进行比较。

情节(t, downlinkStatusWithInterference“- r”t downlinkStatus”——g”、线宽= 2);传奇(“干扰贡献”，“被忽视的干扰”）;包含(“时间”）;ylabel (“下行关闭状态”）;标题(“连接状态与时间的函数”）;ylim (1.2 [0]);网格在

图中包含一个axes对象。标题为Link Status as a Function of Time的axes对象包含两个类型为line的对象。这些对象代表被考虑的干扰，被忽略的干扰。

当选择高斯天线或抛物反射面时

图中显示，在晚上10:54，由于过度干扰，下行链路无法关闭。这是因为卫星10LEO星座的一颗卫星从头顶飞过，它的传输信号被它的主瓣接收。这也可以通过将查看器的当前时间设置为10:54 PM并单击天线主瓣附近的卫星来进行视觉确认。注意，你需要天线工具箱选择抛物面反射器。

如果groundStationAntennaType = =“高斯天线”| | groundStationAntennaType = =“抛物面反射器”v.CurrentTime = datetime(2021 3, 17日,22日,54岁,0);结束

高斯天线

抛物面反射器

当选择具有干扰抵消作用的均匀矩形阵列时

如果您选择了均匀矩形阵列(需要相控阵系统工具箱)，图中显示下行链路可以在场景持续期间关闭，因为阵列指向干扰LEO卫星的方向为零。通过将查看器的当前时间设置为10:54 PM和10:55 PM，也可以从视觉上确认这一点。以便能够手动设置查看器CurrentTime，你必须改变AutoSimulate为true。注意，这将清除模拟数据。此外，还需要重新计算这些时间的权重，并将它们分配给使用的数组pointAt．飞过地面站的卫星是卫星10．点击它可以看到它的名字和轨道。按住鼠标左键或滚动键的同时拖动鼠标，将相机带到所需的位置和方向。旋转滚轮来控制相机的变焦。此外，使辐射图案不透明，以清楚地可视化的位置卫星10对脑叶的影响。你可以看到，这两次，卫星10在叶之间。这是因为阵列指向卫星的是空值，从而消除了来自卫星的干扰。

如果groundStationAntennaType = =“统一矩形数组”将autosimulation设置为true。sc.AutoSimulate = true;将查看器CurrentTime设置为10:54 PM。时间= datetime(2021 3, 17日,22日,54岁,0);v.CurrentTime =时间;计算权重并将它们分配给数组。。currentInterferingSat =干扰ingsat (accessStatus(ac,time) == true);[azd,古人]= aer (rxGs meoSat,时间,CoordinateFrame =“身体”）;[下跌,eln] = aer (rxGs currentInterferingSat,时间,CoordinateFrame =“身体”）;需要相控阵系统工具箱。wd = steervec (getElementPosition (ant) /λ,[wrapTo180 (azd);高龄]);wn = steervec (getElementPosition (ant) /λ,[wrapTo180(下跌);民族解放军的]);rn = (wn * wn) \ (wn * wd);w = rn wd-wn *;pointAt (rxGs、重量= w);使辐射图案不透明。模式(rxGs txMEOFreq,.．.大小= 1000000,.．.透明度= 1);结束

在国际标准时间晚上10:54

您可以在时间设置为10:55 PM的情况下运行上述代码，并观察指向干扰卫星新位置的空值。

在国际标准时间晚上10:55

计算和绘制载波噪声比和载波噪声加干扰比

在解调器输入端计算载波噪声比(CNR)如下所示:

$C / N_{0} ＝ E_{b} / N_{0} + 10 {日志}_{10} B 我 T R 一个 T E + 60$ ，

$C / N ＝ C / N_{0} - 10 {日志}_{10} T x B 一个 n d w 我 d t h$ ，

地点:

$C / N$ 为解调器输入处的载波噪声比，单位为dB。
$E_{b} / N_{0}$ 为解调器输入端载噪功率谱密度比，单位为dB。

计算载波与噪声功率谱密度比。= ebnoDownlink + 10*log10(bitRate) + 60;%计算载波噪声比。cByN = CNoDownlink - 10*log10(txBandwidth);

在解调器输入端计算载波噪声加干涉比(CNIR)如下所示:

$C / （ N + 我）＝ C / （ N_{0} + 我_{0} ） - 10 {日志}_{10} T x B 一个 n d w 我 d t h$

cByNPlusI = CNoPlusInterference - 10*log10(txBandwidth);

绘制CNR和CNIR

情节(t, cByNPlusI“- r”t cByN”——g”、线宽= 2);传奇(“CNIR”，“中国北车”位置=“南”）;包含(“时间”）;ylabel ("CNR或CNIR (dB)"）;标题(CNR和CNIR vs.时间+ groundStationAntennaType);网格在

图中包含一个axes对象。标题为CNR和CNIR随时间变化的高斯天线坐标轴对象包含两个类型为line的对象。这些对象表示CNIR、CNR。

当选择具有MEO卫星跟踪和干扰抵消的均匀矩形阵列(需要相控阵系统工具箱)时，由于没有LEO卫星的干扰，CNR和CNIR会重叠。这是因为阵列指向低轨道卫星的零点。从干扰卫星输入的地面站接收机的最大功率约为-167.3 dBW可以证实这一点，这是非常低的。对于本例中使用的所有其他天线，地面站接收器输入的最大功率要高得多(高斯天线为-125.9 dBW，抛物反射器为-127.3 dBW)。

maxInterferencePowerRxInput = max (interferencePowerRxInput [],“所有”）;disp (在整个场景持续时间内，来自干扰LEO卫星的地面站接收器输入的最大功率为+ maxInterferencePowerRxInput +“瓦分贝”。）;

在整个场景持续时间内，来自干扰LEO卫星的地面站接收机输入的最大功率为-125.9301 dBW。

比较有干扰和无干扰时的链接边距

无干扰链路裕度计算如下:

$米一个 R G 我 N ＝ E_{b} / N_{0} - {（ E_{b} / N_{0} ）}_{R e 问 u 我 r e d}$

marginWithoutInterference = ebnoDownlink - rxGs.RequiredEbNo;

绘制有干扰和无干扰的链接边界。

情节(t, marginWithInterference“- r”t marginWithoutInterference”——g”、线宽= 2);传奇(“干扰”，“无干扰”位置=“南”）;包含(“时间”）;ylabel (“保证金(dB)”）;标题(“链接裕度与时间”+ groundStationAntennaType);网格在

图中包含一个axes对象。标题为Link Margin vs. Time的axes对象包含两个类型为line的对象。这些物体代表有干扰，无干扰。

当链路裕度大于或等于0db时，下行链路关闭。在高斯天线、抛物面反射器(需要天线工具箱)和均匀矩形阵列(需要相控阵系统工具箱)没有干扰抵消的情况下，由于干扰，链路裕度有时会下降到0 dB以下。在这些时候，下行链路是断开的。

进一步的探索

本例演示如何分析卫星通信链路上的干扰。链接关闭时间是这些参数的函数:

卫星的轨道
地面站的位置
发射机和接收机的规格
发射天线和接收天线的规格
如果使用均匀矩形数组，则权重
信号和干扰带宽

修改这些参数，可以观察它们对链路干扰程度的影响。发射机和接收机也可以从“天线工具箱”和“相控阵系统工具箱”中选择不同的天线，观察链路性能。在使用相控阵时，如果您只对主瓣跟踪单个目标而不处理指向空值感兴趣，您可以使用pointAt自动跟踪其他卫星、地面站和地理位置，而无需通过设置手动模拟AutoSimulate为假。呼叫的局限性玩当使用动态操纵相控阵时，你不能在模拟过程中可视化其辐射模式的变化。

辅助函数

该示例使用了helper函数HelperGetNoiseTemperature获取接收天线的噪声温度。

本例还使用该局部函数计算了传输带宽与干扰带宽之间的重叠量。

函数overlapFactor = getOverlapFactor (txFreq txBW、interferenceFreq interferenceBW)% getOverlapFactor提供重叠的干扰带宽量%带传输带宽txFreq_Limits = [txFreq-(txBW/2) txFreq+(txBW/2)];interferenceFreq_Limits = [interferenceFreq - (interferenceBW / 2).．.interferenceFreq + (interferenceBW / 2)];如果(interferenceFreq_Limits(2) < txFreq_Limits(1)) || .输出说明.．.(interferenceFreq_Limits (1) > txFreq_Limits (2))%，表示传输带宽与%干扰带宽，则重叠因子为0overlapFactor = 0;elseif(interferenceFreq_Limits(2) <= txFreq_Limits(2)) &&.．.(interferenceFreq_Limits (1) > = txFreq_Limits (1))%如果干扰带宽完全在传输范围内%带宽，则重叠因子为1overlapFactor = 1;elseif(interferenceFreq_Limits(2) > txFreq_Limits(2)) &&.．.(interferenceFreq_Limits (1) < txFreq_Limits (1))%如果传输带宽完全在干扰范围内%带宽，则重叠系数为传输比%带宽与干扰带宽的比值overlapFactor = txBW / interferenceBW;elseif(interferenceFreq_Limits(2) <= txFreq_Limits(2)) &&.．.(interferenceFreq_Limits (1) < = txFreq_Limits (1))%如果传输带宽的起始边缘在%干扰带宽，则重叠因子为的比值干扰带宽的最后一条边与第一条边的差信号带宽的%，与干扰带宽的%相同overlapFactor = (interferenceFreq_Limits (2) -txFreq_Limits (1)) / interferenceBW;其他的%如果传输带宽的最后边缘位于%干扰带宽，则重叠因子为差值之比%从信号带宽的最后一条边到干扰的第一条边%带宽，与干扰带宽值相同overlapFactor = (-interferenceFreq_Limits txFreq_Limits (1) + (2)) / interferenceBW;结束结束