基于MATLAB对语音信号的处理

%% MATLAB对语音信号处理 %采集音乐信号 Fs = 44100; samples1 = [1,10*Fs]; [audio_data, fs] = audioread('Untitled.flac',samples1); L = length(audio_data); % 听一听该音频 % soundsc(audio_data, fs); %绘制时域图 figure('Name','Original'); tt = (1 : L) * 1/fs; %时间(纵坐标) subplot(4, 1, 1);plot(tt, audio_data);grid; title('时域图');xlabel('时间（s）');ylabel('Amp'); %频域分析 L = 2^nextpow2(L); % 先从原始信号长度确定下一个 2 次幂的新长度,用尾随零填充信号以改善 fft 的性能 Spectrum = fft(audio_data,L); % 快速傅里叶变换，转换为频域 f_x = fs*(0:(L/2))/L; % 定义频域，即变换横坐标量纲为频率 mag_FFT_audio_data = abs(Spectrum / L); % 取幅值，并归一化 subplot(4, 1, 2);plot(f_x,mag_FFT_audio_data(1 : L/2+1)) title('频域图');xlabel('频率（Hz）');ylabel('mag'); % 带通滤波器 (由fdatool生成) Fstop1 = 7200; % First Stopband Frequency Fpass1 = 9600; % First Passband Frequency Fpass2 = 12000; % Second Passband Frequency Fstop2 = 14400; % Second Stopband Frequency Dstop1 = 0.001; % First Stopband Attenuation Dpass = 0.057501127785; % Passband Ripple Dstop2 = 0.0001; % Second Stopband Attenuation dens = 20; % Density Factor % Calculate the order from the parameters using FIRPMORD. [N, Fo, Ao, W] = firpmord([Fstop1 Fpass1 Fpass2 Fstop2]/(Fs/2), [0 1 ... 0], [Dstop1 Dpass Dstop2]); % Calculate the coefficients using the FIRPM function. b = firpm(N, Fo, Ao, W, {dens}); Hd = dfilt.dffir(b); %进行滤波 filtered_audio_data = filter(Hd,audio_data); subplot(4, 1, 3);plot(tt, filtered_audio_data);grid; title('滤波后的时域图');xlabel('时间（s）');ylabel('Amp'); %对滤波之后的信号进行频域分析 Filtered_Spectrum = fft(filtered_audio_data,L); % 快速傅里叶变换，转换为频域 %f_x = fs*(0:(L/2))/L; % 定义频域，即变换横5坐标量纲为频率 Filtered_mag_FFT_audio_data = abs(Filtered_Spectrum / L); % 取幅值，并归一化 subplot(4, 1, 4);plot(f_x,Filtered_mag_FFT_audio_data(1 : L/2+1)) title('滤波之后频域图');xlabel('频率（Hz）');ylabel('mag'); %听一下滤波之后的音乐 %soundsc(filtered_audio_data, fs);