1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
| close all
clear all
clc
c=3e8;
xDim=400;
yDim=400;
tmax=0.8e-6;
sigma=1;
A0=20;
t0=0;
%Coordonée emetteur
X0=100;
Y0=100;
%coordonée réflecteur
Xr=80;
Yr=120;
% Signal emis
%S(t)=A0*exp(-(t-t0)^2/2*sigma^2);
% Distance emetteur-->reflecteur
Ddk=sqrt((Xr-X0)^2+(Yr-Y0)^2);
%Retard emetteur--> reflecteur
tdk=Ddk/c;
%facteur d'attenuation retard propre
Ak=0.8;
tk=6e-8;
%%Coefficient d'attenuation emetteur -->reflecteur
Adk=1/Ddk;
% position du récepteur
X = 150;
Y = 80;
h_recepteur =0.*(0.5e-8:tmax);
figure(1)
subplot(2,2,[2,4])
hs = surf(zeros(400));
% colormap
%colorbar;
caxis([0 0.005])
view(2)
shading interp
axis square
zlim([0 0.1])
xlim([0 xDim])
xlabel('Axe Ox')
ylim([0 yDim])
ylabel('Axe Oy')
axis equal
title(['\fontsize{12}Temps = ',num2str(0),])
subplot(2,2,[1,3])
xlim([0 tmax])
V=axis;
axis([V(1) V(2) 0 0.2])
xlabel('Temps(Seconde)')
hold on
for t=0e-8:0.15e-8:tmax;
for j=1:yDim
for i=1:xDim
r=sqrt((i-X0)^2+(j-Y0)^2);
%Coefficient d'attenuation emetteur-->reflecteur
Ad=1/r;
rk=sqrt((i-Xr)^2+(j-Yr)^2);
%Coefficient d'attenuation reflecteur--> recepteur
Ark=1/rk;
Sr(i,j)= A0*Ad*exp(-(r-c*(t-t0))^2/2*sigma^2)+A0*(Ak*Adk*Ark)*exp(-(rk-c*(t-t0-tdk-tk))^2/2*sigma^2);
end
end
set(hs, 'cdata',Sr);
%set(gca,'ylim',[0 0.1],'ytick',0:0.01:0.1)
hp = plot(t, Sr(X, Y),'.k');
title('Evolution du signal en fonction du temps en un point M(X,Y) ')
drawnow
hold on;
end |
Partager