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;
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);%-pi~pi均匀分布随机变量 %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('time[s]'), ylabel('Magnitude[dB]')
%幅度分布
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('Phase[rad]'), 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);
mp3手表 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(ruu11));
%axis([0 0.015 0 10]);
title(['Ruu, f_d=',num2str(fd),'Hz'])
grid on, xlabel('\tau'), ylabel('Correlation')
legend('8阶','无穷')
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议