Description du signal en dB

**adnino** · 26/03/2014, 15h36

Bonjour,
Pouvez vous m'aider à décrire et visualiser le signal Sr(X,Y) et S(i,j) en dB

merci

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**Newenda** · 26/03/2014, 16h55

Bonjour,

Ta demande n'est pas clair.. comment peut on t'aider à décrire une image qu'on a pas?

As tu un problème de code ? d'unité ? d'affichage des données?

**adnino** · 26/03/2014, 17h01

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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