df_psd.m

来自「source Matlab traites the UWB」· M 代码 · 共 26 行

%
%
%
function y = df_psd(ALPHA,T,Nt,B)

N = ALPHA*Nt;

% $$$ fprintf('ALPHA = %d\n',ALPHA);
% $$$ fprintf('Nt    = %d\n',Nt);
% $$$ fprintf('T     = %.2e\n',T);
% $$$ fprintf('B     = %.2e\n',B);
% $$$ fprintf('2*N = %d, 2*ALPHA*B*T = %d\n',2*N,2*ALPHA*B*T);

y = zeros(1,2*N);

% $$$ y((0:2*ALPHA*B*T-1)+1) = Nt/(2*B*T);
% $$$ y((2*ALPHA*B*T)+1) = Nt/(4*B*T);
% $$$ 
% $$$ y((2*N-2*ALPHA*B*T)+1) = Nt/(4*B*T);
% $$$ y((2*N-2*ALPHA*B*T+1:2*N-1)+1) = Nt/(2*B*T);

y(round([(0:2*ALPHA*B*T-1)+1, (2*N-2*ALPHA*B*T+1:2*N-1)+1])) = Nt/(2*B*T);
y(round([(2*ALPHA*B*T)+1, (2*N-2*ALPHA*B*T)+1])) = Nt/(4*B*T);

return;