数字信号处理.doc

_这个问题问的好!1）用mp3录音，生成文件cricket.wav，把该文件放到matlab文件夹里面（就是你打开matlab后中间顶部的地址）。 2）使用如下程序，做波形显示以及fft变换。 y,Fs,bits=wavread(cricket.wav);%读出信号，采样率和采样位数。 y=y(:,1);%我这里假设你的声音是双声道，我只取单声道作分析，如果你想分析另外一个声道，请改成y=y(:,2) sigLength=length(y); Y = fft(y,sigLength); Pyy = Y.* conj(Y) / sigLength; halflength=floor(sigLength/2); f=Fs*(0:halflength)/sigLength; figure;plot(f,Pyy(1:halflength+1);xlabel(Frequency(Hz); t=(0:sigLength-1)/Fs; figure;plot(t,y);xlabel(Time(s); 3）频率看频谱就有了，声音间隔看声音波形，周期看声音波形。 4）关于去噪声。 a）如果噪声是特定频率的周期噪声（periodic noise），比如说50hz，那么你可以用matlab的filter，作一个低通、高通、带通或者带阻滤波。 b）如果声音是高斯白噪声。那就用自适应滤波（adaptive filter,wiener filter）。这里涉及到对噪声的采样、计算特征值以及决定阶数的问题。 c）幸好我们可以“耍赖”用cool editor。用它打开wav文件，用鼠标把一段噪声圈起来，采样，然后直接选择去噪就可以了。各大网站有介绍。 例子：matlab去除50hz噪声。 我用电脑录了一段声音，里面有50hz的周期噪声（因为受交流电干扰）。而我自己的声音频率最低是90hz。我使用了一个10阶butterworth高通滤波器，边带是70hz(介于50跟90之间)。 问题是，这不能直接用。因为声音文件的采样率是22k，70相对于22k来说太小了。所以我得先把我的声音欠采样，然后再滤波，然后再插值。程序如下。 k,Fs,bits=wavread(mywav.wav); k=k(:,1); y_temp=k(1:90000); dfactor=3; y=decimate(y_temp,dfactor); b,a = butter(10,70/(Fs/(dfactor*2),high); y=filter(b,a,y); y=interp(y,dfactor); sigLength=length(y); Y = fft(y,sigLength); Pyy = Y.* conj(Y) / sigLength; halflength=floor(sigLength/2); f=Fs*(0:halflength)/sigLength; figure;plot(f,Pyy(1:halflength+1);xlabel(Frequency(Hz); sigLength=length(y_temp); Y = fft(y_temp,sigLength); Pyy = Y.* conj(Y) / sigLength; halflength=floor(sigLength/2); f=Fs*(0:halflength)/sigLength; figure;plot(f,Pyy(1:halflength+1);xlabel(Frequency(Hz); t=(0:sigLength-1)/Fs; figure;plot(t,y,t,y_temp);xlabel(Time(s); wavplay(y,Fs); wavplay(y_temp,Fs); 5）回放：使用wavplay函数 wavplay(y,Fs)；谁帮我分析这个matlab语音信号去噪处理程序的算法20标签：matlab,语音信号,算法 帮我在必要的比方注释一下即可，谢谢，多加分给你们clear all; close all;x,Fs,nbits=wavread(F:论文FIR音频滤波器rockroll.wav);figure(10),y1=fft(x);plot(abs(y1)xt1,Fs=wavread(F:论文FIR音频滤波器rockroll.wav);N=max(size(xt1);Ts=1/Fs;t=0:Ts:(N-1)*Ts;F=fft(xt1,8192)/N;Xf1=2*abs(F);SS=Xf1.*conj(Xf1)/8192;Freq=Fs*(0:4096)/8192;figure(1),subplot(211),plot(Freq,xt1(1:4097),r),title(TIME Domain signal,color,r)subplot(212),plot(Freq,abs(SS(1:4097),r),title(Frequency Domain signal,color,r)title(Frequency Domain,color,r)ylabel(Voice FREQUENCY DOMAIN,color,r);xlabel(Frequency/Hz,color,r);grid on;zoom on;fhz1=0,400,400,4000,4000,11025;m=0,0,1,1,0,0;fhz0=fhz1./(Fs/2);b=fir2(100,fhz0,m,chebwin(101,20);tool=0:11025;tool2=freqz(b,1,tool,Fs);figure(2),plot(tool,abs(tool2)xt2=filter(b,1,xt1);F1=fft(xt2,8192)/N;Xf2=2*abs(F1);SS2=Xf2.*conj(Xf2)/8192;soundsc(xt2,Fs);figure(3),subplot(211),plot(Freq,xt2(1:4097),r),title(TIME Domain signal,color,r)subplot(212),plot(Freq,abs(SS2(1:4097),r),title(Frequency Domain signal,color,r)title(Frequency Domain,color,r)ylabel(Voice FREQUENCY DOMAIN,color,r);xlabel(Frequency/Hz,color,r);wavwrite(xt2,Fs,zz_cheb.wav);grid on;zoom on;pause求救MATLAB实现语音信号量化15标签：matlab,语音信号 MATLAB实现语音信号量化，程序，完整版跪谢！ 问题补充2009-12-13 19:09语音信号量化-月见晚樱 回答:1 人气:51 解决时间:2009-12-18 18:31 满意答案好评率：100% 以前做个这个二维矢量量化的给你一个%文件名VQ.m%采用K法算法%码,失真=VQ(原始信源,维数,码率,理想失真)function codebook,distortion=VQ(training_seq,dimension,codebook_size,tolerance)if (nargin=3)tolerance=0.001;endm=round(length(training_seq)/dimension);if (m*dimension-length(training_seq)0updated_codebook(l,:)=partition*training_matrix/sum(partition);endendnewdistortion=newdistortion/m;%继续训练矢量，直到失真没有显著变化while(abs(distortion-newdistortion)/newdistortiontolerance)distortion=newdistortion;newdistortion=0;for i=1:m;training_seq_block=training_seq(i-1)*dimension+1:i*dimension);training_matrix(i,:)=training_seq_block;distortion_block= ;for j=1:codebook_size;distort=sum(training_seq_block-updated_codebook(j,:).2);distortion_block=distortion_block distort;enddistortion_min,ind=min(distortion_block);newdistortion=newdistortion+distortion_min;index(i)=ind;endfor l=1:codebook_size;partition=(index=l);if sum(partition)0updated_codebook(l,:)=partition*training_matrix/sum(partition);endendnewdistortion=newdistor