%輸入為(所有的這個(gè)函數(shù)都要改成本函數(shù))（只適用于單個(gè)系統(tǒng)，對(duì)于求信噪比，單個(gè)系統(tǒng)和多個(gè)系統(tǒng)的求發(fā)有所不同）
%Signal:濾波器的輸入信號(hào)         %Filter_Coef:濾波器系數(shù)
%計(jì)算
%PSL:峰值旁瓣電平                %ISL：積分旁瓣電平
%MSSL：均方旁瓣電平              %LSNR：信噪比損失
%LPG“脈壓增益衰減（loss in procession gain）
function  [PSL,ISL,MSSL,LSNR,LPG]=PSLandISLandMSSLandLSNRandLPG(Signal,Filter_Coef)

%對(duì)于給定一個(gè)輸出序列
%求出主峰值（MainlobePeak）,旁瓣值（Sidelobe），峰值旁瓣值（PeakSidelobe）
%計(jì)算出PSl和ISl，MSSL
OutPut=conv(Signal,Filter_Coef);
listCopy=OutPut;                                                                          %備份數(shù)據(jù)
[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;