jakes模型

clear all, clear all; M=8; N=4*M+2; fd=926;%an=0 的 fd % fd=input(an=0 的 fd:测试值 926：);%an=0 时的最大多普勒频移 wd=2*pi*fd; Ns=50000; % Ns=50000; Ts=1e-6;%样本点数和采样频率 t=0:Ns-1*Ts; %f=-Ns/2:Ns/2-1/(Ns*Ts*fd); t_0=0; %t=t; %ph=unifrnd (-pi,pi,1,M+1);%-pipi 均匀分布随机变量 %set beita_n,an,bn,wn uc=zeros(1,Ns); us=zeros(1,Ns); ruu=zeros(1,Ns); rucuc=zeros(1,Ns); for n=1:M+1 if(n=M+1) beita_n(M+1)=pi/4; an(M+1)=sqrt(2)*cos(beita_n(M+1); bn(M+1)=sqrt(2)*sin(beita_n(M+1); wn(M+1)=wd; else beita_n(n)=pi*n/M; an(n)=2*cos(beita_n(n); bn(n)=2*sin(beita_n(n); wn(n)=wd*cos(2*pi*n/N); end temp1=2/sqrt(N)*(an(n)*cos(wn(n)*t); uc=uc+temp1;%实部 temp2=2/sqrt(N)*(bn(n)*cos(wn(n)*t); us=us+temp2;%虚部 temp3=4/N*(2*cos(wn(n)*t)+cos(wd*t); ruu=ruu+temp3;%自相关 temp4=4/N*(an(n).2/2*cos(wn(n)*t); rucuc=rucuc+temp4; end u=uc+1i*us; figure(1); subplot(311); plot(1:Ns*Ts,10*log10(abs(u) title(Jakes Model, f_d=,num2str(fd),Hz, T_s=,num2str(Ts),s); % axis(0 0.05 -20 10); xlabel(times), ylabel(MagnitudedB) %幅度分布 hold on; subplot(223); hist(abs(u),50); title(Jakes Model, f_d=,num2str(fd),Hz, T_s=,num2str(Ts),s); xlabel(Magnitude), ylabel(Occasions) %相位分布 subplot(224), hist(angle(u),50); title(Jakes Model, f_d=,num2str(fd),Hz, T_s=,num2str(Ts),s); xlabel(Phaserad), ylabel(Occasions) D=10000; % % fc=1000000; f=-D/2:D/2-1/(M*Ts*fd); t1=0:D-1*Ts; for m=1:D ruu_c(m)=ruu(m); rucuc_c(m)=rucuc(m); end %多普勒频移 modelbd=2*besselj(0,wd*t1);%无穷阶自相关 ruu fmodelbd=fft(modelbd); % bd=fft(ruu); bd=fft(ruu_c); % for j=1:1853 % Sf(j)=1.5/(pi*fd*sqrt(1-(f(j)-fc)/fd)2); % end figure(4) % subplot(326); % plot(f,abs(Sf); plot(f,mat2gray(abs(bd); % plot(f,abs(fmodelbd); % plot(fd*t,abs(bd); hold on plot(f,mat2gray(abs(fmodelbd),r); title(Doppler Spectrum,f_d=,num2str(fd),Hz) %axis(-5 5 0 1000)； xlabel(f/f_d), ylabel(Magnitude) legend(8 阶, 无穷) %ruu1=xcorr(uc); %for i=1:Ns %ruu11(i)=ruu1(Ns-1+i); %end figure(2) %自相关 rucuc_if=besselj(0,wd*t1)+besselj(4,wd*t1); % plot(fd*t,rucuc); plot(t1,rucuc_c); hold on; plot(t1,rucuc_if,r); title(Rucuc, f_d=,num2str(fd),Hz) grid on, xlabel(tau), ylabel(Correlation) legend(8 阶, 无穷) figure(3) % plot(fd*t,abs(ruu); plot(t1,abs(ruu_c); hold on; plot(t1,abs(modelbd),r); %subplot(325), plot(t,abs(ruu