校准=真；%执行A部分校准测试systemLevelSimulation=true；%执行B部分系统级模拟showScenarioPlot=true；%显示动态模拟打印更新的步骤

PHYParameters=struct；PHYParameters.TxPower=20；%dBm中的发射机功率PHYParameters.TxGain=0；%发射机天线增益（单位：dBi）PHYParameters.RxGain=-2；%接收机天线增益（单位：dBi）PHYParameters.NoiseFigure=7；%接收机噪声系数（dB）PHYParameters.numtxantenas=1；%发射机天线的数量PHYParameters.numrxanternas=1；%接收天线的数量物理信道带宽=“CBW80”;%系统带宽PHYParameters.TransmitterFrequency=5e9；%发射机频率（单位：Hz）MACParameters=struct；MACParameters.numChannel=3；%非重叠通道的数量MACParameters.CCALevel=-70；%CCA算法中的传输阈值（dBm）

%[x，y，z]方向上的房间数ScenarioParameters=struct；ScenarioParameters.BuildingLayout=[10 2 5]；%每个房间的尺寸（米）[x，y，z]ScenarioParameters.RoomSize=[10 3]；%每个房间的接收器数量。请注意，只能使用一个接收器（STA）%在任何给定时间都处于活动状态。ScenarioParameters.numRxperoom=1；

SimParameters=struct；SimParameters.NumDrops=3；SimParameters.NumTxEventsPerDrop=2；

%变送器总数，假设每个房间有一个变送器（tx）numTx=prod（ScenarioParameters.BuildingLayout）；%使用各向同性天线元件创建发射机站点roomNames=字符串（1，numTx）；对于siteInd=1:numTx房间名（siteInd）=“房间”+siteInd；终止txs=txsite(“笛卡尔”,“姓名”，房间名，...“发射机频率”，PHYParameters.TransmitterFrequency，...“传输功率”，10.^（（PHYParameters.TxPower+PHYParameters.TxGain-30）/10），...“天线”,“各向同性”);

