DL-SCH和PDSCH发送和接收处理链

打开生活的脚本

本例展示了如何使用5G Toolbox™功能建模5G NR物理下行共享通道(PDSCH)链路，包括从传输块生成到接收端比特解码的所有步骤。

介绍

此图显示下行共享信道(DL-SCH)和PDSCH发送和接收处理链。

1. png

这个例子展示了如何对链接级模拟的这些元素建模。

DL-SCH编码
混合ARQ (HARQ)管理
PDSCH编码
多输入多输出(MIMO)预编码
OFDM调制
传播信道和噪声相加
时间同步
OFDM解调
信道估计和均衡
PDSCH解码
DL-SCH解码

关于如何使用链接级模拟来度量吞吐量的例子，请参见NR PDSCH吞吐量。

仿真参数

指定信噪比(SNR)，要模拟的槽数，以及完善的信道估计标志。要了解本例中使用的信噪比定义的更多信息，请参见在链路模拟中使用的信噪比定义。

SNRdB = 10;% SNR in dBtotalNoSlots = 20;%要模拟的槽数perfectEstimation = false;%完美同步和信道估计rng (“默认”);%为可重复性设置默认随机数生成器

载波配置

创建carrier配置对象。该对象控制数字命理，如子载波间距、带宽和循环前缀(CP)长度。本例使用默认的属性集。

载体= nrCarrierConfig

carrier = nrCarrierConfig with properties: NCellID: 1 SubcarrierSpacing: 15 CyclicPrefix: 'normal' NSizeGrid: 52 NStartGrid: 0 NSlot: 0 NFrame: 0 Read-only properties: SymbolsPerSlot: 14 SlotsPerSubframe: 1 SlotsPerFrame: 10

PDSCH和DM-RS配置

创建PDSCH配置对象。指定调制方案(16-QAM)和层数(2)，将所有资源块(RBs)分配给PDSCH(全带分配)。还可以在该对象中指定其他时间分配参数和解调参考信号(DM-RS)设置。

pdsch = nrPDSCHConfig;pdsch。调制=“16 qam”;pdsch。NumLayers = 2; pdsch.PRBSet = 0:carrier.NSizeGrid-1;%全频带分配

显示PDSCH参数。

pdsch

