数字信号处理实验二().doc
文本预览下载声明
实验二 离散时间系统的时域分析
实验室名称： 实验时间：
姓 名： 学号： 专业： 指导教师：
成绩
教师签名： 年 月 日 一、实验目的
1. 通过MATLAB仿真一些简单的离散时间系统，并研究它们的时域特性。
2. 通过MATLAB仿真滤波器，研究它们的性质。
3. 用MATLAB仿真研究离散时间系统的性质。
二、实验内容
Q2.1 M=2运行内容一的程序代码，生成输入x[n]=s1[n]+s2[n]的输出信号。
Q2.2若线性时不变系统由y[n]=0.5(x[n]+x[n-1])变成y[n]=0.5(x[n]-x[n-1])对输入x[n]=s1[n]+s2[n]的影响是什么？
Q2.4 修改程序P2.1，用一个长度为101、最低频率为0、最高频率为0.5的扫频正弦信号作为输入信号（见程序P1.7），计算其输出信号。
Q2.7 运行程序P2.3，对由加权输入得到的y[n]与在相同权系数下输出y1[n]和y2[n]相加得到的yt[n]进行比较，这两个序列是否相等？该系统是线性系统吗？
Q2.12 运行程序P2.4并比较输出序列y[n]和yd[n-10]。
Q2.19 运行程序P2.5，生成式（2.15）所给离散时间系统的冲击响应。
Q2.20 修改程序P2.5产生如下因果线性时不变系统的冲激响应的前45个样本：
y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3]=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3]。
Q2.23 运行程序P2.6，计算输出序列y[n]和y2[n]以及差值信号d[n]。y[n]和y2[n]相等吗？
Q2.28 运行程序P2.7，对序列h[n]和x[n]求卷积，生成y[n],并用FIR滤波器h[n]对输入x[n]滤波，求得y1[n]。
三、实验器材及软件
1. 微型计算机1台
2. MATLAB 7.0软件
四、实验原理
1.离散时间系统中信号运算方法和基本性质；
2.离散时间系统的MATLAB表示和基本运算；
3. MATLAB处理离散时间系统的方法。
五、实验步骤
根据实验题目要求进行分析；
2. 运用所学知识用MATLAB编程实现题目要求；
3. 对结果进行分析总结。
六、实验记录（数据、图表、波形、程序等）
Q2.1 程序代码：
%产生输入信号
n = 0:100;
s1 = cos(2*pi*0.05*n); %一个低频正弦
s2 = cos(2*pi*0.47*n); %一个高频正弦
x = s1+s2;
%滑动平均滤波器的实现
M = input(Desired length of the filter = );
num = ones(1,M);
y = filter(num,1,x)/M;
clf;
%显示输入和输出信号
subplot(2,2,1);
plot(n, s1);
axis([0, 100, -2, 2]);
xlabel(Time index n); ylabel(Amplitude);
title(Signal #1);
subplot(2,2,2);
plot(n, s2);
axis([0, 100, -2, 2]);
xlabel(Time index n); ylabel(Amplitude);
title(Signal #2);
subplot(2,2,3);
plot(n, x);
axis([0, 100, -2, 2]);
xlabel(Time index n); ylabel(Amplitude);
title(Input Signal);
subplot(2,2,4);
plot(n, y);
axis([0, 100, -2, 2]);
xlabel(Time index n); ylabel(Amplitude);
title(Output Signal); axis;
显示结果：
Q2.2 代码（将输出信号的代码做出修改，其余代码与Q2.1相同）：
n=0:100;
s1=cos(2*pi*0.05*n);
s2=cos(2*pi*0.47*n);
x=s1+s2;
M=input( Desired length of the filter = );
num=(-1).^[0:M-1];
y=filter(num,1,x)/M;
clf;
%显示输入和输出信号
subplot(2,2,1);
plot(n, s1);
axis([0, 100, -2, 2]);
xlabel(Time i
显示全部