《用matlab对信号进行频谱分析报告及滤波》由会员分享，可在线阅读，更多相关《用matlab对信号进行频谱分析报告及滤波（11页珍藏版）》请在金锄头文库上搜索。
1、word华东交通大学课程设计论文任务书专业xxx班级xxxxx一、课程设计论文题目应用Matlab对信号进展频谱分析与滤波二、课程设计论文工作：自xxx年x月x日起至xxx年x月x日止。三、课程设计论文的容要求：序号项目等级优秀良好中等与格不与格1课程设计态度评价2出勤情况评价3任务难度评价4工作量饱满评价5任务难度评价6设计中创新性评价7论文书写规化评价8综合应用能力评价综合评定等级学生签名：20年月日课程设计论文评阅意见评阅人职称20 年月日目录l 设计过程步骤52.1 语音信号的采集52.2 语音信号的频谱分析62.3 设计数字滤波器和画出其频谱响应72.4 用滤波器对信号进展滤波9滤波
2、器分析后的语音信号的波形与频谱10心得和经验 (11)参考文献(12)l 设计过程步骤2.1 语音信号的采集我们利用Windows下的录音机，录制了一段开枪发出的声音，时间在1 s。接着在C盘保存为WAV格式，然后在Matlab软件平台下利用函数wavread对语音信号进展采样，并记录下了采样频率和采样点数，在这里我们还通过函数sound引入听到采样后自己所录的一段声音。通过wavread函数和sound的使用，我们完成了本次课程设计的第一步。其程序如下：x,fs,bite=wavread(c:alsndmgr.wav,1000 20000);sound(x,fs,bite);2.2 语音信号
3、的频谱分析首先我们画出语音信号的时域波形；然后对语音信号进展频谱分析，在Matlab中，我们利用函数fft对信号进展快速傅里叶变换，得到信号的频谱特性性。到此，我们完成了课程实际的第二部。其程序如下：n=1024;subplot(2,1,1);y=plot(x(50:n/4);grid on ;title(时域信号)X=fft(x,256);subplot(2,1,2);plot(abs(fft(X);grid on ;title(频域信号);运行程序得到的图形：2.3 设计数字滤波器和画出其频谱响应紧接着着我们做了一个数字滤波器：采样频率10Hz，通带截止频率fp=3Hz，阻带截止频率fs=
4、4Hz通带衰减小于1dB，阻带衰减大于20dB，我们主要使用双线性变换法由模拟滤波器原型设计数字滤波器。程序最后我们还利用Matlab中的函数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); % Prew
5、arp Prototype Passband freq OmegaS = (2/T)*tan(ws/2); % Prewarp Prototype Stopband freq / %设计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*OmegaCN; B = real(poly(z); b0 = k; c
6、s = k*B; ds = real(poly(p); % 双线性变换b,a = bilinear(cs,ds,FS); % 绘制结果freqz(b,a,512,FS); 运行程序得到的图形：2.4 用滤波器对信号进展滤波我们做的是IIR滤波器，所以我们利用函数filter对信号进展滤波。其程序如下：z=filter(b,a,x);subplot(2,1,1);plot(z);title(滤波后信号波形);Z=fft(z,256);subplot(2,1,2);plot(abs(Z);title(滤波后信号频谱);滤波器分析后的语音信号的波形与频谱运行程序得到滤波后的语音信号的波形与频谱图如下：心得和经验 对录音的采集后，需要保存在C盘根目录下，否如此用Matlab读取音频文件读取不了。其余操作只要输入函数按确定键得出图形后，保存就行。参考文献1 谢希仁. 计算机网络第五版M. ：电子工业，2008年2月2 小强计算机网络M ：邮电大学2005年1月致谢
