数字信号处理课程设计报告

课程设计报告

课程名称数字信号处理

课题名称数字滤波器设计及在语音信号分析中的应用

专业通信工程

班级通信工程1101

学号201103020127

姓名皮锋

指导教师胡瑛张细政

2013年9月8日

湖南工程学院

课程设计任务书

课程名称数字信号处理

课题数字滤波器设计及在语音信号分析中的应用

专业班级通信工程1101

学生姓名皮锋

学号201103020127

指导老师胡瑛张细政

审批

任务书下达日期2013 年9月 1 日

任务完成日期2013 年9月8日

《数字信号处理》课程设计任务书

一、设计目的

综合运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应结论，再利用MATLAB 作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。

二、设计要求

1、MATLAB 的使用，掌握MATLAB 的程序设计方法。

2、Windows 环境下语音信号采集的方法。

3、数字信号处理的基本概念、基本理论和基本方法。

4、TLAB 设计FIR 和nR 数字滤波器的方法。

5、用 MATLAB 对信号进行分析和处理

6、计报告4000以上，含程序设计说明，用户使用说明，源程序清单及程序框图。

7、机演示。

8、有详细的文档。文档中包括设计思路、设计仿真程序、仿真结果及相应的分析与结论。

三、进度安排

第一周星期一：课题讲解，查阅资料

星期二: 总体设计，详细设计

星期三：编程，上机调试、修改程序

星期四：上机调试、完善程序

星期五：答辩

星期六-星期天：撰写课程设计报告

附：

课程设计报告装订顺序：封面、任务书、目录、正文、评分、附件（A4大小的图纸及程序清单）。

正文的格式:一级标题用3号黑体,二级标题用四号宋体加粗,正文用小四号宋体;行距为22。

正文的内容:一、课题的主要功能；二、课题的功能模块的划分；三、主要功能的实现；

四、程序调试；五、总结；六、附件（所有程序的原代码，要求对程序写出必要的注释）；

七、评分表。

目录

1. 课题的主要功能 (1)

1.1问题描述 (1)

1.2功能要求 (1)

2. 课题的功能模块的划分 (2)

2.1语音信号采集 (2)

2.3 含噪语音信号合成 (4)

2.4 数字滤波器设计及滤波 (4)

2.5 回放语音信号 (4)

3. 滤波器实现的原理 (5)

3.1 带阻滤波器 (6)

3.2 低通滤波器 (8)

3.3 高通滤波 (10)

4. 程序调试 (12)

5. 总结 (13)

6. 附件 (14)

6.1 带阻滤波 (14)

6.2 低通滤波 (15)

6.3 高通滤波 (16)

7. 评分表 (18)

1.课题的主要功能

1.1问题描述

录制一段课程设计学生的语音信号并保存为文件，要求长度不小于10秒，并对录制的信号进行采样，使用MATLAB绘出采样后的语音信号的时域波形和频谱图。在MATLAB 软件平台下，给原始的语音信号叠加上噪声，噪声类型分为如下几种：（1）白噪声；（2）单频噪色（正弦干扰）；（3）多频噪声（多正弦干扰）；（4）其它干扰，可设置为低频、高频、带限噪声，或Chirp干扰、冲激干扰。绘出叠加噪声后的语音信号时域和频谱图，在视觉上与原始语音信号图形对比，也可通过Windows播放软件从听觉上进行对比，分析并体会含噪语音信号频谱和时域波形的改变。给定滤波器的规一化性能指标（参考指标,实际中依据每个同学所叠加噪声情况而定），对滤波后的语音信号进行回放，感觉滤波前后语音信号的变化。

1.2功能要求

采用IIR滤波器的完全设计函数来设计各型IIR滤波器（Butterworth型与切比雪夫I 型、切比雪夫Ⅱ型和椭圆型的低通、高通、带通、带阻中的至少3种类型），滤掉在MATLAB 软件平台下，给原始的语音信号叠加上噪声，绘出滤波器的频域响应。并用这些数字滤波器对含噪语音信号分别进行滤波处理，比较不同方法下设计出来的数字滤波器的滤波效果，并从理论上进行分析（或解释）。对滤波后的语音信号进行回放，感觉滤波前后语音信号的变化。

2.课题的功能模块的划分

图2-1 课题的功能模块图

2.1语音信号采集

使用Windows自带的录音机，在安静、干扰小的环境下录音。按“开始”—“程序”—“附件”—“娱乐”—“录音机”的顺序打开windows中的录音机软件。如图2-1所示：

图2-2系统录音界面

录制一段自己的话音“信号”，按文件—属性—立即转换—8000KHz，8位，单声道到底顺序，时间控制在10秒以上，然后将音频文件保存“好男人.wav”

2.2语音信号分析

使用MATLAB 绘出采样后的语音信号的时域波形和频谱图，其主要的程序代码和信号时域波形图和频谱图如下： [x1,fs,bits]=wavread('i:/好男人_new'); %sound(x1,fs,bits); y1=fft(x1,32768); figure(1) subplot(2,1,1) plot(x1);

title('原始语音信号时域波形'); subplot(2,1,2) plot(abs(y1));

title('原始语音信号频谱')

0.5

1

1.5

2

2.5

3

3.5x 10

5

-1-0.500.5

1原始语音信号时域波形

0.5

1

1.5

2

2.5

3

3.5x 10

4050100150

200原始语音信号频谱

图2-3 原始语音信号波形

2.3 含噪语音信号合成

（1）在原语音信号上加单频余弦信号，其信号时域波形图和频谱图如下：

0.5

1

1.5

2

2.5

3

3.5x 10

5

-2-101

2加单频余弦信号的语音信号时域波形

00.51 1.52 2.53 3.5x 10

4

5000

10000

15000加单频余弦信号的语音信号频谱

图2-4 加单频余弦信号的波形

2.4 数字滤波器设计及滤波

采用IIR 滤波器的完全设计函数来设计各型IIR 滤波器（Butterworth 型与切比雪夫I 型、切比雪夫Ⅱ型和椭圆型的低通、高通、带通、带阻中的至少3种类型），滤掉在MATLAB 软件平台下，给原始的语音信号叠加上噪声，绘出滤波器的频域响应。并用这些数字滤波器对含噪语音信号分别进行滤波处理，比较不同方法下设计出来的数字滤波器的滤波效果，并从理论上进行分析。

2.5 回放语音信号

对滤波后的语音信号进行回放，感觉滤波前后语音信号的变化。

3. 滤波器实现的原理

MATLAB 信号处理工具箱提供了IIR 滤波器设计的完全工具函数，用户只要调用这些工具函数即可一次性完成设计，而不需要调用那些基本工具函数分步实现。

IIR 滤波器设计的完全工具函数有butter,cheby1,cheby2,ellip 。这些工具函数既可用于设计模拟滤波器，也可用于设计数字滤波器。这里介绍这些函数在IIR 数字滤波器中的应用。在这两类滤波器设计中，这些工具函数调用格式基本相同，只是在频率处理上有所不同。

在MATLAB 滤波器设计工具箱中，数字滤波器采用归一化频率，取值为0~1之间，归一化频率1对应的数字角频率为π，对应的真实频率为采样频率的一半。在应用MATLAB 工具函数设计数字滤波器时应注意这一点。

