【正文】 for nSample = 1:SampleNum H = (randn(Nr,Nt)+sqrt(1)*randn(Nr,Nt))/sqrt(2)。 Cn(nSNR,nSample) = log2(det(1+rho*H*H39。/Nt))。 end Capacity2(nSNR) = mean( Cn(nSNR,:) )。 tempC2 = sort(Cn(nSNR,:))。 end Nt = 5。 % 5 Nr = 1。 SampleNum = 10000。 SNR_dB = 0:2:24。 for nSNR = 1:length(SNR_dB) SNR_dB(nSNR) rho = 10^(SNR_dB(nSNR)/10)。 for nSample = 1:SampleNum H = (randn(Nr,Nt)+sqrt(1)*randn(Nr,Nt))/sqrt(2)。 Cn(nSNR,nSample) = log2(det(1+rho*H*H39。/Nt))。 end Capacity3(nSNR) = mean( Cn(nSNR,:) )。 tempC3 = sort(Cn(nSNR,:))。 end Nt = 7。 % 7 Nr = 1。 SampleNum = 10000。 SNR_dB = 0:2:24。 for nSNR = 1:length(SNR_dB) SNR_dB(nSNR) rho = 10^(SNR_dB(nSNR)/10)。 for nSample = 1:SampleNum H = (randn(Nr,Nt)+sqrt(1)*randn(Nr,Nt))/sqrt(2)。 Cn(nSNR,nSample) = log2(det(1+rho*H*H39。/Nt))。 29 end Capacity4(nSNR) = mean( Cn(nSNR,:) )。 tempC4 = sort(Cn(nSNR,:))。 end Nt = 9。 % 9 Nr = 1。 SampleNum = 10000。 SNR_dB = 0:2:24。 for nSNR = 1:length(SNR_dB) SNR_dB(nSNR) rho = 10^(SNR_dB(nSNR)/10)。 for nSample = 1:SampleNum H = (randn(Nr,Nt)+sqrt(1)*randn(Nr,Nt))/sqrt(2)。 Cn(nSNR,nSample) = log2(det(1+rho*H*H39。/Nt))。 end Capacity5(nSNR) = mean( Cn(nSNR,:) )。 tempC5 = sort(Cn(nSNR,:))。 end figure(1) % nSNR = length(SNR_dB)。 plot(tempC1, (1:SampleNum)/SampleNum, 39。r39。)。 hold on。 plot(tempC2, (1:SampleNum)/SampleNum, 39。b39。)。 hold on。 plot(tempC3, (1:SampleNum)/SampleNum, 39。m39。)。 hold on。 plot(tempC4, (1:SampleNum)/SampleNum, 39。c39。)。 hold on。 plot(tempC5, (1:SampleNum)/SampleNum, 39。g39。)。 hold on。 xlabel(39。信道容量 (bit/s/Hz)39。)。 ylabel(39。概率分布 39。)。 text(,39。從左向右發(fā)射 天線數(shù)依次為 1,3,5,7,939。)。 grid on。 30 SIMO信道容量仿真的 matlab實(shí)現(xiàn) % clear all。clc。 Nt = 1。 % 1 Nr = 1。 % 1 SampleNum = 10000。 % SNR_dB =0:2:24。 % (單位 dB) for nSNR = 1:length(SNR_dB) SNR_dB(nSNR) rho = 10^(SNR_dB(nSNR)/10)。 % for nSample = 1:SampleNum H = (randn(Nr,Nt)+sqrt(1)*randn(Nr,Nt))/sqrt(2)。 % Cn(nSNR,nSample) = log2(real(det(1+rho*H39。*H)))。 end Capacity1(nSNR) = mean( Cn(nSNR,:) )。 % tempC1 = sort(Cn(nSNR,:))。 end Nt = 1。 % 3 Nr = 3。 SampleNum = 10000。 SNR_dB =0:2:24。 for nSNR = 1:length(SNR_dB) SNR_dB(nSNR) rho = 10^(SNR_dB(nSNR)/10)。 for nSample = 1:SampleNum H = (randn(Nr,Nt)+sqrt(1)*randn(Nr,Nt))/sqrt(2)。 Cn(nSNR,nSample) = log2(real(det(1+rho*H39。*H)))。 end Capacity2(nSNR) = mean( Cn(nSNR,:) )。 tempC2 = sort(Cn(nSNR,:))。 end Nt = 1。 % 5 Nr = 5。 31 SampleNum = 10000。 SNR_dB =0:2:24。 for nSNR = 1:length(SNR_dB) SNR_dB(nSNR) rho = 10^(SNR_dB(nSNR)/10)。 for nSample = 1:SampleNum H = (randn(Nr,Nt)+sqrt(1)*randn(Nr,Nt))/sqrt(2)。 Cn(nSNR,nSample) = log2(real(det(1+rho*H39。*H)))。 end Capacity3(nSNR) = mean( Cn(nSNR,:) )。 tempC3 = sort(Cn(nSNR,:))。 end Nt = 1。 % 7 Nr = 7。 SampleNum = 10000。 SNR_dB = 0:2:24。 for nSNR = 1:length(SNR_dB) SNR_dB(nSNR) rho = 10^(SNR_dB(nSNR)/10)。 for nSample = 1:SampleNum H = (randn(Nr,Nt)+sqrt(1)*randn(Nr,Nt))/sqrt(2)。 Cn(nSNR,nSample) = log2(real(det(1+rho*H39。*H)))。 end Capacity4(nSNR) = mean( Cn(nSNR,:) )。 tempC4 = sort(Cn(nSNR,:))。 end Nt = 1。 % 9 Nr = 9。 SampleNum = 10000。 SNR_dB = 0:2:24。 for nSNR = 1:length(SNR_dB) SNR_dB(nSNR) rho = 10^(SNR_dB(nSNR)/10)。 32 for nSample = 1:SampleNum H = (randn(Nr,Nt)+sqrt(1)*randn(Nr,Nt))/sqrt(2)。 Cn(nSNR,nSample) = log2(real(det(1+rho*H39。*H)))。 end Capacity5(nSNR) = mean( Cn(nSNR,:) )。 tempC5 = sort(Cn(nSNR,:))。 end figure(1) % nSNR = length(SNR_dB)。 plot(tempC1, (1:SampleNum)/SampleNum,39。r39。)。 hold on。 plot(tempC2, (1:SampleNum)/SampleNum,39。b39。)。 hold on。 plot(tempC3, (1:SampleNum)/SampleNum,39。m39。)。 hold on。 plot(tempC4, (1:SampleNum)/SampleNum,39。c39。)。 hold on。 plot(tempC5, (1:SampleNum)/SampleNum,39。g39。)。 hold on。 xlabel(39。信道容量 (bit/s/Hz)39。)。 ylabel(39。概率分布 39。)。 text(,39。從左向右接收天線數(shù)依次為 1,3,5,7,939。)。 grid on。 MIMO 信道容量仿真的 matlab 實(shí)現(xiàn) % cle