5g nr downlink vector waveform generation -凯发k8网页登录

waveconfig = nrdlcarrierconfig;    % create a downlink carrier configuration object
waveconfig.label = 'dl carrier 1'; % label for this downlink waveform configuration
waveconfig.ncellid = 0;            % cell identity
waveconfig.channelbandwidth = 40;  % channel bandwidth (mhz)
waveconfig.frequencyrange = 'fr1'; % 'fr1' or 'fr2'
waveconfig.numsubframes = 10;      % number of 1 ms subframes in generated waveform (1, 2, 4, 8 slots per 1 ms subframe, depending on scs)
waveconfig.windowingpercent = 0;   % percentage of windowing relative to fft length
waveconfig.samplerate = [];        % sample rate of the ofdm modulated waveform
waveconfig.carrierfrequency = 0;   % carrier frequency in hz. this property is used for symbol phase
                                   % compensation before ofdm modulation
% define a set of scs specific carriers, using the maximum sizes for a
% 40 mhz nr channel. see ts 38.101-1 for more information on defined
% bandwidths and guardband requirements
scscarriers = {nrscscarrierconfig,nrscscarrierconfig};
scscarriers{1}.subcarrierspacing = 15;
scscarriers{1}.nsizegrid = 216;
scscarriers{1}.nstartgrid = 0;
scscarriers{2}.subcarrierspacing = 30;
scscarriers{2}.nsizegrid = 106;
scscarriers{2}.nstartgrid = 1;

% ss burst configuration
ssburst = nrwavegenssburstconfig;
ssburst.enable = 1;                     % enable ss burst
ssburst.power = 0;                      % power scaling in db
ssburst.blockpattern = 'case b';        % case b (30khz) subcarrier spacing
ssburst.transmittedblocks = [1 1 1 1];  % bitmap indicating blocks transmitted in a 5ms half-frame burst
ssburst.period = 20;                    % ss burst set periodicity in ms (5, 10, 20, 40, 80, 160)
ssburst.ncrbssb = [];                   % frequency offset of ss burst (crb), use [] for the waveform center

% bwp configurations
bwp = {nrwavegenbwpconfig,nrwavegenbwpconfig};
bwp{1}.bandwidthpartid = 1;        % bwp id
bwp{1}.label = 'bwp 1 @ 15 khz';   % label for this bwp
bwp{1}.subcarrierspacing = 15;     % bwp subcarrier spacing
bwp{1}.cyclicprefix = 'normal';    % bwp cyclic prefix for 15 khz
bwp{1}.nsizebwp = 25;              % size of bwp in prbs
bwp{1}.nstartbwp = 12;             % position of bwp, relative to point a, in crbs
bwp{2}.bandwidthpartid = 2;        % bwp id
bwp{2}.label = 'bwp 2 @ 30 khz';   % label for this bwp
bwp{2}.subcarrierspacing = 30;     % bwp subcarrier spacing
bwp{2}.cyclicprefix = 'normal';    % bwp cyclic prefix for 30 khz
bwp{2}.nsizebwp = 50;              % size of bwp in prbs
bwp{2}.nstartbwp = 51;             % position of bwp, relative to point a, in crbs

