cp0602_thr_db_vectors.m

超宽带链路分析主程序书标有example的M文件 能很直观的看到超宽带不同距离的带宽（速率）

function [f_th_dB,PSD_th_dB]=cp0602_thr_dB_vectors(f,PSD,threshold)
N=length(f);
df=f(length(f))/length(f);
%计算带宽
[Ppeak,index]=max(PSD); %PSD的峰值
f_peak=index*df;       %峰值频率
Pth=Ppeak*10^(threshold/10); %由给出的门限得出其对应的PSD
%找出最高频率和最低频率
imax=index;
P0h=PSD(index);
while (P0h>Pth)&(imax<=(N/2))
    imax=imax+1;
    P0h=PSD(imax);
end 
f_high=(imax-1)*df; %高频
imin=index;
P0l=PSD(index);
while (P0l>Pth)&(imin>1)&(index>1)
    imin=imin-1;
    P0l=PSD(imin);
end 
f_low =(min(index,imin)-1)*df; %低频

BW=f_high-f_low;  %信号的带宽
f_th_dB=f(imin:imax); %截断后的频率
PSD_th_dB=PSD(imin:imax); %截断后的功率谱