ebno.zip_Eb/No_ebno _zip

% % % Eb/No Calulation program % % % 2005 4.14 Revision 0 % % no_bit = input (' 1. Modulation mode 1= bpsk, 2 qpsk, 3 8psk?'); fec_mode = input(' 2. FEC mode 1= viterbi + RS, 2 = turbo ?'); code_rate = input(' 3. FEC Code Rate 1=1/2, 2=1/3, 3=2/3, 4=3/4, 5=5/6, 6= 7/8 ?'); es_no = input (' 4. Input measured C+N /N with spectrun analyzer(dB) ? '); atm = input ('5 Dataformat 1= ATM, 2= MPEG-II TS ?' ); c_n = 10*log10((10^(es_no/10))-1); if no_bit == 1 es_qpsk = 10*log10(1/1); elseif no_bit == 2 es_qpsk = 10*log10(2/1); elseif no_bit == 3 es_qpsk = 10*log10(3/1); end; et_no = c_n - es_qpsk; if code_rate == 1 fec = 10*log10(2/1); elseif code_rate == 2 fec = 10*log10(3/1); elseif code_rate == 3 fec = 10*log10(3/2); elseif code_rate == 4 fec = 10*log10(4/3); elseif code_rate == 5 fec = 10*log10(6/5); elseif code_rate == 6 fec = 10*log10(8/7); end; et_no = et_no + fec; if fec_mode == 1 if atm == 1 rs_fec = 10*log10(69/53); elseif atm == 2 rs_fec = 10*log10(204/188); end else rs_fec = 0; end; eb_no = et_no + rs_fec; disp ('Calculated Eb/No = '); disp (eb_no); disp ('dB');