《MATLAB课程设计》报告DOCWord下载.docx
- 文档编号:16543916
- 上传时间:2022-11-24
- 格式:DOCX
- 页数:22
- 大小:248.16KB
《MATLAB课程设计》报告DOCWord下载.docx
《《MATLAB课程设计》报告DOCWord下载.docx》由会员分享,可在线阅读,更多相关《《MATLAB课程设计》报告DOCWord下载.docx(22页珍藏版)》请在冰豆网上搜索。
例如解微分方程、傅里叶正反变换、拉普拉斯正反变换和z正反变换等。
MATLA在信号与系统中的另一主要应用是数值计算与仿真分析,主要包括函数波形绘制、函数运算、冲击响
应与阶跃响应仿真分析、信号的时域分析、信号的频谱分析、系统的S域分析和
零极点图绘制等内容。
数值计算仿真分析可以帮助学生更深入地理解理论知识,并为将来使用MATLABS行信号处理领域的各种分析和实际应用打下基础。
2.本题目的意义
本次课程设计的课题为《基于MATLA的语音信号采集与处理》,学会运用MATLA的信号处理功能,采集语音信号,并对语音信号进行滤波及变换处
理,观察其时域和频域特性,加深对信号处理理论的理解,并为今后熟练使用MATLA进行系统的分析仿真和设计奠定基础。
此次实习课程主要是为了进一步熟悉对matlab软件的使用,以及学会利用matlab对声音信号这种实际问题进行处理,将理论应用于实际,加深对它的理解
、设计任务及指标
利用MATLA对语音信号进行分析和处理,采集语音信号后,利用
软件平台进行频谱分析;
并对所采集的语音信号加入干扰噪声,对加入噪声的信号进行频谱分析,设计合适的滤波器滤除噪声,恢复原信号。
语音信号的“短时谱”对于非平稳信号,它是非周期的,频谱随时间连续变化,因此由傅里叶变换得到的频谱无法获知其在各个时刻的频谱特性。
如果利用加窗的方法从语音流中取出其中一个短断,再进行傅里叶变换,就可以得到该语音的短时
三、设计过程
(1).语音信号的采集及分析
基于声卡进行数字信号的采集。
将话筒插入计算机的语音输入插口上录音机。
按下录音按钮,对话筒说话,说完后停止录音。
要保存文件时
算机上的A/D转换器,把模拟的声音信号变成了离散的量化了的数字信号,放音时,
它又通过D/A转换器,把保存的数字数据恢复为原来的模拟的声音信号。
在
Matlab软件平台下可以利用函数wavread对语音信号进行采样,得到了声音数据
变量x1,同时把x1的采样频率fs=22050Hz和数据位Nbits=16Bit放进了MATALB
的工作空间。
图figure1为原始语音信号的时域图形。
从图中可以看出在时域
环境下,信号呈现出6不规则的信号峰值。
通过freqz函数绘制原始语音信号的频率响应图figure2。
然后对语音信号进行频谱分析,在Matlab中可以利用函数fft对信号行快速傅里叶变换,得到信号的频谱图figure3,从图中可以看出
对各个频点上的随机信号在频域进行抽样,抽样频率为22050Hz。
(2).给原始信号加上一个高频噪声
在Matlab中人为设计一个固定频率5500Hz的噪声干扰信号。
噪声信号通常
为随机序列,在本设计中用正弦序列代替,干扰信号构建命令函数为
d=[Au*sin(2*pi*5500*t)]'
给出的干扰信号为一个正弦信号,针对上面的语音
信号,采集了其中一段。
再对噪音信号进行频谱变换得到其频谱图,从图中可以
看出干扰信号,在4000Hz和6000Hz频点处有一高峰,其中5500Hz正是本设计所要利用的。
(3).设计一个滤波器,滤除高频噪声
由模拟滤波器变换为数字滤波器时,采用的是双线性变换法,它保留的是从模拟到数字域的系统函数表示。
在滤波器的选取时,由于设计方法的侧重点不同作出比较是困难的。
如果FIR滤波器情况下,最优的设计将是椭圆滤波器。
用双线性变换法设计低通滤波器。
的本设计是用双线性变换法设计BW带阻型滤波器。
在MATLAB中,可以利用函数butterworth,设计FIR滤波器,利用MATLAB中的函数freqz画出各滤波器的频率响应。
用设计好的带阻滤波器对含噪语音信号进行滤波,在Matlab中FIR滤波器利用函数fftfilt对信号进行滤波,FIR滤波器
利用函数filter对信号进行滤波。
在一个窗口同时画出滤波前后的波形及频谱。
从图中可以看出,5500Hz看到的高峰消失了,语音信号与开始的一样,滤波器成功的滤除了干扰信号。
利用MATLAB寸语音信号进行分析和处理,采集语音信号后,利用MATLAB软件平台进行频谱分析;
并对所采集的语音信号加入干扰噪声,对加入噪声的信号进行频谱分析,设计合适的滤波器滤除噪声,恢复原信号。
四、结论及分析
1•用MATLA对原始语音信号进行分析,画出它的时域波形和频谱
程序:
fs=22050;
%语音信号采样频率为22050
x1=wavread('
C:
\Documentsand桌面声音文
件.wav'
);
%读取语音信号的数据,赋给变量x1,
sound(x1,22050);
%
y1=fft(x1,1024);
f=fs*(0:
511)/1024;
figure⑴%
plot(x1)%title('
原始语音信号'
xlabel('
timen'
%xylabel('
fuzhin'
%yfigure
(2)
freqz(xl)%
title('
频率响应图'
)figure(3)
subplot(2,1,1);
plot(abs(y1(1:
512)))%
原始语音信号FFT频谱'
)
subplot(2,1,2);
plot(f,abs(y1(1:
512)));
%abs线
原始语音信号频谱'
xlabel('
Hz'
ylabel('
幅值'
播放语音信号
对信号做1024点FFT变换
将0到511,步长为1的序列的值与fs相乘并
除以1024的值,赋值给f
创建图形窗1
做原始语音信号的时域图形
轴的名字是“timen”
轴的名字是“fuzhin”
绘制原始语音信号的频率响应图
创建两行一列绘图区间的第1个绘图区间
做原始语音信号的FFT频谱图
是绝对值,plot是直角坐标下线性刻度曲
时域波形和频谱:
Fignre1
图1原始语音信号
图2语音信号频率响应图
图3原始语音信号FFT与信号频谱
2.给原始的语音信号加上一个高频余弦噪声,频率为5500hz。
对加噪后的语音
进行分析,并画出其信号时域和频谱图。
桌面声音文
读取语音信号的数据,赋给变量x1
将0到511,步长为1的序列的值与fs相乘
并除以1024的值,赋值给f
将0到x1的长度减1后的值除以fs的值,
且步长为1/fs的值,的序列的值,赋予t
\Documentsand
%
t=0:
1/fs:
(length(x1)-1)/fs;
Au=0.03;
d=[Au*sin(2*pi*5500*t)]'
;
x2=x1+d;
sound(x2,22050);
y2=fft(x2,1024);
噪声幅值
所加的噪声是正弦信号
将正弦信号噪声加在语音信号上
对信号y1做1024点FFT变换
对信号y2做1024点FFT变换
figure
(1);
plot(t,x2);
线性图)
加噪后的信号'
%x
%y
figure
(2)%
fuzhi'
%y
plot(f,abs(y2(1:
加噪后的信号频谱'
做加噪后的信号时域图形(绘出t对x2的
创建图形窗2
做原始语音信号的频谱图
轴的名字是“fuzhi”
创建两行一列绘图区间的第2个绘图区间做加噪后的语音信号的频谱图(abs是绝对值,plot是直角坐标下线性刻度曲线)
与原始信号对比,区别:
图4加噪后的信号时域
图
图5原始信号与加噪后信号频谱对比
0到2000有幅值,在4000到6000之间没有幅值,
4000到6000之间出现最大幅值12,超出正常值。
如图5
绘出滤波后的信号频域和时域波形:
先原始信号没加噪音之前
但是在加了噪音之后
原始信号与加噪后信号频谱对比图所示
3.设计合适的滤波器,滤除高频噪声,程序:
%用双线性变换法设计巴特沃思低通滤波器
wp=0.25*pi;
ws=0.3*pi;
Rp=1;
Rs=15;
Fs=22050;
Ts=1/Fs;
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);
[Bap,Aap]=zp2tf(Z,P,K);
[b,a]=lp2lp(Bap,Aap,Wn);
止频率
[bz,az]=bilinear(b,a,Fs);
[H,W]=freqz(bz,az);
figure(9)
plot(W*Fs/(2*pi),abs(H))%
应曲线
grid%
频率/Hz'
频率响应幅度'
通带截止频率
阻带截止频率通带最大衰减(db)阻带最大衰减(db)
将模拟指标转换成数字指标
选择滤波器的最小阶数(估算得到
Butterworth低通滤波器的最小阶数
N和3dB截止频率WC
创建Butterworth低通滤波器原型
将零极点增益转换为普遍分子,分母
将普遍的分子和分母转换为以Wr为截
用双线性变换法实现模拟滤波器到数
字滤波器的转换(模拟转换为数字)
求频率响应
绘制Butterworth低通滤波器频率响
添加图格
Butterworth'
f仁filter(bz,az,z);
%利用上面已做好的滤波器(已bz和
az为特征的滤波器)对x2信号进行滤波,赋值给fl
figure(10)
subplot(2,1,1)
创建两行一列绘图区间的第
1绘图区
间
plot(t,z)
做加噪后的信号时域图形
滤波前的时域波形
'
subplot(2,1,2)
2绘图区
plot(t,f1);
滤波后的时域波形
sound(f1,22050);
p=length(f1);
F0=fft(f1,p);
f=0:
fs/p:
fs*(s-1)/p;
figure(11)
y2=fft(z,1024);
plot(f,abs(y4));
滤波前的频谱
%x
轴的名字是“Hz”
%y
轴的名子是fuzhi
F1=plot(f,abs(F0));
滤波后的频谱
图6低通滤
波器频率响应曲线
时域波形和频谱:
*At'
ifurp1II
匚忙IX
lil.a―
LkiiinJtIwL>
iLelk^v^V.JftduWH-slf?
□Q昌
$1題Qo
◎遲LI目aO
滤液信E勺吋域泱笊
0.4I111
U上-J-
d~™
njIiiiii*iii
■"
nD.10304C.506170.6191
图7滤波前后时域波形图
图8滤波前后频谱图
与原噪音信号对比,区别:
通过对比分析可知,滤波后的输出波形和原始语音加噪声信号的图形发生了一些变化,在3000-6000HZ之间的信号消失,出现0-1000和7000-8000之间的信号。
滤波后的输出波形明显在对应时间幅度比原语音加噪声信号的要小,而且滤波的效果也与滤波器的选择有关,可以看出滤波器的性能差异以及参数的调节会对滤波器产生一定的影响,并且通过回放可以发现滤波前后的声音有变化•低通滤波后,已
很接近原来的声音,人耳几乎辨别不出。
从频谱图中我们还可以看
出声音的能量信号主要集中在低频部分,说明高频语音信号被滤出,滤波器达到既定要求。
五、设计体会
正所谓“纸上得来终觉浅,觉知此事要躬行。
”学习任何知识,仅从理论上去求知,而不去实践、探索是不够的。
因此在学期末来临之际,我们迎来了MATLAB课程设计。
通过为期一周半的MATLA课程设计,我对MATLA这个仿真软件有了更进一步的认识和了解。
在这一周半时间里,我通过自己摸索,查阅资料,并且在老师的指导下完成了:
语音信号的采集及分析;
给原始信号加上一个高频噪声;
设计一个滤波器,滤除高频噪声;
并最终将课程设计报告总结完毕。
在整个设计过程中我懂得了许多东西,也培养了独立思考和设计的能力,树立了对知识应用的信心,相信会对今后的学习工作和生活有非常大的帮助,并且提高了自己的动手实践操作能力,使自己充分体会到了在设计过程中的成功喜
悦。
虽然这个设计做的不怎么好,但是在设计过程中所学到的东西是这次课程设计的最大收获和财富,使我终身受益。
在没有做课程设计以前,觉得课程设计只是对知识的单纯总结,但是通过这次课程设计发现自己的看法有点太片面,课程设计不仅是对前面所学知识的一种
检验,也是对自己能力的一种提高,通过这次课程设计使自己明白了原来的那点知识是非常欠缺的,要学习的东西还很多,通过这次课程设计,明白学习是一个长期积累的过程,在以后的工作和生活中都应该不断的学习,努力提高自己的知识和综合素质。
希望以后像这样的课程设计在多一点。
参考文献:
[1]姚东.MATLAB及在电子信息课程中的应用[M].第二版.北京:
北京电子工业出版社,2000
[2]张智星.MATLAB程序设计与应用•北京:
清华大学出版社,2002
[3]陈怀琛.MATLAB及其在理工课程中的应用指南.西安:
西安电子科技大学出版
社,1999.10
[4]肖伟、刘忠.MATLAB程序设计与应用[M].北京:
清华大学出版社2005
程序清单:
口
■**——******}'
2
\Documents
and桌面
声音文
%读取语音信号的数据,赋给变量
x1,
y仁fft(x1,1024);
将0到511,步长为1的序列的值与
fs相乘并
plot(x1)%
原始语音信号
轴的名子是timen
轴的名子是fuzhin
figure
(2)
figure(3)
512)))
FFT频谱'
%abs
线
ylabel('
%读取语音信号的数据,赋给变量x1
将0到x1的长度减1后的值除以fs的值,
sound(x2,22050);
加噪后的信号
频谱'
%y
创建两行一列绘图区间的第2个绘图区间
做加噪后的语音信号的频谱图(abs是绝对
值,plot是直角坐标下线性刻度曲线)
轴的名字是“
timen”
轴的名字是
“fuzhi”
[Bap,Aap]=zp2tf(Z,P,K);
截止频率
[H,W]=freqz(bz,az);
plot(W*Fs/(2*pi),abs(H))
grid
f1=filter(bz,az,z);
%通带截止频率
%阻带截止频率
%通带最大衰减(db)
%阻带最大衰减(db)
%将模拟指标转换成数字指标
%选择滤波器的最小阶数(估算得到
Butterworth低通滤波器的最小阶数
N和3dB截止频率Wc)
%创建Butterworth低通滤波器原型
%将零极点增益转换为普遍分子,分母
%将普遍的分子和分母转换为以Wn为
%用双线性变换法实现模拟滤波器到数字滤波器的转换(模拟转换为数字)
%求频率响应
%绘制Butterworth低通滤波器频率响
%添加图格
利用上面已做好的滤波器(已bz和az为特征的滤波器)对x2信号进行滤波,赋值给fl
subplot(2,1,1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB课程设计 MATLAB 课程设计 报告 DOC