gongnengfangzhen.m

这是一些自己几个月的劳动成果

clear,clc,clf;
Rmin=1;  Rmax=100;  %km
% Define Rtar (km) from Rmin to Rmax
Rtar=Rmin:Rmax;      %1 km steps
% Radar and Target Parameters
Pt=2*10^5;            %w
Gr=10^4;
sigma=3;             %sq m
w1=0.03;               % m
Lr=10;
Pr1=(Pt*Gr^2*sigma*w1^2)./((4*pi)^3*Lr*(Rtar*10^3).^4);
Pr1dbm=10*log10(Pr1);
% Jammer Related Parameters
Pj=200;               %w
Gj=10^(13/10);
Br=2*10^6;               %Hz
Lj=10^(7/10);
Pr2=(Pj*Gj*Gr*w1^2)./((4*pi)^2*(Rtar*10^3).^2*Lj);
Pr2dbm=10*log10(Pr2);
figure(1);
plot(Rtar,Pr1dbm,'r',Rtar,Pr2dbm,'k');
ylabel('功率增益(dB)');
xlabel('探测距离（单位Km）');grid;
title('单脉冲雷达功率随距离变化关系曲线');
text(50,-130,'脉冲信号噪声功率');
text(70,-80,'干扰机噪声功率');
k=1.38*10^(-23);       
T0=290;             %K
F=10^(5/10);
Prs=Pr1./(k*T0*Br*F);
Prsdbm=10*log10(Prs);
Prsj=Pr1./(Pr2+k*T0*Br*F);
Prsjdbm=10*log10(Prsj);
% Plotting Received Target and Jammer Signals
Prsdbm=Prsdbm.*(Rtar./Rtar);
Prsjdbm=Prsjdbm.*(Rtar./Rtar);
figure(2);
subplot(211);plot(Rtar,Prsdbm,'r');
ylabel('S/N(dB)');grid;
title('不受干扰时热噪声信号比');
subplot(212);plot(Rtar,Prsjdbm,'k');
ylabel(' S/(J+N)(dB)');
xlabel('探测距离（单位Km）');grid;
title('噪声干扰下信噪比');