发送端未知信道状态信息情况下MIMO系统容量仿真

程序功能:在发送端与接收端天线个数不同的情况下，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);