DSP实验报告.doc

实验报告
实验名称____采样、系统性质及滤波;__
数字滤波器的频率响应和
处理算法实现;
___加窗和离散傅式变换_;____
___数字滤波器设计。_________
课程名称数字信号处理A(双语)
班级学号_____B09010533___________
姓名_____赵尚文______________
开课时间 2011/2012学年,第二学期
实验一:采样、系统性质及滤波
,实验内容:
1 实验目的和任务。
(1)熟悉MATLAB的主要操作命令。
(2)学会简单的矩阵输入和数据读写。
(3)掌握简单的绘图命令。
(4)用MATLAB编程并学会创建函数。
(5)观察采样引起的混叠。
(6)判别离散时间系统的时不变性。
(7)卷积计算。
二、观察采样引起的混叠。(来源:p57 习题1-3)
设模拟信号为,t的单位为毫秒(ms)。
设采样频率为3kHz,确定与混叠的采样重建信号。
画出和在范围内的连续波形。(因数字计算机无法真正画出连续波形,可用较密的离散点的连线来近似。)
分别用和在两信号波形上标记出3kHz采样点。
两信号波形是否相同?采样后的两序列是否相同?
三、判别离散时间系统的时不变性。(来源:p105 )
设输入序列为,系统实现对的抽取。
设。取延迟量D(例如D=30)。记,画出、的序列波形。
编程求出系统对的响应以及对的响应
画出、的波形。
该系统是否为时不变的?
四、利用卷积计算信号通过FIR滤波器的输出,并观察输出信号的input-on暂态、input-off
暂态和稳态阶段。(来源:p144 )
考虑两个滤波器,,;输入为周期方波,第一个周期内。
1..分别画出通过两个滤波器的输出、的波形,。
2. 计算图中稳态部分的响应值。

1、实验程序:
% =============
% problem 1
% =============
clear
% estimate x(t) and xa(t) with a much higher sampling freq. 'fs1'
time_period=6; % unit: ms
fs1=50; % unit: kHz
T1=1/fs1; % unit: ms
n1=0:fix(time_period/T1);
x=cos(5*pi*n1*T1)+4*sin(2*pi*n1*T1).*sin(3*pi*n1*T1);
xa=cos(pi*n1*T1);
% obtain x(nT) and xa(nT) with given sampling freq. 'fs'
fs=3;
T=1/fs;
n=0:fix(time_period/T);
x_sample=cos(5*pi*n*T)+4*sin(2*pi*n*T).*sin(3*pi*n*T);
xa_sample=cos(pi*n*T);
figure,plot(n1*T1,x,'r',n1*T1,xa,'b',n*T,x_sample,'ro'),
hold on, stem(n*T,xa_sample,'b:x')
legend('x(t)','xa(t)','x(nT)','xa(nT)'),xlabel('t(ms)')
% =============
% problem 2
% =============
clear
% plot x(n) and x(n-D)
D=30;
N=500;
n=1:N;
x=sin(2*pi/100*n);
for n=1:N+D,
if (n-D)<=0, xD(n)=0;
else xD(n)=x(n-D); end
end
figure,subplot(2,1,1),
plot(1:N,x,'r:',1:length(xD),xD,'b'),
legend('x(n)','xD(n)'),,xlabel('n')
% plot y(n) and yD(n)
for n=1:fix(N/2)
y(n)=x(2*n);
end
for n=1:length(y)+D,
if (n-D)<=0, y_delay(n)=0;
else y_delay(n)=y(n-D); end
end
for n=1:fix(length(xD)/2)
yD(n)=xD(2*n);
end
subplot(2,1,2),
plot(1:length(y)