timeLength = 0.1; % 采样时长，单位秒 samples = timeLength * 44100; % 默认采样率44100，计算采样点数 H = audioDeviceReader(... 'NumChannels', 1,... % 1个通道 'SamplesPerFrame', samples); % 采样点数 audioIn = H(); % 第一次采样 figure('Name','实时频谱','MenuBar'... ,'none','ToolBar','none','NumberTitle','off'); xdata=(1:1:samples/2)/timeLength; axes1= subplot(1,2,1); axes2= subplot(1,2,2); pic= plot(axes1, 1:1:samples, audioIn); % 初始化音频波形图 pic2= bar(axes2,xdata, xdata*0,'r'); % 初始化频谱图 set(axes1,'xlim', [0 samples], 'ylim', ... [-0.15 0.15],'XTick',[],'YTick',[] ); set(axes2,'xlim', [min(xdata) max(xdata)], 'ylim',[0 6] , ... 'xscale','log','XTick',[1 10 100 1e3 1e4],'YTick',[] ); xlabel(axes2,'频率 (Hz)'); xlabel(axes1,'波形'); axes2.Position=[0.040 0.48 00.92 0.48]; % 左，下，宽度，高度 axes1.Position=[0.040 0.06 0.92 0.25]; drawnow; while 3>2 [audioIn,Overrun] = step(H); % 采样 if Overrun > 0 warning(' 数据溢出 %d 位\n',Overrun); end ydata_fft=fft(audioIn); % 傅里叶变换 ydata_abs=abs(ydata_fft(1:samples/2));% 取绝对值 set(pic, 'ydata',audioIn); % 更新波形图数据 set(pic2, 'ydata',log(ydata_abs)); % 更新频谱图数据 drawnow; % 刷新 end

这段代码是用于实时绘制音频信号的波形图和频谱图。首先定义了采样时长和采样点数，然后通过audioDeviceReader对象来读取音频输入。接下来创建了一个图形窗口，分为两个子图，一个用于显示音频波形图，一个用于显示频谱图。初始化了波形图和频谱图的数据，并设置了它们的坐标轴范围和刻度。然后进入一个无限循环，不断采样音频信号并进行傅里叶变换，更新波形图和频谱图的数据，并刷新显示。

注意：此处代码只是一个示例，具体使用时可能需要根据实际情况进行适当修改。