clear all;close all;clc Fs = 22050; [fname,pname] = uigetfile('小提琴.wma');%获取所选文件信息:文件名、路径 file = [pname,fname]; [x1,Fs] = audioread(file);%x1为所读取的音频数据,Fs为采样频率 sound(x1,Fs); %播放音乐 figure(1); subplot(211); plot(x1);%做原始语音信号的时域图形 title('原始语音信号') xlabel('采样点 n'); ylabel('音量 n'); y1=fft(x1); %做length(x1)点的FFT y1=fftshift(y1);%频率分量将会移到坐标中心 subplot(212); plot(abs(y1));%画出原始语音信号的频谱图，这里保证了x轴的点数必须和y轴点数一致 title('原始语音信号的频谱'); grid on; [y,fs]=audioread('小提琴.wma'); %读取音乐信号 size(y) %看音乐信号是双列还是单列 y1=y(:,1); %取单列 fs=fs %获取音乐信号的抽样率 N=length(y1); %音乐信号的长度 Fy1=fft(y1,N); %对信号做快速傅里叶变换 w=2/N*[0:N-1]; figure %画音乐信号的波形和频谱 subplot(2,1,1);plot(y1); grid on title('音乐信号的波形') xlabel('time/s') ylabel('Magnitude') subplot(2,1,2);plot(w,abs(Fy1)); grid on title('音乐信号的频谱') xlabel('Frequency/pi') ylabel('Magnitude') sound(y1,1/2*fs) %以抽样率1/2fs播放音乐信号

这段代码是用MATLAB实现的音频处理程序，主要包括以下功能：

1. 读取音频文件并播放。
2. 绘制原始语音信号的时域图形和频谱图。
3. 读取音乐信号并绘制波形图和频谱图。
4. 使用快速傅里叶变换对音乐信号进行频谱分析。
5. 播放音乐信号。

值得注意的是，该代码只适用于单声道音频文件。如果读取的是双声道音频文件，需要将其转换为单声道才能进行处理。此外，在播放音频文件时，需要根据采样率设置合适的播放速度，否则可能会导致音频变形。