capacity_44.m

计算信道容量的matlab程序

clear all
clc
Nr=4;
Nt=Nr;
t=10;
[corrRT]=mimo_channel(Nr, Nt,t);
[power]=generator_power(Nr);
T=corrRT.*power;
C=chol(T);
N=500;
w1=randn(1,N);
w2=randn(1,N);
w3=randn(1,N);
w4=randn(1,N);
w5=randn(1,N);  
w6=randn(1,N);
w7=randn(1,N);
w8=randn(1,N);
w9=randn(1,N);
w10=randn(1,N);
w11=randn(1,N);
w12=randn(1,N);
w13=randn(1,N);
w14=randn(1,N);
w15=randn(1,N);
w16=randn(1,N);
W=[w1',w2',w3',w4',w5',w6',w7',w8',w9',w10',w11',w12',w13',w14',w15',w16'];
A=C'*W';
alph11=A(1,:);
alph12=A(2,:);
alph13=A(3,:);
alph14=A(4,:);
alph21=A(5,:);
alph22=A(6,:);
alph23=A(7,:);
alph24=A(8,:);
alph31=A(9,:);
alph32=A(10,:);
alph33=A(11,:);
alph34=A(12,:);
alph41=A(13,:);
alph42=A(14,:);
alph43=A(15,:);
alph44=A(16,:);
%%%
a11_real=alph11(10);
a11_imag=alph11(20);
a12_real=alph12(10);
a12_imag=alph12(20);
a13_real=alph13(10);
a13_imag=alph13(20);
a14_real=alph14(10);
a14_imag=alph14(20);
%%%
a21_real=alph21(10);
a21_imag=alph21(20);
a22_real=alph22(10);
a22_imag=alph22(20);
a23_real=alph23(10);
a23_imag=alph23(20);
a24_real=alph24(10);
a24_imag=alph24(20);
%%%
a31_real=alph31(10);
a31_imag=alph31(20);
a32_real=alph32(10);
a32_imag=alph32(20);
a33_real=alph33(10);
a33_imag=alph33(20);
a34_real=alph34(10);
a34_imag=alph34(20);
%%%
a41_real=alph41(10);
a41_imag=alph41(20);
a42_real=alph42(10);
a42_imag=alph42(20);
a43_real=alph43(10);
a43_imag=alph43(20);
a44_real=alph44(10);
a44_imag=alph44(20);
%%%%
  h11=sqrt(1/2)*(a11_real+j*a11_imag); 
  h12=sqrt(1/2)*(a12_real+j*a12_imag); 
  h13=sqrt(1/2)*(a13_real+j*a13_imag); 
  h14=sqrt(1/2)*(a14_real+j*a14_imag);
  %%%
    h21=sqrt(1/2)*(a21_real+j*a21_imag); 
  h22=sqrt(1/2)*(a22_real+j*a22_imag); 
  h23=sqrt(1/2)*(a23_real+j*a23_imag); 
  h24=sqrt(1/2)*(a24_real+j*a24_imag); 

  %%
  h31=sqrt(1/2)*(a31_real+j*a31_imag); 
  h32=sqrt(1/2)*(a32_real+j*a32_imag); 
  h33=sqrt(1/2)*(a33_real+j*a33_imag); 
  h34=sqrt(1/2)*(a34_real+j*a34_imag); 
     %%%
    h41=sqrt(1/2)*(a41_real+j*a41_imag); 
  h42=sqrt(1/2)*(a42_real+j*a42_imag); 
  h43=sqrt(1/2)*(a43_real+j*a43_imag); 
  h44=sqrt(1/2)*(a44_real+j*a44_imag); 

    H=[h11,h12,h13,h14;...
        h21,h22,h23,h24;...
        h31,h32,h33,h34;...
        h41,h42,h43,h44;];
    snr_dB=0:1:20;
    snr=10.^(snr_dB/10);
    for i=1:length(snr)
        C1(i)=log2(det(eye(Nr)+snr(i)*H*H'));
    end

plot(snr_dB,C1,'*');
xlabel('SNR(db)');
ylabel('信道容量(bits/s/Hz)');