rayleigh.m

来自「基频收发机设计以Matlab仿真」· M 代码 · 共 26 行

function [r]=rayleigh(n,na,fc,v)
% clear all
% fc=3000;
% v=10;
% n=1000;
% na=1   

wave_l=1;

a=randn(1,n-2);
b=randn(1,n-2);
fm=v/wave_l;
f=2*fm/(n-1)+(fc-fm):2*fm/(n-1):(fc+fm)-2*fm/(n-1);
SE=real((1.5*((pi*fm*(1-((f-fc)/fm).^2).^(1/2))).^(-1)).^(1/2));

FA=fftshift(abs(fft(a)));
FB=fftshift(abs(fft(b)));

da=(ifft(FA.*SE));
db=(ifft(FB.*SE));

r=na*10*abs(((da(3:n-3)).^2+(db(3:n-3)).^2).^(1/2));
% s_r=hist(r,500);
% 
% subplot(221),plot(f,SE),title('Doppler frequency respond');
% subplot(222),plot(f,FA),title('FA