pdsch = nrPDSCHConfig with properties: NSizeBWP: [] NStartBWP: [] ReservedPRB: {[1x1 nrPDSCHReservedConfig]} ReservedRE: [] Modulation: '16QAM' NumLayers: 2 MappingType: 'A' SymbolAllocation: [0 14] PRBSet:[0 12 3 4 5 6 7 8 9 11 12 13 14 15 16 17…[VRBToPRBInterleaving: 0 VRBBundleSize: 2 NID: [] RNTI: 1 DMRS: [1x1 nrPDSCHDMRSConfig] EnablePTRS: 0 PTRS: [1x1 nrPDSCHPTRSConfig]只读属性:NumCodewords: 1

设置DM-RS参数。为了提高信道估计，增加一个额外的DM-RS位置。

pdsch.DMRS.DMRSAdditionalPosition = 1;

设置DM-RS配置类型和DM-RS长度，决定正交DM-RS序列或DM-RS端口的数量。

DMRSConfigurationType = 1支持最多4个DM-RS端口时DMRSLength = 1。
DMRSConfigurationType = 1时支持最多8个DM-RS端口DMRSLength = 2。
DMRSConfigurationType = 2支持最多6个DM-RS端口时DMRSLength = 1。这是为多用户MIMO (MU-MIMO)而设计的。
DMRSConfigurationType = 2支持多达12个DM-RS端口时DMRSLength = 2。这是为MU-MIMO设计的。

最大层数必须小于等于DM-RS端口数。

pdsch.DMRS.DMRSConfigurationType = 1;pdsch.DMRS.DMRSLength = 2;pdsch。dmr%显示DM-RS属性

ans = nrPDSCHDMRSConfig with properties: DMRSConfigurationType: 1 DMRSReferencePoint: 'CRB0' DMRSAdditionalPosition: 2 DMRSAdditionalPosition: 1 DMRSLength: 2 CustomSymbolSet: [] DMRSPortSet: [] NIDNSCID: 0 NumCDMGroupsWithoutData: 2 DMRSDownlinkR16: 0 Read-only properties: CDMGroups: [0 0] DeltaShifts: [0 0] FrequencyWeights: [2x2 double] TimeWeights: [2x2 double] DMRSSubcarrierLocations: [6x2 double] cdmlength: [2 1]

DL-SCH配置

指定码率，HARQ进程的数量，冗余版本(RV)序列值。该序列在发生错误时控制RV重传。要禁用HARQ重传，您可以设置rvSeq到一个固定值(例如，0)。有关如何用HARQ建模传输通道的更多信息，请参见带有HARQ的5G NR传输通道模型。

NHARQProcesses = 16;%并行HARQ进程的数量rvSeq = [0 2 3 1];

在指定编码率时，要考虑码字的数量。码字数是PDSCH配置对象的一个只读属性，取决于层数。

1码字，最多4层
2码字超过4层

%的编码率如果pdsch。NumCodewords == 1 codeRate = 490/1024;其他的codeRate = [490 490]./1024;结束

创建DL-SCH编码器和解码器对象。若要使用多个进程，请设置MultipleHARQProcesses财产真正的两个对象。不需要指定HARQ进程的数量。DL-SCH编码器和解码器对象最多可以建模16个HARQ进程。在使用DL-SCH编码器和解码器对象执行操作时，要识别活动的HARQ进程，请使用HARQprocessID属性的HARQ实体对象，在下一节中定义。

%创建DL-SCH编码器对象encodeDLSCH = nrDLSCH;encodeDLSCH。MultipleHARQProcesses = true;encodeDLSCH。TargetCodeRate = codeRate;%创建DLSCH解码器对象decodeDLSCH = nrDLSCHDecoder;decodeDLSCH。MultipleHARQProcesses = true;decodeDLSCH。TargetCodeRate = codeRate;decodeDLSCH。LDPCDecodingAlgorithm =“规范化min-sum”;decodeDLSCH。MaximumLDPCIterationCount = 6;

HARQ管理

创建一个HARQ实体对象来管理HARQ进程和DL-SCH编码器和解码器缓冲区。对于每个HARQ进程，一个HARQ实体存储以下元素:

HARQ ID号码。
房车。
传输数，表示某个传输块被传输了多少次。
标志，指示是否需要新数据。当成功接收到传输块或发生序列超时(所有RV传输都失败了)时，需要新数据。
标志，指示是否发生了序列超时(所有RV传输都失败了)。

harqEntity = harqEntity (0: NHARQProcesses-1、rvSeq pdsch.NumCodewords);

通道配置

指定发射和接收天线的数量。

nTxAnts = 8;nRxAnts = 8;%检查层数对于天线数是否有效如果pdsch。NumLayers > min(nTxAnts,nRxAnts) error(“层数(”字符串(pdsch.NumLayers) + +)必须小于min(nTxAnts,nRxAnts) ("+字符串(min (nTxAnts nRxAnts) +“)”)结束

创建一个通道对象。

频道= nrTDLChannel;通道。DelayProfile =“TDL-C”;通道。NumTransmitAntennas = nTxAnts; channel.NumReceiveAntennas = nRxAnts;

设置信道采样率为OFDM信号的采样率。获取OFDM信号的采样率，使用nrOFDMInfo函数。

ofdmInfo = nrOFDMInfo(载体);通道。SampleRate = ofdmInfo.SampleRate;

传输和接收

设置一个循环，模拟槽位的传输和接收。创建一个comm.ConstellationDiagram显示均衡信号的星座。

constPlot = comm.ConstellationDiagram;星座图对象constPlot。ReferenceConstellation = getConstellationRefPoints (pdsch.Modulation);%参考星座值constPlot。EnableMeasurements = 1;%启用EVM测量%初始定时偏移抵消= 0;estChannelGrid = getInitialChannelEstimate(渠道、载体);newPrecodingWeight = getPrecodingMatrix (pdsch.PRBSet pdsch.NumLayers estChannelGrid);为nSlot = 0: totalNoSlots-1%的新位置母舰。NSlot = nSlot;

计算传输块大小

传输块大小是发送到信道编码阶段的比特数。这个值取决于PDSCH的容量。要计算传输块大小，请使用nrTBS函数。

%生成PDSCH索引信息，这是计算传输所需的%的块大小[pdschIndices, pdschInfo] = nrPDSCHIndices(载体,pdsch);%计算传输块大小Xoh_PDSCH = 0;trBlkSizes = nrTBS (pdsch.Modulation、pdsch.NumLayers元素个数(pdsch.PRBSet) pdschInfo.NREPerPRB, codeRate, Xoh_PDSCH);

HARQ处理(缓冲区管理)

本节解释编码器和解码器中的缓冲区管理。

DL-SCH编码器缓冲区:如果活动HARQ进程需要新数据，则生成一个新的传输块。将传输块存储在相应的缓冲区中。如果不需要新的数据，DL-SCH编码器使用它的缓冲位进行重传。
DL-SCH译码器缓冲区:接收器中的软缓冲区存储以前接收到的相同码字的版本。这些缓冲区在成功接收后自动清除(没有CRC错误)。然而，如果RV序列结束时没有成功解码，则通过使用resetSoftBuffer对象的功能。

%根据需要获取新的传输块和刷新解码器软缓冲区为cwIdx = 1: pdsch。NumCodewords如果harqEntity.NewData (cwIdx)创建并存储用于传输的新传输块。trBlk = randi([0 1]，trBlkSizes(cwIdx)，1);setTransportBlock (encodeDLSCH trBlk、cwIdx-1 harqEntity.HARQProcessID);%如果之前的RV序列没有成功结束%解码，刷新软缓冲区如果harqEntity.SequenceTimeout (cwIdx) resetSoftBuffer (decodeDLSCH, cwIdx-1 harqEntity.HARQProcessID);结束结束结束

DL-SCH编码

对传输块进行编码。传输块已经存储在DL-SCH编码器对象的一个内部软缓冲区中。

codedTrBlock = encodeDLSCH (pdsch.Modulation pdsch.NumLayers、pdschInfo.G harqEntity.RedundancyVersion, harqEntity.HARQProcessID);

PDSCH调制和MIMO预编码

从编码的传输块生成PDSCH符号。

pdschSymbols = nrPDSCH(载体,pdsch codedTrBlock);

获取预编码权值。本例假设信道知识用于预编码。(有关如何在接收端使用信道估计来计算下一个插槽中用于传输的权重的示例，请参见NR PDSCH吞吐量。)

precodingWeights = newPrecodingWeight;

预编码PDSCH符号。

pdschSymbolsPrecoded = pdschSymbols * precodingWeights;

PDSCH DM-RS代

生成DM-RS符号和索引。

dmrsSymbols = nrPDSCHDMRS(载体、pdsch);dmrsIndices = nrPDSCHDMRSIndices(载体、pdsch);

映射到资源网格

生成一个空的资源网格。这个网格表示一个槽。

pdschGrid = nrResourceGrid(载体、nTxAnts);

的nrPDSCHIndices函数生成的索引引用的是层，而不是天线。当将PDSCH符号直接映射到层时，这种格式很有用。在这种情况下，得到的资源网格没有预先编码。这张图显示了PDSCH符号映射到与层一样多的资源网格的过程。

2. png

因为本例在将PDSCH符号映射到资源网格之前，将MIMO预编码应用到PDSCH符号，因此MIMO预编码的PDSCH符号指的是天线而不是层。要将层索引转换为天线索引，请使用nrExtractResources函数。该图显示了mimo预编码符号映射到与发射天线一样多的资源网格的过程。

3. png

[~, pdschAntIndices] = nrExtractResources (pdschIndices pdschGrid);pdschGrid (pdschAntIndices) = pdschSymbolsPrecoded;

mimo -预编码并将DM-RS符号映射到资源网格。与PDSCH索引类似，DM-RS索引指的是层。要将这些层指标转换为天线指标，请使用nrExtractResources函数了。

% PDSCH DM-RS预编码和映射为p = 1:size(dmrsSymbols,2) [~，dmrsAntIndices] = nrExtractResources(dmrsIndices(:，p)，pdschGrid);pdschGrid(dmrsAntIndices) = pdschGrid(dmrsAntIndices) + dmrsSymbols(:，p)*precodingWeights(p，:);结束

OFDM调制

ofdm -调制资源网格。

[txWaveform, waveformInfo] = nrOFDMModulate(载体,pdschGrid);

传播渠道

传播通道生成N输入的输出样本N样本。然而，块N输出样本包括信道滤波器瞬态(K样本)。因为同步阶段删除了这个初始瞬态，如果一个槽在通道输出N样品,n - k同步后样本仍保留。n - k样本不足以解码一个槽值的数据。部分槽值样本位于通道滤波器延迟线中，尚未被冲洗。要将所有相关样本从通道滤波器中冲洗出来，可以用零填充输入信号。信道滤波器引入的最大延时影响填充的大小。填充占了所有多径组件引入的延迟和信道滤波器实现延迟。这张图显示了在波形进入信道之前需要零填充。

4. png

用足够的零填充输入信号，以确保生成的信号从通道滤波器中冲洗出来。

chInfo =信息(渠道);maxChDelay = ceil(max(chinfo . pathdelayed *channel.SampleRate)) + chInfo.ChannelFilterDelay;txWaveform = [txWaveform;0 (maxChDelay、尺寸(txWaveform 2)));

通过频道发送信号并添加噪声。

[rxWaveform, pathGains sampleTimes] =通道(txWaveform);噪音= generateAWGN (SNRdB nRxAnts waveformInfo.Nfft,大小(rxWaveform));rx波形= rx波形+噪声;

时间同步

您可以执行完美的或实际的同步。

完美同步假设信道知识(nrPerfectTimingEstimate)。通道返回关于路径增益和路径滤波器脉冲响应的信息。您可以使用该信息来确定与跨所有通道快照和跨所有发射和接收天线的最强多路径组件相关的偏移量。
实际同步将接收到的信号与时域的PDSCH DM-RS符号进行相互关联(nrTimingEstimate)。在一些不利的情况下，由于衰落或噪声，这种相互关系可能很弱，导致错误的时间偏移。这个函数hSkipWeakTimingOffset检查互相关的大小玛格。如果互相关性较弱，则函数忽略当前的定时估计，而使用之前的估计(抵消)。

执行完美或实际的时间估计和同步。

如果perfectEstimation获取路径过滤器以进行完美的时间估计pathFilters = getPathFilters(渠道);(抵消,mag) = nrPerfectTimingEstimate (pathGains pathFilters);其他的(t, mag) = nrTimingEstimate(载体、rxWaveform dmrsIndices, dmrsSymbols);抵消= hSkipWeakTimingOffset(抵消t杂志);结束rxWaveform = rxWaveform(1 +抵消:最终,);

OFDM解调

ofdm -解调同步信号。

rxGrid = nrOFDMDemodulate(载体、rxWaveform);

信道估计

信道估计提供了每个资源元素(RE)的信道效应的表示。均衡器使用这个信息来补偿由信道引入的失真。

可以执行完美的或实用的通道估计。

完美信道估计假设信道知识(nrPerfectChannelEstimate)。完美信道估计表示发射天线和接收天线之间的信道状况。因为均衡器需要发射层和接收天线之间的信道知识，所以必须对完美信道估计应用预编码。
实际信道估计使用PDSCH DM-RS来估计信道条件，并使用噪声平均和插值来获得槽内所有REs的估计。由于DM-RSs是每层指定的，因此得到的实际信道估计表示发射层和接收天线之间的信道条件。实际信道估计包括MIMO预编码操作的影响。

该图显示了下行链路processing链中通道估计的参考点。

5. png

执行完美或实用的信道估计。

如果perfectEstimation%在发射和接收之间执行完美的信道估计%的天线。estChGridAnts = nrPerfectChannelEstimate(载体,pathGains pathFilters,抵消,sampleTimes);%获得完美噪声估计(来自噪声实现)noiseGrid = nrOFDMDemodulate(载波，噪声(1+offset:end，:)));噪音= var (noiseGrid (:));%获取下一个槽位的预编码矩阵newPrecodingWeight = getPrecodingMatrix (pdsch.PRBSet pdsch.NumLayers estChGridAnts);%对estchgri唐应用预编码。得到的估计为%层与接收天线之间的信道估计。estChGridLayers = precodeChannelEstimate (estChGridAnts precodingWeights。');其他的%在层和接收之间执行实际的信道估计%的天线。[estChGridLayers,噪音]= nrChannelEstimate(载体、rxGrid dmrsIndices, dmrsSymbols,“CDMLengths”, pdsch.DMRS.CDMLengths);%在预编码前从estChannelGrid中删除预编码%矩阵计算estChGridAnts = precodeChannelEstimate (estChGridLayers连词(precodingWeights));%获取下一个槽位的预编码矩阵newPrecodingWeight = getPrecodingMatrix (pdsch.PRBSet pdsch.NumLayers estChGridAnts);结束

绘制第一层和第一接收天线之间的信道估计。

网格(abs (estChGridLayers (:,:, 1, 1)));

图中包含一个axes对象。标题为Channel Estimate的axes对象包含一个类型为surface的对象。

标题(信道估计的);包含(OFDM符号的);ylabel (“副载波”);zlabel (“级”);

此时，可以使用信道估计获得下一个槽中传输的预编码矩阵。因为本例假设发射端有信道知识，所以不需要在接收端计算预编码矩阵。关于如何根据接收端的信道估计计算数据传输的预编码矩阵的示例，请参见NR PDSCH吞吐量。

均衡

均衡器使用信道估计来补偿由信道引入的失真。

从接收到的网格和相关信道估计中提取PDSCH符号。的csi输出对每个均衡的PDSCH符号都有信道状态信息(CSI)。CSI是每个PDSCH符号的信道条件的度量。使用CSI对PDSCH解码后解码的软位进行加权，有效地增加了经历更好信道条件的符号的重要性。

[pdschRx, pdschHest] = nrExtractResources (pdschIndices、rxGrid estChGridLayers);[pdschEq, csi] = nrEqualizeMMSE (pdschRx pdschHest,噪音);

绘制均衡符号的星座图。该图包括所有层的星座图。

constPlot。ChannelNames =“层”+ (pdsch.NumLayers: 1:1);constPlot。ShowLegend = true;%的星座为第一层有更高的信噪比%的最后一层。翻转图层，这样星座就不会蒙版%。constPlot (fliplr (pdschEq));

PDSCH解码

解码均衡化PDSCH符号，获得软位码字。

[dlschLLRs, rxSymbols] = nrPDSCHDecode(载体、pdsch pdschEq,噪音);

根据CSI来衡量软位或对数似然比(LLRs)。这种缩放对具有较好信道条件的REs中的符号应用更大的权重。

%根据CSI缩放llrcsi = nrLayerDemap (csi);% CSI层demapping为cwIdx = 1: pdsch。NumCodewordsQm = length(dlschLLRs{cwIdx})/length(rxSymbols{cwIdx});%比特/符号csi {cwIdx} = repmat (csi {cwIdx}。”、Qm 1);%按每个符号的每位展开dlschLLRs{cwIdx} = dlschLLRs{cwIdx} .* csi{cwIdx}(:);%的规模结束

DL-SCH解码

解码llr并检查错误。

decodeDLSCH。TransportBlockLength = trBlkSizes;[decbits, blkerr] = decodeDLSCH (dlschLLRs、pdsch.Modulation pdsch.NumLayers,…harqEntity.RedundancyVersion harqEntity.HARQProcessID);

HARQ进程更新

用得到的块错误状态更新当前HARQ进程，然后推进到下一个进程。此步骤更新HARQ实体中与活动HARQ进程相关的信息。

statusReport = updateAndAdvance (harqEntity blkerr、trBlkSizes pdschInfo.G);

总结当前槽位的HARQ和解码信息。

disp (“槽”+ (nSlot) +“。”+ statusReport);

槽0。HARQ Proc 0: CW0:初始传输通过(RV=0,CR=0.482212)。槽位1。HARQ Proc 1: CW0:初始传输通过(RV=0,CR=0.482212)。槽2。HARQ Proc 2: CW0:初始传输通过(RV=0,CR=0.482212)。位置3。HARQ Proc 3: CW0:初始传输通过(RV=0,CR=0.482212)。槽4。HARQ Proc 4: CW0:初始传输通过(RV=0,CR=0.482212)。 Slot 5. HARQ Proc 5: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 6. HARQ Proc 6: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 7. HARQ Proc 7: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 8. HARQ Proc 8: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 9. HARQ Proc 9: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 10. HARQ Proc 10: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 11. HARQ Proc 11: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 12. HARQ Proc 12: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 13. HARQ Proc 13: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 14. HARQ Proc 14: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 15. HARQ Proc 15: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 16. HARQ Proc 0: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 17. HARQ Proc 1: CW0: Initial transmission passed (RV=0,CR=0.482212). Slot 18. HARQ Proc 2: CW0: Initial transmission passed (RV=0,CR=0.482212).

槽19。HARQ Proc 3: CW0:初始传输通过(RV=0,CR=0.482212)。

结束% for nSlot = 0:totalNoSlots

本地函数

函数噪音= generateAWGN (SNRdB nRxAnts、Nfft sizeRxWaveform)%为给定的SNR in dB (SNRDB)值生成AWGN，即每RE和天线的接收机信噪比%，假设信道是这样的%不影响信号的功率。NRXANTS是接收的数量%的天线。NFFT是OFDM解调中使用的FFT大小。SIZERXWAVEFORM%是接收波形的大小，用于计算的大小%噪声矩阵。由OFDM调制中使用的IFFT大小归一化噪声功率，如OFDM调制器将此归一化应用于传输的%波形。也由接收天线的数量归一化，如%通道模型将此归一化应用到接收的波形通过%违约。信噪比定义为每个接收天线(TS% 38.101 - 4)。信噪比= 10 ^ (SNRdB / 10);%计算线性噪声增益N0 = 1 /√(2.0 * nRxAnts *双(Nfft) *信噪比);噪音= N0 *复杂(randn (sizeRxWaveform) randn (sizeRxWaveform));结束函数wtx = getPrecodingMatrix (PRBSet NLayers hestGrid)给定分配和信道估计，计算预编码矩阵%已分配子载波指数allocSc = (1:12)' + 12*PRBSet(:).';allocSc = allocSc (:);%平均通道估计(~, ~, R, P) = (hestGrid)大小;estAllocGrid = hestGrid (allocSc、::,:);Hest = permute(mean(重塑(estAllocGrid，[]，R,P))，[2 3 1]);%计算分解[~ ~ V] =圣言(命令);wtx = V (:, 1: NLayers)。”;wtx = mtx /√(NLayers);% Normalize by NLayers结束函数estChannelGrid = getInitialChannelEstimate(渠道、载体)%获取用于计算预编码矩阵的初始信道估计。%这个函数假设一个完美的信道估计%通道的克隆chClone = channel.clone ();chClone.release ();%无需过滤即可获得信道路径增益chClone。ChannelFiltering = false;%获取通道路径增益[pathGains, sampleTimes] = chClone ();%完美的定时同步pathFilters = getPathFilters (chClone);抵消= nrPerfectTimingEstimate (pathGains pathFilters);%完美信道估计estChannelGrid = nrPerfectChannelEstimate(载体,pathGains pathFilters,抵消,sampleTimes);结束函数refPoints = getConstellationRefPoints (mod)%计算给定调制的参考星座点%的计划。开关国防部情况下“正交相移编码”《不扩散核武器条约》= 4;情况下“16 qam”《不扩散核武器条约》= 16;情况下“64 qam”《不扩散核武器条约》= 64;情况下“256 qam”《不扩散核武器条约》= 256;结束binaryValues = int2bit (0: nPts-1 log2 (npt));refPoints = nrSymbolModulate (binaryValues (:), mod);结束函数estChannelGrid = precodeChannelEstimate (estChannelGrid W)%对信道估计的最后一个维度应用预编码矩阵W。%将4-D矩阵线性化，并在乘法后重塑K =大小(estChannelGrid, 1);L =大小(estChannelGrid, 2);R =大小(estChannelGrid 3);estChannelGrid =重塑(estChannelGrid K * L * R, []);estChannelGrid = estChannelGrid * W;estChannelGrid =重塑(estChannelGrid K, L, R, []);结束