Энергетический расчет канала дальней тропосферной радиосвязи (Программа для ЭВМ) | страница 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);