log_psd.m

来自「本程序仿真一个从卫星到区域数据中心的下行链路（由一组6个载波构成）」· M 代码 · 共 20 行

function[logpsd,freq,ptotal,pmax]=log_psd(x,n,ts)
y=zeros(1,0);
h=waitbar(0,'For Loop in PSD Calcuation');
for k=1:n
freq(k)=(k-1-(n/2))/(n*ts);
y(k)=x(k)*((-1.0)^k);
waitbar(k/n)
end;
v=fft(y)/n;
psd=abs(v).^2;
pmax=max(psd);
ptotal=sum(psd);
logpsd=10*log10(psd/pmax);
for k=1:n
if(logpsd(k)<-60.0)
logpsd(k)=-60.0;
end
end
close(h)