function?[irwpslrislr]=QP（tfZfSout）
%?tf:輸入信號的長度
%?Zf:輸入的幅度歸一化的信號
%?Sout:輸入的真實信號
%?輸出：irw:沖擊響應寬度；pslr：峰值旁瓣比；islr：積分旁瓣比；

N=length（tf）;
Midle=round（（N+1）/2）;
%IRW
for?s=Midle+1:1:N-1
????if?（?Zf（s）+3<=Zf（Midle）?&&?Zf（s-1）+3>=Zf（Midle））
????????irw=tf（s）;
????????break;
????end
end
irw=irw*2;
%PSLR
pslr=-inf;
for?k=Midle+1:1:N-1
????if?（?Zf（k+1）<=Zf（k）?&&?Zf（k-1）<=Zf（