MATLAB声音信号处理.docVIP

___________________________________________________________ 设计题目：基于MATLAB的语音信号频谱分析 指导老师：岳学东___________________________ 专业：测控技术与仪器 ______ 姓名：张猛___________________ 班级：2班___________________ 学号：20072240228____________ close all; clear all; clc [s,Fs,bit]=wavread(xinyueshenhua.wav); %wavread()函数读取声音信号 a=s(:,1); %做单声道处理 a1=a(10000:89999); %截取80000个信号采样点 siglength1=length(a1); t1=(0:siglength1-1)/Fs; figure(1); subplot(2,3,1) plot(t1,a1); xlabel(Time(s)); title(截取信号序列); A1 = fft(a1); ************************** w1=2/siglength1*[0:(siglength1-1)]; subplot(2,3,4); plot(w1,abs(A1)); 获取截取信号的频谱 xlabel(w); title(截取信号频谱); grid on; %sound(s,Fs); **************************** %4:1简抽样 N1=4; siglength2=siglength1/N1; a2=zeros(siglength2,1); **************************** for(i=1:siglength2) a2(i)=a1(N1*i); end; t2=(0:siglength2-1)/Fs*N1; subplot(2,3,2) plot(t2,a2); xlabel(Time(s)); 对截取信号做4倍的简抽样 title(10:1简抽样信号序列); 并获取其频谱： 可以观察到，当时域信号做4倍的简抽样，其频谱未发生混 A2=fft(a2); 叠。 w2=2/siglength2*(0:siglength2-1) subplot(2,3,5); plot(w2,abs(A2)); xlabel(w); title(10:1简抽样信号频谱); grid on; ***************************** %10:1简抽样 N2=10; ******************************* siglength3=siglength1/N2; a3=zeros(siglength3,1); for(i=1:siglength3) 对信号进行10倍的简抽样并做fft a3(i)=a1(N2*i); 变换获得其频谱，可以观察到幅频响应曲线已经发生混跌 end; t3=(0:siglength3-1)/Fs*N2; subplot(2,3,3) plot(t3,a3); xlabel(Time(s)); title(4:1简抽样信号序列); A3=fft(a3); w3=2/siglength3*(0:siglength3-1) subplot(2,3,6); plot(w3,abs(A3)); xlabel(w); title(4:1简抽样信号频谱); grid on;