数字IIR 滤波器的完全设计函数有： [b,a]=butter(n,wn[,'ftype']) [z,p,k]=butter(n,wn[, 'ftype']) [b,a]=cheby1(n,Rp,wn[,'ftype']) [z,p,k]=cheby1(n, Rp,wn[,'ftype']) [b,a]=cheby2(n,Rs,wn[,'ftype']) [z,p,k]=cheby2(n, Rs,wn[,'ftype']) [b,a]=ellip(n,Rp,Rs,wn[,'ftype']) [z,p,k]=ellip(n, Rp,Rs,wn[,'ftype'])

在上面的调用方式中，n 为滤波器的阶数，wn 为滤波器的截止频率，取值为0~1。

需根据采样频率Fs 来定，如滤波器的截止频率为Fc （Hz ），则wn 的计算公式为： 这样就转换为0~1的归一化频率。其中wp,ws 等边界频率都要根据此公式进行转换。 'ftype'滤波器的类型为：

‘high ’为高通滤波器，截止频率为wn.

‘stop ’为带阻滤波器，截止频率为wn=[w1,w2] (w

a,b 分别为滤波器传递函数分子和分母多项式系数向量；z,p,k 分别为滤波器的零极点和增益。Rp,Rs 分别为所设计滤波器的通带波纹和阻带衰减，单位为dB 。

Fs

Fc

wn *2=

设计好的数字滤波器传递函数具有下面形式：

上述函数采用双线性变换法和频率的预畸变处理将模拟滤波器离散化为数字滤波

器，同时保证模拟滤波器和数字滤波器在wn(或w1,w2)处具有相同的幅频响应。

设计时应注意真实频率和MATLAB 归一化数字频率之间的转换，即6-20式的应用。在进行IIR 数字滤波器设计之前，请大家注意在模拟滤波器中我们用求取最小阶数和截止频率的函数如buttord,cheb1ord,cheb2ord,ellipord ，这些函数完全可以用于IIR 数字滤波器设计中，只不过在模拟滤波器设计中需加可选项's'，在数字滤波器中则不加该项。另外，输出的截止频率也是归一化频率（归一化为0~1）。

3.1 带阻滤波器

Fs=fs; %采样频率

wp=[4500 5500]*2/Fs; %通带边界频率（归一化频率） ws=[4000 6000]*2/Fs; %阻带边界频率（归一化频率）

Rp=1;Rs=30;Nn=128; %通带波纹和阻带衰减以及绘制频率特性的数据点数 [N,Wn]=cheb1ord(wp,ws,Rp,Rs);%求得数字滤波器的最小阶数和归一化截止频率 [b,a]=cheby1(N,Rp,Wn,'stop'); %按最小阶数、通带波纹和截止频率设计数字滤波器 figure(3)

[H,f]=freqz(b,a,Nn,Fs); %求得滤波器的频率特性 subplot(2,1,1),plot(f,20*log10(abs(H))); xlabel('频率/Hz');ylabel('振幅/dB');grid on; subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel('频率/Hz');ylabel('相位/^o');grid on; figure(4)

y=filtfilt(b,a,x3); %对输入信号进行滤波 y4=fft(y,32768); subplot(2,1,1)

m

n

z m a z a a z n b z b b z A z B z H ----++++++++=

=)1()2()1()1()2()1()()()(11

plot(y);

title('加单频余弦信号滤波后的语音信号时域波形');

subplot(2,1,2) plot(abs(y4));

title('加单频余弦信号滤波后的语音信号频谱')% %sound(y,fs,bits);

图3-1 带阻滤波器性能指标

0200040006000

800010000120001400016000

-150

-100

-50

频率/Hz

振幅/d B

0200040006000

800010000120001400016000

-800

-600-400-200

0频率/Hz

相位/o

0.5

1

1.5

2

2.5

3

3.5x 10

5

-1-0.500.5

1加单频余弦信号滤波后的语音信号时域波形

00.51 1.52 2.53 3.5x 10

4

50

100

150加单频余弦信号滤波后的语音信号频谱

图3-2 加单频余弦信号滤波后波形

3.2 低通滤波器

Fs=fs; %采样频率

wp=4000*2/Fs;ws=4300*2/Fs; %根据采样频率将滤波器边界频率进行转换 Rp=1;Rs=20; %通带波纹和阻带衰减

Nn=128; %显示滤波器频率特性的数据长度

[N,Wn]=buttord(wp,ws,Rp,Rs); %求得数字滤波器的最小阶数和截止频率（归一化频率）

[b,a]=butter(N,Wn); %设计Butterworth 低通数字滤波器 figure(3)

[H,f]=freqz(b,a,Nn,Fs); %用Nn 点绘出频率特性 subplot(2,1,1),plot(f,20*log10(abs(H))); xlabel('频率/Hz');ylabel('振幅/dB');grid on; subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel('频率/Hz');ylabel('相位/^o');grid on; figure(4)

y=filter(b,a,x3); %对输入信号进行滤波 y4=fft(y,32768); subplot(2,1,1) plot(y);

title('加单频余弦信号滤波后的语音信号时域波形'); subplot(2,1,2) plot(abs(y4));

title('加单频余弦信号滤波后的语音信号频谱') sound(y,fs,bits);

0200040006000

800010000120001400016000

-1000

-500

500

频率/Hz

振幅/d B

0200040006000

800010000120001400016000

-4000

-3000-2000-1000

0频率/Hz

相位/o

图3-3 低通滤波器性能指标

0.5

1

1.5

2

2.5

3

3.5x 10

5

-1-0.500.5

1加单频余弦信号滤波后的语音信号时域波形

0.5

1

1.5

2

2.5

3

3.5x 10

4050100150

200加单频余弦信号滤波后的语音信号频谱

图3-4 加单频余弦信号滤波后波形

3.3 高通滤波

Fs=1000; %采样频率

wp=50*2/Fs;ws=20*2/Fs; %根据采样频率将滤波器边界频率进行转换 Rp=1;Rs=20; %通带波纹和阻带衰减

Nn=128; %显示滤波器频率特性的数据长度

[N,Wn]=buttord(wp,ws,Rp,Rs); %求得数字滤波器的最小阶数和截止频率（归一化频率）

[b,a]=butter(N,Wn,'high'); %设计Butterworth 高通数字滤波器 figure(3)

[H,f]=freqz(b,a,Nn,Fs); %用Nn 点绘出频率特性 subplot(2,1,1),plot(f,20*log10(abs(H))); xlabel('频率/Hz');ylabel('振幅/dB');grid on; subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel('频率/Hz');ylabel('相位/^o');grid on; figure(4)

y=filter(b,a,x3); %对输入信号进行滤波 y4=fft(y,32768); subplot(2,1,1) plot(y);

title('加单频余弦信号滤波后的语音信号时域波形'); subplot(2,1,2) plot(abs(y4));

title('加单频余弦信号滤波后的语音信号频谱')% sound(y,fs,bits);

050100150200

250300350400450500

-80

-60-40-20

0频率/Hz

振幅/d B

050100150200

250300350400450500

-400

-300-200-100

0频率/Hz

相位/o

图3-5 高通滤波器性能指标

0.5

1

1.5

2

2.5

3

3.5x 10

5

-1-0.500.5

1加单频余弦信号滤波后的语音信号时域波形

0.5

1

1.5

2

2.5

3

3.5x 10

4050100150

200加单频余弦信号滤波后的语音信号频谱

图3-6 加单频余弦信号滤波后波形

4. 程序调试

在设计低通滤波器时，首先通带边界频率wp 和阻带边界频率ws 取值不当，分别取3500*2/Fs 和4000*2/Fs ，使滤波效果太差，后来经过不断调试才取值为4000*2/Fs 和4300*2/Fs ，使滤波效果很理想。

对输入信号进行滤波时，函数总是书写错误，写成y=filter(b,a,x3)，y4=fft(x3,32768)，因为没有理解这两个函数，所以不会修改，后来问老师才知道，要把y4=fft(x3,32768)中的x3改成y ，进过这样一改，结果就正确了。

在设计带阻滤波器时，语句[b,a]=cheby1(N,Rp,Wn,'stop')中有个stop 表示带阻，但首先一直没加，后来看以前做过的上机实验才知道写错了。

在运行程序的过程中，总是发生错误终止。如下图所示：

图4-1程序调试错误图示

后来经过调试才知道程序中有多处应用了sound 函数，才导致这样的错误，在注释掉多余的sound函数，只剩下最后一个后，程序正常工作，而且听到了滤波后干净的歌声。

5.总结

经过为期一周的数字信号处理课程设计，我重新学到了许多新的知识，而且复习了以前学过但是却逐渐淡忘的是知识。我这次做的题目是采用IIR滤波器的完全设计函数来设计各型IIR滤波器（Butterworth型与切比雪夫I型、切比雪夫Ⅱ型和椭圆型的低通、高通、带通、带阻中的至少3种类型），绘出滤波器的频域响应；并用这些数字滤波器对含噪语音信号分别进行滤波处理，比较不同方法下设计出来的数字滤波器的滤波效果，并从理论上进行分析。我耐心的坚持把课程设计完成。有互联网这个帮手，对于不懂的可以从网上搜索，也可以查找很多有用的函数、源代码及其它有用的资料。我参照课本知识对我找到的那些源代码加以修改和整合，在用上老师提供的那些加噪的源代码做出程序。虽然中间遇到很多问题，但通过自己的仔细排查，老师的帮助，同学之间的互相讨论基本上解决了。

总之，通过这次数字信号课程设计，我受益匪浅。

6.附件

程序源代码：

6.1 带阻滤波

[x1,fs,bits]=wavread('i:/好男人_new');

sound(x1,fs,bits);

y1=fft(x1,32768);

figure(1)

subplot(2,1,1)

plot(x1);

title('原始语音信号时域波形');

subplot(2,1,2)

plot(abs(y1));

title('原始语音信号频谱')

N=length(x1)-1;

t=0:1/fs:N/fs;

d=[0.8*cos(2*pi*5000*t)]';

x3=x1+d;

%sound(x3,fs,bits);

figure(2)

subplot(2,1,1)

plot(x3)

title('加单频余弦信号的语音信号时域波形')

y3=fft(x3,32768);

subplot(2,1,2)

plot(abs(y3));

title('加单频余弦信号的语音信号频谱')

Fs=fs; %采样频率

wp=[4500 5500]*2/Fs; %通带边界频率（归一化频率）

ws=[4000 6000]*2/Fs; %阻带边界频率（归一化频率）

Rp=1;Rs=30;Nn=128; %通带波纹和阻带衰减以及绘制频率特性的数据点数

[N,Wn]=cheb1ord(wp,ws,Rp,Rs);%求得数字滤波器的最小阶数和归一化截止频率[b,a]=cheby1(N,Rp,Wn,'stop'); %按最小阶数、通带波纹和截止频率设计数字滤波器

figure(3)

[H,f]=freqz(b,a,Nn,Fs); %求得滤波器的频率特性

subplot(2,1,1),plot(f,20*log10(abs(H)));

xlabel('频率/Hz');ylabel('振幅/dB');grid on;

subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))

