信号课程设计 伯晓乐 029110028.docx
- 文档编号:10568317
- 上传时间:2023-05-26
- 格式:DOCX
- 页数:17
- 大小:175.93KB
信号课程设计 伯晓乐 029110028.docx
《信号课程设计 伯晓乐 029110028.docx》由会员分享,可在线阅读,更多相关《信号课程设计 伯晓乐 029110028.docx(17页珍藏版)》请在冰点文库上搜索。
信号课程设计伯晓乐029110028
信号与系统课程设计(第二题)
姓名:
伯晓乐
学号:
2012029110028
教师:
孔斌
题目:
语音信号处理与识别
时间:
2014.05.23
第一部分理论自学内容阐述
1采样定理
时域采样定理
频带为F的连续信号f(t)可用一系列离散的采样值f(t1),f(t1±Δt),f(t1±2Δt),...来表示,只要这些采样点的时间间隔Δt≤1/(2F),便可根据各采样值完全恢复原来的信号f(t)。
这是时域采样定理的一种表述方式。
时域采样定理的另一种表述方式是:
当时间信号函数f(t)的最高频率分量为fM时,f(t)的值可由一系列采样间隔小于或等于1/(2fM)的采样值来确定,即采样点的重复频率f≥(2fM)。
图为模拟信号和采样样本的示意图。
时域采样定理是采样误差理论、随机变量采样理论和多变量采样理论的基础。
频域采样定理 对于时间上受限制的连续信号f(t)(即当│t│>T时,f(t)=0,这里T=T2-T1是信号的持续时间),若其频谱为F(ω),则可在频域上用一系列离散的采样值来表示,只要这些采样点的频率间隔ω≦π/tm。
2理想低通滤波器
对于理想低通滤波器而言,在物理上是不可能实现的,但是我们能设计出接近理想特性的滤波器。
虽然理想低通滤波器在实际中是不能实现的,但是我们希望找到一种区分可实现性与不可实现性的标准,这就是佩利-维纳(Paley-Wiener)准则。
由佩利-维纳准则给出了,一个网络幅度函数物理可实现的必要条件,但绝非充要条件。
物理可实现性在文献中有不同定义方法,这里采用最低限度的定义把物理可实现性系统和不可实现系统区分开来。
我们可以直观地看到,一个物理可实现系统在激励加入之前是不可能有响应输出的,这称为因果条件。
这个条件在时域里的表述为:
物理可实现系统的单位冲激响应必须是有起因的,即
。
从频域来看,如果幅度函数满足平方可积条件,即
,佩利和维纳证明了对于幅度函数物理可实现的必要条件是
,它被称为佩利-维纳准则。
第二部分案例或课题陈述及实现方案
【设计思路】
滤波器是一种用来消除干扰杂讯的器件,将输入或输出经过过滤而得到纯净的交流电。
采用MATLAB设计滤波器,使原来非常繁琐复杂的程序设计变成了简单的函数调用,为滤波器的设和实现开辟了广阔的天地,尤其是Matlab工具箱使各个领域的研究人员可以直观方便地进行科学研究与工程应用。
MATLAB 信号处理工具箱为滤波器设计及分析提供了非常优秀的辅助设计工具, 在设计数字滤波器时, 善于应用MATLAB进行辅助设计, 能够大大提高设计效率。
我们都知道,语音信号经过调制在大自然中传播,必然会夹杂着噪声一起传过来,那么我们该如何设计电路,从而达到对语音信号的提取,滤除噪音信号呢。
根据资料显示,人类的发生频率范围在300hz-3400hz
2.语音信号处理与识别
目的:
理解时域和频域尺度变换基本概念,掌握信号时频域分析方法,正确理解采样定理,准确理解滤波器的概念。
内容:
(1)使用Matlab中wavrecord命令录制一段3秒的语音信号,使用wavplay命令播放,录制命令和播放命令中的采样频率设置成相同和不同两种情况,对观察到的现象进行分析并结合课本中的知识对该现象进行解释;
(2)使用不同的采样频率录制一段3秒的语音信号,画出信号的时域波形和频谱;找到语音信号的主要频谱成分所在的带宽;观察并分析不同采样频率对波形和频谱的影响;寻找声音信号不出现明显失真的最低采样频率;
(3)录制一段男生的语音信号和一段女生的语音信号,对两段音频信号进行混合,设计滤波器将混合的语音信号分开成单独的男声和女声信号,如果分离效果不好,对原因进行解释。
答:
(1)采样频率8000hz,播放频率分别为8000hz和6000hz和10000hz,可明显听出声音的不同,采样频率与播放频率相同的声音听起来正常,播放频率低的时候,声音变低沉;播放频率高的时候,声音变尖锐。
结合已学知识可略作解释,声音信号以较采样频率低的频率播放,频域被压缩,则时域上被拉伸,即听起来较“低沉”声音信号以较采样频率高的频率播放,频域被拉伸,则时域上被压缩,即听起来较“尖锐”。
(2)实验所用matlab程序代码以及信号的时域波形和频谱如下:
可以找到语音信号的主要频谱成分所在的带宽在600Hz——1300Hz,不同采样频率对波形和频谱的影响,幅值有所不同,采样频率较低的幅值较小,波形形状相同。
据图可以看出,声音信号不出现明显失真的最低采样频率大约在1800HZ。
(3)实验所用matlab程序代码以及信号的时域波形和频谱如下:
录制了一段男声信号和一段女生信号,混合后在利用低通和高通滤波器分别滤得较低频和较高频率的两断信号,播放两段信号,可以较清晰的听出单独的男声和女声,但似乎分离不彻底,信号中有所混杂,原因分析是男声和女声毕竟有共同频率,,无法彻底分离。
(1)fs=10000;
>>x=wavrecord(3*fs,fs);
>>wavplay(x,fs)
>>wavplay(x,6000)
>>wavplay(x,8000)
(2)>>Fs1=6000;
>>y=wavrecord(3*Fs1,Fs1);
>Inwavrecordat42
>>Fs=11025;
>>y1=wavrecord(3*Fs,Fs);
>Inwavrecordat42
>>wavplay(y1)
>Inwavplayat41
>>wavplay(y)
>Inwavplayat41
>>figure
(1)
>>subplot(2,1,1)
>>N=length(y);
>>N1=length(y1);
>>t=(0:
N-1)/Fs1;
>>t1=(0:
N1-1)/Fs;
>>w=Fs1*linspace(0,1,N)-Fs1/2;
>>w1=Fs*linspace(0,1,N1)-Fs/2;
>>Y=fft(y,N);
>>Y1=fft(y1,N1);
>>Y=fftshift(Y);
>>Y1=fftshift(Y1);
>>figure
(1)
>>figure
(1)
>>subplot(2,1,1)
>>plot(t,y);
>>title('6000采样频率时域图');
>>subplot(2,1,2)
>>plot(t1,y1);
>>title('11025采样频率时域图');
>>figure
(2)
>>figure
(2)
>>subplot(2,1,1)
>>plot(w,abs(Y))
>>title('6000采样频率频域图');
>>subplot(2,1,2)
>>plot(w1,abs(Y1))
>>title('11025采样频率频域图');
>>Fs2=1300;
>>y2=wavrecord(3*Fs2,Fs2);
>Inwavrecordat42
>>wavplay(y2,Fs2)
>Inwavplayat41
>>wavplay(y,Fs1)
>Inwavplayat41
>>N2=length(y2);
>>t2=(0:
N2-1)/Fs2;
>>Y2=fft(y2,N2);
>>Y2=fftshift(Y2);
>>figure(3)
>>subplot(2,1,1)
>>plot(t2,y2)
>>title(‘1300采样频率时域图’)
>>subplot(2,1,2)
>>w2=Fs2*linspace(0,1,N2)-Fs2/2;
>>plot(w2,abs(Y2))
>>subplot(2,1,1)
>>title(‘1300采样频率频域图’)
(3)>>y3=wavrecord(5*Fs3,Fs3);
Warning:
WAVRECORDwillberemovedinafuturerelease.UseAUDIORECORDERinstead.
>Inwavrecordat42
>>wavplay(y3)
Warning:
WAVPLAYwillberemovedinafuturerelease.UseAUDIOPLAYERinstead.
>Inwavplayat41
>>y5=y3+y4;
>>wavplay(y5)
Warning:
WAVPLAYwillberemovedinafuturerelease.UseAUDIOPLAYERinstead.
>Inwavplayat41
>>y4=y4/2;
>>y5=y3+y4;
>>wavplay(y5)
Warning:
WAVPLAYwillberemovedinafuturerelease.UseAUDIOPLAYERinstead.
>Inwavplayat41
>>y4=y4/1.5;
>>y5=y3+y4;
>>wavplay(y5)
Warning:
WAVPLAYwillberemovedinafuturerelease.UseAUDIOPLAYERinstead.
>Inwavplayat41
>>N5=length(y5);
>>t=(0:
N5-1)/11025;
>>figure(4)
>>figure(4)
>>subplot(2,1,1)
>>w5=11025*linspace(0,1,N5)-11025/2;
>>Y5=fft(y5,N5);
>>Y5=fftshift(Y5);
>>t5=(0:
N5-1)/11025;
>>t=(0:
N-1)/Fs1;
>>Fs5=11025;
>>plot(t5,Y5);
Warning:
ImaginarypartsofcomplexXand/orYargumentsignored
>>plot(t5,abs(Y5));
>>subplot(2,1,2)
>>plot(w5,Y5)
Warning:
ImaginarypartsofcomplexXand/orYargumentsignored
>>plot(w5,abs(Y5))
>>title('混合信号频域图');
title('混合信号频域图');
|
Error:
TheinputcharacterisnotvalidinMATLABstatementsorexpressions.
>>subplot(2,1,1)
>>plot(t5.y5);
Attempttoreferencefieldofnon-structurearray.
>>plot(t5,abs(y5));
>>title('混合信号时域')
>>subplot(2,1,2)
>>title('混合信号频域')
>>N3=length(y3);
>>N4=length(y4);
>>t3=(0:
N3-1)/Fs3;
>>t4=(0:
N4-1)/Fs4;
>>w3=Fs3*linspace(0,1,N3)-Fs3/2;
>>w4=Fs4*linspace(0,1,N4)-Fs4/2;
>>Y3=fft(y3,N3);Y4=fft(y4,N4);
>>Y3=fftshift(Y3);Y4=fftshift(Y4);
>>figure(5)
>>figure(5)
>>subplot(2,1,1)
>>plot(t3,abs(y3))
>>title('男生时域')
>>subplot(2,1,2)
>>plot(t4,abs(y4))
>>title('女生时域')
>>
figure(6)
>>figure(6)
>>subplot(2,1,1)
>>plot(w3,abs(Y3))
>>title('男生频域')
>>subplot(2,1,2)
>>plot(w4,abs(Y4))
>>title('女生频域')
>>
>>b=fir1(6,0.1179,'low');
>>filter(b,1,y5);
>>wavplay(y5)
Warning:
WAVPLAYwillberemovedinafuturerelease.UseAUDIOPLAYERinstead.
>Inwavplayat41
>>b=fir1(60,0.1179,'low');
>>freqz(b,1);
>>y6=filter(b,1,y5);
>>wavplay(y6)
Warning:
WAVPLAYwillberemovedinafuturerelease.UseAUDIOPLAYERinstead.
>Inwavplayat41
>>Y6=fft(y6,N5);
>>figure(7)
>>Y6=fftshift(Y6);
>>plot(w5,abs(Y6))
>>Y6=fftshift(Y6);
>>plot(w5,abs(Y6))
>>b=fir1(60,0.1179,'high');
>>freqz(b,1);
>>y6=filter(b,1,y5);
>>wavplay(y6)
Warning:
WAVPLAYwillberemovedinafuturerelease.UseAUDIOPLAYERinstead.
>Inwavplayat41
>>freqz(b,1);
>>Y6=fft(y6,N5);
>>Y6=fftshift(Y6);
>>plot(w5,abs(Y6))
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号课程设计 伯晓乐 029110028 信号 课程设计
![提示](https://static.bingdoc.com/images/bang_tan.gif)