基于matlab的音频信号低通滤波实验_低通滤波处理例题-程序员宅基地

技术标签： matlab 音频低通滤波

基于matlab的音频信号低通滤波实验

一、实验任务

录制语音信号，查看其时域频谱图。
加入高频噪声信号，并查看其时域频谱图。
设计低通滤波器，将带有高频噪声的语音信号经过低通滤波器，将噪声滤除，查看其输出时域频谱图。

二、仿真软件

matlab2013a

三、实验原理

语音信号：这次大作业录制的语音信号是WAV格式，通过audioread（）函数读取了语音数字信号的采样频率和数据。数据是两列N行的矩阵，两列表示双通道立体声，我只取第一列进行处理。
高频噪声信号：
模拟信号噪声：xn(t)=Acos(2πft),f为噪声频率，设定8kHz。
语音信号的采样频率为fs，周期为T，则：
以语音信号采样频率对模拟噪声采样得xn(nT)=Acos(2πfnT)，xn(n)=Acos(2πfn/fs)，得数字噪声信号矩阵（n取值范围是1到语音信号长度N）。噪声信号矩阵与语音信号矩阵相加得到加噪信号矩阵。
低通滤波器设计：
这次采用IIR巴特沃斯低通滤波器。低通滤波器得技术指标为:通带边界频率fp=3000Hz，阻带边界频率Fs=4000Hz,通带最大衰减1db，阻带最小衰减60db。按照IIR设计出滤波器。因为人声在300Hz到3.4KHz,噪声得频率在8kHz，因此噪声信号被滤除，语音信号稍微缺失一点。
四、实验代码



clear

clc

[audio,fs]=audioread('test.wav');%声音读取

audio = audio(:,1); %双通道变单通道

n=length(audio);

 

 

T = 1/fs;%采样间隔

t = (0:n-1)*T;%时间轴

f = (0:n-1)/n*fs;%频率轴

 

 %快速傅里叶变换

audio_fft=fft(audio,n)*T; 

 

%加噪声

tt =(1:n);

