Description du signal en dB
Bonjour,
Pouvez vous m'aider à décrire et visualiser le signal Sr(X,Y) et S(i,j) en dB
merci
Code:
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
| for t=0e-8:60e-10:tmax;
for i=1:Ymax
for j=1:Xmax
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;
Ddr=sqrt((Xr-X0).^2+(Yr-Y0).^2);
tdr=Ddr./c;
Adr=1./Ddr;
Sr(i,j)= A0*Ad*exp(-(r-c*(t-t0))^2/2*sigma^2)+ sum(A0*(A*Adr.*Ark).*exp(-(rk-c*(t-t0-tdr-Thau)).^2/2*sigma^2));
end
end
title(['\fontsize{12}Temps = ',num2str(t),])
set(hs, 'cdata',Sr);
hp = plot(t, Sr(X, Y),'.k');
title('Evolution du signal en fonction du temps en un point M(X,Y) ')
title(['\fontsize{12}Temps = ',num2str(t),])
drawnow
hold on;
end |
Description du signal en dB
oui, je m'excuse j’étais pas clair,
enfaite, l'amplitude du signal est très petite sur le graph, du coup je me suis dit de visualiser le signal en décibel (dB)
Code:
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 |
Description du signal en dB
Bonjour,
j'ai essayé avec cette fonction y=mag2dby; mais ça ne marche pas, le programme se bloque
Avez-vous d'autres propositions ?
merci
Description du signal en dB
Description du signal en dB
c'est bien en ajoutant le 20 log10 que ca marche
Merci beaucoup
