的应用matlab对信号进行频谱分析报告及滤波(数字信号处理)

《的应用matlab对信号进行频谱分析报告及滤波(数字信号处理)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案存档资料成绩：　华东交通大学理工学院课程设计报告书所属课程名称数字信号处理题目应用Matlab对信号进行频谱分析及滤波分院　　电信分院　　　专业班级电子信息工程*班学　　号**************　　　　　　　　学生姓名邓群　　　　　　　　指导教师****　　　2013年6月27日精彩文档实用标准文案序号项目等　　　级优秀良好中等及格不及格1课程设计态度评价2出勤情况评价3任务难度评价4工作量饱满评价5任务难度评价6设计中创新性评价7论文书写规范化评价8综合应用能力评价综合评定等级

2、课程设计（论文）评阅意见评阅人**　职称讲师2013年6月27日一、目录一、目录2二、设计内容及要求3三、设计方法与步骤33.1语音信号的采集33.2语音信号的频谱分析43.3设计数字滤波器和画出其频谱响应53.4用滤波器对信号进行滤波73.5滤波器分析后的语音信号的波形及频谱8四、课程设计心得8五、参考文献9精彩文档实用标准文案二、设计内容及要求利用Windows下的录音机，录制了一段个人自己的语音。在C盘保存为WAV格式，然后在Matlab软件平台下．利用函数wavread对语音信号进行采样，

3、并记录下了采样频率和采样点数，在这里我们还通过函数sound引入听到采样后自己所录的一段声音。然后画出语音信号的时域波形；然后对语音信号进行频谱分析，在Matlab中，我们利用函数fft对信号进行快速傅里叶变换，得到信号的频谱特性性。紧接着着我们做了一个数字滤波器：采样频率10Hz，通带截止频率fp=3Hz，阻带截止频率fs=4Hz通带衰减小于1dB，阻带衰减大于20dB，我们主要使用双线性变换法由模拟滤波器原型设计数字滤波器。程序最后我们还利用Matlab中的函数freqz画出了该滤波器的频率响

4、应。课设报告要求：1、字体必须是四号字体，宋体字，行间距设为1.25；2、报告书内容要涵盖课程设计题目、详细设计过程、频谱图的分析、心得体会和参考文献。3、给出可运行的源程序代码。三、设计方法与步骤3.1语音信号的采集我们利用Windows下的录音机，录制了一段开枪发出的声音，时间在1s内。接着在C盘保存为WAV格式，然后在Matlab软件平台下．利用函数wavread对语音信号进行采样，并记录下了采样频率和采样点数，在这里我们还通过函数sound引入听到采样后自己所录的一段声音。通过wavrea

5、d函数和sound的使用，我们完成了本次课程设计的第一步。其程序如下：精彩文档实用标准文案[x,fs,bite]=wavread('c:alsndmgr.wav',[100020000]);sound(x,fs,bite);3.2语音信号的频谱分析首先我们画出语音信号的时域波形；然后对语音信号进行频谱分析，在Matlab中，我们利用函数fft对信号进行快速傅里叶变换，得到信号的频谱特性性。到此，我们完成了课程实际的第二部。其程序如下：n=1024;subplot(2,1,1);y=plot(x(

6、50:n/4));gridon;title('时域信号')X=fft(x,256);subplot(2,1,2);plot(abs(fft(X)));gridon;title('频域信号');运行程序得到的图形：精彩文档实用标准文案3.3设计数字滤波器和画出其频谱响应紧接着着我们做了一个数字滤波器：采样频率10Hz，通带截止频率fp=3Hz，阻带截止频率fs=4Hz通带衰减小于1dB，阻带衰减大于20dB，我们主要使用双线性变换法由模拟滤波器原型设计数字滤波器。程序最后我们还利用Matlab中的函

7、数freqz画出了该滤波器的频率响应。其程序如下：%采样频率10Hz，通带截止频率fp=3Hz，阻带截止频率fs=4Hz%通带衰减小于1dB，阻带衰减大于20dB%使用双线性变换法由模拟滤波器原型设计数字滤波器精彩文档实用标准文案T=0.1;FS=1/T;fp=3;fs=4;wp=fp/FS*2*pi;ws=fs/FS*2*pi;Rp=1;%通带衰减As=20;%阻带衰减%频率预畸OmegaP=(2/T)*tan(wp/2);%PrewarpPrototypePassbandfreqOmegaS=

8、(2/T)*tan(ws/2);%PrewarpPrototypeStopbandfreq精彩文档实用标准文案%设计butterworth低通滤波器原型N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(OmegaP/OmegaS)));OmegaC=OmegaP/((10^(Rp/10)-1)^(1/(2*N)));[z,p,k]=buttap(N);%获取零极点参数p=p*OmegaC;k=k*OmegaC^N;B=re