雷达公式实现

代码来源：《雷达系统分析与设计(MATLAB版)(第三版)》

function [snr] = radar_eq(pt,freq,g,sigma,b,nf,loss,range)
% This program implements Eq.(1.63)
%% Inputs:% pt——峰值功率，W% freq——雷达中心频率，Hz% g——天线增益，dB% sigma——目标截面积，m^2% b——带宽，Hz% nf——噪声系数，dB% loss——雷达损耗，dB% range——目标距离（可以是一个单值或者一个向量），m%% Outputs:% snr——SNR（单值或者向量，取决于目标距离），dB%% 实现
c = 3.0e+8; % speed of light
lambda = c / freq; % wavelength
p_peak = 10*log10(pt); % convert peak power to dB
lambda_sqdb = 10*log10(lambda^2); % compute wavelength square in dB
sigmadb = 10*log10(sigma); % convert sigma to dB
four_pi_cub = 10*log10((4.0*pi)^3); % (4pi)^3 in dB
k_db = 10*log10(1.38e-23); % Boltzmann's constant in dB
to_db = 10*log10(290); % noise temp. in dB
b_db = 10*log10(b); % bandwidth in dB
range_pwr4_db = 10*log10(range.^4); % vector of target range^4 in dB
% implement Equation (1.63)
num = p_peak + g + g + lambda_sqdb + sigmadb; % 分子
den = four_pi_cub + k_db + to_db + b_db + nf + loss + range_pwr4_db; % 分母
snr = num - den;
return 
end