% coreset and search space configurations
coresets = {nrcoresetconfig};
coresets{1}.coresetid = 1;                      % coreset id
coresets{1}.duration = 3;                       % coreset symbol duration (1,2,3)
coresets{1}.frequencyresources = [1 1 0 1];     % bitmap indicating blocks of 6 prb for coreset (rrc - frequencydomainresources)
coresets{1}.cceregmapping = 'noninterleaved';   % mapping: 'interleaved' or 'noninterleaved'
coresets{1}.regbundlesize = 3;                  % l (2,6) or (3,6)
coresets{1}.interleaversize = 2;                % r (2,3,6)
coresets{1}.shiftindex = waveconfig.ncellid;    % set to ncellid
searchspaces = {nrsearchspaceconfig};
searchspaces{1}.searchspaceid = 1;              % search space id
searchspaces{1}.coresetid = 1;                  % coreset associated with this search space
searchspaces{1}.searchspacetype = 'ue';         % search space type, 'ue' or 'common'
searchspaces{1}.slotperiodandoffset = [5 0];    % allocated slot period and slot offset of search space pattern
searchspaces{1}.duration = 2;                   % number of slots in the block of slots in pattern period
searchspaces{1}.startsymbolwithinslot = 0;      % first symbol of each coreset monitoring opportunity in a slot
searchspaces{1}.numcandidates = [8 8 4 2 0];    % number of candidates at each al (set to 0 if the al doesn't fit in coreset)

pdcch = {nrwavegenpdcchconfig};
pdcch{1}.enable = 1 ;                     % enable pdcch sequence
pdcch{1}.label = 'ue 1 - pdcch @ 15 khz'; % label for this pdcch sequence
pdcch{1}.bandwidthpartid = 1;             % bandwidth part of pdcch transmission
pdcch{1}.power = 1.1;                     % power scaling in db
pdcch{1}.coding = 1;                      % enable dci coding
pdcch{1}.searchspaceid = 1;               % search space
pdcch{1}.slotallocation = 0;              % allocated slots indices for pdcch sequence
pdcch{1}.period = 5;                      % allocation period in slots
pdcch{1}.aggregationlevel = 8;            % aggregation level (1,2,4,8,16 cces)
pdcch{1}.allocatedcandidate = 1;          % pdcch candidate in search space (1 based)
pdcch{1}.rnti = 11;                       % rnti
pdcch{1}.dmrsscramblingid = 1;            % pdcch and dm-rs scrambling nid
pdcch{1}.dmrspower = 0;                   % additional dm-rs power boosting in db
pdcch{1}.datablocksize = 20;              % dci payload size
pdcch{1}.datasource = 'pn9';              % dci data source

pdsch = {nrwavegenpdschconfig};           % create a pdsch configuration object for the first ue
pdsch{1}.enable = 1;                      % enable pdsch sequence
pdsch{1}.label = 'ue 1 - pdsch @ 15 khz'; % label for this pdsch sequence
pdsch{1}.bandwidthpartid = 1;             % bandwidth part of pdsch transmission
pdsch{1}.power  = 0;                      % power scaling in db
pdsch{1}.coding = 1;                      % enable the dl-sch transport channel coding
pdsch{1}.datasource = 'pn9';              % channel data source
pdsch{1}.targetcoderate = 0.4785;         % code rate used to calculate transport block sizes
pdsch{1}.xoverhead = 0;                   % rate matching overhead
pdsch{1}.modulation = 'qpsk';             % 'qpsk', '16qam', '64qam', '256qam'
pdsch{1}.numlayers = 2;                   % number of pdsch layers
pdsch{1}.rvsequence = [0 2 3 1];          % rv sequence to be applied cyclically across the pdsch allocation sequence
pdsch{1}.vrbtoprbinterleaving = 0;        % disable interleaved resource mapping
pdsch{1}.vrbbundlesize = 2;               % vrb-toprb-interleaver parameter

pdsch{1}.symbolallocation = [2 9];    % first symbol and length
pdsch{1}.slotallocation = 0:9;        % allocated slot indices for pdsch sequence
pdsch{1}.period = 15;                 % allocation period in slots
pdsch{1}.prbset = [0:5, 10:20];       % prb allocation
pdsch{1}.rnti = 11;                   % rnti for the first ue
pdsch{1}.nid = 1;                     % scrambling for data part

pdsch{1}.reservedcoreset = 1;           % rate matching pattern, defined by coreset ids
pdsch{1}.reservedprb{1}.prbset = [];    % rate matching pattern, defined by set of prb (rrc 'bitmaps')
pdsch{1}.reservedprb{1}.symbolset = [];
pdsch{1}.reservedprb{1}.period = [];

% antenna port and dm-rs configuration (ts 38.211 section 7.4.1.1)
pdsch{1}.mappingtype = 'a';                % pdsch mapping type ('a'(slot-wise),'b'(non slot-wise))
pdsch{1}.dmrspower = 0;                    % additional power boosting in db
pdsch{1}.dmrs.dmrsconfigurationtype = 2;   % dm-rs configuration type (1,2)
pdsch{1}.dmrs.numcdmgroupswithoutdata = 1; % number of dm-rs cdm groups without data. the value can be one of the set {1,2,3}
pdsch{1}.dmrs.dmrsportset = [];            % dm-rs antenna ports used ([] gives port numbers 0:numlayers-1)
pdsch{1}.dmrs.dmrstypeaposition = 2;       % mapping type a only. first dm-rs symbol position (2,3)
pdsch{1}.dmrs.dmrslength = 1;              % number of front-loaded dm-rs symbols (1(single symbol),2(double symbol))
pdsch{1}.dmrs.dmrsadditionalposition = 0;  % additional dm-rs symbol positions (max range 0...3)
pdsch{1}.dmrs.nidnscid = 1;                % scrambling identity (0...65535)
pdsch{1}.dmrs.nscid = 0;                   % scrambling initialization (0,1)

% pt-rs configuration (ts 38.211 section 7.4.1.2)
pdsch{1}.enableptrs = 0;             % enable or disable the pt-rs (1 or 0)
pdsch{1}.ptrspower = 0;              % additional pt-rs power boosting in db
pdsch{1}.ptrs.timedensity = 1;       % time density (l_pt-rs) of pt-rs (1,2,4)
pdsch{1}.ptrs.frequencydensity = 2;  % frequency density (k_pt-rs) of pt-rs (2,4)
pdsch{1}.ptrs.reoffset = '00';       % pt-rs resource element offset ('00','01','10','11')
pdsch{1}.ptrs.ptrsportset = 0;       % pt-rs antenna ports must be a subset of dm-rs ports

pdsch{2} = pdsch{1};                   % create a pdsch configuration object for the second ue
pdsch{2}.enable = 1;
pdsch{2}.label = 'ue 2 - pdsch @ 30 khz';
pdsch{2}.bandwidthpartid = 2;          % pdsch mapped to the second bwp
pdsch{2}.rnti = 12;                    % rnti for the second ue
pdsch{2}.symbolallocation = [0 12];
pdsch{2}.slotallocation = [2:4, 6:20];
pdsch{2}.prbset  = [25:30, 35:38];     % prb allocation, relative to bwp

csirs = {nrwavegencsirsconfig};
csirs{1}.enable = 0;
csirs{1}.label = 'csi-rs @ 15 khz';
csirs{1}.bandwidthpartid = 1;
csirs{1}.power = 3; % power scaling in db

csirs{1}.csirstype = {'nzp','zp'};
csirs{1}.rownumber = [3 5];
csirs{1}.density = {'one','one'};
csirs{1}.subcarrierlocations = {6 4};
csirs{1}.numrb = 25;
csirs{1}.rboffset = 12;
csirs{1}.symbollocations = {13 9};
csirs{1}.csirsperiod = {[5 0], [5 0]};
csirs{1}.nid = 5;

csirs{2} = nrwavegencsirsconfig;
csirs{2}.enable = 0;
csirs{2}.label = 'csi-rs @ 30 khz';
csirs{2}.bandwidthpartid = 2;
csirs{2}.power = 3; % power scaling in db
csirs{2}.csirstype = {'nzp','nzp'};
csirs{2}.rownumber = [1 1];
csirs{2}.density = {'three','three'};
csirs{2}.subcarrierlocations = {0 0};
csirs{2}.numrb = 50;
csirs{2}.rboffset = 50;
csirs{2}.symbollocations = {6 10};
csirs{2}.csirsperiod = {[10 1], [10 1]};
csirs{2}.nid = 0;

waveconfig.ssburst = ssburst;
waveconfig.scscarriers = scscarriers;
waveconfig.bandwidthparts = bwp;
waveconfig.coreset = coresets;
waveconfig.searchspaces = searchspaces;
waveconfig.pdcch = pdcch;
waveconfig.pdsch = pdsch;
waveconfig.csirs = csirs;
% generate complex baseband waveform
[waveform,info] = nrwaveformgenerator(waveconfig);

figure;
plot(abs(waveform));
title('magnitude of 5g downlink baseband waveform');
xlabel('sample index');
ylabel('magnitude');

samplerate = info.resourcegrids(1).info.samplerate;
nfft = info.resourcegrids(1).info.nfft;
figure;
spectrogram(waveform(:,1),ones(nfft,1),0,nfft,'centered',samplerate,'yaxis','minthreshold',-130);
title('spectrogram of 5g downlink baseband waveform');

disp('modulation information associated with bwp 1:')
disp(info.resourcegrids(1).info)