在matlab中编程（如何实时处理）

Question

我想在matlab中制作光谱图， 这是我的代码：

% Record your voice for 100 seconds.
recObj = audiorecorder;
disp('Start speaking.')
recordblocking(recObj, 100);

% Store data in double-precision array.
my= getaudiodata(recObj);
figure;
specgram(my,512);

问题是，当我说话时我想要显示频谱图，所以它应该在我说话时更新。当音频来自麦克风时，如何绘制频谱图？所以我应该能够实时看到频谱图

我也试过这个

% Record your voice for 100 seconds.
recObj = audiorecorder;
disp('Start speaking.')
a=0;
figure;
while a<60
    recordblocking(recObj, 100);

    % Store data in double-precision array.
    my= getaudiodata(recObj);

    specgram(my,512);
    a=a+1;
end

但它只会在while循环捕获时显示频谱图（因此在运行60次后）

Answer 1

MATLAB本质上是单线程的。一次只能发生一件事。这使得实时任务有些困难。如您所述，recordblocking在100秒过去之前不会将控制权返回给您的脚本。关键在于阻止。

解决这个问题的方法是使用callbacks and non-blocking functions。 audiorecorder对象有一些方法和属性可以实现这种行为。

audiorecorder properties：
- StartFcn：设置一个在启动异步记录时执行的功能
- StopFcn：停止录制时要执行的功能
- TimerFcn：录制期间每TimerPeriod秒执行的功能。

然后record method将在后台开始录制和处理，按指示调用上述功能。

通过定期更新回调函数中的数据，您可以更新您的绘图。不幸的是，以允许实时更新的有效方式做同样不容易。但这应该让你开始。

Answer 2

这是一种可能的实现方式。主要问题是你忘记在每个循环结束时调用DRAWNOW：

Fs = 8000;                    %# sampling frequency in Hz
T = 1;                        %# length of one interval signal in sec
t = 0:1/Fs:T-1/Fs;            %# time vector
nfft = 2^nextpow2(Fs);        %# n-point DFT
numUniq = ceil((nfft+1)/2);   %# half point
f = (0:numUniq-1)'*Fs/nfft;   %'# frequency vector (one sided)

%# prepare plots
figure
hAx(1) = subplot(211);
hLine(1) = line('XData',t, 'YData',nan(size(t)), 'Color','b', 'Parent',hAx(1));
xlabel('Time (s)'), ylabel('Amplitude')
hAx(2) = subplot(212);
hLine(2) = line('XData',f, 'YData',nan(size(f)), 'Color','b', 'Parent',hAx(2));
xlabel('Frequency (Hz)'), ylabel('Magnitude (dB)')
set(hAx, 'Box','on', 'XGrid','on', 'YGrid','on')
%#specgram(sig, nfft, Fs);

%# prepare audio recording
recObj = audiorecorder(Fs,8,1);

%# Record for 10 intervals of 1sec each
disp('Start speaking...')
for i=1:10
    recordblocking(recObj, T);

    %# get data and compute FFT
    sig = getaudiodata(recObj);
    fftMag = 20*log10( abs(fft(sig,nfft)) );

    %# update plots
    set(hLine(1), 'YData',sig)
    set(hLine(2), 'YData',fftMag(1:numUniq))
    title(hAx(1), num2str(i,'Interval = %d'))
    drawnow                   %# force MATLAB to flush any queued displays
end
disp('Done.')

我只是在每次迭代中显示频率成分。您应该能够修改它以显示光谱图，如果您想...

Answer 3

最明显的事情是将代码放在循环中以继续更新数字。但是请注意，Matlab并不是专门为这类任务设计的，所以我不知道你会有多少成功。您是否尝试使用谷歌搜索免费软件为您执行此操作？如果没有任何东西没有这样做，我会感到惊讶。