%接收器总数，假设每个房间有一个发射器（tx）numRx=numTx*ScenarioParameters.numrxperoom；%创建接收者站点roomNames=字符串（1，numRx）；对于siteInd=1:numRx房间名（siteInd）=“房间”+（国防部（siteInd-1，numTx）+1）+"-"+ceil（siteInd/numTx）；终止rxs=rxsite(“笛卡尔”,“姓名”，房间名，“天线”,“各向同性”);%dBm中的接收机噪声功率T=290；%温度（开尔文）k=物理常数(“玻尔兹曼”);%玻尔兹曼常数%采样率（Hz）fs=wlanSampleRate（wlanHESUConfig(“信道带宽”、物理参数、信道带宽）；rxNoisePower=10*log10（k*T*fs）+30+PHYParameters.NoiseFigure；

种子=rng（6）；%种子随机数发生器和存储状态如果showScenarioPlot hGrid=tgaxBuildResidentialGrid（ScenarioParameters.RoomSize、ScenarioParameters.BuildingLayout、，...numTx、numRx、MACParameters.NumChannels）；终止如果校准fprintf('正在运行校准…\n');%预分配输出output=struct；output.sinr=零（SimParameters.NumDrops，numTx）；%用于存储SINR值对于drop=1:SimParameters.NumDrops%在每个房间放置接收器[关联，txChannels，rxChannels，txPositions，rxPositions]=tgaxDropNodes(...txs、rxs、ScenarioParameters、MACParameters.numChannel）；%所有发射机均激活activeTx=true（numTx，1）；%每个房间只选择一个接收器rxAlloc=randi（[1 ScenarioParameters.NumRxPerRoom]，numTx，1）；activeRx=重塑（rxAlloc==1:ScenarioParameters.numrxperoom，[]，1）；%生成传播模型propModel=TGaxResidential(“房间大小”、场景参数、房间大小）；%获取每个接收器的发射器索引tnum=repmat（（1:numTx），1，numRx/numTx）；%SINR计算-每个非重叠信道上的环路numChannels=numel（唯一（txChannels））；对于k=1:numChannels%使用第k个非重叠通道tind=txChannels==k；外皮=false（大小（activeRx））；外皮（activeRx）=rxChannels（activeRx）==k；%获取每个活动接收器感兴趣的发射器的索引tsigind=tnum（外皮）；%计算信噪比输出sinr（下降，tind）=sinr（rxs（外皮），txs（tind），...“破产管理人收益”，PHYParameters.RxGain，...“接收者噪音功率”，rxNoisePower，...“传播模型”，propModel，...“信号源”，txs（tsigind））；终止%打印节点和链接如果展示场景图%更新绘图数据掩码=txChannels==rxChannels'；tgaxUpdatePlot（hGrid、txPositions、rxPositions、activeTx、activeRx、mask、TXChannel、RxChannel、，...斯普林特('方框1测试2“仅下行链路”校准，下降#%d/%d'，drop，SimParameters.NumDrops））；终止终止%绘制SINR的CDF，并与校准曲线进行比较tgaxCalibrationCDF（output.sinr，“SS1Box1Test2”,“长期无线电特性”)；fprintf('校准完成\n')终止

如果系统级模拟%预先分配产出output.numPkts=zeros（numRx，1）；output.numPktErrors=zeros（numRx，1）；output.sinrMeas=nan（numRx，SimParameters.nummtxeventsperdrop，SimParameters.NumDrops）；output.sinrEst=nan（numRx，SimParameters.numdtxeventsperdrop，SimParameters.NumDrops）；output.pktErrorRate=0；%为每个可能的发射机创建一个波形配置%示例：所有节点的链路和无线电参数都相同。cfgHEBase=wlanHESUConfig；cfgHEBase.ChannelBandwidth=PHYParameters.ChannelBandwidth；%信道带宽cfgHEBase.NumTransmitAntennas=PHYParameters.NumTxAntennas；%发射天线的数量cfgHEBase.SpatialMapping=“傅里叶”;%空间映射矩阵cfgHEBase.NumSpaceTimeStreams=1；%时空流的数量cfgHEBase.GuardInterval=0.8；%保护间隔持续时间cfgHEBase.helttype=4；%HE-LTF压缩模式cfgHEBase.APEPLength=1e3；%有效负载长度（字节）cfgHEBase.ChannelCoding=“自民党”;%信道编码cfgHEBase.MCS=3；%调制编码方案cfgHE=单元（numTx，1）；对于txidx=1:numTx-cfgHE{txidx}=cfgHEBase；终止fprintf('正在运行系统级模拟…\n')对于drop=1:SimParameters.NumDrops fprintf('正在运行的drop#%d/%d…\n'，drop，SimParameters.NumDrops）；%在每个房间放置接收器[关联，txChannels，rxChannels，txPositions，rxPositions]=tgaxDropNodes(...txs、rxs、ScenarioParameters、MACParameters.numChannel）；%生成传播模型propModel=TGaxResidential(“房间大小”、场景参数、房间大小）；%计算所有链路的信号强度信号强度=信号强度（rxs、txs、propModel、，“类型”,“权力”,...“破产管理人收益”，PHYParameters.RxGain）；%dBm中的所有信号强度%阈值信号低于噪声级以减少模拟时间信号强度（信号强度%不在同一非重叠信道上的阈值信号信号强度（~（txChannels==rxChannels'））=-Inf；%屏蔽不可忽略的发射机-接收机链路%模拟得到线性指标不可忽略的MAK=信号强度>-Inf；%重置不可忽略的通道，为%电流降NonNegibleIdx=查找（NonNegibleMask）'；%为放置中的每个可能的活动链接创建一个通道%配置。在本例中，链路和无线电参数为%对所有节点都一样。tgaxChan=细胞（numel（非忽略指数），1）；对于i=1:numel（不可忽略）%给定链路的发射机索引txIdx=mod（不可忽略的eIdx（i）-1，numTx）+1；%通道配置。每个链路的信道实现如下所示：%由于使用了全局随机流，因此不同。tgaxhanbase=wlanTGaxChannel；tgaxhanbase.DelayProfile=“D型”；tgaxhanbase.numtransmitantenas=cfgHE{txix}.numtransmitantenas；tgaxhanbase.NumReceiveAntennas=PHYParameters.numrxantenas；tgaxhanbase.TransmitReceiveDistance=10；%非直瞄距离（米）tgaxhanbase.ChannelBandwidth=cfgHE{txix}.ChannelBandwidth；tgaxhanbase.LargeScaleFadingEffect=“没有”；tgaxChanBase.EnvironmentalSpeed=0；%米/秒，静止tgaxhanbase.SampleRate=fs；tgaxChanBase.NormalizeChannelOutputs=false；%存储在单元阵列中，并重置通道以生成新的%回应tgaxhan{i}=tgaxhanbase；复位（tgaxChan{i}）；终止对于txevent=1:SimParameters.NumTxEventsPerDrop fprintf('正在运行传输事件#%d/%d…\n'，txevent，SimParameters.numtxeventsperdop）；%通过清晰通道评估确定有源发射机和接收机[activeTx，activeRx]=tgaxCCA（信号强度，MACParameters.CCALevel）；%情节场景和链接如果showScenarioPlot tgaxUpdatePlot（hGrid、txPositions、rxPositions、activeTx、activeRx、不可忽略的掩码、txChannels、rxChannels、，...斯普林特('物理层系统级模拟，丢弃#%d/%d，传输事件#%d/%d',...drop，SimParameters.NumDrops，txevent，SimParameters.numtxeventsperdop））；终止%使用activeTx和activeRx提取活动链接的元素cfgHEActive=cfgHE（activeTx）；associationActive=association（activeTx，activeRx）；nonnegiblemaskactive=nonnegiblemask（activeTx，activeRx）；signalStrengthActive=signalStrength（activeTx，activeRx）；%创建包含活动通道的阵列tgaxhanactive=cell（size（associationActive））；matchIdx=nonnegibleidx==find（activeTx&activeRx'）；tgaxhanactive（nonnegiblemaskactive）=tgaxChan（any（matchIdx，1））；%为每个不可忽略的活动链路生成波形，并%组合每个接收器的波形[rxWavs、TXPSDAUActive、signalPower、Interpower]=TGAXGeneratorX波形(...CFG活性、TGA活性、不可忽略的活性、信号强度活性、联想活性）；%对每个链路运行PHY链路模拟，并确定数据包%已成功解码。估计干扰功率%被传送到接收器以代替干扰电源%测量算法。numActiveRxs=sum（activeRx）；pktError=false（numActiveRxs，1）；sinrMeas=nan（numActiveRxs，1）；对于rxIdx=1:numActiveRxs[pktError（rxIdx），sinrMeas（rxIdx）]=tgaxModelPHYLink(...rxWavs{rxIdx}、cfgHEActive{rxIdx}、rxNoisePower、interpower（rxIdx）、txpsduative{rxIdx}）；终止%为有源接收机存储输出output.numPktErrors（activeRx）=output.numPktErrors（activeRx）+PKERROR；output.numPkts（activeRx）=output.numPkts（activeRx）+1；output.sinrMeas（activeRx，txevent，drop）=sinrMeas；%计算每个接收机的预期SINRsinrEst=10*log10（信号功率/（干扰功率+10^（（rxNoisePower-30）/10））；output.sinrEst（activeRx、txevent、drop）=sinrEst；终止终止%计算平均数据包错误率output.pktErrorRate=sum（output.numpkterrorres）/sum（output.numPkts）；disp(“模拟完成”)disp(['发射机的平均数据包错误率：'num2str（output.pkterrorate）]；终止rng（种子）；%恢复随机状态