说明：  单天线相比 多天线的信道容量仿真比较 SNR = 10 ^ (SNR/10); n_t = length(R_t); % number of Tx antennas n_r = length(R_r); % number of Rx antennas Pout = 0.1; % outage probability nsample = 1000; % number of channel realizations % nsample = 100 / Pout; R_r_sqrt = R_r^0.5; R_t_sqrt = R_t^0.5; cwf = zeros(nsample, 1); cep = zeros(nsample, 1); Hw = (randn(n_r, n_t, nsample) + sqrt(-1) * randn(n_r, n_t, nsample)) / sqrt(2); for i = 1:nsample H = R_r_sqrt * Hw(:,:,i) * R_t_sqrt; lambda = eig(H * H'); [mu pl] = WFL(lambda, SNR); cwf(i) = sum(log2(1 + pl .* lambda)); cep(i) = sum(log2(1 + SNR/n_t * lambda)); end Ce = mean([cwf cep]); C = sort([cwf cep]); Co = C(round(Pout * nsample),:); Cwf.ergodic = Ce(1); Cep.ergodic = Ce(2); Cwf.outage = Co(1); Cep.outage = Co(2);
(Simulation and comparison of channel capacity between single antenna and multi antenna SNR = 10 ^ (SNR/10); n_t = length(R_t); % number of Tx antennas n_r = length(R_r); % number of Rx antennas Pout = 0.1; % outage probability nsample = 1000; % number of channel realizations % nsample = 100 / Pout; R_r_sqrt = R_r^0.5; R_t_sqrt = R_t^0.5; cwf = zeros(nsample, 1); cep = zeros(nsample, 1); Hw = (randn(n_r, n_t, nsample) + sqrt(-1) * randn(n_r, n_t, nsample)) / sqrt(2); for i = 1:nsample H = R_r_sqrt * Hw(:,:,i) * R_t_sqrt; lambda = eig(H * H'); [mu pl] = WFL(lambda, SNR); cwf(i) = sum(log2(1 + pl .* lambda)); cep(i) = sum(log2(1 + SNR/n_t * lambda)); end Ce = mean([cwf cep]); C = sort([cwf cep]); Co = C(round(Pout * nsample),:); Cwf.ergodic = Ce(1); Cep.ergodic = Ce(2); Cwf.outage = Co(1); Cep.outage = Co(2);)