n = 2^11; delta = 2*pi/n; t = 0:delta:2*pi-delta; c = sin(2*pi*65.41*t); d = sin(2*pi*73.42*t); e = sin(2*pi*82.41*t); f = sin(2*pi*87.3*t); g = sin(2*pi*98*t); a = sin(2*pi*110*t); b = sin(2*pi*123.47*t); cc = sin(2*pi*65.41/2*t); gg = sin(2*pi*98/2*t); ggnoise = gg.*(1.5*rand(1,n)); ce3k = [d e c cc gg]; ce3knoise = ce3k.*(1.5*rand(1,n*5));