matlab入门 绘制mp3文件时频波形

一、读取mp3文件，绘制时域和频域波形

%读取音乐文件
%moon.mp3音频文件长度280s
%返回结果fs的长度为12391536*2
%采样率44100
%12391536/44100=280和实际音频长度一致
[y,Fs]=audioread('moon.mp3');
N = length(y);

%计算音频时域绘图数据
t = (0:N-1)/Fs;

%计算频域绘图数据
Y = fft(y(:,1));
f = (0:N-1)*Fs/N;

%==========================================================================
%截取mp3音频前一半
halfN = fix(N/2);
halfy = y(1:halfN,1);

%计算音频时域绘图数据
halft = (0:halfN-1)/Fs;

%计算频域绘图数据
halfY = fft(halfy);
halff = (0:halfN-1)*Fs/halfN;

%==========================================================================
%画图
%完整音频
subplot(221);plot(t,y(:,1));%时域波形
subplot(222);plot(f,abs(Y));%频域波形

%前一半音频
subplot(223);plot(halft,halfy);%时域波形
subplot(224);plot(halff,abs(halfY));%频域波形

二、播放PCM音频

%读取音乐文件
%moon.mp3音频文件长度280s
%返回结果fs的长度为12391536*2
%采样率44100
%12391536/44100=280和实际音频长度一致
[y,Fs]=audioread('moon.mp3');
N = length(y);

%==========================================================================
%播放音频
player = audioplayer(y, Fs, 16);
%开始播放
play(player);
%停止
stop(player);