capacity_11.m

计算信道容量的matlab程序

%clear all
%clc
hold on
Nr=1;
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);
W=w1';
A=C'*W';

a11_real=alph11(10);
a11_imag=alph11(20);


  h11=sqrt(1/2)*(a11_real+j*a11_imag); 
    H=[h11]
    snr_dB=0:1:20;
    snr=10.^(snr_dB/10);
    for i=1:length(snr)
        C4(i)=log2(det(eye(Nr)+snr(i)*H*H'));
    end

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