capacity_33.m

计算信道容量的matlab程序

%clear all
%clc
hold on
Nr=3;
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(4,:);
alph22=A(5,:);
alph23=A(6,:);
%alph24=A(8,:);
alph31=A(7,:);
alph32=A(8,:);
alph33=A(9,:);
%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;...
        h21,h22,h23;...
        h31,h32,h33];
    snr_dB=0:1:20;
    snr=10.^(snr_dB/10);
    for i=1:length(snr)
        C3(i)=log2(det(eye(Nr)+snr(i)*H*H'));
    end

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