版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、1、设计原理设计数字滤波器的任务就是寻求一个因果稳定的线性时不变系统,并使系统函数H(z)具有指定的频率特性。 2、设计内容:以Matlab实现语音信号的低通滤波器设计为例:(1)语音信号的采集 利用Windows下的录音机,录制一段自己的话音,时间在1s内。然后在Matlab软件平台下,利用函数wavread对语音信号进行采样,记住采样频率和采样点数。通过wavread函数的使用,我们很快理解了采样频率、采样位数等概念。这里我直接采用了一段现成的.wav格式的语音信号。(2)语音信号的频谱分析 首先画出语音信号
2、的时域波形;然后对语音号进行快速傅里叶变换,得到信号的频谱特性。 程序如下: clc;clear;close all; fs=32768; %语音信号采样频率为32768,即215 x1=wavread('qq.wav'); %读取语音信号的数据,赋给变量x1 sound(x1,32768); %播放语音信号 y1=fft(x1,1024); %对信号做1024点FFT变换 f=fs*(0:511)/1024; figure(1); plot(x1)
3、 %做原始语音信号的时域图形 title('原始语音信号');xlabel('time n');ylabel('fuzhi n'); figure(2);freqz(x1) %绘制原始语音信号的频率响应图 title('频率响应图') figure(3);subplot(2,1,1); plot(abs(y1(1:512) %做原始语音信号的FFT频谱图 title('原始语音信号FFT频谱'); subplot(2,1
4、,2); plot(f,abs(y1(1:512); title('原始语音信号频谱') xlabel('Hz');ylabel('fuzhi'); 波形如下: 注意:原始语言信号FFT频谱和原始语言信号频谱的区别是:前者是频率为1递增的频谱,而后者是以f=fs*(0:511)/1024递增;另外,后者是在“不小于原始信号的频率(采样定理)”上完全展开的频谱。(3)给原始的语音信号加上一个高频余弦噪声,频率为5kHz。画出加噪后的语音信号时域和频谱图,与原
5、始信号对比,可以很明显的看出区别。 程序如下:%给原始的语音信号加上一个高频余弦噪声,频率为5kHz。画出加噪后的语音信号时域和频谱图,与原始信号对比,可以很明显的看出区别。fs=32768;x1=wavread('qq.wav');f=fs*(0:511)/1024;t=0:1/32768:(size(x1)-1)/32768; %将所加噪声信号的点数调整到与原始信号相同Au=0.03;d=Au*cos(2*pi*5000*t)' %噪声为5kHz的余弦信号x2=x1+d;sound(x2,327
6、68); %播放加噪声后的语音信号y2=fft(x2,1024);figure(1)plot(t,x2)title('加噪后的信号');xlabel('time n');ylabel('fuzhi n');figure(2)subplot(2,1,1);plot(f,abs(y1(1:512);title('原始语音信号频谱');xlabel('Hz');ylabel('fuzhi');subplot(2,1,2);plot(f,abs(y2(1:512);tit
7、le('加噪后的信号频谱');xlabel('Hz');ylabel('fuzhi'); 波形如下: (4)双线性变换法设计Butterworth滤波器 程序如下: fs=32768;x1=wavread('qq.wav');t=0:1/32768:(size(x1)-1)/32768;Au=0.03;d=Au*cos(2*pi*10000*t)'x2=x1+d;wp=0.25*pi;ws=0.3*pi;Rp=1;Rs=15;Fs=32768;Ts=1/Fs;
8、wp1=2/Ts*tan(wp/2); %将模拟指标转换成数字指标ws1=2/Ts*tan(ws/2); N,Wn=buttord(wp1,ws1,Rp,Rs,'s');%选择滤波器的最小阶数Z,P,K=buttap(N); %创建butterwo
9、rth模拟滤波器Bap,Aap=zp2tf(Z,P,K);b,a=lp2lp(Bap,Aap,Wn); bz,az=bilinear(b,a,Fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换H,W=freqz(bz,az); %绘制频率响应曲线figure(1)plot(W*F
10、s/(2*pi),abs(H)gridxlabel('频率Hz')ylabel('频率响应幅度')title('Butterworth')f1=filter(bz,az,x2);figure(2)subplot(2,1,1)plot(t,x2) %画出滤波前的时域图title('滤波前的时域波形');subplot(2,1,2)plot(t,f1);
11、160; %画出滤波后的时域图title('滤波后的时域波形');sound(f1,32768); %播放滤波后的信号F0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x2,1024);subplot(2,1,1);plot(f,abs(y2(1:512);
12、 %画出滤波前的频谱图title('滤波前的频谱')xlabel('Hz');ylabel('fuzhi');subplot(2,1,2)F1=plot(f,abs(F0(1:512); %画出滤波后的频谱图title('滤波后的频谱')xlabel('Hz');ylabel('fuzhi'); 波形如下: (5)窗函数法设计滤波器:
13、160;程序如下:fs=32768;x1=wavread('qq.wav');t=0:1/32768:(size(x1)-1)/32768;Au=0.03;d=Au*cos(2*pi*10000*t)'x2=x1+d;wp=0.25*pi;ws=0.3*pi;wdelta=ws-wp;N=ceil(6.6*pi/wdelta); %取整wn=(0.2+0.3)*pi/2
14、;b=fir1(N,wn/pi,hamming(N+1); %选择窗函数,并归一化截止频率figure(1)freqz(b,1,512)f2=filter(bz,az,x2)figure(2)subplot(2,1,1)plot(t,x2)title('滤波前的时域波形');subplot(2,1,2)plot(t,f2);title('滤波后的时域波形');sound(f2,32768);
15、60; %播放滤波后的语音信号F0=fft(f2,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x2,1024);subplot(2,1,1);plot(f,abs(y2(1:512);title('滤波前的频谱')xlabel('Hz');ylabel('fuzhi');subplot(2,1,2)F2=plot(f,abs(F0(1:512);title('滤波后的频谱')xlabel('Hz');ylabel('fuzhi'); 波形如下: (6)回放语音信号在Matlab中,函数sound可以对声音进行回放。其调用格式:sound(x,fs,bits),x为要播放的音频信号,fs采样频率,bits采样位。在运行Matlab程序时,可以听到回放的声音,经过比较,明显感觉滤波前后的声音有变化。 3.结语遇到的问题:出现错误提示:? Error using => wavreadError using =>
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 供应链管理实战课程设计
- 湖北省十堰市2024年七年级上学期语文期中考试试卷【附答案】
- 《几分之几》(教学设计)-2023-2024学年三年级上册数学人教版
- 湘科版(2017秋) 一年级上册2.1、初识大自然 教案
- 油气资源长途运输合同
- 茶楼装修居间协议模板
- 家居装修监理链家居间服务
- 酒店建设土石方合同样本
- 油田泥浆处理运输合同
- 游泳馆装修合同正规样本
- 2024年四川能投兴文电力限公司招聘【重点基础提升】模拟试题(共500题)附带答案详解
- 2024年中国社会科学院民族文学研究所第二批专业技术人员招聘1人历年【重点基础提升】模拟试题(共500题)附带答案详解
- 新闻记者采编报导人员岗位从业资格考试题含答案
- 汽车维修技术创新案例分析报告
- 车辆维修服务方案投标(2篇)
- 【小说双文本】2024届高考名校模拟小说双文本阅读
- 2024中考语文名著导读通关训练-《朝花夕拾》(原卷版)
- 高中信息技术学考复习19信息系统的组成与功能训练含答案
- (高清版)JGJT 178-2009 补偿收缩混凝土应用技术规程
- 中医诊疗技术规范标准
- 25题战略规划岗位常见面试问题含HR问题考察点及参考回答
评论
0/150
提交评论