MIMO-OFDM系统容量汇编.ppt

下载文档 降价啦

55
0
约6.27千字
约 38页
2016-09-24 发布于湖北
举报
版权申诉
保障服务

MIMO-OFDM系统容量汇编.ppt

1、本文档共38页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

MIMO-OFDM系统容量汇编

* %画收、发天线个数为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