%
%
% 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');
ebno.zip_Eb/No_ebno _zip
版权申诉
5星 · 超过95%的资源 22 浏览量
2022-07-14
13:43:42
上传
评论
收藏 613B ZIP 举报
我虽横行却不霸道
- 粉丝: 75
- 资源: 1万+
评论1