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

waveconfig = nrulcarrierconfig;    % create an uplink carrier configuration object
waveconfig.label = 'ul carrier 1'; % label for this uplink 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;

% 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 = 10;             % 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 = 51;              % size of bwp in prbs
bwp{2}.nstartbwp = 40;             % position of bwp, relative to point a, in crbs

pusch = {nrwavegenpuschconfig};           % create a pusch configuration object for the first ue
pusch{1}.enable = 1;                      % enable pusch sequence
pusch{1}.label = 'ue 1 - pusch @ 15 khz'; % label for this pusch sequence
pusch{1}.bandwidthpartid = 1;             % bwp of pusch transmission
pusch{1}.power = 0;                       % power scaling in db
pusch{1}.coding = 1;                      % enable the ul-sch transport channel coding
pusch{1}.nid = 1;                         % scrambling for data part
pusch{1}.rnti = 11;                       % rnti for the first ue
pusch{1}.transformprecoding = false;      % transform precoding
pusch{1}.targetcoderate = 0.47;           % code rate used to calculate transport block sizes
pusch{1}.xoverhead = 0;                   % rate matching overhead
% transmission settings
pusch{1}.transmissionscheme = 'codebook'; % 'codebook','noncodebook'
pusch{1}.modulation = 'qpsk';             % 'pi/2-bpsk','qpsk','16qam','64qam','256qam'
pusch{1}.numlayers = 2;                   % number of pusch layers
pusch{1}.numantennaports = 4;             % number of antenna ports
pusch{1}.tpmi = 0;                        % transmitted precoding matrix indicator (0...27)
pusch{1}.rvsequence = [0 2 3 1];          % rv sequence to be applied cyclically across the pusch allocation sequence
pusch{1}.frequencyhopping = 'interslot';  % frequency hopping configuration
pusch{1}.secondhopstartprb = 10;          % resource block offset for second hop
% data source
pusch{1}.datasource = 'pn9';              % channel data source

pusch{1}.mappingtype = 'a';         % pusch mapping type ('a'(slot-wise),'b'(non slot-wise))
pusch{1}.symbolallocation = [0 14]; % first symbol and length
pusch{1}.slotallocation = [0 1];    % allocated slots indices for pusch sequence
pusch{1}.period = 5;                % allocation period in slots
pusch{1}.prbset = 0:10;             % prb allocation

% antenna port and dm-rs configuration (ts 38.211 section 6.4.1.1)
pusch{1}.dmrspower = 0;                    % additional power boosting in db
pusch{1}.dmrs.dmrsconfigurationtype = 1;   % dm-rs configuration type (1,2)
pusch{1}.dmrs.numcdmgroupswithoutdata = 2; % number of dm-rs cdm groups without data. the value can be one of the set {1,2,3}
pusch{1}.dmrs.dmrsportset = [0 2];         % dm-rs antenna ports used ([] gives port numbers 0:numlayers-1)
pusch{1}.dmrs.dmrstypeaposition = 2;       % mapping type a only. first dm-rs symbol position (2,3)
pusch{1}.dmrs.dmrslength = 1;              % number of front-loaded dm-rs symbols (1(single symbol),2(double symbol))
pusch{1}.dmrs.dmrsadditionalposition = 2;  % additional dm-rs symbol positions (max range 0...3)
pusch{1}.dmrs.nidnscid = 1;                % scrambling identity for cp-ofdm (0...65535). use empty ([]) to use physical layer cell identity
pusch{1}.dmrs.nscid = 0;                   % scrambling initialization for cp-ofdm (0,1)
pusch{1}.dmrs.nrsid = 0;                   % scrambling identity for dft-s-ofdm dm-rs (0...1007). use empty ([]) to use physical layer cell identity
pusch{1}.dmrs.grouphopping = true;         % group hopping configuration. this property is used only when transform precoding is enabled
pusch{1}.dmrs.sequencehopping = false;     % sequence hopping configuration. this property is used only when transform precoding is enabled

% pt-rs configuration (ts 38.211 section 6.4.1.2)
pusch{1}.enableptrs = 0;            % enable or disable the pt-rs (1 or 0)
pusch{1}.ptrspower = 0;             % additional pt-rs power boosting in db for cp-ofdm
pusch{1}.ptrs.timedensity = 1;      % time density (l_pt-rs) of pt-rs (1,2,4)
pusch{1}.ptrs.frequencydensity = 2; % frequency density (k_pt-rs) of pt-rs for cp-ofdm (2,4)
pusch{1}.ptrs.numptrssamples = 2;   % number of pt-rs samples (ngroupsamp) for dft-s-ofdm (2,4)
pusch{1}.ptrs.numptrsgroups = 2;    % number of pt-rs groups (nptrsgroup) for dft-s-ofdm (2,4,8)
pusch{1}.ptrs.reoffset = '00';      % pt-rs resource element offset for cp-ofdm ('00','01','10','11')
pusch{1}.ptrs.ptrsportset = 0;      % pt-rs antenna ports must be a subset of dm-rs ports for cp-ofdm
pusch{1}.ptrs.nid = 0;              % pt-rs scrambling identity for dft-s-ofdm (0...1007)

