MIMO典型算法MATLAB程序

clc; clear; NUM=1; M=4; %接收天线数 N=8; %发送天线数 Lt=4; %发送端选择天线数 anglespread=5; d=3; angle=30; % anglespread 角度扩展在 0 －－10 变化，越小，相关性越高 % angle 角度在0 －－60 变化 越大， 相关性越高 % d 距离在1 －－5 变化 越小，相关性越高 I =eye(M); Lsnr=0; % low noise power Hsnr=20; % high noise power i=1; SNR=(Lsnr:Hsnr); for snr=SNR ga(i)=10^(snr/10); i=i+1; end Ccsa1=zeros(1,Hsnr-Lsnr+1); CLcsa1=zeros(1,Hsnr-Lsnr+1); Cexh1=zeros(1,Hsnr-Lsnr+1); Crand1=zeros(1,Hsnr-Lsnr+1); Cnorm1=zeros(1,Hsnr-Lsnr+1); Cgcicle1=zeros(1,Hsnr-Lsnr+1); CLow_gcicle1=zeros(1,Hsnr-Lsnr+1); Rr=mimocor(anglespread,angle,d,M); %correlation at Rx Rt=mimocor(anglespread,angle,d,N); [h]=mimoCN(Lt,M); [h1]=mimoCN(N,M); H=mimocorchannel(N,M,h1,Rr,Rt); G=H'*H; for num=1:NUM %%%%%%%%%%%%%%%穷举法 [Hexh]=mimoexh(G,N,Lt); %%%%%%%%范数法 [Hnorm]=mimonorm(G,N,Lt); %%%%%%%%%%%%%%%%%随机选取法 [Hrand]=mimorand(G,N,Lt); %%% Gcircle算法 [Hgcircle]=mimogcicle(G,N,Lt); %%% 改进Gcircle算法 [HLow_gcicle]=mimoLow_gcicle(G,N,Lt); %%%%%%%%%%%%%%%%CSA算法 [Rcsa]=mimoCSA(Rt,N,Lt); %correlation at Tx [Hcsa]=mimocorchannel(Lt,M,h,Rr,Rcsa); %%%%%%%%%%%%%%%%%L-CSA [~,RLcsa]=mimoLowCSA(Rt,N,Lt); [HLcsa]=mimocorchannel(Lt,M,h,Rr,RLcsa); %容量计算 for i=1:Hsnr-Lsnr+1 Ccsa(i)=log2(real(det(I+ga(i)*Hcsa*Hcsa'))); CLcsa(i)=log2(real(det(I+ga(i)*HLcsa*HLcsa'))); Cexh(i)=log2(real(det(I+ga(i)*Hexh*Hexh'))); Crand(i)=log2(real(det(I+ga(i)*Hrand*Hrand'))); Cnorm(i)=log2(real(det(I+ga(i)*Hnorm*Hnorm'))); Cgcicle(i)=log2(real(det(I+ga(i)*Hgcircle*Hgcircle'))); CLow_gcicle(i)=log2(real(det(I+ga(i)*HLow_gcicle*HLow_gcicle'))); end Ccsa1=Ccsa1+Ccsa; CLcsa1= CLcsa1+CLcsa; Cexh1=Cexh1+Cexh; Crand1=Crand1+Crand; Cnorm1=Cnorm1+Cnorm; Cgcicle1=Cgcicle1+Cgcicle; CLow_gcicle1=CLow_gcicle1+CLow_gcicle; end Ccsa=Ccsa1/NUM; CLcsa=CLcsa1/NUM; Cexh=Cexh1/NUM; Crand=Crand1/NUM; Cnorm=Cnorm1/NUM; Cgcicle=Cgcicle1/NUM; CLow_gcicle=CLow_gcicle1/NUM; plot(SNR,Cexh,'md'); hold on; plot(SNR,Cnorm,'r+'); hold on; plot(SNR,Ccsa,'r*'); hold on; plot(SNR,CLcsa,'b.'); hold on; plot(SNR,Crand,'cx'); hold on; plot(SNR,Cgcicle,'ks'); hold on; plot(SNR,CLow_gcicle,'gh'); xlabel('SNR[dB]'); ylabel('Capacity'); legend('Cexh','Cnorm','Ccsa','CLcsa','Crand','Cgcicle','CLow_gcicle'); title('Capacity Plots '); grid on;