音频降噪matlab仿真.pdf

下载文档 降价啦

39
0
约1.61千字
约 6页
2020-10-09 发布于江苏
举报
版权申诉
保障服务

音频降噪matlab仿真.pdf

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

数字信号处理大作业班级： 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); % 保存输出信号