xlabel('频率/Hz');ylabel('相位/^o');grid on;

figure(4)

y=filtfilt(b,a,x3); %对输入信号进行滤波

y4=fft(y,32768);

subplot(2,1,1)

plot(y);

title('加单频余弦信号滤波后的语音信号时域波形');

subplot(2,1,2)

plot(abs(y4));

title('加单频余弦信号滤波后的语音信号频谱')%

%sound(y,fs,bits);

6.2 低通滤波

[x1,fs,bits]=wavread('i:/好男人_new');

%sound(x1,fs,bits);

y1=fft(x1,32768);

figure(1)

subplot(2,1,1)

plot(x1);

title('原始语音信号时域波形');

subplot(2,1,2)

plot(abs(y1));

title('原始语音信号频谱')

N=length(x1)-1;

t=0:1/fs:N/fs;

d=[0.8*cos(2*pi*5000*t)]';

x3=x1+d;

%sound(x3,fs,bits);

figure(2)

subplot(2,1,1)

plot(x3)

title('加单频余弦信号的语音信号时域波形')

y3=fft(x3,32768);

subplot(2,1,2)

plot(abs(y3));

title('加单频余弦信号的语音信号频谱')

Fs=fs; %采样频率

wp=4000*2/Fs;ws=4300*2/Fs; %根据采样频率将滤波器边界频率进行转换

Rp=1;Rs=20; %通带波纹和阻带衰减

Nn=128; %显示滤波器频率特性的数据长度

[N,Wn]=buttord(wp,ws,Rp,Rs); %求得数字滤波器的最小阶数和截止频率（归一化频率）[b,a]=butter(N,Wn); %设计Butterworth低通数字滤波器

figure(3)

[H,f]=freqz(b,a,Nn,Fs); %用Nn点绘出频率特性

subplot(2,1,1),plot(f,20*log10(abs(H)));

xlabel('频率/Hz');ylabel('振幅/dB');grid on;

subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))

xlabel('频率/Hz');ylabel('相位/^o');grid on;

figure(4)

y=filter(b,a,x3); %对输入信号进行滤波

y4=fft(y,32768);

subplot(2,1,1)

plot(y);

title('加单频余弦信号滤波后的语音信号时域波形');

subplot(2,1,2)

plot(abs(y4));

title('加单频余弦信号滤波后的语音信号频谱')

sound(y,fs,bits);

6.3 高通滤波

[x1,fs,bits]=wavread('i:/好男人_new');

%sound(x1,fs,bits);

y1=fft(x1,32768);

figure(1)

subplot(2,1,1)

plot(x1);

title('原始语音信号时域波形');

subplot(2,1,2)

plot(abs(y1));

title('原始语音信号频谱')

N=length(x1)-1;

t=0:1/fs:N/fs;

d=[0.8*cos(2*pi*40*t)]';

x3=x1+d;

%sound(x3,fs,bits);

figure(2)

subplot(2,1,1)

plot(x3)

title('加单频余弦信号的语音信号时域波形')

y3=fft(x3,32768);

subplot(2,1,2)

plot(abs(y3));

