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%252Fpage4.html%3

ここに私のフィルターコードがあります。

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 に答える