Энергетический расчет канала дальней тропосферной радиосвязи (Программа для ЭВМ) | страница 12
x=sqrt(h01);
Kr1=kramp(x);
Pk1=(1-Kr1)/2;
x=sqrt(h02);
Kr2=kramp(x);
Pk2=(1-Kr2)/2;
disp(['Poh1=', num2str(Pk1),' когерентный прием']);
disp(['Poh2=', num2str(Pk2),' когерентный прием']);
Pn1=exp(-h01/2)/2;
Pn2=exp(-h02/2)/2;
disp(['Poh1=', num2str(Pn1),' некогерентный прием']);
disp(['Poh2=', num2str(Pn2),' некогерентный прием']);
elseif Tm==3
disp('Фазовая манипуляция');
x=sqrt(2*h01);
Kr1=kramp(x);
Pk1=(1-Kr1)/2;
x=sqrt(2*h02);
Kr2=kramp(x);
Pk2=(1-Kr2)/2;
disp(['Poh1=', num2str(Pk1),' когерентный прием']);
disp(['Poh2=', num2str(Pk2),' когерентный прием']);
elseif Tm==4
disp('Относительная фазовая манипуляция');
x=sqrt(2*h01);
Kr1=kramp(x);
Pk1=(1-Kr1)^2/2;
x=sqrt(2*h02);
Kr2=kramp(x);
Pk2=(1-Kr2)^2/2;
disp(['Poh1=', num2str(Pk1),' сравнение полярностей']);
disp(['Poh2=', num2str(Pk2),' сравнение полярностей']);
Pn1=exp(-h01)/2;
Pn2=exp(-h02)/2;
disp(['Poh1=', num2str(Pn1),' сравнение фаз']);
disp(['Poh2=', num2str(Pn2),' сравнение фаз']);
end
end
disp('Цифровая обработка сигнала');
h01c=h01*log2(m);
h02c=h02*log2(m);
if Tm==1
disp('Амплитудная модуляция');
disp('Однополярные прямоугольные импульсы');
x=sqrt(h01c/(m-1)^2);
Kr1=kramp(x);
Pk1=(m-1)/m*(1-Kr1);
x=sqrt(h02c/(m-1)^2);
Kr2=kramp(x);
Pk2=(m-1)/m*(1-Kr2);
set(handles.edit33, 'String', Pk1);
set(handles.edit37, 'String', Pk2);
disp(['Poh1=', num2str(Pk1),' когерентный прием']);
disp(['Poh2=', num2str(Pk2),' когерентный прием']);
x=sqrt(h01c/(2*(m-1)^2));
Kr1=kramp(x);
Pk1=(m-1)/m*(1-Kr1);
x=sqrt(h02c/(2*(m-1)^2));
Kr2=kramp(x);
Pk2=(m-1)/m*(1-Kr2);
set(handles.edit32, 'String', Pk1);
set(handles.edit36, 'String', Pk2);
disp(['Poh1=', num2str(Pk1),' некогерентный прием']);
disp(['Poh2=', num2str(Pk2),' некогерентный прием']);
disp('Биполярные прямоугольные импульсы');
x=sqrt(2*h01c/(m-1)^2);
Kr1=kramp(x);
Pk1=(m-1)/m*(1-Kr1);
x=sqrt(2*h02c/(m-1)^2);
Kr2=kramp(x);
Pk2=(m-1)/m*(1-Kr2);
set(handles.edit35, 'String', Pk1);
set(handles.edit39, 'String', Pk2);
disp(['Poh1=', num2str(Pk1),' когерентный прием']);
disp(['Poh2=', num2str(Pk2),' когерентный прием']);
x=sqrt(h01c/(m-1)^2);
Kr1=kramp(x);
Pk1=(m-1)/m*(1-Kr1);
x=sqrt(h02c/(m-1)^2);
Kr2=kramp(x);
Pk2=(m-1)/m*(1-Kr2);
set(handles.edit34, 'String', Pk1);
set(handles.edit38, 'String', Pk2);
disp(['Poh1=', num2str(Pk1),' некогерентный прием']);
disp(['Poh2=', num2str(Pk2),' некогерентный прием']);
elseif Tm==2
disp('Частотная модуляция');
disp('FSK с разрывом фазы c ФНЧ');
x=sqrt(h01c/2);
Kr1=kramp(x);
Pk1=(m-1)/m*(1-Kr1);
x=sqrt(h02c/2);
Kr2=kramp(x);