5G蜂窝网络基础仿真

(1) 蜂窝网络拓扑生成

% 定义蜂窝网络参数
numCells = 7;          % 7个六边形小区
numUEsPerCell = 10;    % 每个小区10个用户
cellRadius = 500;      % 小区半径（米）

% 生成六边形蜂窝布局
txsites = txsite('Name','HexGrid', ...
    'AntennaHeight', 25, ...
    'TransmitterPower', 40);  % 基站参数

% 生成基站位置（六边形布局）
[bsPositions, cellLabels] = hHexGrid(numCells, cellRadius);

% 随机生成用户位置（均匀分布）
uePositions = hGenerateRandomUEs(bsPositions, cellRadius, numUEsPerCell);

% 可视化网络拓扑
show(bsPositions); 
hold on;
scatter(uePositions(:,1), uePositions(:,2), 'filled');
title('5G蜂窝网络布局');

 (2) 5G NR信号生成与链路级仿真

% 配置5G NR参数
carrier = nrCarrierConfig('NSizeGrid', 52, 'SubcarrierSpacing', 15); % 52 RB, 15 kHz子载波间隔
pdsch = nrPDSCHConfig('Modulation', 'QPSK', 'NumLayers', 1);         % QPSK调制

% 生成5G NR波形
waveform = nrWaveformGenerator(carrier, pdsch);

% 添加信道模型（TDL-C信道）
channel = nrTDLChannel('DelayProfile', 'TDL-C', 'SampleRate', carrier.SampleRate);
rxWaveform = channel(waveform);

% 解调与误码率计算
[~, ber] = nrPDSCHDecode(carrier, pdsch, rxWaveform);
disp(['误码率 (BER): ', num2str(ber)]);

2. 大规模MIMO波束成形仿真

(1) 波束成形权重计算

% 配置大规模MIMO天线阵列
numTxAntennas = 64;   % 64天线基站
numRxAntennas = 2;    % 用户设备2天线

% 生成信道矩阵H（基于空间信道模型）
H = (randn(numRxAntennas, numTxAntennas) + 1i*randn(numRxAntennas, numTxAntennas))/sqrt(2);

% 计算SVD波束成形权重
[U, S, V] = svd(H);
precodingWeights = V(:,1:numRxAntennas);  % 主特征模式

% 应用预编码
txSignal = precodingWeights * dataSymbols;

(2) 波束方向图可视化 

% 定义天线阵列参数
antennaArray = phased.URA('Size', [8 8], 'ElementSpacing', 0.5); % 8x8 URA

% 计算波束方向图
freq = 28e9;    % 毫米波频段28 GHz
beamformer = phased.PhaseShiftBeamformer('SensorArray', antennaArray, ...
    'OperatingFrequency', freq, 'Direction', [30; 0]); % 指向30度方位角

% 可视化波束
pattern(antennaArray, freq, 'PropagationSpeed', physconst('LightSpeed'), ...
    'CoordinateSystem', 'polar', 'Type', 'powerdb');

3. 5G信道建模

(1) CDL信道模型配置

% 配置CDL信道参数（3GPP TR 38.901）
cdl = nrCDLChannel('DelayProfile', 'CDL-D', ...  % 密集城区场景
    'CarrierFrequency', 3.5e9, ...              % 3.5 GHz频段
    'SampleRate', 30.72e6, ...                  % 采样率
    'TransmitAntennaArray', struct('Size', [8 8]), ...  % 8x8 URA
    'ReceiveAntennaArray', struct('Size', [2 2]));      % 2x2 URA

% 生成信道冲激响应
[pathGains, sampleTimes] = cdl();

(2) 信道容量分析

% 计算MIMO信道容量（基于香农公式）
SNR_dB = 20;           % 信噪比（dB）
SNR = 10^(SNR_dB/10);
H = (randn(numRxAntennas, numTxAntennas) + 1i*randn(numRxAntennas, numTxAntennas))/sqrt(2);
capacity = log2(det(eye(numRxAntennas) + SNR/numTxAntennas * (H * H')));

disp(['信道容量: ', num2str(capacity), ' bps/Hz']);

4. 网络级性能评估

(1) 小区间干扰分析

% 计算用户接收SINR
noisePower = 1e-12;     % 噪声功率
interferencePower = calculateInterference(uePositions, bsPositions); % 自定义干扰计算函数
SINR = signalPower ./ (interferencePower + noisePower);

% 绘制SINR热力图
heatmap(SINR, 'Title', '用户SINR分布');

(2) 吞吐量仿真 

% 基于CQI（信道质量指示）的吞吐量计算
cqiTable = nrCQITable();                     % 5G CQI表
sinrValues = 0:2:30;                         % SINR范围（dB）
throughput = zeros(size(sinrValues));

for i = 1:length(sinrValues)
    cqi = getCQIFromSINR(sinrValues(i));      % 自定义CQI映射函数
    [modulation, codeRate] = getMCSFromCQI(cqi, cqiTable); % 获取MCS参数
    throughput(i) = calculateThroughput(modulation, codeRate); % 计算吞吐量
end

plot(sinrValues, throughput);
xlabel('SINR (dB)'); ylabel('Throughput (Mbps)');

5. 关键工具箱与资源

• MATLAB工具箱:

  • 5G Toolbox: 5G NR协议栈仿真

  • Communications Toolbox: 信号处理与调制

  • Phased Array System Toolbox: 天线阵列与波束成形

  • Wireless Network Simulator: 网络级仿真