5g nr downlink aclr measurement -凯发k8网页登录

% select the nr-tm waveform parameters
nrtm = "nr-fr1-tm1.2";  % nr-tm name and properties
bw   = "20mhz";  % channel bandwidth
scs  = "15khz";  % subcarrier spacing
dm   = "fdd";  % duplexing mode
% create generator object for the above nr-tm
tmwavegen = hnrreferencewaveformgenerator(nrtm,bw,scs,dm);
% ensure no windowing to highlight impact of filtering on aclr
tmwavegen = makeconfigwritable(tmwavegen);
tmwavegen.config.windowingpercent = 0;    
% generate waveform
[tmwaveform,tmwaveinfo] = generatewaveform(tmwavegen);
samplingrate = tmwaveinfo.info.samplingrate;   % waveform sampling rate (hz)
% visualize the associated prb and subcarrier resource grids
displayresourcegrid(tmwavegen);

aclrparameters = haclrparametersnr(tmwavegen.config);
disp(aclrparameters);

% design filter
lpfilt = designfilt('lowpassfir',...
    'passbandfrequency',aclrparameters.measurementbandwidth/2,...
    'stopbandfrequency',aclrparameters.bandwidth/2,...
    'passbandripple',0.1,...
    'stopbandattenuation',80,...
    'samplerate',samplingrate);
% apply filter
filttmwaveform = filter(lpfilt,tmwaveform);

% apply required oversampling
resampled = resample(tmwaveform,aclrparameters.osr,1);           % not filtered
filtresampled = resample(filttmwaveform,aclrparameters.osr,1);   % filtered
% create and configure a memoryless nonlinearity to model the amplifier
nonlinearity = comm.memorylessnonlinearity;
nonlinearity.method = 'rapp model';
nonlinearity.smoothness = 3;              % p parameter          
nonlinearity.lineargain = 0.5;            % db
nonlinearity.outputsaturationlevel = 2;   % it limits the output signal level  
% signal conditioning to control the hpa input back-off level
resampled = resampled/max(abs(resampled));                % not filtered
filtresampled = filtresampled/max(abs(filtresampled));    % filtered
% apply the amplifier model to the nr waveforms
txwaveform = nonlinearity(resampled);           % not filtered
txfiltwaveform = nonlinearity(filtresampled);   % filtered

% calculate nr aclr
aclr = haclrmeasurementnr(aclrparameters,txwaveform);               % not filtered
filtaclr = haclrmeasurementnr(aclrparameters,txfiltwaveform);   % filtered

% evm configuration parameters  
evmcfg.plotevm = false;
evmcfg.samplerate = aclrparameters.samplingrate;
evmcfg.label = tmwavegen.configuredmodel{1};
% measure the evm related statistics for the transmitted waveform without filtering
evminfo = hnrdownlinkevm(tmwavegen.config,txwaveform,evmcfg)

% measure the evm related statistics for the transmitted waveform with filtering
evminfofilt = hnrdownlinkevm(tmwavegen.config,txfiltwaveform,evmcfg);

haclrresultsnr(aclr,aclrparameters,txwaveform,'(not filtered)');

haclrresultsnr(filtaclr,aclrparameters,txfiltwaveform,'(filtered)');

function haclrresultsnr(aclr, aclrparameters, waveform, arg)
    minaclr = 45;
    % 4th input argument is plot title qualifier
    if nargin > 3
        titletext = [' ' arg];
    elseif nargin > 1
        titletext = [];
    else
        titletext = [];
        waveform = [];
    end
    
    % aclr values and ticks for bar chart
    values = round([aclr(1:end/2) 0 aclr(end/2 1:end)],1);
    tick = 1:numel(values);
    ticklabel = tick-ceil(numel(tick)/2);
    labelvec = tick;
    labelvec(ceil(end/2)) = []; % do not plot label for 0db aclr on channel
    
    % plot nr spectrum
    if ~isempty(waveform)
        figure;
        [spectrum, frequency] = pwelch(waveform, kaiser(8192*4,19), [], [], ...
            aclrparameters.samplingrate, 'centered', 'power');
        frequency = frequency * 10^(-6); % mhz
        spectrum = 10*log10(spectrum / max(spectrum));
        adjacentchannellabel = [ticklabel(1:floor(length(ticklabel)/2)) ...
        ticklabel(floor(length(ticklabel)/2) 2:end)];
        % select 'x' and 'y' limits to show the adjacent channels in the plot
        xlimitright = aclrparameters.adjacentchanneloffset   (aclrparameters.measurementbandwidth/2);
        xlimitright = xlimitright * 10^(-6); % mhz
        xlimitleft = aclrparameters.adjacentchanneloffset - (aclrparameters.measurementbandwidth/2);
        xlimitleft = xlimitleft * 10^(-6); % mhz
        ylimits = [min(spectrum)-20 max(spectrum) 10];
        ylim(ylimits);
        xlim([min(frequency) max(frequency)])
        hold on; 
        for i = 1:length(aclrparameters.adjacentchanneloffset)
            patch('xdata',[xlimitright(i) xlimitright(i) xlimitleft(i) ...
                xlimitleft(i)],'ydata', [ylimits fliplr(ylimits)], ...
                'facecolor','y','facealpha',0.2) % plot adjacent channels
            text(aclrparameters.adjacentchanneloffset(i)*10^(-6), i, sprintf('%d', ...
                adjacentchannellabel(i)), 'horizontalalignment', 'center', ...
                'verticalalignment', 'top'); % plot adjacent channel labels
        end
        plot(frequency, spectrum);
        hold off;
        xlabel('frequency (mhz)');
        ylabel('normalized power (db)');
        title(strcat ('nr spectrum', titletext));
        legend('adjacent channels', 'location', 'southeast')
    end
    
    % plot nr aclr    
    figure;
    hold on;
    yline(minaclr,'r'); 
    bar(values, 'basevalue', 0, 'facecolor', 'yellow');
    hold off;
    set(gca, 'xtick', tick, 'xticklabel', ticklabel, 'ylim', ...
        [0 0.2*max(values) max(values)]);
    for i = labelvec
        text(i, values(i), sprintf('%0.1f db',values(i)), ...
            'horizontalalignment', 'center', 'verticalalignment', 'top');
    end
    title(strcat('nr adjacent channel leakage ratio', titletext));
    xlabel('adjacent channel offset');
    ylabel('adjacent channel leakage ratio (db)');
    legend('minimum required aclr');
    
end