title('加单频余弦信号的语音信号频谱')

Fs=1000; %采样频率

wp=50*2/Fs;ws=20*2/Fs; %根据采样频率将滤波器边界频率进行转换

Rp=1;Rs=20; %通带波纹和阻带衰减

Nn=128; %显示滤波器频率特性的数据长度

[N,Wn]=buttord(wp,ws,Rp,Rs); %求得数字滤波器的最小阶数和截止频率（归一化频率）

数字信号处理课程设计报告

抽样定理的应用 摘要 抽样定理表示为若频带宽度有限的，要从抽样信号中无失真地恢复原信号，抽样频率应大于2倍信号最高频率。抽样频率小于2倍频谱最高频率时，信号的频谱有混叠。抽样频率大于2倍频谱最高频率时，信号的频谱无混叠。 语音信号处理是研究用数字信号处理技术和语音学知识对语音 信号进行处理的新兴学科，是目前发展最为迅速的学科之一，通过语音传递信息是人类最重要，最有效，最常用和最方便的交换信息手段，所以对其的研究更显得尤为重要。 Matlab语言是一种数据分析和处理功能十分强大的计算机应用 软件，它可以将声音文件变换成离散的数据文件，然后用起强大的矩阵运算能力处理数据。这为我们的本次设计提供了强大并良好的环境！ 本设计要求通过利用matlab对模拟信号和语音信号进行抽样，通过傅里叶变换转换到频域，观察波形并进行分析。 关键词：抽样Matlab

目录 一、设计目的： (2) 二、设计原理： (2) 1、抽样定理 (2) 2、MATLAB简介 (2) 3、语音信号 (3) 4、Stem函数绘图 (3) 三、设计内容： (4) 1、已知g1(t)=cos(6πt),g2(t)=cos(14πt),g3(t)=cos(26πt),以抽样频率 fsam=10Hz对上述三个信号进行抽样。在同一张图上画出g1(t)，g2(t)，g3(t)及其抽样点,对所得结果进行讨论。 (4) 2、选取三段不同的语音信号，并选取适合的同一抽样频率对其进 行抽样，画出抽样前后的图形，并进行比较，播放抽样前后的语音。 (6) 3、选取合适的点数，对抽样后的三段语音信号分别做DFT，画图 并比较。 (10) 四、总结 (12) 五、参考文献 (13)

dsp课程设计实验报告

