基于Matlab实现cell-free Massive MIMO生成用户平均误码率随平均接收信噪比变化曲线（源码）.rar

clear all; close all; %% System parameter settings par.squareLength = 2; %in kilometer par.L = 64; % number of APs par.K = 32; % user terminals (set not larger than MR!) par.N = 1; % number of antennas per AP par.mod = 'BPSK'; % modulation type: 'BPSK','QPSK','16QAM','64QAM' par.nbrOfRealizations = 1000; % number of channel realizations per setup par.nbrOfSetups = 1; % Number of Monte-Carlo setups, re-generate AP and UE locations per setup par.SNRdB_list = 5:5:20; % list of SNR [dB] values to be simulated par.tau_c = 200; % Length of coherence block par.tau_p = 30;% Length of pilot sequences % -- initialization: set up Gray-mapped constellation alphabet (according to IEEE 802.1switch (par.mod) switch (par.mod) case 'BPSK' par.symbols = [ -1 1 ]; case 'QPSK' par.symbols = [ -1-1i,-1+1i, ... +1-1i,+1+1i ]; case '16QAM' par.symbols = [ -3-3i,-3-1i,-3+3i,-3+1i, ... -1-3i,-1-1i,-1+3i,-1+1i, ... +3-3i,+3-1i,+3+3i,+3+1i, ... +1-3i,+1-1i,+1+3i,+1+1i ]; case '64QAM' par.symbols = [ -7-7i,-7-5i,-7-1i,-7-3i,-7+7i,-7+5i,-7+1i,-7+3i, ... -5-7i,-5-5i,-5-1i,-5-3i,-5+7i,-5+5i,-5+1i,-5+3i, ... -1-7i,-1-5i,-1-1i,-1-3i,-1+7i,-1+5i,-1+1i,-1+3i, ... -3-7i,-3-5i,-3-1i,-3-3i,-3+7i,-3+5i,-3+1i,-3+3i, ... +7-7i,+7-5i,+7-1i,+7-3i,+7+7i,+7+5i,+7+1i,+7+3i, ... +5-7i,+5-5i,+5-1i,+5-3i,+5+7i,+5+5i,+5+1i,+5+3i, ... +1-7i,+1-5i,+1-1i,+1-3i,+1+7i,+1+5i,+1+1i,+1+3i, ... +3-7i,+3-5i,+3-1i,+3-3i,+3+7i,+3+5i,+3+1i,+3+3i ]; end par.Es = mean(abs(par.symbols).^2);% extract average symbol energy par.Q = log2(length(par.symbols)); % number of bits per symbol par.bits = de2bi(0:length(par.symbols)-1,par.Q,'left-msb'); % precompute bit labels %% start simulation time_elapsed = 0; % track simulation time bits = randi([0 1],par.K,par.Q,par.nbrOfSetups,par.nbrOfRealizations);% generate random bit stream par.detector = {'CMMSE','CMMSE-SIC','CMF'}; % initialize result arrays (detector x SNR) res.BER = zeros(length(par.detector),length(par.SNRdB_list)); % bit error rate % Prepare to save simulation results %% Go through all setups: 每次Setup下AP位置和用户位置随机生成，由此导致不同的AP-用户关联、导频分配、信道矩阵等 for n = 1:par.nbrOfSetups %Generate one setup with UEs and APs at random locations [BETTA,R,pilotIndex,D,APpositions,UEpositions,distances] = generateSetup(par.L,par.K,par.N,par.tau_p,par.squareLength); figure(1) scatter(real(APpositions),imag(APpositions),'s');hold on; scatter(real(UEpositions),imag(UEpositions),'c'); grid on; xlabel('X Location (km)'); ylabel('Y Location (km)'); legend('Access Point','User Point'); for snr=1:length(par.SNRdB_list) % SNR loop tic; fprintf('Number of Setup: %3d; SNR: %3d\n',n,par.SNRdB_list(snr)); SNR_dB = par.SNRdB_list(snr); % Current SNR point in dBs SNR_lin = 10.^(SNR_dB/10); % Linear SNR N0 = par.Es*sum(sum(BETTA*par.N))/(par.L*par.N)/SNR_lin; [Hhat_tot,H_tot,B_tot,C_tot] = functionChannelEstimates(R,par.nbrOfRealizations,par.L,par.K,par.N,par.tau_p,pilotIndex,par.Es,N0,BETTA); noise = sqrt(N0)*sqrt(0.5)*(randn(par.N*par.L,par.nbrOfRealizations)+1i*randn(par.N*par.L,par.nbrOfRealizations)); for t = 1 : par.nbrOfRealizations fprintf('Number of Setup: %3d; SNR: %3d; nbrOfH: %3d\n',n,par.SNRdB_list(snr),t); idx = bi2de(bits(:,:,n,t),'left-msb')+1; % generate transmit symbol s = par.symbols(idx).'; H = reshape(H_tot(:,t,:),[par.N*par.L par.K]); Hhat = reshape(Hhat_tot(:,t,:),[par.N*par.L par.K]); y = H*s + noise(:,t); for d=1:length(par.detector) switch (par.detector{d}) % select algorithms case 'CMMSE' bithat = CMMSE_EstH(noise(:,t),par,H,Hhat,N0,s,C_tot); res.BER(d,snr) = res.BER(d,snr) + sum(sum(bits(:,:,n,t)~=bithat))/(par.K*par.Q); case 'CMMSE-SIC' bithat = CMMSE_SIC_EstH(noise(:,t),par,H,Hhat,N0,s,C_tot); res.BER(d,snr) = res.BER(d,snr) + sum(sum(bits(:,:,n,t)~=bithat))/(par.K*par.Q); case 'CMF' bithat = CMF_EstH(noise(:,t),par,H,Hhat,N0,s,C_tot); res.BER(d,snr) = res.BER(d,snr) + sum(sum(bits(:,:,n,t)~=bithat))/(par.K*par.Q); end end end time=toc; time_elapsed = time_elapsed + time; fprintf('estimated remaining simulation time: %3.0f min.\n',time*((par.nbrOfSetups-n+1)*(length(par.SNRdB_list)-snr+1)/60)); end end % normalize results res.BER = res.BER/(par.nbrOfRealizations*par.nbrOfSetups); res.time_elapsed = time_elapsed; % -- show results (generates fairly nice Matlab plot) marker_style = {'bo-','rs-','mv-','kp-','bo-.','rs-.','mv-.','kp-.','bo:','rs:','mv:','kp:'}; figure(2) for d=1:(length(par.detector)) if d==1 semilogy(par.SNRdB_list,res.BER(d,:),marker_style{d},'LineWidth',1) hold on else semilogy(par.SNRdB_list,res.BER(d,:),marker_style{d},'LineWidth',1) end end grid on hold off; xlabel('Average SNR per receive antenna [dB]'); ylabel('Bit error rate (BER)'); legend([par.detector],'location','southwest');