doppler_psd.m

多径衰落模型的基本建模

clear all

N = 100;
fm = 50;  %unit: 1Hz

deltaf = 2*fm/(N-1);
T = 1/deltaf;
sf0 = 1.5/(pi*fm);
for n = 1:(N-2)/2
    sf(n) = 1.5/(pi*fm*sqrt(1-(n*deltaf/fm)^2));
end

classicf = [fliplr(sf),sf0,sf];

% figure;
t=-49:1:49; 
plot(t,classicf);
xlabel('f/ Hz');       
ylabel('s(f)');