音频降噪Matlab仿真.docxVIP

下载本文档

19
0
约1.26千字
约 5页
2017-01-31 发布于重庆
举报
版权申诉

音频降噪Matlab仿真.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

音频降噪Matlab仿真

数字信号处理大作业班级： 1401012_姓名：齐翔奡 _ 学号： 14010120082 输入信号的时域波形及其功率谱密度：叠加噪声后的音频信号的时域图形及功率谱密度：经过带通滤波器的音频信号的时域和功率谱密度：程序解读：clc;clear all;close all; [wav,fs]=audioread(GDGvoice8000.wav);t_end=1/fs *length(wav); % 计算声音的时间长度Fs=50000; % 仿真系统采样率t=1/Fs:1/Fs:t_end; % 仿真系统采样时间点% 利用插值函数将音频信号的采样率提升为Fs=50KHz wav=interp1([1/fs:1/fs:t_end],wav,t,spline);% 设计300Hz～3400Hz的带通预滤波器H(z) [fenzi,fenmu]=butter(6,[300 3400]/(Fs/2));nt = wgn(1,length(t),0.1); % 噪声nt=nt/(max(abs(nt)));%归一化噪声wav_noise = wav + nt;% 对音频信号进行滤波 wav_after = filter(fenzi,fenmu,wav_noise);figure(1);subplot(2,1,1); plot(wav(53550:53750));title(语音信号时域波形);axis([0 200 -0.3 0.3]);subplot(2,1,2); psd(wav, 10000, Fs);title(语音信号功率谱密度);axis([0 25000 -20 10]);figure(2);subplot(2,1,1); plot(wav_noise(53550:53750));title(加噪声后的语音信号时域波形);axis([0 200 -0.3 0.3]);subplot(2,1,2); psd(wav_noise, 10000, Fs);title(加噪声后的语音信号功率谱密度);axis([0 25000 -20 10]);figure(3);subplot(2,1,1); plot(wav_after(53550:53750));title(滤波后的语音信号时域波形);axis([0 200 -0.3 0.3]);subplot(2,1,2); psd(wav_after, 10000, Fs);title(滤波后的语音信号功率谱密度);axis([0 25000 -20 10]);sound(wav/max(wav), Fs); % 播放解调音频sound(wav_after/max(wav_after), Fs); % 播放解调音频wavwrite(demod_out,Fs,SSBDemod_OUT.wav); % 保存输出信号