数字信号处理实验讲义
??(n?p)?eq,0?n?15xa(n)??
??0,else
2n=0:15; %定义序列的长度是15
p=8; q=2; x=exp(-1*(n-p).^2/q); %利用fft 函数实现富氏变换 close all; subplot(3,1,1); stem(abs(fft(x)))
p=8; q=4; x=exp(-1*(n-p).^2/q); %改变信号参数，重新计算 subplot(3,1,2); stem(abs(fft(x))) p=8; q=8; x=exp(-1*(n-p).^2/q); subplot(3,1,3); stem(abs(fft(x)))
4、衰减正弦序列
?e??nsin2?fn,0?n?15x(n)??b
0,else?n=0:15; %定义序列的长度是15
a=0.1; f=0.0625; x=exp(-a*n).*sin(2*pi*f*n); close all; subplot(2,1,1); stem(x); subplot(2,1,2); stem(abs(fft(x)))
5、三角波序列
?n?1,0?n?3?x(n)??8?n,4?n?7 c?0,else?for i=1:4 %设置信号前4 个点的数值
x(i)=i; %注意：MATLAB 中数组下标从1 开始 end
for i=5:8 %设置信号后4 个点的数值 x(i)=9-i; end
close all; subplot(2,1,1); stem(x); %绘制信号图形 subplot(2,1,2); stem(abs(fft(x,16))) %绘制信号的频谱
6、反三角序列
?4?n,0?n?3?xd(n)??n?3,4?n?7
?0,else?for i=1:4 %设置信号前4 个点的数值
x(i)=5-i; %注意：MATLAB 中数组下标从1 开始
22
数字信号处理实验讲义
end
for i=5:8 %设置信号后4 个点的数值
x(i)=i-4; end
close all;
subplot(2,1,1); stem(x); %绘制信号图形
subplot(2,1,2); stem(abs(fft(x,16))) %绘制信号的频谱
23
数字信号处理实验讲义
实验五 用双线性变换法设计IIR数字滤波器
一、实验目的
1、掌握双线性变换法及脉冲相应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。
2、观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。
3、熟悉Butterworth滤波器、Chebyshev滤波器和椭圆滤波器的频率特性。
二、实验原理与方法
1、脉冲响应不变法
用数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应ha(t),让h(n)正好等于ha(t)的采样值，即
h(n)?ha(nT)
其中T为采样间隔，如果以Ha(S)及H(z)分别表示ha(t)的拉式变换及h(n)的Z变换，则
1?2?H(Z)|Z?esT??Ha(S?jm)
Tm???T 2、双线性变换法
S平面与z平面之间满足以下映射关系：
TS21?z?12,S???j?;z?rej?S?,z????1 TST1?z1?21?s平面的虚轴单值地映射于z平面