mimo.rar_in_mimo_one ring_one-ring_one-ring model

% plot_CCDF.m % Plot the CCDF curves of Fig. 7.3. %MIMO-OFDM Wireless Communications with MATLAB�� Yong Soo Cho, Jaekwon Kim, Won Young Yang and Chung G. Kang %2010 John Wiley & Sons (Asia) Pte Ltd clear all; clc; clf Ns = 2.^[6:10]; b=2; M=2^b; Nblk = 1e3; %mod_object = modem.qammod('M',M, 'SymbolOrder','gray'); %Es=1; A=sqrt(3/2/(M-1)*Es); zdBs = [4:0.1:10]; N_zdBs = length(zdBs); %Ray_fnc = inline('z/s2*exp(-z^2/(2*s2))','s2','z'); CCDF_formula=inline('1-((1-exp(-z.^2/(2*s2))).^N)','N','s2','z'); % Eq.(7.9) for n = 1:length(Ns) N=Ns(n); x = zeros(Nblk,N); sqN=sqrt(N); for k = 1:Nblk %msgint=randint(1,N,M); X=A*modulate(mod_object,msgint); X = mapper(b,N); x(k,:) = ifft(X,N)*sqN; CFx(k) = PAPR(x(k,:)); end s2 = mean(mean(abs(x)))^2/(pi/2); CCDF_theoretical=CCDF_formula(N,s2,10.^(zdBs/20)); for i = 1:N_zdBs %zdB=zdBs(i); %z=10^(zdB/20); %CCDF_theoretical(i)=CCDF_formula(N,s2,z); CCDF_simulated(i) = sum(CFx>zdBs(i))/Nblk; end semilogy(zdBs,CCDF_theoretical,'k-'); hold on; grid on; semilogy(zdBs(1:3:end),CCDF_simulated(1:3:end),'k:*'); end axis([zdBs([1 end]) 1e-2 1]); title('OFDM system with N-point FFT'); xlabel('PAPR0[dB]'); ylabel('CCDF=Probability(PAPR>PAPR0)'); legend('Theoretical','Simulated');