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

pucch = {nrWavegenPUCCH3Config};                   % Create a PUCCH format 3 configuration object for the first UE
pucch{1}.Enable = 1;                               % Enable PUCCH sequence
pucch{1}.Label = 'UE 1 - PUCCH Format 3 @ 15 kHz'; % Label for this PUCCH sequence
pucch{1}.BandwidthPartID = 1;                      % BWP of PUCCH transmission
pucch{1}.Power = 0;                                % Power scaling in dB

pucch{1}.Modulation = 'QPSK';            % 'pi/2-BPSK','QPSK'
pucch{1}.FrequencyHopping = 'intraSlot'; % Frequency hopping configuration
pucch{1}.SecondHopStartPRB = 10;         % Resource block offset for second hop
pucch{1}.GroupHopping = 'enable';        % Group hopping configuration
pucch{1}.HoppingID = 1;                  % Hopping identity
pucch{1}.RNTI = 11;                      % RNTI for the first UE
pucch{1}.NID = 0;                        % Scrambling identity

pucch{1}.SymbolAllocation = [3 11]; % First symbol and length
pucch{1}.SlotAllocation = [3 4];    % Allocated slots indices for PUCCH sequence
pucch{1}.Period = 6;                % Allocation period in slots
pucch{1}.PRBSet = 0:9;              % PRB allocation

pucch{1}.AdditionalDMRS = 1;        % Additional DM-RS
pucch{1}.DMRSPower = 1;             % Additional power boosting for DM-RS in dB

pucch{1}.Coding = 1;
pucch{1}.TargetCodeRate = 0.15;
pucch{1}.NumUCIBits = 20;
pucch{1}.NumUCI2Bits = 10;
pucch{1}.DataSourceUCI = 'PN9';
pucch{1}.DataSourceUCI2 = 'PN9';

pucch{2} = nrWavegenPUCCH2Config;                  % Create a PUCCH format 2 configuration object for the second UE
pucch{2}.Label = 'UE 2 - PUCCH Format 2 @ 30 kHz'; % Label for this PUCCH sequence
pucch{2}.BandwidthPartID = 2;                      % PUCCH mapped to 2nd BWP
pucch{2}.SymbolAllocation = [10 2];                % Symbol allocation
pucch{2}.SlotAllocation = 0:2;                     % Slot allocation
pucch{2}.Period = [];                              % Specify no repetitions of the slot pattern for this PUCCH
pucch{2}.RNTI = 12;                                % RNTI for the second UE
pucch{2}.NID0 = 0;                                 % DM-RS scrambling identity

pucch{3} = nrWavegenPUCCH0Config;                  % Create a PUCCH format 0 configuration object for the third UE
pucch{3}.Label = 'UE 3 - PUCCH Format 0 @ 30 kHz'; % Label for this PUCCH sequence
pucch{3}.BandwidthPartID = 2;                      % PUCCH mapped to 2nd BWP
pucch{3}.SymbolAllocation = [1 2];                 % Symbol allocation
pucch{3}.PRBSet = 40;                              % PRB allocation
pucch{3}.FrequencyHopping = 'interSlot';           % Frequency hopping
pucch{3}.SecondHopStartPRB = 30;                   % Resource block offset for second hop
pucch{3}.InitialCyclicShift = 3;                   % Initial cyclic shift
pucch{3}.NumUCIBits = 2;                           % Number of UCI bits containing HARQ-ACK
pucch{3}.DataSourceSR = 1;                         % SR data source

pusch = {nrWavegenPUSCHConfig};
pusch{1}.Enable = 0;
pusch{1}.Label = 'PUSCH @ 15 kHz';
pusch{1}.BandwidthPartID = 1;

srs = {nrWavegenSRSConfig};
srs{1}.Enable = 0;
srs{1}.Label = 'SRS @ 15 kHz';
srs{1}.BandwidthPartID = 1;

waveconfig.SCSCarriers = scscarriers;
waveconfig.BandwidthParts = bwp;
waveconfig.PUCCH = pucch;
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('Information associated with the first PUCCH:')
disp(info.WaveformResources.PUCCH(1))

disp('Modulation information associated with BWP 1:')
disp(info.ResourceGrids(1).Info)