clear N5=4; Nr=16; SNR=6; SNR= 10~(SNR/10); simulation=100; capacity0fAuer=[]; for Lr=1:16 capacity0fSum=6; for sim=1:simulation H=sqrt(1/2)*(randn(Nr,Ns)+1jxrandn(Nr,Ns));fullAntenna=[1:Hr]; B=eye(Ns,Ns); Alpha=[]; H_sel=[]; for j=1:Nr f=H(j，:)； h=f'; alpha=h'*h; Alpha=[Alpha alpha]; end for n=1:Lr [max0fAlpha,index]=max(Alpha)； fullAntenna(index)=[]; H_sel=[H_sel;H(index,:)]; if (n<Lr) f=H(index,:); h=f'; alpha=Alpha(index); a=(Bxh)/sqrt((Ns/SNR)+alpha); B=B-a a'; Alpha(index)=[]; for k=1:length( fullAntenna) Alpha(k)=Alpha(k)-(abs(a'*h))^2; for k=1:length( fullAntenna) Alpha(k)=Alpha(k)-(abs(a'*h))^2; end end end capacity0fSelected=1og2(det(eye(Ns)+SNR/NS*(H_sel'*H_sel)))；capacity0fSum=capacity0fSum+capacity0fSelected; I end capacity0fAver=[capacity0fAuer,capacity0fSum/sinulation];end plot(1:Lr,capacity0fAver); xlabel('Lr'); ylabel('capacity(bit/s/Hz)'); grid on; hold on写出这串matlab代码的注释

这段 Matlab 代码是一个无线通信系统的模拟程序，用于计算基于选择性反馈的无线系统中不同天线选择数量下的容量。以下是代码的注释：

% 清除变量
clear N5=4; Nr=16; SNR=6;

% 将 SNR 从 dB 转换为线性值
SNR= 10^(SNR/10);

% 模拟次数
simulation=100;

% 初始化变量
capacity0fAuer=[];

% 循环枚举天线选择数量
for Lr=1:16
    capacity0fSum=6;
    
    % 多次模拟求平均容量
    for sim=1:simulation
        % 生成随机信道矩阵 H
        H=sqrt(1/2)*(randn(Nr,Ns)+1jxrandn(Nr,Ns));
        
        % 选取完整天线集
        fullAntenna=[1:Hr];
        
        % 初始化选择后的天线集合、选择后的反馈矩阵 B、以及每个天线的贡献量 Alpha
        H_sel=[];
        B=eye(Ns,Ns);
        Alpha=[];
        for j=1:Nr
            f=H(j,:)； h=f';
            alpha=h'*h;
            Alpha=[Alpha alpha];
        end
        
        % 逐步选择天线
        for n=1:Lr
            % 选择贡献最大的天线
            [max0fAlpha,index]=max(Alpha)；
            fullAntenna(index)=[];
            H_sel=[H_sel;H(index,:)];
            
            % 更新反馈矩阵 B 和 Alpha
            if (n<Lr)
                f=H(index,:);
                h=f';
                alpha=Alpha(index);
                a=(Bxh)/sqrt((Ns/SNR)+alpha);
                B=B-a a';
                Alpha(index)=[];
                for k=1:length(fullAntenna)
                    Alpha(k)=Alpha(k)-(abs(a'*h))^2;
                    for k=1:length(fullAntenna)
                        Alpha(k)=Alpha(k)-(abs(a'*h))^2;
                    end
                end
            end
        end
        
        % 计算选择后的容量
        capacity0fSelected=1og2(det(eye(Ns)+SNR/NS*(H_sel'*H_sel)))；
        capacity0fSum=capacity0fSum+capacity0fSelected;
    end
    
    % 计算平均容量并保存结果
    capacity0fAver=[capacity0fAuer,capacity0fSum/sinulation];
end

% 绘制容量曲线
plot(1:Lr,capacity0fAver);
xlabel('Lr');
ylabel('capacity(bit/s/Hz)');
grid on;
hold on;