Question

我想在Matlab中创建一个连续的声音（40赫兹）并用节拍调制声音（参见下面的简单代码）。我的问题是我想在EEG实验中使用这种声音作为刺激，但是节拍和持续的40 Hz频率会产生很多频率。如何减少频率？我还用fft来检查我用这种刺激产生的频率（太多）。有没有更好的方法来创造一个轻微的节拍？我做错了吗？我希望它足够清楚，谢谢

clear all

Accent = 3; % amplitude for the beat
Freq = 40; % frequency
Dur = 4; % Duration
Samp = 48000; % sampling points

%% signal
Sine = sin((1:Dur*Samp)*2*pi*Freq/Samp);
Beat = repmat([ones(1,Samp/8)*Accent ones(1,Samp/8) ones(1,Samp/8)*Accent  ones(1,Samp/8) repmat(ones(1,Samp/8),1,4)],1,Dur);

S = Sine.*Beat;
S = S/max(abs(S));
% plot(S)
% sound(S,Samp)
FT = abs(fft(S)); 

%% plots
subplot(2,1,1)
plot(S); title(['Signal ' num2str(Freq) ' Hz + Beat'])
subplot(2,1,2)
plot(((0:length(S)-1)/length(S))*Samp,FT);
title('Signal FFT')
xlabel('Frequency')
axis([0 45 0 1])

Answer 1

我认为你的问题在于数学。

您正在将正弦与sqare波相乘（请参阅此处：https://en.wikipedia.org/wiki/Square_wave）+偏移量（但这并不重要，因为它具有频率0）。

当您使用方波（具有无限频率）的变换对峰值（对于正弦）进行卷积时，频率空间中的情况相同。

我认为最简单的解决方案是使用维基百科中的公式并且只使用前几个术语（取决于你想要达到方形的接近程度）而不是方波来创建波浪。

在Matlab中用节拍编程声音

1 个答案: