音频信号叠加噪声及滤波

一、前言
二、信号分析及加噪
三、滤波去噪
四、总结

很抱歉大家，最近经常有朋友私信问我关于这篇信号处理的一些问题，因为最近比较忙所以没能一一回复，给大家说句抱歉，希望那些给我私信的人可以看到。
大家问的问题大多是运行了我笔记上的代码但是结果不太好，可能是因为之前编辑错误所以笔记上的代码出现问题，在此向大家提供源码与音频，有需要的同学可以自取

链接如下：
https://download.csdn.net/download/zsisinterested/87222376

一、前言

之前一直对硬件上的内容比较关注，但是可能是因为硬件方面的东西可能真的是比较杂，而且需要渗透的东西太多了，所以学习进展比较缓慢。

因为也很少有单纯的硬件学习研究，总是会伴随着各种理论需要硬件做支撑，所以还是想要慢慢接触理论学习。但是之前总找不到切入点，不知道从哪里开始，就一直拖着。最近稍微接触了一点信号处理，就用这个当作切入点，开始接触理论学习。

二、信号分析及加噪

信号处理选用了matlab做工具，选了一个最简单的语音信号处理方式，用MATLAB先对原语音信号添加一个正弦噪声，在通过合适的滤波器去除噪声。

原语音信号需要自己录制一段或者找一段简短的音频，然后先对原语音信号进行分析。原语音信号尽量频率成分少一点，因为是初学，所以目的在于对信号处理有一个大概的了解，复杂的可以留待以后深入学习在做研究。

原语音信号做成.mp4或者.wav格式都可以，尽量做单声道的信号方便分析。如果录制时候是双声道，为了方便，就在处理时候选取单声道处理。

clc;clear
[x,Fs]=audioread('E:\test\tests.wav'); %读文件
N=length(x);
X=fft(x,N); %做傅里叶变换
k=0:N-1;
D=Fs/N;%计算频率分辨率
stem(k*D,abs(X),'Marker','none');
axis([0,4000,0,150])

用傅里叶变换可以看一下语音信号的频谱，由于人声的频率比较低，所以频率大多都集中在低频段，得到这个信息就得考虑一下需要如果有噪声的话，需要用什么样的滤波器。

下面对这个语音信号加上噪音，噪音采用正弦波 0.01* sin(2* pi* f) 函数，频率为1800HZ，加上噪声以后再对信号频谱观察

f=1800;
noise=0.01*sin(2*pi*f*k/Fs);
noise=noise'；
xa=x+noise;
Xa=fft(xa,N);
stem(k*D,abs(Xa),'Marker','none');
axis([0,4000,0,150])

对信号进行分析，得到频谱：

可以很明显的看到在1800HZ左右出现一个信号幅度，这个就是我们添加的噪音。但是真实的信号处理应该的有用信号的频带比较宽，而噪声的频带也会和有用信号更加的接近，或者说更复杂。

三、滤波去噪

对带有噪声的信号进行处理，一般选用合适的滤波器对信号进行过滤，根据频率滤除噪声。上面可以看到有用信号大部分都在1500HZ频段以下，而噪声在1800，所以就可以选用带阻滤波器对信号进行处理。
低通滤波器可以只过滤掉1800HZ的信号，留下其他信号。

先设计一个巴特沃斯带阻滤波器

Fs=8000; %采样频率
fp1=1500;fp2=2100;
fs1=1750;fs2=1850;
wp1=fp1/Fs*2*pi;   wp2=fp2/Fs*2*pi;   %通带截止频率
wp=[wp1,wp2];
ws1=fs1/Fs*2*pi;   ws2=fs2/Fs*2*pi;   %阻带截止频率
ws=[ws1,ws2];
Rp=1;As=30;
[n,wc]=buttord(wp/pi,ws/pi,Rp,As)
[b,a]=butter(n,wc,'stop') %求数字带阻滤波器系数
[H,w]=freqz(b,a);
dbH=20*log10(abs(H)/max(abs(H)));
plot(w/2/pi*Fs,dbH,'r');
axis([0,4000,-40,2]);

