音乐信号处理1Word格式.docx
- 文档编号:20522659
- 上传时间:2023-01-23
- 格式:DOCX
- 页数:14
- 大小:99.37KB
音乐信号处理1Word格式.docx
《音乐信号处理1Word格式.docx》由会员分享,可在线阅读,更多相关《音乐信号处理1Word格式.docx(14页珍藏版)》请在冰豆网上搜索。
3、通过信息处理实践的课程设计,掌握设计信息处理系统的思维方法和基本开发过程。
2.设计内容和要求(包括原始数据、技术参数、条件、设计要求等):
一、学习Matlab软件
1、熟悉Matlab的运行环境;
2、学会并掌握m语言编程;
3、根据所设计系统的需要会合理设定需完成系统的各项参数;
4、根据所设计系统的需要,优化程序设计最优系统。
二、实践设计要求:
1、根据所选题目,设计实现系统的原理框图。
2、编写m语言程序,给出系统不同节点输出波形。
3、每人一组,写出设计报告。
三、参考题目
题目1:
音乐信号的滤波处理
(1)
1、在windows下录制自己的一段语音信号或选取一段MP3文件;
2、在语音信号上分别叠加均匀白噪声和高斯白噪声,使信噪比为(学号)dB;
3、对叠加噪声前后的信号进行频谱分析,确定降噪的滤波器指标;
4、根据滤波器指标利用双线性变换法设计滤波器,在Matlab平台下编写程序,并观察分析滤波器的幅频特性、相频特性和群延时,以及滤波前后信号时域特性和频域特性。
题目2:
音乐信号的滤波处理
(2)
2、在语音信号上分别叠加均匀白噪声和高斯白噪声,使信噪比为(学号)dB;
4、根据滤波器指标利用冲激响应不变法设计滤波器,在Matlab平台下编写程序,并观察分析滤波器的幅频特性、相频特性和群延时,以及滤波前后信号时域特性和频域特性。
题目3:
数字基带传输系统眼图分析仿真
1、根据数字基带传输系统原理,建立数字基带传输系统原理框图;
2、生成任意二进制信号作为数字基带信号的信源,设信道为高斯白噪声,观察并分析(学号)dB信噪比条件下系统的输出结果,并利用眼图分析方法给出眼图分析结果;
3、对数字基带信号进行M序列加扰,在接收端对其进行解扰,其余条件同2。
根据要求,在Matlab平台下编写程序。
题目4:
音乐信号的滤波处理(3)
4、根据滤波器指标利用窗函数法设计滤波器,在Matlab平台下编写程序,并观察分析滤波器的幅频特性、相频特性和群延时,以及滤波前后信号时域特性和频域特性。
题目5:
音乐信号的滤波处理(4)
4、根据滤波器指标利用频率抽样法设计滤波器,在Matlab平台下编写程序,并观察分析滤波器的幅频特性、相频特性和群延时,以及滤波前后信号时域特性和频域特性。
题目6:
平顶抽样与自然抽样仿真
1、生成任一高斯白噪声信号f(t);
2、将高斯白噪声信号的带宽限制在(学号)kHz,形成信号x(t);
3、根据平顶抽样和自然抽样的原理框图,对x(t)信号分别进行平顶抽样和自然抽样。
在Matlab平台下编写程序,观察并分析抽样前后信号的时域与频域结果;
4、对抽样后的信号进行恢复,比较分析抽样前与恢复后的信号。
注:
以上6个题目中学号为班级学号的后两位,如0705094123,学号为23
选题要求:
学号对6取余,余数为1,选择题目1;
余数为2,选择题目2;
余数为3,选择题目3;
余数为4,选择题目4;
余数为5,选择题目5;
余数为0,选择题目6。
3.设计工作任务及工作量的要求〔包括课程设计计算说明书(论文)、图纸、实物样品等〕:
每个同学独立完成自己的任务,每人写一份设计报告,在课程设计论文中写明自己设计的部分,给出设计的matlab的程序与仿真波形图。
4.主要参考文献:
1.南利平等.通信原理简明教程(第2版).北京:
清华大学出版社.2009
2.郑君里等.信号与系统(第二版).北京:
高等教育出版社.2000
3.matlab相关书籍
5.设计成果形式及要求:
设计的程序;
课程设计说明书。
6.工作计划及进度:
6月12日~6月13日相关资料查阅;
6月17日~6月18日相关软件的学习;
6月21日~6月23日方案设计,建立模型并实现;
6月24日书写课程设计说明书;
6月25日答辩、成绩考核。
系主任审查意见:
签字:
年月日
目录
一、设计目的……………………………………………………………
二、设计要求……………………………………………………………
三、设计原理……………………………………………………………
四、设计步骤……………………………………………………………
1、设计方案……………………………………………………………………
2、设计过程……………………………………………………………………
3、设计结果……………………………………………………………………
五、设计总结……………………………………………………………
六、参考资料……………………………………………………………
一、设计目的
二、设计要求
三、设计原理
双线性变换法原理:
为了克服脉冲响应不变法所存在的缺点,提出采用双线性变换法。
该法的基本思想是首先按给定的指标设计一个模拟滤波器,其次将这个模拟滤波器的系统函数H(s),通过适当的数学变换方法把无限宽的频带,变换成频带受限的系统函数H(s~)。
最后再将H(s~)进行常规z变换,求得数字滤波器的系统函数H(z)。
这样由于在数字化以前已经对频带进行了压缩,所以数字化以后的频响可以做到无混叠效应。
显然,这里寻找压缩频带而又能满足上述映射条件的变换式是个关键。
设将s平面映射到s~平面存在下列的关系式
(1)
式中s~=σ~+jω~,C为变换常数,在式的右边是以ω~表示的周期函数,其周期为2π/T。
如果考虑频率特性则分别以s=jω,s~=jω~代入式
(1)故得:
(2)
现以ω为纵坐标,Ω~=ω~T为横坐标,则模拟与数字频率变量间的关系将如图所示。
结合式
(1)不难看出;
s的左半平面与s~的左半平面相对应;
s的右半平面与s~的右半平面相对应;
s平面的虚轴与s~平面的虚轴相对应。
它们之间主要的区别在于s平面-∞<
ω<
∞的无限频率范围,被映射到s~平面的主值范围内,即:
也就是说,通过变换式
(1)把整个
s平面映射到s~平面以+-ωs/2为边界的水平窄区内。
这避免数字化后可能出现的频谱混叠提供了必要条件。
为了求出数字滤波器的系统函数,最后还得通过常规z变换将s~平面变换到z平面上来,其关系式为:
(3)
显见,这时在s~左半平面的窄区就被映射到z平面的单位圆内。
现将式(3代入式(1最后求得:
或
(4)
该式是两个线性函数之比,称为线性分式变换,若把它展开求z,则得:
(5)
可见,其反变换也是线性分式函数,所以这种变换是双向的,因此叫做双线性变换。
双线性变换仍然具有将s的左半平面映射到z平面单位圆内;
jω轴映射到单位圆上的基本性质。
因为当σ=0时,│z│=1,说明s平面jω轴映射到z平面单位圆上。
当σ<
0时,上式中的分母大于分子│z│<
1,说明s左半平面映射到z平面单位圆内。
因而它们一一对应,有着单值关系,是一种保角交换。
因此一个稳定的模拟滤波器,通过双线性变换只能得到一个(唯一)稳定的数字滤波器。
它不象脉冲响应不变法那样,由于在z平面与s平面之间的映射存在着多值关系,以至在s平面上许多不同的ω值,映射到z平面后都重叠在一个点子上,造成频谱混叠。
由此可见,双线性变换法将;
频带严格限制在ωs/2范围内,从根本上消除了频谱混叠。
四、设计步骤
(一)、设计方案
2、在语音信号上分别叠加均匀白噪声和高斯白噪声,使信噪比为13dB;
(二)、设计过程
所选语音信号为aplacenearby.wav
1、读取声音信号并画出频谱图
[y,fs,bits]=wavread('
aplacenearby.wav'
);
%读取声音信号
sound(y,fs,bits);
%播放
N=length(y);
%求出声音信号长度
M=length(fs);
Y=fft(y,N);
%进行傅里叶变换
figure
(1);
subplot(311);
plot(y)
title('
声音信号的波形'
gridon
subplot(312);
plot(abs(Y));
%axis([-5000,185000,-1.2,1.2]);
声音信号的幅频响应曲线'
xlabel('
\omega/\pi'
ylabel('
|H(e^j^\omega)|'
subplot(313)
plot(angle(Y));
%axis([-200,4500,-10,600]);
声音信号的相频响应曲线'
\phi(\omega)'
2、在音频信号上加入均匀白噪声和高斯白噪声,画出加入白噪声后的频谱图
I=1;
%噪声强度
y1=y+1*rand(M,1);
%加均匀白噪声
sound(y1,fs,bits);
y2=awgn(y,13);
%加高斯白噪声,信噪比为学号13
sound(y2,fs,bits);
Y1=fft(y1,N);
Y2=fft(y2,N);
figure
(2);
plot(y1)
加均匀白噪声后的声音信号波形'
plot(abs(Y1));
y1声音信号的幅频响应曲线'
plot(angle(Y1));
y1声音信号的相频响应曲线'
plot(y2)
加高斯白噪声后的声音信号波形'
plot(abs(Y2));
y2声音信号的幅频响应曲线'
plot(angle(Y2));
y2声音信号的相频响应曲线'
3、对比信号的频谱
figure(3);
原始声音信号的幅频响应曲线'
加均匀白噪声后的幅频响应曲线'
subplot(313);
加高斯白噪声后的幅频响应曲线'
4、确定滤波器的设计参数并对加噪声信号去噪处理
Ft=8000;
Fp=1000;
Fs=1200;
T=1;
Ap=1;
As=50;
%衰减系数
wp=(2/T)*pi*Fp/Ft;
ws=(2/T)*pi*Fs/Ft;
fp=2*Ft*tan(wp/2);
fs=2*Fs*tan(wp/2);
[bLP,aLP]=buttord(wp,ws,Ap,As,'
s'
%求带通滤波器的阶数和截止频率
[bBP,aBP]=butter(bLP,aLP,'
%求s域的频率响应参数
[b,a]=bilinear(bBP,aBP,0.5);
%利用双线性变换法实现s域到z域的变换
figure(5);
subplot(211);
plot(abs(b));
%axis([0,7,-0.2,1.2]);
gridon
xlabel('
w(rad)'
|H(jw)|'
title('
Æ
µ
×
º
¯
Ê
ý
'
subplot(212);
plot(20*log10(abs(b)));
%axis([0,2*pi,-120,20]);
20*lg|H(jw)|(db)'
20*lg|H(jw)|--w'
yy1=filter(b,a,y1);
%滤波函数
yy2=filter(b,a,y2);
%%%%%%%%%%%%%%%%
YY1=fft(yy1);
%求滤波后的信号
YY2=fft(yy2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(4);
subplot(221);
plot(yy1);
yy1Â
Ë
²
¨
ó
Ð
Î
subplot(222);
plot(abs(YY1),'
r'
subplot(223);
plot(yy2);
yy2Â
subplot(224);
plot(abs(YY2),'
(三)、设计结果
五、设计总结
本次课程设计我是用双线性变换法设计滤波器,用matlab进行仿真,这次可程设计是数字信号处理和matlab软件的一次结合,让我更深的理解了理论和实践相结合的重要性,我将会继续努力,多多实践,理论和实践最好的联系起来,将自己的实力进行提高。
六、参考资料
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 音乐 信号 处理
![提示](https://static.bdocx.com/images/bang_tan.gif)