0

我的项目降噪。

1.这是我从麦克风输入的代码并保存到.wav

读入文件

clear all;

close all;

mic1= dsp.AudioRecorder;
hmfw = dsp.AudioFileWriter('myspeech.wav','FileFormat','WAV');
disp('Speak into microphone now');
time_end = 10;
tic;
while toc <= time_end
    step(hmfw, step(mic1));
end

release(mic1);
release(hmfw);`

disp('Recording complete');
[f,fs] = audioread('C:\Users\Admin\Documents\MATLAB\myspeech.wav');`

在记录之前我如何绘制频谱图

时间 = 10 频率 = 0 - 8000

2.如何在频率 500 - 2000 Hz 之间绘制降噪后的频谱图

像这样图谱图和光谱

这里链接

https://www.google.co.th/search?q=spectrogram+matlab&newwindow=1&rlz=1C1JPGB_enTH637TH637&espv=2&source=lnms&tbm=isch&sa=X&ei=PZliVf3GMcuIuATCroIg&ved=0CAcQ_AUoAQ&biw=1920&bih=979#imgrc=s05zemtY2IFy7M%253A%3BBgWyAcO6UJomJM%3Bhttp% 253A%252F%252Fwww.aquaphoenix.com%252Flecture%252Fmatlab10%252Fimages-large%252Fmatlab_audio_funky_plot_spectrogram.jpg%3Bhttp%253A%252F%252Fwww.aquaphoenix.com%252Flecture%252Fmatlab10%602Fpage478%3B

这是我的过滤器代码。

n = 7;
beginFreq = 500 / (fs/2);
endFreq = 2000 / (fs/2);
[b,a] = butter(n, [beginFreq, endFreq], 'bandpass');

谢谢你。

4

1 回答 1

0

试试这个：

**** The previous part of your code goes here ***
disp('Recording complete');
[y,fs] = audioread('myspeech.wav');

y = y(:,1); % Take only one channel from the recording


n = 7;
beginFreq = 500 / (fs/2);
endFreq = 2000 / (fs/2);
[b,a] = butter(n, [beginFreq, endFreq], 'bandpass');

%
dt = 1/fs;
L = size(y,1);                    % Length of signal
t = (0:L-1)'*dt;                  % Time vector

figure(1)

plot(t,y(:,1))
title('Recording')
xlabel('time (milliseconds)')

NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y(:,1),NFFT)/L;
f = fs/2*linspace(0,1,NFFT/2+1);

% Plot spectrum of channel 1 from the original signal
figure()
plot(f,2*abs(Y(1:NFFT/2+1))); hold on
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')

% Filter the two channels
filt_y = filter(b,a,y);

filt_Y = fft(filt_y,NFFT)/L;

% Plot the filtered spectrum of channel 1 from the original signal
plot(f,2*abs(filt_Y(1:NFFT/2+1))) 
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
于 2015-05-24T10:38:29.873 回答