np=7;
pt=165.8e3;
freq=3e+9;
g=34.1539;
sigmam=0.5;
sigmaa=4;
te=290.0;
b=1.0e+6;
nf=6.0;
loss=8.0;
SNR_1=(10^1.3)/(2*7)+sqrt((((10^1.3)^2)/(4*7*7))+((10^1.3)/7));
LNCI=10*log10((1+SNR_1)/SNR_1);
loss_total=loss+LNCI;
range=linspace(15e3,100e3,1000);
snrmnci=radar_eq(np*pt,freq,g,sigmam,te,b,nf,loss_total,range);
snrm=radar_eq(pt,freq,g,sigmam,te,b,nf,loss_total,range);
snranci=radar_eq(np*pt,freq,g,sigmaa,te,b,nf,loss_total,range);
snra=radar_eq(pt,freq,g,sigmaa,te,b,nf,loss,range);
rangekm=range./1000;
figure(1)
subplot(2,1,1)
plot(rangekm,snrmnci,'k',rangekm,snrm,'k-.')
grid
legend('With 7-pulse NCI','Single pluse')
ylabel('SNR - dB');
title('Missile case')
subplot(2,1,2)
plot(rangekm,snranci,'k',rangekm,snra,'k-.')
grid
legend('With 7-pulseNCI','Sigle pulse')
ylabel('SNR - dB');
title('Aircraft case')
xlabel('Detection range - Km')
评论7