0

我想在 Matlab 中创建一个连续的声音（40 Hz）并用节拍调制声音（参见下面的简单代码）。我的问题是我想在脑电图实验中使用这种声音作为刺激，但是节拍和持续的 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])
4

1 回答 1

0

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

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

当您将峰值（正弦波）与方波（具有无限频率）的变换进行卷积时，频率空间中的情况也是如此。

我认为最简单的解决方案不是使用方波，而是使用维基百科的公式和前几个术语创建自己的波（取决于您希望接近方形的程度）

于 2012-06-19T09:11:50.093 回答