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
| %%%%%%%%%%%%%%%%
clear all
close all
N0=1;
Z=1/(2*pi*N0);
delta_P0=1;
delta_Q0=N0/delta_P0;
moyenne_Q0=1;
moyenne_P0=1;
%%%%%%%%%%%%%%%%%
delta_P1=1;
delta_Q1=N0/delta_P0;
moyenne_Q1=2;
moyenne_P1=2;
%%%%%%%%%%%%%%%
i=0
for q=0:0.1:5, i=i+1; j=0;
for p=0:0.1:5,
j=j+1;
X01=(q-moyenne_Q0).^2/(2*delta_Q0^2);
X00=(p-moyenne_P0).^2/(2*delta_P0^2);
w0(i,j)=Z.*exp(-X00-X01);
end,
end,
i=0;
for q=0:0.1:5, i=i+1; j=0;
for p=0:0.1:5,
j=j+1;
X11=(q-moyenne_Q1).^2/(2*delta_Q1^2);
X10=(p-moyenne_P1).^2/(2*delta_P1^2);
w1(i,j)=Z.*exp(-X10-X11);
end,
end,
q=0:0.1:5;
p=0:0.1:5;
%figure(1)
hold on
imagesc(q,p,w1)
figure(2)
hold off
imagesc(q,p,w0) |
Partager