rmeddis@38: frequencies=[1000 1250];
rmeddis@38: amplitudesdB=[20 23];
rmeddis@38: nFrequencies=length(frequencies);
rmeddis@38: 
rmeddis@38: dt=0.0001;
rmeddis@38: 
rmeddis@38: toneDuration=.010;
rmeddis@38: time=dt:dt:toneDuration;
rmeddis@38: 
rmeddis@38: % fixed ramp, silenceDuration, toneDuration
rmeddis@38: rampDuration=0.005;
rmeddis@38: rampTime=dt:dt:rampDuration;
rmeddis@38: ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ...
rmeddis@38:     ones(1,length(time)-length(rampTime))];
rmeddis@38: ramp=ramp.*fliplr(ramp);
rmeddis@38: 
rmeddis@38: silenceDuration=0.010;
rmeddis@38: silenceDurationLength=round(silenceDuration/dt);
rmeddis@38: initialSilence=zeros(1,silenceDurationLength);
rmeddis@38: 
rmeddis@38: silenceToneDuration=toneDuration + silenceDuration;
rmeddis@38: silenceToneDurationLength=round(silenceToneDuration/dt);
rmeddis@38: 
rmeddis@38: totalDuration=silenceToneDuration*nFrequencies;
rmeddis@38: totalDurationLength=round(totalDuration/dt);
rmeddis@38: stimulus=zeros(1,totalDurationLength);
rmeddis@38: toneBeginPTR=1;
rmeddis@38: 
rmeddis@38: for i=1:nFrequencies
rmeddis@38:     frequency=frequencies(i);
rmeddis@38:     dBSPL=amplitudesdB(i);
rmeddis@38:     amplitude=28e-6* 10.^(dBSPL/20);
rmeddis@38:     tone=amplitude*sin(2*pi*frequency*time);
rmeddis@38:     tone=tone.*ramp;
rmeddis@38:     stimulus(toneBeginPTR:toneBeginPTR+silenceToneDurationLength-1)=...
rmeddis@38:         [initialSilence tone];    
rmeddis@38:     toneBeginPTR=toneBeginPTR+silenceToneDurationLength;
rmeddis@38: end
rmeddis@38: figure(2), plot( stimulus')