% Simulation des franges de Young
close all; clear
% -------------------------------------------------------
lambda = 632e-9

a = 0.0004   % séparation des fentes
s = 1        % distance à l'écran
l = 0.00005  % largeur de chaque fente

y = -0.03:0.0003:0.03;  % variable le long de l'écran

Ii = 4* (cos(a*pi*y / (s*lambda))).^2;
clear I
figure
plot (y, Ii)
grid on
title ('Interférence par deux fentes a=0.4 mm, l=50um')

Id = (sinc1(pi*l/(lambda*s)*y)).^2;
figure;
plot (y, Id);  grid on;   title ('Effet de la diffraction')

Is = 1./ (s^2 + y.^2);
figure;
plot (y, Is);  grid on;   title ('Effet de la sphéricité')

I = Ii .* Id .* Is;
figure;
plot (y, I);  grid on;   title ('Franges de Young')

RGB = spectrumRGB (lambda*1e9);

rgb = I'/max(I) * RGB;
franges_lin (y, I, rgb)
title ('Intensité des franges de Young')

%-----------------------------

lam = (450:50:650)*1e-9;
color = 'bcgyr';
figure;   hold on
for i = 1:length(lam)
    Ii = 4* (cos(a*pi*y ./ (s*lam(i)))).^2;
    Id = (sinc1(pi*l./(lam(i)*s).*y)).^2;
    Is = 1./ (s^2 + y.^2);
    
    I(i,:) = Ii .* Id .* Is;
    plot (y, I(i,:), color(i));
end
hold off
grid on
legend ('bleu', 'turquoise', 'vert', 'jaune', 'rouge')

I_tot = sum(I,1) / length(lam);
figure; plot(y, I_tot)
grid on; title ('Intensité moyenne des franges multicolor')

% calcul des RGB pour les divers lambda
for i = 1:length(lam)
    rgbi(i,:) = spectrumRGB(lam(i)*1e9);
end
% calibration des RGB pour que la somme des 6 couleurs donne [1 1 1] = blanc
sum_RGB = sum(rgbi,1)
for i=1:3
    RGBi(:,i) = rgbi(:,i) ./ sum_RGB(i);
end

clear rgb_lam
for i = 1:length(lam)
    max_I = max(I(i,:));
    rgb_lam (i,:,:) =  I(i,:)'/max_I * RGBi(i,:);
end
rgb_tot = sum (rgb_lam,1);
franges_lin (y, I_tot, rgb_tot);
xlabel('x (m)'); ylabel('y (m)')
title ('Franges en lumière blanche')