整个程序可以画出来滤波器的样子，但实际用到的时候，只需要用到里面的b、a参数，也就是滤波器系统函数的系数，只需要用filter( )函数设定滤波器系数就可以了。

可以看到滤波器的样子：

使用该滤波器对信号进行滤波，xa为叠加过噪声的信号，进行滤波以后，可以在通过傅里叶变换观察一下频谱

y=filter(b,a,xa);
Y=fft(y,N);
stem(k*D,abs(Y),'Marker','none');
axis([0,4000,0,150])

频谱如下：
可以明显看到1800HZ的频率被抑制

最后，通过对比频谱图观察一下整个滤波过程：

四、总结

很久没有写过博客了，之前总想着切入信号处理方面或者其他理论方面，结果总是找不到切入点，硬件方面因为太复杂进度又比较慢，所以这么长时间几乎没什么进展。只是趁着最近学习数字信号处理正好找到一个切入点，以后可能也会多多的注意这个方面。当然，硬件学习也会一直持续下去。

【注】个人学习笔记，请不吝赐教！

信号处理——MATLAB音频信号加噪、滤波相关推荐

MATLAB绘制加噪和去噪图像
wname = 'db5'; %设置小波名 [Lo_D,Hi_D,Lo_R,Hi_R] = wfilters(wname); % 计算与给定小波名相关的四个滤波器 subplot(221); stem ...
MATLAB绘制加噪和去噪图像2
load sinsin % 打开原始图像 init = 2055615866; randn('seed',init); x = X/18 + randn(size(X)); % 产生噪声图像 % 用w ...
音频信号调制，解调，加噪，去噪，滤波，matlab实现
题目分析调制解调通信系统实现. 录制一段语音信号,对其进行幅度调制解调,对比两端语音信号的频谱,并播放看其是否有偏差.在此基础上,将调制后的语音信号加上一个噪声,再滤波,观察恢复语音信号的失真现象. ...
【语音处理】基于matlab GUI低通滤波器语音信号加噪与去噪【含Matlab源码 1708期】
⛄一.语音处理简介(附课程作业报告) 1 语音信号的特点通过对大量语音信号的观察和分析发现,语音信号主要有下面两个特点: ①在频域内,语音信号的频谱分量主要集中在300-3400Hz的范围内.利用这 ...
【滤波器】基于低通滤波器语音信号加噪与去噪含Matlab源码
1 简介 1.1 课题的背景与意义通过语音传递信息是人类最重要.最有效.最常用和最方便的交换信息的形式.让计算机能听懂人类的语言,是人类自计算机诞生以来梦寐以求的想法,用现代手段研究语音信号,使人们 ...
基于matlab的音频信号低通滤波实验
基于matlab的音频信号低通滤波实验一.实验任务录制语音信号,查看其时域频谱图. 加入高频噪声信号,并查看其时域频谱图. 设计低通滤波器,将带有高频噪声的语音信号经过低通滤波器,将噪声滤除,查看 ...
matlab语音去除白噪声_基于的MATLAB的语音加噪去噪处理
龙源期刊网 http://www.qikan.com.cn 基于的 MATLAB 的语音加噪去噪处理作者:张大林何威李瑶瑶来源:<中国科技博览> 2019 年第 01 期 [ 摘 ...
【滤波器】基于汉宁窗FIR滤波器实现语音信号加噪去噪含Matlab源码
1 简介结合数字滤波器的理论基础和设计方法,在MATLAB程序语言环境下,设计出有限长单位脉冲响应(FIR)数字滤波器,同时利用GUI界面设计FIR数字滤波器人机交互平台,该系统平台界面直观.操作简 ...
干货分享基于MATLAB的带噪图像的高斯滤波论文
基于MATLAB的带噪图像的高斯滤波摘要:图像常常被强度随机信号(也称为噪声)所污染．一些常见的噪声有椒盐(Salt & Pepper)噪声.脉冲噪声.高斯噪声等．椒盐噪声含有随机出现的黑白 ...

信号处理——MATLAB音频信号加噪、滤波