复苏的IEEE 802.15.4z超宽频信号

rng (7);psduTx =兰迪([0,1],1016,1);cfgTx = = lrwpanHRPConfig(模式“HPRF”PSDULength =长度(psduTx) / 8, SamplesPerPulse = 2);waveTx = lrwpanWaveformGenerator (psduTx cfgTx);

负载uwbCapture.mat情节(真实(capturedHPRF),“这”)

散点图(capturedHPRF)

rxSPC = 2;%样本/芯片/脉冲,鉴于receive-side A / DFs = 499.2 e6 * rxSPC;coarseFC = comm.CoarseFrequencyCompensator(调制=“帕姆”SampleRate = Fs, FrequencyResolution = 1 e3);[coarseFCsignal,频率]= coarseFC (capturedHPRF);流(“粗频率偏移:%。2 f MHz \ n”1、频率/ e6)

散点图(coarseFCsignal)

fineFC = comm.CarrierSynchronizer(调制=“帕姆”SamplesPerSymbol = rxSPC);[fineFCsignal,阶段]= fineFC (coarseFCsignal);散点图(fineFCsignal)

fineFCsignalReal =实际(fineFCsignal);symbolSync = comm.SymbolSynchronizer(调制=PAM /相移键控/ QAM的SamplesPerSymbol = rxSPC TimingErrorDetector =早期的后期(non-data-aided)”);recoveredHPRF = symbolSync (fineFCsignalReal);

散点图(recoveredHPRF)

