Приложение б текст программы расчета помехоустойчивости приема сигналов с фм

Расчет вероятности ошибочного приема сигналов с ФМ и ОФМ

clc, clear all; close all;

%Размер шрифта

FS=12;

%Толщина линий

LW=1;

%Размеры графика

Lt=4; Bm=3; Wh=10; Ht=12.5;

%Параметры графика

ns=51;

snrmi=0;snrma=12;

ds=(snrma-snrmi)/(ns-1);

%Расчет вероятности ошибки

snr_db=snrmi:ds:snrma;

snr=10.^(0.1*snr_db);

p1=qfunc(sqrt(2*snr));

p2=2*p1.*(1-p1);

p3=.5*exp(-snr);

%Построение графиков (логарифмическая шкала)

figure;

ymin=1.e-06; ymax=1;

hold on;

semilogy(snr_db,p1,'k-');

semilogy(snr_db,p2,'k--');

semilogy(snr_db,p3,'k-.');

hold off;

grid on;

xlim([snrmi snrma])

ylim([ymin ymax])

line([snrma snrma],[ymin ymax])

line([snrmi snrma],[ymax ymax])

%Параметры графиков

set(gcf,'Color' , [1,1,1])

set(gcf,'Units','centimeters')

set(gcf,'Position',[Lt Bm Wh Ht])

%Position = [left bottom width height]

%Параметры осей

set(gca,'YScale','log')

set(gca,'Units','centimeters')

set(gca,'FontName','Times New Roman','FontSize',FS)

set(findobj('Type','line'),'Color',[0 0 0],'LineWidth',LW)

set(findobj('Type','text'),'FontName','Times New Roman','FontSize',FS)

legend('ФМ','ОФМ (когер.)','ОФМ (некогер.)')

ylabel('Вероятность ошибки, {\itP_{b}}')

xlabel('Отношение сигнал/шум, {\itP_{b}} /{\itN}_{0}, дБ')

line([20 20],[ymin ymax],'Color',[0,0,0])

line([0 20],[ymax ymax],'Color',[0,0,0])

Приложение в текст программы расчета спектральной плотности мощности сигналов с фм

%Расчет СПМ сигналов с ФМ

%Нормировка по амплитуде

clc; clear all; close all;

%Размер шрифта

FS=12;

%Толщина линий

LW=1;

%Размеры графика

Lt=4; Bm=3; Wh=12; Ht=8.5;

%Параметры сигналов

E0=1;

Tb=1;

%Параметры графика

lf_max=5.;

nf=501;

fmi=0;fma=lf_max/Tb;

df=(fma-fmi)/(nf-1);

%Расчет СПМ

f=fmi:df:fma;

y=E0*Tb*sinc(f*Tb).^2;

%Нормировка

y= y/max(y);

%Построение графика (абсолютная шкала)

plot(f,y)

%Параметры графика

set(gcf,'Color' , [1,1,1])

set(gcf,'Units','centimeters')

set(gcf,'Position',[Lt Bm Wh Ht])

%Position = [left bottom width height]

%Параметры осей

set(gca,'Units','centimeters')

set(gca,'FontName','Times New Roman','FontSize',FS)

set(findobj('Type','line'),'Color',[0 0 0],'LineWidth',LW)

set(findobj('Type','text'),'FontName','Times New Roman','FontSize',FS)

ylabel('Нормированная СПМ, {\it\Psi_{s}}({\itf})')

xlabel('Нормированная частота, {\itfT_{b}}')

%Построение графика (логарифмическая шкала)

figure;

umi=-50; uma=0;

plot(f,10*log10(y));

axis([fmi fma umi uma])

%Параметры графика

set(gcf,'Color' , [1,1,1])

set(gcf,'Units','centimeters')

set(gcf,'Position',[Lt Bm Wh Ht])

%Position = [left bottom width height]

%Параметры осей

set(gca,'Units','centimeters')

set(gca,'FontName','Times New Roman','FontSize',FS)

set(findobj('Type','line'),'Color',[0 0 0],'LineWidth',LW)

set(findobj('Type','text'),'FontName','Times New Roman','FontSize',FS)

ylabel('Нормированная СПМ, {\it\Psi_{s}}({\itf}), дБ')

xlabel('Нормированная частота, {\itfT_{b}}')