minimum-norm-Doppler-frequency.zip_Doppler frequency_doppler

%LPNM %参考自相关 %f_max = 333; clear; f_max = 333; thera = 1; t1 = 0:0.0001:0.5 Y= BesselJ(0,2*pi*f_max*t1); r1= thera^2.*Y; %plot(t1,r1,'r');hold on; axis([0 0.015 -1 1]); % tau_max sigma_0= 1;sigma_0_2=1;N=1000;p=2;s_o=1; Ni=9; ts=1/(5*10^5); [f_i_n,c_i_n]=LPNM_opt_Jakes(N,f_max,sigma_0,p,Ni,s_o); [f_i_n2,c_i_n2]=LPNM_opt_Jakes(N,f_max,sigma_0,p,Ni+1,s_o); %k = round(f_i_n*ts*2^21) th=rand(Ni,1)*2*pi; kmax=500000; k = (1:kmax); th2=rand(Ni+1,1)*2*pi; u = 0; for i=1:Ni u= u+ c_i_n(i)*cos(2*pi*f_i_n(i)*k*ts+th(i)); end u2 = 0; for i=1:Ni+1 u2= u2+ c_i_n2(i)*cos(2*pi*f_i_n2(i)*k*ts+th2(i)); end r = 0.01:0.01:6; norm_u = abs (u+j*u2); for l=1:length(r) F(l)=length(find(norm_u<=r(l))); end F=F/length(norm_u); figure; plot(r,F);hold on; rf = 1- exp(-r.^2/2); %取theta=1 plot(r,rf,'k'); legend('仿真模型','参考模型'); figure; step = 50; c1= 0:step:10000; to = c1*ts; B=zeros(1,length(c1)); qrelati = zeros(1,length(c1)); for i = 1:length(c1) for j = 1:kmax-(i-1)*step qrelati(i) = qrelati(i) + u(j)*u((i-1)*step+j); end qrelati(i) = qrelati(i) / kmax; B(i)=qrelati(i); end; % plot(i*ts,qrelati(i),'r');hold on; % A=ts.*[1:10000]; plot(to,B);hold on; %高斯分布 thera = 1; Y= BesselJ(0,2*pi*f_max*to); r= thera^2.*Y'; plot(to,r,'k');