rise/fall time
显示 更早的评论
Hi, can anyone help me to generate 2 sounds of 500 and 1000 Hz with 25 ms duration with ISI of 30 ms and rise/fall time of 10 ms with Fs of 44100 and I want to 500 Hz be 10 dB less than 1000 Hz
采纳的回答
Shae Morgan
2020-7-31
编辑：Shae Morgan
2020-7-31
%set-up parameters
fs=44100; %CD quality - also conveniently divisible by 30 and 25
stim_dur=.025; %duration in seconds
ISI_dur=.030; %ISI duration in seconds
ramp_dur=.010; %ramp duration in seconds
t=0:1/fs:stim_dur-1/fs;
f1=500;
f2=1000;
%create signals
xt1=sin(2*pi*f1*t);
xt2=sin(2*pi*f2*t);
%setup ramp
rampSamps = floor(fs*ramp_dur);
window=hanning(2*rampSamps)'; %hanning window is cosine^2 this will change depending on the kind of ramp you want
w1=window(1:ceil((length(window))/2)); %use the first half of hanning function for onramp
w2=window(ceil((length(window))/2)+1:end); %use second half of hanning function of off ramp
w1 = [w1 ones(1,length(xt1)-length(w1))];
w2 = [ones(1,length(xt1)-length(w2)) w2];
%ramp stimuli
xt1_ramped = xt1.*w1.*w2;
xt2_ramped = xt2.*w1.*w2;
%generate ISI
ISI=zeros(ISI_dur*fs,1);
%final stimulus
out=[xt1_ramped';ISI;xt2_ramped'];
plot(out)
9 个评论
Shae Morgan
2020-7-31
Edited for transposition to allow concatenation*
I hope this is what you're looking for!
Parisa
2020-7-31
Shae Morgan
2020-7-31
See below
sound(out,fs) %play the sound
audiowrite('filename.wav',out,fs) %save it
If this answer works for you, please accept! :)
Parisa
2020-7-31
Parisa
2020-8-2
Shae Morgan
2020-8-3
编辑：Shae Morgan
2020-8-3
sure,
so the amplitude of each sine wave is a multiplier added to the sin function
amp = 1 by default, but since we added a ramp to these stimuli, their total average loudness will vary - so the question is, do you want the peak (after the ramp on but before the ramp off) to be 10 dB louder? or do you want the overall ramped stimulus to be 10 dB louder?
%% option 1
% %scale stimuli and then ramp
amp1=1;
amp2=10^(-10/20);
%create signals
xt1=amp1*sin(2*pi*f1*t);
xt2=amp2*sin(2*pi*f2*t);
%check
rms1=rms(xt1);
rms2=rms(xt2);
dB_diff=20*log10(rms1/rms2);
%proceed with th rest of the code and perform the ramping
%% option 2
%ramp stimuli and then scale
rms_xt1_ramped=rms(xt1_ramped); %check the average intensity of the ramped 500 Hz signal
amp2_10dB=10^(-10/20)*rms_xt1_ramped; %find the new amplitude that is 10 dB lower
%normalize the xt2_ramped stimulus and then scale to -10 dB
xt2_ramped_norm=xt2_ramped./rms(xt2_ramped);
xt2_ramped_10dB=xt2_ramped_norm.*amp2_10dB;
%check
rms_xt2_ramped=rms(xt2_ramped_10dB);
20*log10(rms_xt1_ramped/rms_xt2_ramped)
Parisa
2020-8-3
Shae Morgan
2020-8-7
I looked at your question and submitted an answer - there you asked for a 20 ms ISI, so I answered the question you asked there, not here. Please update the question if it isn't what you need and I'll be happy to submit an answer.
If the solution works, please choose to accept my answer! It helps me a lot as I try to help you!
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Logical 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