DSP 课程设计实验 一、语音信号的频谱分析： 要求首先画出语音信号的时域波形，然后对语音信号进行频谱分析。在MATLAB 中，可以利用函数fft 对信号进行快速傅立叶变换，得到信号的频谱特性，从而加深对频谱特性的理解。 其程序为： >> [y,fs,bits]=wavread('I:\',[1024 5120]); >> sound(y,fs,bits); >> Y=fft(y,4096); >> subplot(221);plot(y);title('原始信号波形'); | >> subplot(212);plot(abs(Y));title('原始信号频谱'); 程序运行结果为： 二、设计数字滤波器和画出频率响应： 根据语音信号的特点给出有关滤波器的性能指标： 低通滤波器性能指标，p f =1000Hz ，c f =1200Hz ，s A =100dB ，p A =1dB ； 高通滤波器性能指标，c f =4800Hz ，p f =5000Hz ，s A =100dB ，p A =1dB ； 带通滤波器性能指标，1p f =1200Hz ，2p f =3000Hz ，1c f =1000Hz ，2c f =3200Hz ，s A =100dB ， p A =1dB ；

】 要求学生首先用窗函数法设计上面要求的三种滤波器，在MATLAB中，可以利用函数firl 设计FIR滤波器；然后再用双线性变换法设计上面要求的三种滤波器，在MATLAB中，可以利用函数butte、cheby1和ellip设计IIR滤波器；最后，利用MATLAB中的函数freqz画出各种滤波器的频率响应，这里以低通滤波器为例来说明设计过程。 低通： 用窗函数法设计的低通滤波器的程序如下： >> fp=1000;fc=1200;As=100;Ap=1;fs=22050; >> wc=2*fc/fs;wp=2*fp/fs; >> N=ceil(/*(wc-wp)/2))+1; >> beta=*; >> Win=Kaiser(N+1,beta); 、 >>b=firl(N,wc,Win); >>freqz(b,1,512,fs); 程序运行结果： 这里选用凯泽窗设计，滤波器的幅度和相位响应满足设计指标，但滤波器长度（N=708）太长，实现起来很困难，主要原因是滤波器指标太苛刻，因此，一般不用窗函数法设计这种类型的滤波器。 用双线性变换法设计的低通滤波器的程序如下： >> fp=1000;fc=1200;As=100;Ap=1;fs=22050; >> wc=2*fc/fs;wp=2*fp/fs; 》 >> [n,wn]=ellipord(wp,wc,Ap,As); >> [b,a]=ellip(n,Ap,As,wn); >> freqz(b,a,512,fs); ^

DSP课程设计总结报告

课程设计总结报告课程名称DSP控制器及其应用 设计题目万年历设计 业专电子信息工程 班级 姓名 学号

指导教师 报告成绩 信息工程学院 年六月十三日二〇一四 录目 言前 (3) 设计要求第一章4.....................................................................................基本要求1.14.....................................................................................

系统的组成和工作原理第二章5............................................................. 芯片的工作原理VC5509APGE2.1DSPTMS3205.............................. 液晶显示器的工作原理2.2LCD16026..............................................主电路图及程序流程图第三章.. (7) 主电路图3.17...................................................................................... 程序总流程图3.27.............................................................................. 程序分块流程图3.38..........................................................................软件程序设计第四章9.............................................................................

数字信号处理课设+语音信号的数字滤波

语音信号的数字滤波 ——利用双线性变换法实现IIR数字滤波器的设计一．课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法，掌握利用数字滤波器对信号进行滤波的方法，掌握数字滤波器的计算机仿真方法，并能够对设计结果加以分析。 二．设计方案论证 1.IIR数字滤波器设计方法 IIR数字滤波器是一种离散时间系统，其系统函数为 假设M≤N，当M＞N时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数和，它 是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。 2.用双线性变换法设计IIR数字滤波器 脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到Ｚ平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用z=e sT转换 平面的-π/T～π到Z平面上。也就是说，第一步先将整个S平面压缩映射到S 1 /T一条横带里；第二步再通过标准变换关系z=e s1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图1所示。 图1双线性变换的映射关系 为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的-π/T到π/T段上，可以通过以下的正切变换实现

dsp课程设计实验报告总结

DSP课程设计总结（2013-2014学年第2学期） 题目： 专业班级：电子1103 学生姓名：万蒙 学号：11052304 指导教师： 设计成绩： 2014 年6 月

目录 一设计目的----------------------------------------------------------------------3 二系统分析----------------------------------------------------------------------3 三硬件设计 3.1 硬件总体结构-----------------------------------------------------------3 3.2 DSP模块设计-----------------------------------------------------------4 3.3 电源模块设计----------------------------------------------------------4 3.4 时钟模块设计----------------------------------------------------------5 3.5 存储器模块设计--------------------------------------------------------6 3.6 复位模块设计----------------------------------------------------------6 3.7 JTAG模块设计--------------------------------------------------------7 四软件设计 4.1 软件总体流程-----------------------------------------------------7 4.2 核心模块及实现代码---------------------------------------8 五课程设计总结-----------------------------------------------------14

DSP课程设计 (3)

深圳大学考试答题纸 (以论文、报告等形式考核专用) 二○一四～二○一五学年度第1 学期 课程编号 课程 名称 单片机/ARM /DSP技术实践 主讲 教师 评 分 学号姓名 专业年级 题 目： 基于DSP2812的课程设计

一、实验要求 由外接的信号发生器产生一正弦信号（电压范围：0~3V），通过DSP的AD功能对此正弦信号进行采集，通过DSP的SCI功能与PC机之间进行通信，把所采集的AD信号发送至PC机端，在超级终端上进行实时显示。 二、实验原理 2.1 ADC概述 ADC，即模/数转换器，将模拟量转换成数字量，提供给控制器使用。TMS320F2812片上有一个12位分辨率、具有流水线结构的模/数转换器，其机构框图如图1所示。其前端为2个8选1多路切换器和2路同时采样/保持器，构成16个模拟输入通道，模拟通道的切换由硬件自动控制，并将各模拟通道的转换结果顺序存入16个结果寄存器中。 图1 ADC机构框图 2.2 ADC模块特点 （1）带2个8选1多路切换器和双采样/保持器的12位的ADC，共有16个模拟输入通道； （2）模拟量输入范围：0.0V-3.0V；

（3）转换率：在25MHZ的ADC时钟下为80ns; (4)转换结果存储在16个结果存储器中； （5）转换结果=4095*（输入的模拟信号-ADCLO）/3； （6）多种A/D触发方式：软件启动、EVA和EVB； （7）灵活中断方式：可以在每次转换结束或每隔一次转换结束触发中断； 3．AD C转换步骤 （1）初始化DSP系统； （2）设置PIE中断矢量表； （3）初始化ADC模块； （4）将ADC中断的入口地址装入PIE中断矢量表中，开中断； （5）软件启动ADC转换； （6）等待ADC中断； （7）在ADC中断中读取ADC转换结果，软件启动下一次ADC中断。 三、实验实现 3.1硬件方案设计 本实验以TMS320F2812为核心控制部件，利用软件编程，通过ADC模块对试验箱上的信号发生器发出的正弦信号进行采集，由于试验箱上的信号发生器只能调节到2V，所以此次实验只针对2V的正弦信号，再通过串口线与PC机连接，将采集转换的数字信号传送到PC机端的串口助手，并还原成采集时的电压值。硬件框架图如图2所示。本次ADC采用SEED-DEC2812的AD接口的ADCINA6通道。 图2 硬件框架图

数字信号处理课设共18页文档

数字信号处理课程设计 姓名：刘倩 学号：201014407 专业：信息与计算科学 实验一：常见离散信号产生和实现 一、实验目的： 1、加深对常用离散信号的理解； 2、掌握matlab 中一些基本函数的建立方法。 二、实验原理： 1.单位抽样序列 在MATLAB 中可以利用zeros()函数实现。 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： 2．单位阶越序列 在MATLAB 中可以利用ones()函数实现。 3．正弦序列 在MATLAB 中 4．复指数序列 在MATLAB 中 5．指数序列 在MATLAB 中

实验内容：由周期为10的正弦函数生成周期为20的余弦函数。 实验代码： n=0:30; y=sin(0.2*pi*n+pi/2); y1=sin(0.1*pi*n+pi/2); subplot(121) stem(n,y); xlabel ('时间序列n');ylabel('振幅');title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); subplot(122) stem(n,y1); xlabel ('时间序列n');ylabel('振幅'); title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); 实验结果： 实验二：离散系统的时域分析 实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。实验原理：离散系统 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号， 记系统单位冲激响应 则系统响应为如下的卷积计算式：

当N k d k ,...2,1,0==时，h[n]是有限长度的（n ：[0，M]），称系统为FIR 系统；反之，称系统为IIR 系统。 在MATLAB 中，可以用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应。 实验内容：用MATLAB 计算全解 当n>=0时，求用系数差分方程y[n]+y[n-1]-6y[n-2]=x[n]描述的一个离散时间系统对阶跃输入x[n]=8μ[n]的全解。 实验代码： n=0:7; >> [y,sf]=filter(1,[1 1 -6],8*ones(1,8),[-7 6]); >> y1(n+1)=-1.8*(-3).^n+4.8*(2).^n-2; >> subplot(121) >> stem(n,y); >> title('由fliter 函数计算结果'); >> subplot(122) >> stem(n,y1); >> title('准确结果'); 实验结果： 结果分析：有图可得由fliter 函数得出的结果与计算出的准确结果完全一致。 实验三FFT 算法的应用

DSP实验报告

电气信息工程学院 D S P技术与综合训练 实验报告 班级 08通信1W 姓名丁安华 学号 08313115 指导老师倪福银刘舒淇 2011年09 月

目录 实验一 LED演示 1.1.实验目的 -------------------------------------------------P2 1. 2.实验设备-------------------------------------------------P2 1. 3.实验原理-------------------------------------------------P2 1. 4.实验程序设计流程------------------------------------------P3 1. 5.实验程序编写----------------------------------------------P4 1. 6.实验步骤-------------------------------------------------P7 1. 7.实验结果与分析--------------------------------------------P7实验二键盘输入 2.1.实验目的 -------------------------------------------------P8 2.2.实验设备-------------------------------------------------P8 2. 3.实验原理-------------------------------------------------P8 2. 4.实验程序设计流程------------------------------------------P9 2. 5.实验程序编写----------------------------------------------P10 2. 6.实验步骤-------------------------------------------------P14 2. 7.实验结果与分析--------------------------------------------P14实验三液晶显示器控制显示 3.1.实验目的 -------------------------------------------------P15 3.2.实验设备-------------------------------------------------P15 3.3.实验原理-------------------------------------------------P15 3. 4.实验程序设计流程------------------------------------------P17 3. 5.实验程序编写----------------------------------------------P18 3. 6.实验步骤-------------------------------------------------P22 3. 7.实验结果与分析--------------------------------------------P23实验四有限冲激响应滤波器（FIR）算法 4.1.实验目的 -------------------------------------------------P23 4.2.实验设备-------------------------------------------------P23 4.3.实验原理-------------------------------------------------P24 4.4.实验程序设计流程------------------------------------------P25 4. 5.实验程序编写----------------------------------------------P25 4. 6.实验步骤-------------------------------------------------P27 4. 7.实验结果与分析--------------------------------------------P28

数字信号处理课程规划报告

数字信号处理课程设计报告《应用Matlab对信号进行频谱分析及滤波》 专业： 班级： 姓名： 指导老师： 二0 0五年一月一日

目录 设计过程步骤（） 2.1 语音信号的采集（） 2.2 语音信号的频谱分析（） 2.3 设计数字滤波器和画出其频谱响应（） 2.4 用滤波器对信号进行滤波（） 2.5滤波器分析后的语音信号的波形及频谱（） ●心得和经验（）

设计过程步骤 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 语音信号的频谱分析 首先我们画出语音信号的时域波形；然后对语音信号进行频谱分析，在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('频域信号'); 运行程序得到的图形：

DSP实验报告

实验一 程序的控制与转移 一、实验目的 1、掌握条件算符的使用。 2、掌握循环操作指令（BNAZ ）和比较操作指令（CMPR ） 二、实验设备 计算机、ZY13DSP12BD 实验箱、5402EVM 板。 三、实验原理 程序控制指令主要包括分支转移、子程序调用、子程序返回、条件操作及循环操作等。通过传送控制到程序存储器的其他位置，转移会中断连续的指令流。转移会影响在PC 中产生和保护的程序地址。其中转移可以分为两种形式的，一种是有条件的，另一种是无条件的。 四、实验内容 编写程序，实现计算y= ∑=5 1 i i x 的值。 五、实验步骤 1、用仿真机将计算机与ZY13DSP12BD 实验箱连接好，并依次打开实验箱电源、仿真机电源，然后运行CCS 软件。 2、新建一个项目：点击Project －New ，将项目命名为example2，并将项目保存在自己定义的文件夹下。 3、新建一个源文件example2.asm 。将该文件添加到工程example2.pjt 中。 4、在工程管理器中双击example2.asm ，编写源程序： .tiltle ”example2.asm ” .mmregs STACK .usect ”STACK ”,10H ;堆栈的设置 .bss x,5 ;为变量分配6个字的存储空间 .bss y,1 .def start .data table: .word 10,20,3,4,5 ;x1,x2,x3,x4,x5 .text Start: STM #0,SWWWSR ;插入0个等待状态 STM #STACK+10H,sp ;设置堆栈指针 STM #x,AR1 ;AR1指向x RPT #4 ;下一条被重复执行5遍 MVPD table,*AR1+ ;把程序存储器中的数据传送到数据存储器 LD #0,A ;A 清零 CALL SUM ;调用求和函数 end: B end SUM: STM #x,AR3 ;AR3指向x STM #4,AR2 ;AR2=4 loop: ADD *AR3+,A ;*AR3+A-->A,然后AR3+ BANZ loop,*AR2- ;如果AR2的值不为0,则跳到loop 处;否则执行下一条指令 STL A,*(y) ;把A 的低16位赋给变量y

DSP课程设计报告

共享知识分享快乐 盛年不重来，一日难再晨。及时宜自勉，岁月不待人。 数据采集处理和控制系统设计 一课程设计要求 1.基本DSP硬件系统设计要求 ①基本DSP硬件系统以TMS320C54x系列为核心处理器，包括最小系统、存储器扩展、显示器、键盘、AD、DA等电路模块； ②硬件设计画出主要芯片及电路模块之间的连接即可，重点考查电路模块方案设计与系统地址分配； ③设计方案以电路示意图为主，辅以必要的文字说明。 2.基本软件设计要求 ①看懂所给例程，画出例程输出波形示意图； ②修改例程程序，使之输出其它波形，如方波、三角波、锯齿波等均可； ③设计方案以程序实现为主，辅以必要的文字说明。 3.课程设计报告要求 ①硬件系统设计：设计思路、设计系统功能、主要芯片选型及使用方法、设计方案说明、电路示意图 ②软件系统设计：示例程序功能解读及输出波形示意图、设计软件功能、设计思路、实现源码（带程序注释） ③报告总结 二系统分析 利用实验箱的模拟信号产生单元产生不同频率的信号，或者产生两个频率的信号的叠加。在DSP 中采集信号，并且对信号进行频谱分析，滤波等。通过键盘或者串口命令选择算法的功能，将计算的信号频率或者滤波后信号的频率在LCD 上显示。主要功能如下： (1)对外部输入的模拟信号采集到DSP 内存，会用CCS 软件显示采集的数据波形。 (2)对采集的数据进行如下算法分析： ①频谱分析：使用fft 算法计算信号的频率。 ②对信号进行IIR 滤波或FIR 滤波，并且计算滤波前后信号的频率。 ③外部键盘或者从计算机来的串口命令选择算法功能，并且将结果在 LCD 上显示。 绘制出DSP系统的功能框图、使用AD（Altium Designer）绘制出系统的原理图和PCB 版图。 在 DSP 中采集信号，用CCS 软件显示采集的数据波形，以及对采集的数据进行算法分析。 三硬件设计 3.1 硬件总体结构

数字信号处理课程设计

数字信号处理 课 程 设 计 院系：电子信息与电气工程学院 专业：电子信息工程专业 班级：电信班 姓名： 学号： 组员：

摘要 滤波器设计在数字信号处理中占有极其重要的地位，FIR数字滤波器和IIR 滤波器是滤波器设计的重要组成部分。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR和IIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。 关键词数字滤波器 MATLAB 窗函数法巴特沃斯

目录 摘要 (1) 1 引言 (1) 1.1课程设计目的 (1) 1.2 课程设计内容及要求 (1) 1.3课程设计设备及平台 (1) 1.3.1 数字滤波器的简介及发展 (1) 1.3.2 MATLAB软件简介 (2) 2 课程设计原理及流程 (4) 3.课程设计原理过程 (4) 3.1 语音信号的采集 (4) 3.2 语音信号的时频分析 (5) 3.3合成后语音加噪声处理 (7) 3.3.1 噪声信号的时频分析 (7) 3.3.2 混合信号的时频分析 (8) 3.4滤波器设计及消噪处理 (10) 3.4.1 设计IIR和FIR数字滤波器 (10) 3.4.2 合成后语音信号的消噪处理 (13) 3.4.3 比较滤波前后语音信号的波形及频谱 (13) 3.4.4回放语音信号 (15) 3.5结果分析 (15) 4 结束语 (15) 5 参考文献 (16)

DSP技术与课程设计实验报告二(精)

东南大学自动化学院 实验报告 课程名称： D SP 原理及C 程序开发 第二次实验 实验名称：基于DSP 系统的实验——指示灯、拨码开关和定时器院（系）：自动化专业：自动化 姓名：学号： 实验室：实验组别： 同组人员：实验时间：2012 年 4 月 18日 评定成绩：审阅教师： 第一部分实验：基于DSP 系统的实验——指示灯和拨码开关 一．实验目的 1. 了解ICETEK –F28335-A 评估板在TMS320F28335DSP 外部扩展存储空间上的扩展。 2. 了解ICETEK –F28335-A 评估板上指示灯和拨码开关扩展原理。 3. 学习在C 语言中使用扩展的控制寄存器的方法。 二．实验设备 计算机，ICETEK –F28335-A 实验箱（或ICETEK 仿真器+ICETEK–F28335-A 评估板+相关连线及电源）。 三．实验原理

1．TMS320F28335DSP 的存储器扩展接口 存储器扩展接口是DSP 扩展片外资源的主要接口，它提供了一组控制信号和地址、数据线，可以扩展各类存储器和存储器、寄存器映射的外设。 -ICETEK –F28335-A 评估板在扩展接口上除了扩展了片外SRAM 外，还扩展了指示灯、DIP 开关和D/A 设备。具体扩展地址如下： 0x180004- 0x180005：D/A 转换控制寄存器 0x180001：板上DIP 开关控制寄存器 0x180000：板上指示灯控制寄存器 -与ICETEK –F28335-A 评估板连接的ICETEK-CTR 显示控制模块也使用扩展空间控制主要设备： 208000-208004h ：读-键盘扫描值，写-液晶控制寄存器 208002-208002h ：液晶辅助控制寄存器 208003-208004h ：液晶显示数据寄存器 2．指示灯与拨码开关扩展原理

dsp课程设计报告(2)分析

华北水利水电大学North China University of Water Resources and Electric Power DSP课程设计 题目： FIR数字低通滤波器 学院信息工程学院 专业电子信息工程 姓名 学号 指导教师

摘要 (1) 一. 绪论 (1) 1.1设计背景 (1) 1.2设计方法分析 (1) 二．FIR滤波器设计过程 (2) 2.1 FIR滤波器原理 (2) 2.2 FIR滤波器的实现方法 (3) 2.3 FIR滤波器的MATLAB实现 (4) 2.4 设计流程图 (6) 三．MATLAB和 CCS操作步骤及仿真结果 (7) 3.1 matlab中的.M文件的编写 (7) 3.2 工程文件的建立 (12) 3.3 仿真结果及分析 (12) 四．心得与总结 (12)

摘要 当前,数字信号处理技术受到了人们的广泛关注,其理论及算法随着计算机技术和微电子技术的发展得到了飞速地发展,并被广泛应用于语音和图象处理、数字通信、谱分析、模式识别和自动控制等领域。数字滤波器是数字信号处理中最重要的组成部分之一,几乎出现在所有的数字信号处理系统中。设计中通过MATLAB环境中图形化的方式建立数字信号处理的模型进行DSP的设计和仿真验证，将设计的图形文件.mdl直接转换成C语言程序在CCS中运行。利用MATLAB 软件开发产品加速了开发周期，比直接在CCs中编程方便快捷了很多，对于任何复杂功能的DSP系统，只需要进行少量的添加和修改就能完成功能正确的C语言程序设计。 一. 绪论 1.1设计背景 一个实际的应用系统中，由于设备或者是外界环境的原因，总存在各种干扰，使信号中混入噪声，譬如音频信号中高频成分的噪声使得音乐听起来刺耳，失去了原有悦耳的音质。为了提高信号质量，可以对信号进行滤波，从噪声中提取信号，即对一个具有噪声和信号的混合源进行采样，然后经过一个数字滤波器，滤除噪声，提取有用信号。DSP（数字信号处理器）与一般的微处理器相比有很大的区别，它所特有的系统结构、指令集合、数据流程方式为解决复杂的数字信号处理问题提供了便利，本文选用TMS320C54X作为DSP处理芯片，通过对其编程来实现FIR滤波器。对数字滤波器而言,从实现方法上,有FIR滤波器和无限冲激响应(IIR)滤波器之分。由于FIR滤波器只有零点,因此这一类系统不像IIR系统那样易取得比较好的通带与阻带衰减特性。但是FIR系统有自己突出的优点:①系统总是稳定的;②易实现线性相位;③允许设计多通带(阻带)滤波器。其中后两项是IIR系统不易实现的。 1.2设计方法分析 FIR滤波器的设计方法分析 数字滤波器依据冲激响应的宽度划分为有限冲激响应(FIR)滤波器和无限冲激响应滤波(IIR)。FIR 滤波器是有限长单位冲激响应滤波器，在结构上是非递归型的，有限冲激响应滤波器(FIR)，具有以下的优点：（1）可以在幅度特性随意设计的同时，保证精确、严格的线性相位；（2）由于FIR滤波器的单位脉冲响应h（n）是有限长序列，因此F I R 滤波器没有不稳定的问题；（3）由于FIR 滤

dsp实验报告 哈工大实验三 液晶显示器控制显示实验

实验三液晶显示器控制显示实验 一. 实验目的 通过实验学习使用2407ADSP 的扩展I/O 端口控制外围设备的方法，了解液晶显示器的显示控制原理及编程方法。 二. 实验设备 计算机，ICETEK-LF2407-EDU 实验箱。 三.实验原理 ICETEK-LF2407-A 是一块以TMS320LF2407ADSP 为核心的DSP 扩展评估板，它通过扩展接口与实验箱的显示/控制模块连接，可以控制其各种外围设备。 液晶显示模块的访问、控制是由2407ADSP 对扩展I/O 接口的操作完成。 控制I/O 口的寻址：命令控制I/O 接口的地址为0x8001，数据控制I/O 接口的地址为0x8003 和0x8004，辅助控制I/O 接口的地址为0x8002。 显示控制方法： ◆液晶显示模块中有两片显示缓冲存储器，分别对应屏幕显示的象素，向其中写入数 值将改变显示，写入“1”则显示一点，写入“0”则不显示。其地址与象素的对应 方式如下： ◆发送控制命令：向液晶显示模块发送控制命令的方法是通过向命令控制I/O 接口 写入命令控制字，然后再向辅助控制接口写入0。下面给出的是基本命令字、解释 和 C 语言控制语句举例。 ?显示开关：0x3f 打开显示；0x3e 关闭显示； ?设置显示起始行：0x0c0+起始行取值，其中起始行取值为0 至63； ?设置操作页：0x0b8+页号，其中页号取值为0-7； ?设置操作列：0x40+列号，其中列号为取值为0-63； ◆写显示数据：在使用命令控制字选择操作位置(页数、列数)之后，可以将待显示的 数据写入液晶显示模块的缓存。将数据发送到相应数据控制I/O 接口即可。

数字信号处理课程设计 1

（一）用窗函数法设计FIR数字滤波器 一、设计题目用窗函数法设计FIR数字低通滤波器 二、设计目的 1. 熟悉设计线性相位数字滤波器的一般步骤。 2. 掌握用窗函数法设计FIR数字滤波器的原理和方法。 3. 熟悉各种窗函数的作用以及各种窗函数对滤波特性的影响。 4. 学会根据指标要求选取合适的窗函数。 三、设计原理 窗函数法又称为傅里叶级数法，FIR数字滤波器的设计问题就是要所设计的FIR数字滤波器的响应H(ejw)去逼近所要求的理想滤波器的响应Hd(ejw)。从单位取样响应序列来看，就是使所设计的滤波器的h(n)逼近理想单位取样响应序列hd(n)。而且Hd(ejw)=逐段恒定的，且在频带边界处有不连续点，因此序列hd(n)是无限长的，通过直接截取无限长序列以得到有限长序列的办法，可以形象的比喻为h(n)通过一个窗口所看到的一段hd(n)。因此，h(n)也可以表达为hd(n)和一个窗函数w(n)的乘积，h(n)=w(n)hd(n)。这里的窗函数就是矩形序列RN(n)。 四、实现方法 用MATLAB编程实现给定指标要求的滤波器设计 五、设计内容及要求 1、各窗函数图（假设N=67;） N=67;

n=0:N-1; wn1=ones(1,N); stem(n, wn1);矩形窗 figure; wn2=hamming(N); stem(n, wn2);海明窗 figure; wn3=BARTLETT(N); stem(n, wn3);巴特列特 figure; wn4= Hanning(N); stem(n, wn4);汉宁窗 将窗函数分别画出来 2、计算理想低通滤波器单位冲激响应的源程序function[hd]=ideal(wc,N) q=(N-1)/2; n=0:N-1; m=n-q+eps; hd=sin(wc*m)./(pi*m); 3、计算频率响应的源程序 function[H]=fr(b,a,w); m=0:length(b)-1; l=0:length(a)-1; num=b*exp(-j*m'*w); den=a*exp(-j*l'*w); H=num./den;

DSP实验报告

学校代码学号分类号密级 DSP实验报告 院系名称 专业名称 年级 学生姓名 指导老师 年月日

实验一数据存储实验 一、实验目的 1. 掌握 TMS320C54X 程序空间的分配； 2. 掌握 TMS320C54X 数据空间的分配； 3. 能够熟练运用TMS320C54X 数据空间的指令。 二、实验设备 计算机，CCS 3.1版软件，DSP仿真器，E300实验箱，DSP-54XP CPU板。 三、实验系统相关资源介绍 1. 本实验指导书是以TMS320VC5416为例，介绍其相关的内部和外部存储器资源。对于其他类型的CPU请参考查阅相关的数据手册。） 下面给出TMS320VC5416的存储器分配表： 对于数据存储空间而言，映射表相对固定。值得注意的是内部寄存器都映射到数据存储器空间内。因此在编程时这些特定的空间不能作其他用途。 对于程序空间而言，其映射表和CPU 的工作模式有关。当MP/MC 引脚为高电平时，CPU 工作在微处理器模式；当MP/MC引脚为低电平时，CPU工作在微计算机模式。具体的MP和MC模式下的程序和数据映射关系如上图所示。 2. 样例程序实验操作简单说明： 本实验程序将对0x1000 开始的8 个地址空间，填写入0xAAAA 的数据，然后读出，并存储到以0x1008开始的8个地址空间，在CCS中可以观察DATA存储器空间地址0x1000～0x100F 值的变化。 四、实验步骤与内容 1. 在进行 DSP实验之前，需先连接好仿真器、实验箱及计算机，连接方法如下所示： 2. E300 底板的开关SW4 的第1位置ON，其余位置OFF，SW5全部置ON，其余开关不做设置要求。 3. 上电复位 在硬件安装完成后，确认安装正确、各实验部件及电源连接无误后，启动计算机，接通仿真器电源，此时，仿真器上的“红色指示灯”应点亮，否则DSP开发系统与计算机连接存

DSP课程设计报告—刘雅琪

一、课程设计的目的和要求 1.1课程设计目的： 本课程是DSP技术类课程配套的课程设计，要求学生通过高级语言或汇编语言编程实现较复杂的功能。通过课程设计，使学生加深对DSP芯片 TMS320C54的结构、工作原理的理解，获得DSP应用技术的实际训练，掌握设计较复杂DSP系统的基本方法。 1.2课程设计要求 1、认真查阅资料 2、课程设计前认真预习 3、遵守课程设计时间安排 4、认真保质保量完成设计要求 5、认真书写报告 二、系统功能介绍及总体设计方案 2．1 功能介绍 随着信息技术和计算机技术的飞速发展, 数字信号处理技术在众多领域得到广泛应用。数字滤波器由于其精度高、稳定性好、使用灵活等优点, 广泛应用在各种数字信号处理领域。数字滤波器根据冲击响应函数的时域特性, 可以分为FIR （有限长冲激响应滤波器）和IIR（无限长冲激响应滤波器）。FIR滤波器与IIR 滤波器相比,具有严格的线性相位,幅度特性可任意等优点。而且, FIR 滤波器的单位抽样响应是有限长的, 故一定是稳定的, 他又可以用快速傅里叶变换（FFT）算法来实现过滤信号，可大大提高运算效率。 本课程设计的是一个等波纹FIR 低通滤波器，其具体参数为：采样频率 F s=1000Hz，通带频率F pass=150Hz截止频率F sg=250Hz，通带衰减A pass =0.5dB 阻带衰减A stop=80dB。 2．2 总体设计方案： 先进行Matlab 程序设计产生待滤波数据（借助设计工具FDATOOL 产生设计系数），将其导入CCS在CCS!进行仿真调试运行，得到了输入和输出的波形及

其频谱。 三、主要设计内容和步骤 3?1 FIR 数字滤波器的原理分析 3.1.1FIR 数字滤波器 数字滤波器原理一般具有如下差分方程 N A N -1 y(n) f a k X( n- k) 、b k y( n-k) k z 0 k =0 式中x(n)为输入序列，y(n)为输出序列，兔和b k 为滤波器系数，N 是滤波器阶 N J 数。当所有的b k 均为零，则有 y (n )八 时(n-k) ⑵ k=0 (2) 式是FIR 滤波器的差分方程，其一般形式为 N -1 y(n)「h k x( n-k) (3) k ￡ 对(3)式进行z 变换，整理后可得FIR 滤波器的传递函数 FIR 的直接型结构: x(n) ——B —— h(0) n "(2) ' b ——■—— 1 h(N ￡3) ?_m h(N ￡2) I ---- ■—1 ? ■ ” ￡ ■ 1 h(N ￡l) z ￡i z ￡i y(n) H(z) = 迩 2 X o -k h)k)z

数字信号处理课程设计(对音乐信号的各种处理)

实验1 1、音乐信号的音谱和频谱观察 ○1使用wavread语句读取音乐信号，获取抽样率； ○2输出音乐信号的波形和频谱，观察现象； ○3使用sound语句播放音乐信号，注意不同抽样率下的音调变化，解释现象。 clear all;close all;clc [a,fs,bit]=wavread('c:\MATLAB6p5\work\陪你一起看草原.wav'); size(a); y1=a(:,1); a1=y1(10000:60000) figure; subplot(2,1,1),plot(a); subplot(2,1,2),plot(a1); x1=resample(a1,2,1); %y=resample(x,p,q)返回量的长度是向量x的p/q倍sound(x1,fs); %sound(a,fs); N1=length(a1); F1=fft(a1,N1); w=2/N1*[0:N1-1]; %频谱图横坐标设置 figure; plot(w,abs(F1)); N2=length(a1); t=0:1/N2:1/N2*(N2-1); title('傅利叶变换'); %傅利叶变换; figure; plot(a1); title('时域波形'); %时域波形;

1，以二倍的抽样率听声音信号时，音乐播放的特别快，像被压缩了，播放的时间比原信号短。 2，以二分之一的抽样率听声音信号时，音乐播放的特别慢，像被拉长了，播放的时间比原信号长。 3，原信号频谱截止频率为0.5*pi 实验2 2、音乐信号的抽取（减抽样） ○1观察音乐信号频率上限，选择适当的抽取间隔对信号进行减抽样（给出两种抽取间隔， 代表混叠和非混叠）； ○2输出减抽样音乐信号的波形和频谱，观察现象，给出理论解释； ○3播放减抽样音乐信号，注意抽样率的改变，比较不同抽取间隔下的声音，解释现象。 clear all;close all;clc [a,fs,bit]=wavread('c:\MATLAB6p5\work\陪你一起看草原.wav'); size(a); y1=a(:,1); a1=y1(10000:60000)