MIMO-OFDM系统容量绪论.ppt

* %画收、发天线个数为2×2、4×4、6×6、8×8时的容量 for n = 1:length(Nt), N = Nt(n); for m = 1:M, H = (randn(N,N)+1i*randn(N,N))/sqrt(2); rho = 10^(SNR/10); CU(m,n) = log2(real(det(eye(N)+rho*H*H/N))); end [cdf,c] = hist(CU(:,n),n_bins); plot(c,1-(cumsum(cdf))/M) end grid on; * 发送端未知CSI时MIMO系统容量CCDF仿真图 * clear all close all clc M = 5000; R = [0 0.95];相关矩阵 SNR = [0:2:20]; figure; xlabel(SNR [dB], FontSize,18); ylabel(容量（bit/s/Hz）, FontSize,18); title(相关信道的容量比较 ,FontSize,18) grid on hold on; 实例4_2 相关对MIMO系统容量影响仿真 * for l = 1:length(R), R_t = eye(2); R_r = [1 R(l);R(l) 1] for snr_idx = 1:length(SNR), snr = 10^(SNR(snr_idx)/10); for m = 1:M, Hw = (randn(2,2)+1i*randn(2,2))/sqrt(2); H = R_r^(.5)*Hw*R_t^(.5); C(m,snr_idx)=log2((det(eye(2)+snr*H*H/2))); end Capacity(snr_idx,l) = mean(C(:,snr_idx)); end end * plot(SNR,Capacity(:,1),k-,SNR,Capacity(:,2),kx-); for l = 1:length(R), R_t = eye(2); R_r = [1 R(l);R(l) 1] for snr_idx = 1:length(SNR), snr = 10^(SNR(snr_idx)/10); for m = 1:M, Hw = (randn(2,2)+1i*randn(2,2))/sqrt(2); H = R_r^(.5)*Hw*R_t^(.5); [gamma,eigs] = pwr_modes(H,snr); C(m,snr_idx) = sum(log2(real(1+eigs.*gamma*snr/2))); end Capacity(snr_idx,l) = mean(C(:,snr_idx)); * end end plot(SNR,Capacity(:,1),k-.,SNR,Capacity(:,2),kd-); legend(R=0 CSI未知,R=0.95 CSI未知, R=0 CSI已知,R=0.95 CSI已知,2); * %计算gamma和奇异值子程序，子程序名称：pwr_modes.m function [g,l] = pwr_modes(H,rho) %这部分程序参看注水算法。 N_tilde = size(H,1); l_tilde = real(eig(H*H)); l = l_tilde(find(l_tilde~=0)); N = length(l); mu = (N_tilde+sum(N_tilde./(rho*l)))/N; g = mu-N_tilde./(rho*l); while (length(find(g = 0)) ~= 0) l = l(find(g 0)); N = length(l); mu = (N_tilde+sum(N_t