Энергетический расчет канала дальней тропосферной радиосвязи (Программа для ЭВМ) | страница 19
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);
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);
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);
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);
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);
Pk2=(m-1)/m*(1-Kr2);
set(handles.edit33, 'String', Pk1);
set(handles.edit37, 'String', Pk2);
disp(['Poh1=', num2str(Pk1),' FSK']);
disp(['Poh2=', num2str(Pk2),' FSK']);
disp('MSK без разрыва фазы c ФНЧ');
x=sqrt(h01c/1.5);
Kr1=kramp(x);
Pk1=(m-1)/m*(1-Kr1);
x=sqrt(h02c/1.5);
Kr2=kramp(x);
Pk2=(m-1)/m*(1-Kr2);
set(handles.edit32, 'String', Pk1);
set(handles.edit36, 'String', Pk2);
disp(['Poh1=', num2str(Pk1),' MSK']);
disp(['Poh2=', num2str(Pk2),' MSK']);
disp('MSK с минимальным отклонением фазы, без разрыва фазы c ФНЧ');
x=sqrt(h01c);
Kr1=kramp(x);
Pk1=(m-1)/m*(1-Kr1);
x=sqrt(h02c);
Kr2=kramp(x);
Pk2=(m-1)/m*(1-Kr2);
set(handles.edit35, 'String', Pk1);
set(handles.edit52, 'String', Pk2);
disp(['Poh1=', num2str(Pk1),' MSK min']);
disp(['Poh2=', num2str(Pk2),' MSK min']);
elseif Tm==3
disp('Фазовая модуляция');
x=sqrt(2*h01c*(sin(pi/m))^2);