程序功能:在发送端与接收端天线个数不同的情况下,MIMO系统的容量(SNR=2*(1+序号mod 5) )
首先,我们定义了发送端和接收端天线的个数,并计算了对应的SNR值。然后,在每个SNR和天线配置的情况下,运行多次仿真来计算容量。最后,计算平均容量并绘制MIMO系统容量SNR曲线。
% 定义 MIMO 系统参数
M = 2; % 发送天线数
N = 2; % 接收天线数
SN=55;%序号
% 定义信道参数
numIterations = 1000; % 仿真迭代次数
SNR_dB = 0:20; % 信噪比范围(dB)
% 初始化容量数组
capacity = zeros(size(SNR_dB));
for snrIdx = 1:length(SNR_dB)
snr = 10^(SNR_dB(snrIdx)/10); % dB转换为线性单位
% 初始化每个 SNR 下的容量累加器
capacitySum = 0;
for iter = 1:numIterations
% 生成复高斯随机信道
H = sqrt(0.5)*(randn(N,M) + 1i*randn(N,M));
% 计算每个接收天线的信噪比
rxSNR = snr / N;
% 计算每个接收天线上的信号功率
rxPower = rxSNR;
% 计算每个接收天线的接收信号
rxSignal = sqrt(rxPower) * H;
% 计算接收信号的容量
capacitySum = capacitySum + log2(det(eye(N) + rxSignal' * rxSignal / rxPower));
end
% 取平均容量
capacity(snrIdx) = capacitySum / numIterations;
end
% 绘制容量随信噪比变化的仿真曲线
figure;
plot(SNR_dB, capacity, 'b-o');
xlabel('SNR (dB)');
ylabel('平均容量');
title('MIMO系统平均容量随信噪比的仿真曲线');
% 计算 SNR=2*(SN mod 5)dB 时的平均容量
targetSNR_dB = 2*mod(SN,5); % 根据题目要求计算 SNR 值
targetIdx = find(SNR_dB == targetSNR_dB); % 找到对应的索引
targetCapacity = capacity(targetIdx); % 取出对应的容量值
fprintf('SNR = %d dB 时的平均容量为: %.4f\n', targetSNR_dB, targetCapacity);