Energy_efficient_hybridbeamforming.rar_Energy-Efficient_Precodin

%% simulation clear all clc close all NM = 64; % transmit antenna K = 16; % receive antenna N = 8; % RF chain number M = 8; % each RF chain connected number antenna L = 3; % the number of cluster S = 5; % the iterations times SNR_db = -30:5:0; SNR = 10.^(SNR_db./10); realization = 1000; N0 = 1; j = sqrt(-1); B = 1; for index = 1:length(SNR) ea=1; %每个天线发射的功率，也即信号向量中每个元素的功率 sigma_2=ea/SNR(index); %噪声功率 for itime = 1:realization H = generate_channel(NM,K,L); [U1 S1 V1] = svd(H); [D A] = SIC_based_precoding(H,N,M,K,sigma_2,S); %利用SIC方法其中特征值和特征向量没有用到SVD分解 [D1 A1] = optimal_precoding(H,N,M,K,sigma_2); %利用SIC方法其中特征值和特征向量用到SVD分解 P1 = A*D; P2 = A1*D1; W1 = SIC_based_precoding(H*P1,N,M,K,sigma_2,S); W2 = U(:,1:N); [U V S] = svd(H); Wopt = U(:,1:N); Fopt = S(:,1:N); Rn = inv(sigma_2*Wopt'*Wopt); capasity_C1(index,itime) = log2(real(det(eye(N)+1*Rn*W1'*H*P1*P1'*H'*W1))); %利用SIC方法其中特征值和特征向量没有用到SVD分解 capasity_C2(index,itime) = log2(real(det(eye(N)+1*Rn*W2'*H*P2*P2'*H'*W2))); %利用SIC方法其中特征值和特征向量用到SVD分解 capasity_C(index,itime) = log2(real(det(eye(N)+1*Rn*Wopt'*H*Fopt*Fopt'*H'*Wopt))); % 全数字波束成形 end end plot(SNR_db,mean(capasity_C'),'-k') hold on; plot(SNR_db,mean(capasity_C2'),'-r') plot(SNR_db,mean(capasity_C1'),'-b') grid on; axis([-30 0 0 25]); xlabel('SNR(dB)'); ylabel('容量/(bit/s/Hz)');