function y = cloche(f0,t,Fe)
%        y = cloche(f0,t,Fe)
% cloche (meme enveloppe pour partiels)
%--------------------------------------

h = 1/Fe; th = 0:h:t;

%f = [117 234 278 358 468 578 610 620 700  788 968 1050 1130 1260 1372 1570 1890 2000 2110 2230 2580 2900]/234*f0;
f = [0.5  1 1.188 1.530 2.0000 2.470 2.607 2.650 2.991 3.367 4.137 4.487 4.829 5.385 5.863 6.709 8.077 8.547 9.017 9.530 11.026 12.393]*f0;
am = [350 950 500 150 700 100 250 370 1000 180 300 100  150  300  100   100  50   20   10   35   5   15];
n = length(f);

y = zeros(1,length(th));
for k = 1:n
    z = sin(2*pi*f(k)*th);
    y = y + am(k)*z;    
end

tr = [0 .005 .2 1.3 t*.8 t];
yr = [0  .6    1  .4 .1 0];
y = y.*envelop(tr,yr,Fe);
y = .98*y/(max(abs(y)));