pusch{1}.enableack = true;        % enable or disable harq-ack
pusch{1}.numackbits = 5;          % number of harq-ack bits
pusch{1}.betaoffsetack = 1;       % power factor of harq-ack
pusch{1}.datasourceack = 'pn9';   % harq-ack data source
pusch{1}.enablecsi1 = true;       % enable or disable csi part 1
pusch{1}.numcsi1bits = 10;        % number of csi part 1 bits
pusch{1}.betaoffsetcsi1 = 2;      % power factor of csi part 1
pusch{1}.datasourcecsi1 = 'pn9';  % csi part 1 data source
pusch{1}.enablecsi2 = true;       % enable or disable csi part 2
pusch{1}.numcsi2bits = 10;        % number of csi part 2 bits
pusch{1}.betaoffsetcsi2 = 2;      % power factor of csi part 2
pusch{1}.datasourcecsi2 = 'pn9';  % csi part 2 data source
pusch{1}.enablecguci = false;     % enable or disable cg-uci
pusch{1}.numcgucibits = 10;       % number of cg-uci bits
pusch{1}.betaoffsetcguci = 2;     % power factor of cg-uci
pusch{1}.datasourcecguci = 'pn9'; % cg-uci data source
pusch{1}.enableulsch = true;      % enable or disable ul-sch when there is uci transmission on pusch
pusch{1}.uciscaling = 1;          % scaling factor (0.5, 0.65, 0.8, 1)

pusch{2} = pusch{1};                 % create a pusch configuration object for the second ue
pusch{2}.enable = 1;
pusch{2}.label = 'ue 2 - pusch @ 30 khz';
pusch{2}.bandwidthpartid = 2;        % pusch mapped to the second bwp
pusch{2}.rnti = 12;                  % rnti for the second ue
pusch{2}.symbolallocation = [0 12];
pusch{2}.slotallocation = [5 6 7 8];
pusch{2}.prbset = 5:10;              % prb allocation, relative to bwp
pusch{2}.period = 10;
pusch{2}.transformprecoding = 1;
pusch{2}.frequencyhopping = 'interslot';
pusch{2}.numlayers = 1;
pusch{2}.rnti = 1;
pusch{2}.dmrs.grouphopping = false;
pusch{2}.dmrs.dmrsportset = 1;

srs = {nrwavegensrsconfig};
srs{1}.enable = 0;
srs{1}.label = 'srs @ 15 khz';
srs{1}.bandwidthpartid = 1;
srs{1}.power = 3; % power scaling in db

srs{1}.numsrsports = 1;             % number of srs ports (1,2,4)
srs{1}.numsrssymbols = 4;           % number of srs symbols in a slot (1,2,4)
srs{1}.symbolstart = 10;            % time-domain position of the srs in the slot. (8...13) for normal cp and (6...11) for extended cp
srs{1}.slotallocation = 2;          % allocated slots indices
srs{1}.period = 5;                  % allocation period in slots
srs{1}.frequencystart = 0;          % frequency position of the srs in bwp in rbs
srs{1}.nrrc = 0;                    % additional offset from freqstart specified in blocks of 4 prbs (0...67)
srs{1}.csrs = 13;                   % bandwidth configuration c_srs (0...63). it controls the allocated bandwidth to the srs
srs{1}.bsrs = 2;                    % bandwidth configuration b_srs (0...3). it controls the allocated bandwidth to the srs
srs{1}.bhop = 1;                    % frequency hopping configuration (0...3). set bhop < bsrs to enable frequency hopping
srs{1}.ktc = 2;                     % comb number (2,4). it indicates the allocation of the srs every ktc subcarriers
srs{1}.kbartc = 0;                  % subcarrier offset of the srs sequence (0...ktc-1)
srs{1}.cyclicshift = 0;             % cyclic shift number (0...ncsmax-1). ncsmax = 8 for ktc = 2 and ncsmax = 12 for ktc = 4.
srs{1}.repetition = 1;              % repetition factor (1,2,4). it indicates the number of equal consecutive srs symbols in a slot
srs{1}.groupseqhopping = 'neither'; % group or sequence hopping ('neither', 'grouphopping', 'sequencehopping')
srs{1}.nsrsid = 0;                  % scrambling identity (0...1023)
srs{1}.srspositioning = false;      % enable srs for user positioning

srs{2} = srs{1};
srs{2}.enable = 0;
srs{2}.label = 'srs @ 30 khz';
srs{2}.bandwidthpartid = 2;
srs{2}.numsrssymbols = 2;
srs{2}.symbolstart = 12;
srs{2}.slotallocation = [5 6 7 8];
srs{2}.period = 10;
srs{2}.bsrs = 0;
srs{2}.bhop = 0;

waveconfig.scscarriers = scscarriers;
waveconfig.bandwidthparts = bwp;
waveconfig.pusch = pusch;
waveconfig.srs = srs;
% generate complex baseband waveform
[waveform,info] = nrwaveformgenerator(waveconfig);

figure;
plot(abs(waveform));
title('magnitude of 5g uplink 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 uplink baseband waveform');

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