重庆大学光电学院数字信号处理实验

五、实验过程原始记录（实验源程序、图标、计算等） 1、用matlab 编制信号产生子程序，产生以下典型信号序列。

()100处有一单位脉冲，1在,)(单位脉冲序列1≤≤=n n

n δ

()()100，，1处及以后为4，在0前为4在,单位阶跃序列2≤≤==n n n

n u

()()n

R

8

矩形序列

3

()()9

3

5

sin

3

正弦序列

4≤

≤

⎪⎪

⎭

⎫

⎝

⎛

+

•

=n

n

n

x

π

π

()()90)4(5)3(4)2(3)1(2)(序列5≤≤-+-+-+-+=n n n n n

n n x δδδδδ

()()90)3(2)2()1(2)(序列6≤≤-+-+-+=n n n

n n n h δδδδ

2、利用MATLAB 编程完成序列(5)和(6)的移位、反转、加法、乘法等运算，并绘制运算后序列的波形。

()())4(5)3(4)2(3)1(2)(序列5-+-+-+-+=n n n n n n x δδδδδ ()())3(2)2()1(2)(序列6-+-+-+=n n n

n n h δδδδ

3、利用MATLAB编制一个计算两个序列线性卷积的通用程序，计算上述两序列，并绘制卷积后序列的波形。

六、实验结果分析

1．如何产生方波信号序列和锯齿波信号序列？

答：方波信号由square函数可得，锯齿波信号由循环语句产生。

%方波

T=0:0.001:2*pi; %方波的时间向量

y=square(6*T); %方波w=6,周期为T=2*pi/6

subplot(2,1,1);

plot(T,y); %画方波

title('方波');

axis([T(1)-1 T(end)+1 -2 2]); %坐标轴区域

hold on

一个序列

x (n )的离散时间傅里叶变换就是它的频谱函数。

()()

()1

20

N j nk N k

n N

X k X e

x n W ω

πω-====∑

四、使用仪器、器材

计算机，MATLAB

五、实验过程原始记录（实验源程序、图标、计算等）

（1）用matlab 编程产生并画出信号x1(n)、x2(n)、x3(n)、x4(n)、x5(n)。

⎪⎩

⎪

⎨⎧≤≤-≤≤+==n

n n n n n x n R n x 其它,074,

830,1)()

()(241⎪⎩⎪

⎨⎧≤≤-≤≤-=n

n n n n n x 其它,

074,

330,

4)(3 n n x 4

cos

)(4π

=n n x 8

sin

)(5π

=

（2）用matlab 编制FFT 函数对上述信号进行频谱分析，并画出上述信号谱图。

3、设计一个BW型低通DF，满足：通带截频fp=100Hz，衰耗满足Apmax=3dB，阻带截频fT=400Hz,衰耗ATmin=15dB,取样频率fs=2000Hz。写出其模拟滤波

器方程及数字滤波器方程。

4、设计一个BW型高通DF，满足：通带截频fp=400Hz，衰耗满足Apmax=3dB，

阻带截频fT=350Hz,衰耗ATmin=15dB,取样频率fs=1000Hz。写出其模拟滤波器方程及数字滤波器方程。

5、设计一个CB型带通DF，满足：通带边界频率为100Hz～500Hz，通带衰耗小

于3dB，过渡带宽20Hz，阻带衰耗大于15dB，取样频率fs=2000Hz。写出其模拟滤波器方程及数字滤波器方程。

六、实验结果分析

1、与脉冲响应不变设计法相比较，简述双线性变换设计法的优缺点。

五、实验过程原始记录（实验源程序、图标、计算等）

1.语音载入和播放

2.频谱

3、混音器

4、均衡器

5、滤波器

6、变声器

7、效果图