pslandislandmsslandlsnrandlpg.m

来自「相位编码脉冲压缩雷达的多普勒补偿算法」· M 代码 · 共 24 行

%输入为(所有的这个函数都要改成本函数)（只适用于单个系统，对于求信噪比，单个系统和多个系统的求发有所不同）
%Signal:滤波器的输入信号         %Filter_Coef:滤波器系数
%计算
%PSL:峰值旁瓣电平                %ISL：积分旁瓣电平
%MSSL：均方旁瓣电平              %LSNR：信噪比损失
%LPG“脉压增益衰减（loss in procession gain）
function  [PSL,ISL,MSSL,LSNR,LPG]=PSLandISLandMSSLandLSNRandLPG(Signal,Filter_Coef)

%对于给定一个输出序列
%求出主峰值（MainlobePeak）,旁瓣值（Sidelobe），峰值旁瓣值（PeakSidelobe）
%计算出PSl和ISl，MSSL
OutPut=conv(Signal,Filter_Coef);
listCopy=OutPut;                                                                          %备份数据
[Maxlist,MaxLocation]=max(abs(OutPut));
MainlobePeak=Maxlist;
OutPut(MaxLocation)=0;
Sidelobe=OutPut;
PeakSidelobe=max(abs(Sidelobe));                                                           %峰值旁瓣
ISL=10*log10(sum(abs(Sidelobe).^2)/MainlobePeak.^2);                                       %积分旁瓣电平ISL
PSL=10*log10(abs(PeakSidelobe).^2/abs(MainlobePeak).^2);                                   %峰值旁瓣电平PSL
MSSL=10*log10(sum(abs(Sidelobe).^2)/(length(OutPut)*abs(MainlobePeak).^2));                %均方旁瓣电平MSSL
LSNR=10*log10(max(abs(listCopy).^2)/(sum(abs(Signal).^2)*sum(abs(Filter_Coef).^2)));       %信噪比损失LSNR
LPG=20*log10(MainlobePeak/length(Signal));                                                 %脉压增益衰减LPG
OutPut=listCopy;