Question

数字信号处理实验fft进行谱分析，可以应用在什么方面

Answer 1

数字信号处理综合设计 

一、实验目的

1．学会MATLAB的使用，掌握MATLAB的程序设计方法；

2．掌握在Windows环境下语音信号采集的方法；

3．掌握数字信号处理的基本概念、基本理论和基本方法；

4．掌握MATLAB设计FIR和IIR数字滤波器的方法；

5．学会用MATLAB对信号进行分析和处理。

二、实验原理

参考《数字信号处理》教材。

三、主要实验仪器及材料

微型计算机、Matlab6.5教学版、TC编程环境。

四、实验内容

1．语音信号的采集

要求利用windows下的录音机（开始—程序—附件—娱乐—录音机，文件—属性—立即转换—8000KHz，8位，单声道）或其他，录制一段自己的话音，时间控制在1秒左右。然后在MATLAB下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。通过wavread函数的使用，要求理解采样频率、采样位数等概念。

wavread函数调用格式：

y=wavread(file)，读取file所规定的wav文件，返回采样值放在向量y中。

[y,fs,nbits]=wavread(file)，采样值放在向量y中，fs表示采样频率（Hz），nbits表示采样位数。

y=wavread(file,N)，读取前N点的采样值放在向量y中。

y=wavread(file,[N1,N2])，读取从N1点到N2点的采样值放在向量y中。

2．语音信号的频谱分析

要求首先画出语音信号的时域波形；然后对语音信号进行频谱分析，在MATLAB中，可以利用函数fft对信号进行快速付立叶变换，得到信号的频谱特性；从而加深对频谱特性的理解。

3．设计数字滤波器和画出频率响应

根据语音信号的特点给出有关滤波器的性能指标：1）低通滤波器性能指标，fp=1000Hz，fc=1200 Hz， As=100dB，Ap=1dB；2）高通滤波器性能指标，fc=2800 Hz，fp=3000 Hz As=100dB，Ap=1dB；3）带通滤波器性能指标，fp1=1200 Hz，fp2=3000 Hz，fc1=1000 Hz，fc2=3200 Hz，As=100dB，Ap=1dB。要求学生首先用窗函数法设计上面要求的三种滤波器，在MATLAB中，可以利用函数fir1设计FIR滤波器，然后在用双线性变换法设计上面要求的三种滤波器；之后再利用函数butter和cheby1设计上面要求的三种IIR滤波器。最后，利用MATLAB中的函数freqz画出各滤波器的频率响应。

4．用滤波器对信号进行滤波

比较FIR和IIR两种滤波器的性能，然后用性能好的各滤波器分别对采集的信号进行滤波，在MATLAB中，FIR滤波器利用函数fftfilt对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。

5．比较滤波前后语音信号的波形及频谱

要求在一个窗口同时画出滤波前后的波形及频谱。

6．回放语音信号

在MATLAB中，函数sound可以对声音进行回放。其调用格式：

sound(x,fs,bits)；

可以感觉滤波前后的声音有变化。

五、实验思考

1．双线性变换法中Ω和ω之间的关系是非线性的，在实验中你注意到这种非线性关系了吗？从哪几种数字滤波器的幅频特性曲线中可以观察到这种非线性关系？ 

2．能否利用公式完成脉冲响应不变法的数字滤波器设计？为什么？

六、实验报告要求 

1. 简述实验原理及目的。

2. 2.按照实验步骤及要求，比较各种情况下的滤波性能。

3. 3.总结实验所得主要结论。

4. 4.简要回答思考题。