谱减法的Matlab程序.doc_第1页
谱减法的Matlab程序.doc_第2页
谱减法的Matlab程序.doc_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

精品文档谱减法的Matlab程序:-clear;xx,fs=wavread(I:aa0a.wav); % 读语音信号team,row=size(xx); if row=2 x=(xx(:,1)+xx(:,2)/2; % 将多通道求平均值 yy=x;else x=xx; yy=x;endx=x-mean(x)+0.1*rand(length(x),1); % 去直流分量并加噪N=length(x);n=220; % 每一帧长n1=160; %每两帧重合的长度frame=floor(N-n)/(n-n1); %帧数for i=1:frame y1=x(i-1)*(n-n1)+1:(i-1)*(n-n1)+n).*hamming(n); % 每一帧对应的语音信号 fy=fft(y1,n); %求fft nen(i,:)=abs(fy).2; %求能量 ang(i,:)=angle(fy); % 求角度endyuzhi=sum(sum(nen(2:5,:)/(4*n); % 求阈值，信号开始阶段的能量平均值for i=1:frame nen(i,:)=nen(i,:)-yuzhi; % 减谱 nen(i,find(nen(i,:)0)=0; % 将小于0的部分赋值为0endfor i=1:frame nen(i,:)=sqrt(nen(i,:); jie=nen(i,:).*exp(j*ang(i,:); %求频域函数 out(i,:)=real(ifft(jie)./hamming(n); %求逆fftendzong=out(1,:); % 对out求导，以便于原信号比较jiewei=n;for i=2:frame %将一系列的帧组合还原 zong(jiewei-n1+1:jiewei)=(zong(jiewei-n1+1:jiewei)+out(i,1:n1)/2; jiewei=jiewei+n-n1; zong=zong;out(i,n1+1:end);endfigure(1); subplot(211);plot(x); %画加噪的原始语音信号axis(1,(n-n1)*frame+n,min(x),max(x);subplot(212);specgram(x,fs,1024,n,n1); % 对应的鱼谱图figure(2);subplot(211);plot(zong); %画增强的语音信号axis(1,(n-n1)*frame+n,min(zong),max(zong);subplot(212);specgram(zong,fs,1024,n,n1);wavplay(x,fs); % 输出音频wavplay(zong,fs
人人文库> 全部分类> 应用文书 > 事务文书

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论