Iscrtavanje OFDM BER versus EB / Ne

[grozdov]

Pozdrav;
Ja sam novi u Matlab, ja je dobio ovaj Matlab koda iz jedne knjige ...Želim zemljište BER versus EB / Ne ...može netko pomoć mene u ovom?
Hvala

% ****** Ovaj je simultion ostvariti prijenos OFDM

para = 128;%
broj paralelnih kanala
fftlen = 128;% FFT length
noć = 128;% Broj prijevoznici
nd = 6;%
No informacija OFDM simbol za jednu petlju
ml = 2;% Modulacija razini
sr = 250000% Symbol rate
br = sr * ml% Bit rate per carrier
gilen = 32;% Dužina branič interval bodova
ebno = 1:5;% EB / Ne probajte i 100

% ****** Glavna petlja *************************

nloop = 100;% Broj petlje simulacije
noe = 0;% Broj pogreške podataka
kimati = 0;%
No Br. prenose podaci
eop = 0;% error Broj paketa
nop = 0;% NE.prenose u paketima

za i = 1: nloop
% ********** Odašiljač
% Podaci generatora
za k = 1: length (ebno)
Rand seldata = (1,
st. * nd * ml)> 0,5;
% Serijsko-paralelna pretvorba
paradata = preoblikovati (seldata, parasailing, nd * ml);

% QPSK modulaciju
[ich, qch] = qpskmod (paradata, parasailing, nd, ml);

% I i Q su podaci povećao kmod puta za normalizovati podataka
kmod = 1/sqrt (2);
ich1 = ich kmod .*;
qch1 = qch .* kmod;

IFFT%
x = ich1 .* qch1 I;% Parellel se pretvaraju u serijskoj
y = ifft (x);
ich2 = real ;
qch2 = Imag ;

% Guard Interval umetanje
[ich3, qch3] = giins (ich2, qch2, fftlen, gilen, ND);
fftlen2 = fftlen gilen

% Prigušenja
spow = suma (ich3. ^ 2 qch3. ^ 2) / nd. / para;
Attn = 0.5 * spow sr * / * 10 br. ^ (-ebno/10);
Attn = sqrt (Attn);

Primatelj% plus buke
[ich4, qch4] = češalj (ich3, qch3, Attn (k));

% Guard uklanjanje
[ich5, qch5] = girem (ich4, qch4, fftlen2, gilen, ND);

% FFT

RX = ich5 .* qch I;
ry = fft (RX);
ich6 = real (ry);
qch6 = Imag (ry);

% Demodulator
ich7 = ich6./kmod;
qch7 = qch6./kmod;
[demodata] = qpskdemod (ich7, qch7, parasailing, nd, ml);

% Paralelno-serijska pretvorba
demodata1 = preoblikovati (demodata, 1,
st. * nd * ml);

% BER izračuni
noe2 (k) = suma (ABS (demodata1-seldata));
nod2 (k) = duljina (seldata);

% Cummulative na br.pogrešaka i podataka u noe i kimati
noe = noe (k) noe2;
klimanje glavom kimati = (k) nod2;Izračun% PER
ako noe2 ~ = 0;
eop = eop
kraj
nop = nop 1;
BER = noe / drijemanje
po = eop / nop;
BER (k) = noe / drijemanje
kraj

kraj

Izlazni rezultati ***********%

 

[cipher_crypto]

možeš objasniti gdje točno trebate pomoć?Mislim da sam vidio ovaj kod.morate dobiti BER rezultat.parcela, a zatim protiv vrijednosti odgovaraju oni.