noise=0.02*cos(8000*2*pi/fs*tt');%加噪声

s=audio+noise;

s_fft=fft(s,n)*T; 

 

%设计IIR低通滤波器

rp = 1;

rs=60;

Ft=fs;

Fp=2000;

Fs=4000;                                        

 

wp=2*pi*Fp/Ft;

ws = 2*pi*Fs/Ft ;   %求出待设计的模拟滤波器的边界频率

 

[N,wn]=buttord(wp,ws,rp,rs,'s');    %低通滤波器的阶数和截止频率

[b,a]=butter(N,wn,'s');             %S域频率响应的参数即：滤波器的传输函数

[bz,az]=bilinear(b,a,0.5);          %利用双线性变换实现频率响应S域到Z域的变换

 

figure(2);%低通滤波器特性

[h,w]=freqz(bz,az);

title('IIR低通滤波器');

plot(w*fs/(2*pi),abs(h));

grid;

 

 

%滤波

z=filter(bz,az,s);

z_fft=fft(z);     %滤波后的信号频谱

 

figure(1); 

 

%绘出原始音频时域波

subplot(2,3,1);

plot(t,audio);   

xlabel('时间/s');

ylabel('幅度');

title('初始信号波形');  

grid;

 

 %绘出原始音频频域频谱

subplot(2,3,4); 

audiof = abs(audio_fft);

plot(f(1:(n-1)/2),audiof(1:(n-1)/2));

title('初始信号频谱');

xlabel('频率/Hz');

ylabel('幅度');

grid;

 

%绘出加噪音频时域波

subplot(2,3,2)

plot(t,s);

title('加噪声后信号波形')

xlabel('时间/s');

ylabel('幅度');

grid;

 

 %绘出加噪音频频域频谱

subplot(2,3,5)

sf = abs(s_fft);

plot(f(1:(n-1)/2),sf(1:(n-1)/2));

xlabel('频率/Hz');

ylabel('幅度');

title('加噪声后信号信号频谱');

grid;

 

%绘出滤波音频时域波

subplot(2,3,3);

plot(t,z);

title('低通滤波后的信号波形');

xlabel('时间/s');

ylabel('幅度');

grid;

 

%绘出滤波音频频域波

subplot(2,3,6);

zf = abs(z_fft);

plot(f(1:(n-1)/2),zf(1:(n-1)/2));

title('低通滤波后信号的频谱');

xlabel('频率/Hz');

ylabel('幅度');

grid;

 

audio_final = [audio;s;z]；%原始语音，加噪语音，滤波语音的合成音频矩阵

sound(audio_final,fs); %播放语音

五、实验结果

滤波器幅频特性

在这里插入图片描述
各个时频域图

在这里插入图片描述

本文链接：https://blog.csdn.net/qq_43160569/article/details/86544542

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

FTP命令字和返回码_ftp 登录返回230-程序员宅基地

文章浏览阅读3.5k次，点赞2次，收藏13次。为了从FTP服务器下载文件，需要要实现一个简单的FTP客户端。FTP（文件传输协议）是 TCP/IP 协议组中的应用层协议。FTP协议使用字符串格式命令字，每条命令都是一行字符串，以“\r\n”结尾。客户端发送格式是：命令+空格+参数+"\r\n"的格式服务器返回格式是以：状态码+空格+提示字符串+"\r\n"的格式，代码只要解析状态码就可以了。读写文件需要登陆服务器，特殊用..._ftp 登录返回230

centos7安装rabbitmq3.6.5_centos7 安装rabbitmq3.6.5-程序员宅基地

文章浏览阅读648次。前提：systemctl stop firewalld 关闭防火墙关闭selinux查看getenforce临时关闭setenforce 0永久关闭sed-i'/SELINUX/s/enforcing/disabled/'/etc/selinux/configselinux的三种模式enforcing：强制模式，SELinux 运作中，且已经正确的开始限制..._centos7 安装rabbitmq3.6.5

idea导入android工程,idea怎样导入Android studio 项目？-程序员宅基地

文章浏览阅读5.8k次。满意答案s55f2avsx2017.09.05采纳率：46%等级：12已帮助：5646人新版Android Studio/IntelliJ IDEA可以直接导入eclipse项目,不再推荐使用eclipse导出gradle的方式2启动Android Studio/IntelliJ IDEA，选择 import project3选择eclipse 项目4选择 create project f..._android studio 项目导入idea 看不懂安卓项目

浅谈AI大模型技术：概念、发展和应用_ai大模型应用开发-程序员宅基地

文章浏览阅读860次，点赞2次，收藏6次。AI大模型技术已经在自然语言处理、计算机视觉、多模态交互等领域取得了显著的进展和成果，同时也引发了一系列新的挑战和问题，如数据质量、计算效率、知识可解释性、安全可靠性等。城市运维涉及到多个方面，如交通管理、环境监测、公共安全、社会治理等，它们需要处理和分析大量的多模态数据，如图像、视频、语音、文本等，并根据不同的场景和需求，提供合适的决策和响应。知识搜索有多种形式，如语义搜索、对话搜索、图像搜索、视频搜索等，它们可以根据用户的输入和意图，从海量的数据源中检索出最相关的信息，并以友好的方式呈现给用户。_ai大模型应用开发

非常详细的阻抗测试基础知识_阻抗实部和虚部-程序员宅基地

文章浏览阅读8.2k次，点赞12次，收藏121次。为什么要测量阻抗呢？阻抗能代表什么？阻抗测量的注意事项... ...很多人可能会带着一系列的问题来阅读本文。不管是数字电路工程师还是射频工程师，都在关注各类器件的阻抗，本文非常值得一读。全文13000多字，认真读完大概需要2小时。一、阻抗测试基本概念阻抗定义：阻抗是元器件或电路对周期的交流信号的总的反作用。AC 交流测试信号 (幅度和频率)。包括实部和虚部。图1 阻抗的定义阻抗是评测电路、元件以及制作元件材料的重要参数。那么什么是阻抗呢？让我们先来看一下阻抗的定义。首先阻抗是一个矢量。通常，阻抗是_阻抗实部和虚部

小学生python游戏编程arcade----基本知识1_arcade语言 like-程序员宅基地

文章浏览阅读955次。前面章节分享试用了pyzero,pygame但随着想增加更丰富的游戏内容，好多还要进行自己编写类，从今天开始解绍一个新的python游戏库arcade模块。通过此次的《连连看》游戏实现，让我对swing的相关知识有了进一步的了解，对java这门语言也有了比以前更深刻的认识。java的一些基本语法，比如数据类型、运算符、程序流程控制和数组等，理解更加透彻。java最核心的核心就是面向对象思想，对于这一个概念，终于悟到了一些。_arcade语言 like

随便推点

【增强版短视频去水印源码】去水印微信小程序+去水印软件源码_去水印机要增强版-程序员宅基地

文章浏览阅读1.1k次。源码简介与安装说明：2021增强版短视频去水印源码去水印微信小程序源码网站去水印软件源码安装环境（需要材料）：备案域名–服务器安装宝塔-安装 Nginx 或者 Apachephp5.6 以上-安装 sg11 插件小程序已自带解析接口，支持全网主流短视频平台，搭建好了就能用注：接口是公益的，那么多人用解析慢是肯定的，前段和后端源码已经打包，上传服务器之后在配置文件修改数据库密码。然后输入自己的域名，进入后台，创建小程序，输入自己的小程序配置即可安装说明：上传源码，修改data/_去水印机要增强版