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)