A=sqrt(2);
E=1;
N=100;
L=500;
fc=1e9;
fs=5e9;
vc=3e8;
M = 1e6;
x_o=0; y_o=0;
A_r = zeros(M,1);
for i=1:M
x_or=rand(N,1)*2000-1000;
y_or=rand(N,1)*2000-1000;
d=sqrt(x_or.^2+y_or.^2)+sqrt((x_or-500).^2+(y_or).^2);
temp=d/(vc/fc);
po=(temp-floor(temp))*2*pi;
A_r(i)=abs(sum(A*exp(j*po)));
end
[f_est,xi]=ksdensity(A_r);
f_th = raylpdf(xi, sqrt(N*E));
plot(xi, f_th, 'r', xi, f_est, 'b');
legend('Rayleigh PDF', 'Statistics');
grid on;
x_ob=rand(N,1)*2000-1000;
y_ob=rand(N,1)*2000-1000;
Dt = 500e-5;
t = (0:L-1)*Dt;
tt=repmat(t,N,1);
vv = [1.0, 10.0];
A_r_v = zeros(length(vv), L);
for ii = 1:length(vv)
v = vv(ii);
xx_r = 500 - v*tt;
yy_r = 0;
xx_ob = repmat(x_ob, 1, L);
yy_ob = repmat(y_ob, 1, L);
d1=sqrt(xx_ob.^2+yy_ob.^2);
d2=sqrt((xx_ob-xx_r).^2+(yy_ob-yy_r).^2);
d = d1+d2;
temp=d/(vc/fc);
psi=(temp-floor(temp))*2*pi;
A_r_v(ii,:)=abs(sum(A*exp(j*psi)));
end
figure;
plot(t, A_r_v(1,:), 'm', t, A_r_v(2,:), 'b');
grid on;
legend('1m/s ', '10m/s');
xlabel('t'); ylabel('A');
评论0