Приложение 2

Программа исследования влияния отношение сигнал/шум на вероятность обнаружения первичного сигнала.

inp_t=input('do u want to add noise: ','s');

if(inp_t=='y'|inp_t=='Y')

d = input('Enter the SNR in dB: ');

figure

Y = awgn(y,d);

Pxx1 = periodogram(Y);

Hpsd = dspdata.psd(Pxx1,'Fs',Fs);

plot(Hpsd);

disp('adding noise');

chek1 = Pxx1(25).*10000;

chek2 = Pxx1(49).*10000;

chek3 = Pxx1(62).*10000;

chek4 = Pxx1(89).*10000;

chek5 = Pxx1(105).*10000;

if(chek1 < 8000)

disp('User 1 is not present.');

else

disp('User 1 is present.');

end

if(chek2 < 8000)

disp('User 2 is not present.');

else

disp('User 2 is present.');

end

if(chek3 < 8000)

disp('User 3 is not present.');

else

disp('User 3 is present.');

end

if(chek4 < 8000)

disp('User 4 is not present.');

else

disp('User 4 is present.');

end

if(chek5 < 8000)

disp('User 5 is not present.');

else

disp('User 5 is present.');

end

end

Приложение 3

Программа исследования влияния замирание на вероятность обнаружения первичного сигнала.

temp = input('Do you want to attenuate the signals? [Y/N]: ','s');

if(temp == 'Y' | temp == 'y')

Yff=0;

aF = input('Enter the percentage to attenuate the signal: ');

figure

tem = aF/100;

tm = 1-tem;

Z = y.*tm;

disp('attenuating');

grid on

plot(Z);

figure

Pxx4 = periodogram(Z);

Hpsd = dspdata.psd(Pxx4,'Fs',Fs);

plot(Hpsd);

chek1 = Pxx4(25).*10000;

chek2 = Pxx4(49).*10000;

chek3 = Pxx4(62).*10000;

chek4 = Pxx4(89).*10000;

chek5 = Pxx4(105).*10000;

if(chek1 < 8000)

disp('User 1 is not present.');

else

disp('User 1 is present.');

end

if(chek2 < 8000)

disp('User 2 is not present.');

else

disp('User 2 is present.');

end

if(chek3 < 8000)

disp('User 3 is not present.');

else

disp('User 3 is present.');

end

if(chek4 < 8000)

disp('User 4 is not present.');

else

disp('User 4 is present.');

end

if(chek5 < 8000)

disp('User 5 is not present.');

else

disp('User 5 is present.');

end

end

end

%if rerun the program

temp = input('Do you want to re-run the program? [Y/N]: ','s');

if(temp == 'Y' | temp == 'y')

disp('\n\nEnter the users again.\n\n');

else

break;

end

end