MATLAB声音信号处理.docx

MATLAB声音信号处理.docx《数字信号处理》课程设计设计题目:基于MATLAB的语音信号频谱分析指导老师:岳学东 专业:测控技术与仪器 姓名:M 班级:2> 学号:20072240228 日期: closeall;clearall;clc[s/Fs,bit]=wavread(,');a=s(:,l);al=a(10000:89999);siglengthl=length(al);tl=(O:siglengthl-l)/Fs;figure(l);subplot(2z3,l)plot(tl,al);xlabelfTimefs)');titled截取信号序列J;%wavread()函数读収声音信号%做单声道处理%-- 510:-- 54:-- 5Time(s)Al=fft(al);wl=2/siglengthl*[0:(siglengthl-l)];subplot(2,3,4);plot(wl,abs(Al));xlabel('w');title?截取信号频谱〕;gridon;%sound(s,Fs);**************************获取截取佶号的频谱截取信号频谱4003002001000Time(s)w****************************012%4:1简抽样Nl=4;siglength2=siglengthl/Nl;Time(s)wa2=zeros(siglength2J);***********************************************************for(i=l:siglength2)siglength3=siglengthl/N2;a2(i)=al(Nl*i);a3=zeros(siglength3,l);end;for(i=l:siglength3)对信号进行10倍的简抽样并做fftt2=(0:siglength2-l)/Fs*Nl;a3(i)=al(N2*i);变换获得其频谱,可以观察到幅频subplot(2,3,2)响应曲线己经发生混跌plot(t2za2);end;xlabel('Tjme(s)');对截取信号做4倍的简抽样t3=(O:siglength3-l)/Fs*N2;title('10:l简抽样信号序列J;并获収其频谱:subplot(2,3,3)可以观察到,当时域信号做plot(t3,a3);4倍的简抽样,其频谱未发生混xlabelfTime(s)');A2=fft(a2);柱。title('4:l简抽样信号序列J;w2=^siglength2*(0:siglength2-l)subplot(2,3,5);A3=fft(a3);plot(w2,abs(A2));w3=2;/siglength3*(0:siglength3-l)xlabel('w');subplot(2,3,6);titleC10:l简抽样信号频谱J;plot(w3,abs(A3));gridon;*****************************xlabel('w');titleC4:l简抽样信号频谱J;%10:l简抽样gridon;*******************************N2=10;nl=O:(siglengthl-l)bl=cos(*pi*nl)**********************************调制载波,其角频率依据截取信号频谱通带截止频率获得**********************************cl=bl'.*alfigure(2)nl=O:(siglengthl-l)bl=cos(*pi*nl)cl=bl'.*alfigure(2)subplot(2,2,l)plot(tl,cl);xlabel('Time(s)');title/调制后信号J;Cl=fft(cl);subplot(2,2,3)plot(wl,abs(Cl));xlabelCw');title(调制后信号的频谐);gridon;************************************Time(s)用载波对信号进行调制,并对其做fft变换c2=blf.*clsubplot(2z2,2)plot(tl,c2);xlabelfTimefs)1);********************************************信号解调处理,观察解调后信号的频谱可知,-- 2 4 6