OFDM.zip_OFDM ，Matlab_OFDM_MATLAB_matlab_matlab ofdm_ofdm matlab

%------------------------------------------------------- % TP de Communications Numériques : transmission OFDM %------------------------------------------------------- clear all; close all; N = 2^15; M = 4; Ns = N/M; Np = 256; L = 64; Ntirage = 1; s = zeros(1, Ns); sest = zeros(1, Ns); aEst = zeros(1,N); x = zeros(1,Ns + L*Ns/Np); r = zeros(1,Ns); Teb = zeros(1,10); % Filtre modélisant le canal %--------------------------- h = [1 -0.7]; H = fft(h,Np) /sqrt(Np); % DFT de h pour Zero Forcing % Chaine de communication %------------------------ for EbN0db = 0:9 %data = fopen('data.txt','a'); for k = 1:Ntirage % Génération d'une trame binaire uniformément distribuée a = round(rand(1,N)); % Mise sur constellation : modulation 16-QAM for j = 1:Ns word = a(1+(j-1)*4:j*4); s(j) = MQAM(word); end sigmas2 = var(s); Es = sigmas2; sigmab2 = Es/log2(2^M) * 10^(-EbN0db /10); %----------------------------------------------- % Emetteur OFDM j = 1:Np; index = 1:Np+L; while j(end) <= Ns ofdmSymb = s(j); % S/P xk = sqrt(Np)*ifft(ofdmSymb,Np); % IDFT prefix = xk(end-L+1:end); % préfixe cyclique x(index) = [prefix xk]; % P/S j = j+Np; index = index +Np+L; end %---------------------------------------------- % Passage dans le canal xf = canal(x, h); % Ajout de Bruit blanc gaussien complexe b = bcgaussien (sigmab2, length(xf)); y = xf + b; %---------------------------------------------- % Récepteur OFDM j = 1:Np; p = 1:Np+L; while j(end) <= Ns yk = y(p); % S/P yk = yk(L+1:end); % suppression du préfixe rk = fft(yk,Np) /sqrt(Np); % DFT rkt = rk ./ H; % ZF r(j) = rkt; % P/S j = j+Np; p = p+Np+L; end %----------------------------------------------- % Décision for j=1:Ns sest(j) = seuil(r(j)); end % Demapping for j = 1:Ns word = unmapping(sest(j)); aEst(1+(j-1)*4:j*4) = word; end % Calcul des erreurs bit diff = mod(aEst - a, 2); Teb(EbN0db +1) = Teb(EbN0db +1) + sum(diff); end % Teb moyen Teb(EbN0db +1) = Teb(EbN0db +1) / (Ntirage * N); EbN0db Teb(EbN0db +1) % écriture des Teb dans un fichier %fprintf(data,'%i %f\n',EbN0db,Teb(EbN0db +1)); %fclose(data); end semilogy([0:9],Teb);