Tpream = 5;%手册,迭代规范;阈值可以增加直到1代码检测Tnoise = 0.01;%用于能量检测,对规范序言探测器的阈值Tzero = 25/100;% 0符号振幅阈值,最高比例的振幅,应该高于噪音cfgRcv = = lrwpanHRPConfig(模式“HPRF”SamplesPerPulse = rxSPC);%创建一个对象跟踪received-frame属性preambleCodesHPRF = 25:32;preambleDet = comm.PreambleDetector;syncFound = false;codeIdx = preambleCodesHPRF (1) 1;hprfPreambleDurations = [16、24、32、48、64、96、128、256];而~ syncFound & & codeIdx <马克斯(preambleCodesHPRF)%考虑所有代码,直到检测codeIdx = codeIdx + 1;为极性= [1]%考虑所有极性,直到检测%从传播创建同步代码:代码=极性* lrwpan.internal.HRPCodes (codeIdx);spreadingFactorHPRF = cfgRcv.PreambleSpreadingFactor;%扩散因子,L(总是4 HPRF模式)序言= 0(长度(代码)* spreadingFactorHPRF, 1);序言(1:spreadingFactorHPRF:结束)=代码;释放(preambleDet);preambleDet。序言=序言;meanCorr =意味着(abs(过滤器(flipud(序言),1,recoveredHPRF (recoveredHPRF > Tnoise))));preambleDet。阈值= Tpream * meanCorr;[preamPos, corrMat] = preambleDet (recoveredHPRF);durationIdx =找到(元素个数(preamPos) > hprfPreambleDurations, 1,“最后一次”);如果~ isempty (durationIdx) cfgRcv。PreambleDuration = hprfPreambleDurations (durationIdx);syncFound = true;打破;%中发现,没有必要去探索其他的极性结束结束结束如果syncFound流(“发现同步代码# % d。”,codeIdx);cfgRcv。CodeIndex = codeIdx;%跟踪确认帧特征其他的错误(“不同步字段中输入数据。);结束

%片剩余获取与集成和转储(帧长度是未知,现在,在PHR)recoveredHPRF = recoveredHPRF / max (abs (recoveredHPRF));%正常化振幅frameStart = preamPos(1)长度(序言)+ 1;ternarySymbols =极性* recoveredHPRF (frameStart:结束);ternarySymbols (abs (ternarySymbols) < Tzero) = 0;ternarySymbols (ternarySymbols > Tzero) = 1;ternarySymbols (ternarySymbols < -Tzero) = 1;

%检测陕西林业局通过识别匹配chip-by-chip基础上的95%berThr = 1 - 95/100;sfdFound = false;sfdStart = 1 +(序言)* cfgRcv.PreambleDuration长度;为sfdNum = 0:4从表15-7c %陕西林业局值cfgRcv。SFDNumber = sfdNum;陕西林业局= lrwpan.internal.getSFD (cfgRcv);spreadedSFD =陕西林业局。*中*极性;spreadedSFD = spreadedSFD (:);bitErrs =总和(spreadedSFD ~ = ternarySymbols (sfdStart: sfdStart +长度(spreadedSFD) 1));sfdBER = bitErrs /长度(spreadedSFD);如果sfdBER < berThr sfdFound = true;打破;结束结束如果sfdFound流(“发现陕西林业局# % d。”,sfdNum);其他的警告(“同步后没有发现陕西林业局。”)结束

numChipsPerGap = 512;%部分15.2.9 [2]。sfdEnd = sfdStart +长度(spreadedSFD) - 1;gapStarts = [];idx = sfdEnd + 1;而idx < = (ternarySymbols) -numChipsPerGap长度+ 1如果~任何(ternarySymbols (idx: idx + numChipsPerGap-1))%长512芯片的差距如果~ isempty (gapStarts) & & idx = = gapStarts(结束)+ 512%的框架发现,不视为STS差距gapStarts(结束)= [];ternarySymbols = ternarySymbols (1: idx);打破;结束gapStarts = [gapStarts idx];% #好< AGROW >idx = idx + numChipsPerGap;%的差距继续;结束%到最后一个非零,避免毫无意义的比较nextIncr =找到(ternarySymbols (idx + 1: idx + numChipsPerGap), 1,“最后一次”);如果isempty (nextIncr) nextIncr = numChipsPerGap;结束idx = idx + nextIncr + 1;结束stsAfterSFD = ~ isempty (gapStarts) & & gapStarts (1) = = sfdEnd + 1;如果cfgRcv isempty (gapStarts)。STSPacketConfiguration = 0;其他的cfgRcv。NumSTSSegments = max(1,元素个数(gapStarts) 1);% 1段如果帧结束的差距cfgRcv。STSPacketConfiguration = 0 + (cfgRcv.NumSTSSegments>0) + ~stsAfterSFD;% STSPacketConfiguration = 3不支持万博1manbetx结束%最后差距不开始在结束的象征,因为过去的l - 1符号无足轻重如果stsAfterSFD cfgRcv。STSSegmentLength = (diff(gapStarts(1:2))-numChipsPerGap + spreadingFactorHPRF-1)/numChipsPerGap;% 512年单位芯片elseif~ isempty (gapStarts) cfgRcv。STSSegmentLength = (length(ternarySymbols)-gapStarts -numChipsPerGap + spreadingFactorHPRF-1)/numChipsPerGap;% 512年单位芯片结束如果cfgRcv。STSPacketConfiguration fprintf(“发现% d STS段长度% d (s);STS % d。\ n”包配置,cfgRcv。NumSTSSegments cfgRcv。STSSegmentLength, cfgRcv.STSPacketConfiguration) stsLen = numChipsPerGap*(cfgRcv.NumSTSSegments+1) + cfgRcv.NumSTSSegments * cfgRcv.STSSegmentLength * 512;其他的流(“没有发现STS。\ n”);stsLen = 0;结束

如果~ stsAfterSFD phrStart = sfdEnd + 1;其他的phrStart = sfdEnd + 1 + stsLen;结束%需要找到调制类型(124.8 vs 249.6 MHz)% Guarbands区分因素(分别为75%和50%的符号)numPHRBits = 19;minSymbolDuration = 16;%在芯片中,所有的数据速率、编码的组合% PHR将跨度至少19 x16芯片。检查范围内% guarband统计:numZeros =总和(0 = = ternarySymbols (phrStart: phrStart + numPHRBits * minSymbolDuration-1));zerosPct = numZeros / (numPHRBits * minSymbolDuration);mode249MHz = zerosPct < (50 + (75 - 50) / 2) / 100;如果mode249MHz cfgRcv。MeanPRF = 249.6;其他的cfgRcv。MeanPRF = 124.8;结束%,确定卷积编码使用约束% 3或7长度,使用不同的符号映射(表15-10c-f [2])。为CL = 7 [3] cfgRcv。ConstraintLength = CL;isPHR = true;[cwPHR, phrEnd] = helperUWBHPRFDemod (isPHR、ternarySymbols phrStart, cfgRcv);[secdedPass, PSDULength] = helperUWBPHRDecode (cwPHR cfgRcv, CL);如果secdedPass cfgRcv。PSDULength = PSDULength;打破;结束结束

payloadStart = phrEnd + 1;decodedPSDU = helperUWBPayloadDecode (ternarySymbols payloadStart、cwPHR cfgRcv);[~,误码率]= biterr (psduTx decodedPSDU);流(“比特误码率:% 0.2 f \ n ',ber)