Matlab音频信号的基本处理与分析

Matlab音频信号的基本处理与分析

1 音频信号的读取与时域分析

原文链接：Matlab音频信号的基本处理与分析

close all
clear all
clc
 
%读取音频文件
info =audioinfo('sample_orig.mp3');%获取音频文件的信息
[audio,Fs] = audioread('sample_orig.mp3');%读取音频文件
sound(audio,Fs);%播放音频文件
audiolength = length(audio);%获取音频文件的数据长度
t = 1:1:audiolength;
 
figure(1),
plot(t,audio(1:audiolength));
xlabel('Time');
ylabel('Audio Signal');
title('原始音频文件信号幅度图');

其中info为此音频的详细信息包括数据位宽，采样率等。

2 音频的频域分析

%音频的频域
y =fft(audio);
y=20*log(abs(y))/log(10);        %换算成dBW单位
f=[0:(Fs/audiolength):Fs/2];     %转换横坐标以Hz为单位
y=y(1:length(f));
 
figure(2),
plot(f,y);
xlabel('频率(Hz)','fontsize',8,'position',[180,-34.38,0]); 
ylabel('功率(dBW)');
title('原始音频文件信号频谱图');legend('Audio');

由此可见此音频的主要频率集中在低频段。

3 音频的叠加与写出

%音频信号的叠加
[audio1,Fs1] = audioread(‘test.wav’);%读取音频文件

audio2 = audio1(1:audiolength);
audio3 = audio2+0.1*audio;
sound(audio3,Fs);%播放音频文件
audiowrite(‘test1.wav’,audio3,Fs)%写出新的音频文件

4 音频信号的拼接

%音频信号的拼接
audio4  = [audio;audio1];
sound(audio4,Fs);%播放音频文件
audiowrite('test2.wav',audio3,Fs)%写出新的音频文件