语音信号的录制和处理
题目要求:
1.录制一段自己的语音信号,并对录制的信号进行采样;
2.画出采样后的语音信号的时域波形和频谱图;
3.给定滤波器的性能指标,采用窗函数法和双线性变换法设计滤波器,
并划出滤波器的频域响应;
4.用该滤波器对采集的信号进行滤波,画出滤波后信号的时域波形和频谱,
并对滤波前后的信号进行对比,分析信号的变化;
5.回放语音信号;
6.设计一个信号处理系统界面。
报告要求:
根据自己录制的语音信号选择合适的滤波器,并按照题目要求完成报告
程序1:用MATLAB对原始语音信号进行分析,画出它的时域波形和频谱
%语音信号采样频率为22050
x1=wavread(‘f:/s1.wav’); %读取语音信号的数据,赋给变量x1
fs=22050;
sound(x1,22050); %播放语音信号
y1=fft(x1,1024); %对信号做1024点FFT变换
f=fs*(0:511)/1024;
figure(1)
plot(x1) %做原始语音信号的时域图形
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,2);
plot(f,abs(y1(1:512)));
title(‘原始语音信号频谱’)
xlabel(‘Hz’);
ylabel(‘fuzhi’);
程序2:给原始的语音信号加上一个高频余弦噪声,频率为5kHz。画出加噪后的语音信号时域和频谱图,与原始信号对比,可以很明显的看出区别。
x1=wavread(‘f:/s1.wav’);
fs=22050;
f=fs*(0:511)/1024;
t=0:1/22050:(size(x1)-1)/22050;
Au=0.03;
d=[Aucos(2pi5000t)]’;
x2=x1+d;
sound(x2,22050)%播放加噪声后的语音信号
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)));
title(‘加噪后的信号频谱’);
xlabel(‘Hz’);
ylabel(‘fuzhi’);
程序3:双线性变换法设计Butterworth滤波器
x1=wavread(‘f:/s1.wav’);
t=0:1/22050:(size(x1)-1)/22050;
Au=0.03;
d=[Aucos(2pi5000t)]’;
x2=x1+d;
wp=0.25pi;
ws=0.3pi;
Rp=1;
Rs=15;
Fs=22050;
Ts=1/Fs;
wp1=2/Tstan(wp/2); %将模拟指标转换成数字指标
ws1=2/Tstan(ws/2);
[N,Wn]=buttord(wp1,ws1,Rp,Rs,‘s’); %选择滤波器的最小阶数
[Z,P,K]=buttap(N); %创建butterworth模拟滤波器
[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(WFs/(2pi),abs(H))
grid
xlabel(‘频率/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); %画出滤波后的时域图
title(‘滤波后的时域波形’);
sound(f1,22050); %播放滤波后的信号
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))); %画出滤波前的频谱图
title(‘滤波前的频谱’)
xlabel(‘Hz’);
ylabel(‘fuzhi’);
subplot(2,1,2)
F1=plot(f,abs(F0(1:512))); %画出滤波后的频谱图
title(‘滤波后的频谱’)
xlabel(‘Hz’);
ylabel(‘fuzhi’);
程序4:窗函数法设计滤波器
x1=wavread(‘f:/s1.wav’);
t=0:1/22050:(size(x1)-1)/22050;
Au=0.03;
d=[Aucos(2pi5000t)]’;
x2=x1+d;
wp=0.25pi;
ws=0.3pi;
wdelta=ws-wp;
N=ceil(6.6*pi/wdelta); %取整
wn=(0.2+0.3)pi/2;
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,22050); %播放滤波后的语音信号
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’);
语音信号的录制和处理相关推荐
- 数字信号处理的MATLAB实践(一)语音信号的录制和读取
萌生写这个系列的想法的原因是电脑里的代码太多,但删掉又怪可惜的,我的老师说学习和思考需要留下痕迹,这样后来者能有迹可循,不至于在浩如烟海的资料中无所适从,我是赞同这个说法的. 这个系列主要是我在学习数 ...
- 用计算机的声音编辑工具录制一段语音信号,语音信号的频谱分析实验指导
信号与系统中语音信号频谱分析实验指导 1.实验目的 (1)掌握傅里叶变换的物理意义,深刻理解傅里叶变换的内涵: (2)了解MATLAB对声音信号的处理指令: (3)了解计算机存储信号的方式及语音信号的 ...
- 在Matlab平台上实现对语音信号的去噪研究和仿真
在Matlab平台上实现对语音信号的去噪研究和仿真 摘要 语音信号在数字信号处理中占有极其重要的地位,因此选择通过对语音信号的研究来巩固和掌握数字信号处理的基本能力十分具有代表性.对数字信号处理离不开 ...
- Matlab语音采集与读写程序,基于MATLAB的语音信号录制采集和分析的程序设计
理 论广 角 ● I 基于 MATLAB的语音信号录制采集和分析的程序设计 刘 晓炯 (西北民族大学电气工程学院 甘肃 兰州I 730030) [摘 要]语音信号处理技术是语音处理领域中新近发展起来的 ...
- 基于Matlab的声波信号处理,基于声卡和Matlab平台的语音信号增强处理系统
第29卷第6期 V01.29 No.6 企业技术开发 TECHNOLOGICAL DEVELOPMENT 0F ENTERPRISE 2010年3月 Mar.2010 基于声卡和Matlab平台的语 ...
- matlab语音信号分析实验报告,DSP实验二 语音信号分析与处理,南京理工大学紫金学院实验报告,信号与系统...
实验二语音信号分析与处理 学号姓名 注:1)此次实验作为<数字信号处理>课程实验成绩的重要依据,请同学们认真.独立完成,不得抄袭. 2)请在授课教师规定的时间内完成: 3)完成作业后,请以 ...
- matlab编写信号采集程序,MATLAB语音信号采集课程设计
MATLAB语音信号采集课程设计 数字信号处理 课程设计报告 课设题目:语音信号的采集与处理 学 院:信息科学与工程学院 专 业:电子信息工程 班 级:电子0903 姓 名: 学 号: 指导教师:张海 ...
- 语音识别(20201030)——麦克风阵列(采集语音信号)
目录 麦克风阵列的原理 麦克风阵列的功能 麦克风阵列解决噪声干扰 存在的干扰 算法 1.多通道降混响算法(室内) 2.多通道点干扰抑制算法(室外) 3.多通道目标说话人提取算法(声纹,不适用) 麦克风 ...
- 语音信号的基音参数提取
文章目录 一.目的及原理 1.短时自相关法 2.倒谱法 二.两种方法实现 1.短时自相关法 1.程序代码 2.运行结果 2.倒谱法 1.程序代码 2.仿真结果 三.分析对比 1.自相关法运行结果 2. ...
最新文章
- Spring 为啥默认把 bean 设计成单例的?
- 4 种大学生兼职,助你走向经济独立!
- 安装好MongoDB,但服务中没有MongoDB服务的解决办法
- CCKS 2019 | 百度 CTO 王海峰详解知识图谱与语义理解
- SDN精华问答 | SDN可以做什么?
- mybatis报错:Could not find result map java.lang.Integer
- python输出24进制时间_python-如何在给定时间原点的情况下将日期时间从十进制转换为“%y-%m-%d%H:%M:%S”?...
- iOS 评论APP撰写评论
- #《机器学习》_周志华(西瓜书)南瓜书_第4章 决策树
- 如何在iconfont.cn 下载ttf格式的图标
- 如何查看电脑ip和端口
- 55、nginx rewrite
- Windows10设置暗色主题
- jQuery动态绑定事件或者原生js动态绑定事件
- A类博士后年薪30万起!北邮人工智能学院郭军教授实验室招聘博士后研究人员...
- [python爬虫]爬取电影,电视剧
- Metal 框架之创建纹理及纹理采样
- 微信第一个小程序----------EGO微商小程序部署搭建
- Matlab App Designer 中表格的操作
- PDN connection